Capture–Annihilation Equilibrium in Astrophysical Systems

• Capture–annihilation equilibrium is defined as the steady-state balance where particle injection (capture) is exactly offset by removal (annihilation), typically modeled by dN/dt = C - A N².
• In astrophysical systems, this equilibrium underpins dark matter accumulation in stars and governs key observables like annihilation luminosity and thermal heating through measurable capture and annihilation coefficients.
• The concept also applies to non-astrophysical contexts such as exciton capture in semiconductors, highlighting its broad relevance in understanding driven dissipative systems.

Searching arXiv for the cited papers to ground the article in current literature. arXiv search query: (Bell et al., 2023) Capture–annihilation equilibrium is the steady-state regime of an open many-body system in which a capture or injection process is balanced by annihilation or removal. In the astrophysical literature, the canonical example is the accumulation of dark matter in a gravitating body, where the total number $N(t)$ obeys

$\frac{dN}{dt}=C-A\,N^2,$

with $C$ the capture rate and $A\,N^2$ twice the annihilation rate, $\Gamma_{\rm ann}=\tfrac12 A N^2$. At late times, the equilibrium condition $dN/dt=0$ gives

$$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$

This same balance structure reappears in stars, planets, neutron stars, and in formally analogous non-astrophysical settings such as exciton capture at defects and kinetic annihilation models (Bell et al., 2023).

1. Formal definition and governing equations

In dark-matter capture problems, the equilibrium variable is the accumulated population inside the host object. Bell–Busoni–Robles–Virgato formulate the neutron-star problem through $dN/dt=C-A N^2$, with $C$ determined by dark-matter scattering on neutron-star constituents and $A$ set by the annihilation cross section together with the spatial distribution of the captured population (Bell et al., 2023). Croon & Sakstein use the same structure for massive stars undergoing pair-instability, writing

$\frac{dN}{dt}=C-A\,N^2,$0

where $\frac{dN}{dt}=C-A\,N^2,$1 is defined so that $\frac{dN}{dt}=C-A\,N^2,$2 equals twice the total annihilation rate (Croon et al., 2023).

When evaporation is relevant, the minimal equation acquires a linear sink term. In low-mass red-giant-branch stars, Lopes & Lopes write

$\frac{dN}{dt}=C-A\,N^2,$3

with $\frac{dN}{dt}=C-A\,N^2,$4 the evaporation coefficient and $\frac{dN}{dt}=C-A\,N^2,$5 twice the annihilation rate per pair. In most stellar contexts considered there, and in particular on the RGB for $\frac{dN}{dt}=C-A\,N^2,$6–$\frac{dN}{dt}=C-A\,N^2,$7\,GeV, evaporation is negligible, so the standard quadratic balance is recovered (Lopes et al., 2021).

An analogous structure appears in condensed-matter defect physics. Wang et al. describe the exciton density $\frac{dN}{dt}=C-A\,N^2,$8 in a monolayer transition-metal dichalcogenide through

$\frac{dN}{dt}=C-A\,N^2,$9

where $C$0 is generation, $C$1 is a single-body Auger-capture loss, and $C$2 is a two-body Hartree–Fock contribution (Wang et al., 2014). The common feature across these systems is that equilibrium is defined by the vanishing of the net population derivative, even though the microscopic origin of the sink term differs.

2. Capture and annihilation coefficients

The capture coefficient encodes scattering kinematics, the phase-space distribution of incident particles, and the internal structure of the host. In neutron stars, the optically thin capture rate is written as an integral over the incident speed distribution $C$3, the stellar radius, the Schwarzschild metric function $C$4, and the local down-scattering rate $C$5. In the optically thick or geometric limit, the star acts as a gravitationally focused hard sphere, giving a geometric capture rate $C$6 (Bell et al., 2023). The microscopic input is the DM–nucleon cross section

$C$7

with $C$8 potentially carrying momentum dependence $C$9 (Bell et al., 2023).

For stellar applications outside the neutron-star regime, the capture rate is typically expressed as a radial integral over the object combined with a velocity integral over the halo distribution. Lopes & Lopes write

$A\,N^2$0

where $A\,N^2$1 and $A\,N^2$2 is the per-target scattering rate into bound orbits (Lopes et al., 2021). In the Sun, the same logic appears with hydrogen as the dominant SD target and an effective capture cross section $A\,N^2$3 built from operator-dependent form factors (Liang et al., 2013).

The annihilation coefficient converts microscopic annihilation into a macroscopic depletion rate. For an isothermal neutron-star distribution of radius $A\,N^2$4,

$A\,N^2$5

with

$A\,N^2$6

for the fiducial scaling quoted by Bell–Busoni–Robles–Virgato (Bell et al., 2023). In Croon & Sakstein, the same quantity is written as

$A\,N^2$7

where $A\,N^2$8 is defined from the quasi-thermal dark-matter profile (Croon et al., 2023). Catena likewise writes for the Earth

$A\,N^2$9

with effective volumes $\Gamma_{\rm ann}=\tfrac12 A N^2$0 determined by the thermalized central distribution (Catena, 2016).

3. Equilibrium solution and equilibration timescale

The simplest capture–annihilation equation admits a closed analytic solution. DarkCapPy gives

$\Gamma_{\rm ann}=\tfrac12 A N^2$1

and therefore

$\Gamma_{\rm ann}=\tfrac12 A N^2$2

For $\Gamma_{\rm ann}=\tfrac12 A N^2$3, $\Gamma_{\rm ann}=\tfrac12 A N^2$4, so equilibrium is effectively reached (Green et al., 2018). The Sun formulation is identical: $\Gamma_{\rm ann}=\tfrac12 A N^2$5 with equilibrium when $\Gamma_{\rm ann}=\tfrac12 A N^2$6 (Liang et al., 2013).

The physical interpretation is consistent across applications: in equilibrium, the annihilation rate of pairs is $\Gamma_{\rm ann}=\tfrac12 A N^2$7, because each annihilation event removes two captured particles. This point is stated explicitly in the Earth and Sun treatments [(Green et al., 2018); (Liang et al., 2013)].

The timescale $\Gamma_{\rm ann}=\tfrac12 A N^2$8 determines whether steady state is astrophysically relevant. In neutron stars, for scattering cross sections that saturate the capture rate, Bell–Busoni–Robles–Virgato find that capture–annihilation equilibrium is typically reached on a timescale of less than $\Gamma_{\rm ann}=\tfrac12 A N^2$9 year for vector interactions and $dN/dt=0$0 years for scalar interactions (Bell et al., 2023). In low-mass RGB stars, Lopes & Lopes quote

$dN/dt=0$1

which is much shorter than the RGB lifetime $dN/dt=0$2–$dN/dt=0$3\,yr, so equilibrium holds throughout ascent to the TRGB (Lopes et al., 2021). By contrast, Earth capture can fail to equilibrate within the age of the Earth unless annihilation is enhanced (Green et al., 2018).

4. Suppression effects, thermalization, and common misconceptions

A central technical point is that equilibrium is controlled by $dN/dt=0$4 and $dN/dt=0$5, not by a requirement that all microscopic processes be unsuppressed. In neutron stars, momentum-suppressed scattering of the form $dN/dt=0$6 reduces capture through $dN/dt=0$7, and $dN/dt=0$8-wave annihilation introduces an additional factor $dN/dt=0$9 in $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$0. Nevertheless, the algebraic equilibrium condition

$$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$1

still holds; only the numerical value of $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$2 and the equilibration time shift because $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$3 and/or $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$4 are reduced (Bell et al., 2023).

This is closely related to operator dependence in solar and terrestrial WIMP capture. In the Sun, the six spin-dependent operators $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$5 and $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$6 yield different capture rates and hence different $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$7. The paper on direct detection and solar capture emphasizes that $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$8 has no $$N_{\rm eq}=\sqrt{\frac{C}{A}}, \qquad \Gamma_{\rm ann}=\frac{C}{2}.$$9-suppression and typically gives the largest $dN/dt=C-A N^2$0, while $dN/dt=C-A N^2$1, $dN/dt=C-A N^2$2, and $dN/dt=C-A N^2$3 carry extra powers of $dN/dt=C-A N^2$4, and $dN/dt=C-A N^2$5 is suppressed by $dN/dt=C-A N^2$6 (Liang et al., 2013). In the Earth, DarkCapPy shows that a light mediator can multiply the annihilation coefficient by a thermal average of the Sommerfeld factor, $dN/dt=C-A N^2$7, with boosts of $dN/dt=C-A N^2$8–$dN/dt=C-A N^2$9 in favorable regions of parameter space (Green et al., 2018).

A frequent misconception is that capture–annihilation equilibrium requires complete thermalization of the captured population. Bell–Busoni–Robles–Virgato explicitly show that this is not necessary in neutron stars: the annihilation zone is only of order metres in size, much smaller than the kilometers-scale orbits of partially thermalised dark matter, and in most cases $C$0. Their conclusion is that maximal annihilation heating can be achieved without complete thermalization of the captured dark matter (Bell et al., 2023). This suggests that the relevant criterion is localization into the annihilation region rather than full kinetic equilibration with the surrounding medium.

5. Astrophysical realizations

In old, isolated neutron stars, capture and subsequent annihilation can heat the star. Kinetic heating requires sufficient scattering to deposit the infalling dark matter’s kinetic energy, and appreciable annihilation heating additionally requires capture–annihilation equilibrium. Bell–Busoni–Robles–Virgato find that both can be achieved for all types of dark matter–baryon interactions considered, including momentum- or velocity-suppressed cases (Bell et al., 2023).

In massive stars subject to pair-instability, Croon & Sakstein assume capture–annihilation equilibrium and model local energy injection in a local thermal equilibrium approximation. They find significant changes to the masses of astrophysical black holes formed by (pulsational) pair-instability supernovae when the ambient dark matter density satisfies $C$1. For $C$2 GeV, dark matter is primarily confined to the core, prolongs core helium burning, increases oxygen production, and drives stronger pulsations, leading to lighter black holes. For $C$3 GeV, dark matter is significant in the envelope, can support the star through annihilation energy, and allows heavier black holes to form by avoiding pair-instability (Croon et al., 2023).

In low-mass RGB stars, capture–annihilation equilibrium makes the dark-matter population approximately constant during the RGB phase and mostly independent of stellar mass. Lopes & Lopes find that the equilibrium number scales as

$C$4

for the models quoted, and that the annihilation luminosity

$C$5

can reduce the helium-core mass at flash and dim the tip of the red giant branch. For $C$6 GeV, deviations from the standard TRGB luminosity of $C$7 can be achieved under conditions that can be realistic in the inner parsec of the Milky Way (Lopes et al., 2021).

The Sun and Earth provide complementary realizations. In the Sun, numerical evaluation shows $C$8 through much of the parameter space within current bounds, so $C$9 is generally valid for the SD operators studied (Liang et al., 2013). In the Earth, Catena and the DarkCapPy framework emphasize that equilibration is more difficult; resonances in capture occur near $A$0 for terrestrial elements, and Sommerfeld enhancement can be decisive for reaching $A$1 (Catena, 2016, Green et al., 2018).

6. Analogous formulations beyond dark-matter astrophysics

The phrase “capture–annihilation equilibrium” also describes mathematically similar balance laws outside gravitational dark-matter systems. In monolayer metal dichalcogenides, Wang et al. analyze Auger-mediated exciton capture at defects. Because excitons are tightly bound, $A$2 is large, and the correlated electron–hole structure enhances defect-assisted capture. The resulting rates can be $A$3–$A$4 times larger than for uncorrelated carriers at carrier densities in the $A$5–$A$6 range, with capture times in the sub-picosecond to a few picoseconds range (Wang et al., 2014). The steady-state condition is

$A$7

with positive-root solution

$A$8

Here the linear and quadratic loss channels play roles analogous to one-body escape and two-body annihilation in astrophysical problems (Wang et al., 2014).

A different analogue arises in the ideal Rayleigh gas with annihilation. Nota–Winter–Lods derive a linear Boltzmann equation in the Boltzmann–Grad limit in which the tagged particle is removed upon collision with annihilating obstacles and elastically scattered by non-annihilating ones. In the spatially homogeneous steady state, the balance condition becomes

$A$9

or equivalently a linear Fredholm equation for $\frac{dN}{dt}=C-A\,N^2,$00 (Nota et al., 2019). This is not a quadratic population law, but it retains the defining feature of equilibrium as an exact balance between inflow and annihilative loss. A plausible implication is that capture–annihilation equilibrium is best understood as a structural principle of driven dissipative systems rather than as a concept restricted to WIMP phenomenology alone.