Two-Centered Black Hole Index

• The two-centered black hole index is a protected contribution from a bound state of two BPS centers in 4D supersymmetric string compactifications, encoding wall-crossing phenomena.
• It is computed via localization on the compact S² phase space and geometric quantization, incorporating both bosonic and fermionic corrections along with error-function completions.
• This index is pivotal in understanding primitive wall-crossing formulae and modularity, unifying microscopic string theory descriptions with macroscopic supergravity analyses.

Searching arXiv for the cited papers to ground the article in the published literature. arXiv search query: (Pioline et al., 11 Jul 2025) Multi-centered Black Hole Quantum Mechanics and Generalized Error Functions The two-centered black hole index is the protected contribution of a bound state of two BPS centers to the total BPS index in four-dimensional supersymmetric string compactifications. In the formulations relevant to type II Calabi–Yau compactifications, $\mathcal N=2$ supergravity, heterotic string theory on $T^6$, and $\mathbb Z_N$ CHL models, it is determined by the relative dynamics of two mutually non-local charges, by localization on the classical phase space of two-center solutions, or by extracting the polar part of meromorphic Siegel or Jacobi generating functions. Its central roles are to encode the jump of the index across walls of marginal stability, to supply the two-body factor in primitive wall-crossing formulae, and to provide the non-holomorphic completion needed for modularity in mock Jacobi and related automorphic structures (Pioline et al., 11 Jul 2025, Manschot et al., 2011, Manschot et al., 2010).

1. Physical set-up and classical two-center geometry

For two mutually non-local BPS dyons of charges $\gamma_1,\gamma_2$ in four dimensions, the center-of-mass decouples, and the relative motion is described by an $\mathcal N=4$ supersymmetric quantum mechanics on $\mathbb R^3$ with coordinate $\vec r=\vec x_1-\vec x_2$. The effective Lagrangian can be written schematically as

$$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$

where $m=m_1m_2/(m_1+m_2)$ is the reduced mass,

$$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$

$T^6$0 is the real FI-parameter, $T^6$1 is the Dirac monopole potential, and $T^6$2 are two complex fermions. Classical BPS bound states solve $T^6$3, so

$T^6$4

and their moduli space is the two-sphere $T^6$5 of radius $T^6$6 (Pioline et al., 11 Jul 2025).

In four-dimensional $T^6$7 supergravity, the same structure appears through Denef’s equilibrium condition. For two centers of charges $T^6$8, the relative separation is fixed,

$T^6$9

so after removing the overall center-of-mass motion the allowed relative orientations form a two-dimensional symplectic manifold $\mathbb Z_N$0. The pull-back of the canonical two-form gives

$\mathbb Z_N$1

and the moment map for rotations about the $\mathbb Z_N$2-axis is

$\mathbb Z_N$3

Equivalently, the classical angular momentum vector is

$\mathbb Z_N$4

These formulations identify the two-centered problem as the quantization of a compact $\mathbb Z_N$5 phase space in the absence of scaling solutions (Manschot et al., 2011, 0807.4556).

2. Quantization, localization, and the configurational factor

Quantizing the two-center phase space produces the universal configurational factor multiplying the single-center indices carried by the two constituents. In geometric quantization, one views $\mathbb Z_N$6 as a compact Kähler manifold, with stereographic coordinate

$\mathbb Z_N$7

Kähler potential

$\mathbb Z_N$8

and prequantum line bundle $\mathbb Z_N$9. The holomorphic basis is

$\gamma_1,\gamma_2$0

Without spin one gets $\gamma_1,\gamma_2$1 states, while including the fermionic $\gamma_1,\gamma_2$2 correction gives exactly $\gamma_1,\gamma_2$3 states. Since $\gamma_1,\gamma_2$4, the two-center configurational index is

$\gamma_1,\gamma_2$5

Including the sign from fermion zero modes and the horizon contributions $\gamma_1,\gamma_2$6, the full two-center BPS index becomes

$\gamma_1,\gamma_2$7

(0807.4556).

The same factor follows from equivariant localization on $\gamma_1,\gamma_2$8. The only fixed points of the $\gamma_1,\gamma_2$9-action are the north and south poles, $\mathcal N=4$0, corresponding to the two collinear configurations. Localization gives the refined two-center factor

$\mathcal N=4$1

and in the unrefined limit

$\mathcal N=4$2

Accordingly, for two distinct centers,

$\mathcal N=4$3

The factor

$\mathcal N=4$4

is the character of the spin-$\mathcal N=4$5 representation carried by the relative motion (Manschot et al., 2011, Manschot et al., 2010).

3. Supersymmetric quantum mechanics and the error-function completion

The two-body supersymmetric quantum mechanics admits a direct localization computation of the refined Witten index

$\mathcal N=4$6

By the standard argument for supersymmetric $\mathcal N=4$7-models with a potential, $\mathcal N=4$8 localizes onto constant modes. One finds the finite-dimensional integral

$\mathcal N=4$9

The fermion integral yields

$\mathbb R^3$0

and after passing to spherical coordinates and performing the Gaussian $\mathbb R^3$1-integral one obtains

$\mathbb R^3$2

In the unrefined limit,

$\mathbb R^3$3

Here

$\mathbb R^3$4

(Pioline et al., 11 Jul 2025).

The purely holomorphic or step-function index counts normalizable bound states and is obtained by replacing $\mathbb R^3$5,

$\mathbb R^3$6

The difference $\mathbb R^3$7 is the continuum or spectral-asymmetry piece,

$\mathbb R^3$8

with standard integral representation

$\mathbb R^3$9

Equivalently,

$\vec r=\vec x_1-\vec x_2$0

Introducing $\vec r=\vec x_1-\vec x_2$1 and the stability ratio

$\vec r=\vec x_1-\vec x_2$2

the completed index can be written as

$\vec r=\vec x_1-\vec x_2$3

or equivalently

$\vec r=\vec x_1-\vec x_2$4

The non-holomorphic term exactly cancels the anomaly of the holomorphic term under $\vec r=\vec x_1-\vec x_2$5-duality, so the completion transforms as an ordinary non-holomorphic Jacobi form of weight $\vec r=\vec x_1-\vec x_2$6, with $\vec r=\vec x_1-\vec x_2$7 in the conventions of the cited treatment. Since only one error-function is needed, the result is a depth-one mock Jacobi form (Pioline et al., 11 Jul 2025).

4. Wall-crossing, rational invariants, and continuum scattering states

Across the two-center wall $\vec r=\vec x_1-\vec x_2$8, the total index jumps by the index of the two-center solution, with the centers carrying rational invariants

$\vec r=\vec x_1-\vec x_2$9

The primitive wall-crossing formula is

$$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$0

For the refined index,

$$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$1

one introduces

$$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$2

and the refined jump is

$$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$3

The physical rationale is that Bose–Fermi combinatorics can be traded for Maxwell–Boltzmann statistics once $$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$4 is replaced by $$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$5 (Manschot et al., 2010).

A complementary macroscopic description comes from the supersymmetric quantum mechanics on Taub–NUT for the relative dynamics of two $$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$6-BPS centers. In that setting the localized index contains a temperature-independent contribution from discrete BPS bound states and a temperature-dependent contribution from a spectral asymmetry in the continuum. For each charge sector $$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$7, in the limit $$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$8, $$L=\frac m2(\dot{\vec x}^{\,2}+D^2+2\bar\lambda \dot\lambda)+(-U D+\vec A\cdot \dot{\vec x})+(\nabla U\cdot \bar\lambda \sigma \lambda),$$9,

$m=m_1m_2/(m_1+m_2)$0

The first two terms arise from discrete BPS bound states, while the complementary-error-function and Gaussian term come from a spectral asymmetry in the continuum. Upon the identifications $m=m_1m_2/(m_1+m_2)$1, $m=m_1m_2/(m_1+m_2)$2, and $m=m_1m_2/(m_1+m_2)$3, this agrees with the Fourier coefficients of the non-holomorphic completion $m=m_1m_2/(m_1+m_2)$4. This provides a physical derivation of the completion term $m=m_1m_2/(m_1+m_2)$5 required for modularity (Murthy et al., 2018).

5. Negative discriminant states and bound-state metamorphosis

In heterotic string theory on $m=m_1m_2/(m_1+m_2)$6, and similarly in CHL orbifolds, one can define the T-duality invariants

$m=m_1m_2/(m_1+m_2)$7

and the discriminant

$m=m_1m_2/(m_1+m_2)$8

Classical single-centered $m=m_1m_2/(m_1+m_2)$9-BPS black holes exist only when $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$0; in four-dimensional $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$1 supergravity their horizon area is proportional to $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$2, and no regular attractor flow can be sustained for $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$3. Accordingly, any microscopic index for states with $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$4 must arise from genuinely multi-center configurations (Sen, 2011).

The exact dyon index $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$5 is extracted from Fourier coefficients of $$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$6, and one expects a schematic decomposition

$$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$7

For a bound state of two half-BPS centers,

$$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$8

with

$$U(r)=-\frac12\left(\frac{\kappa}{|r|}-c\right),\qquad \kappa\equiv \langle \gamma_1,\gamma_2\rangle\in \mathbb Z,$$9

where $T^6$00 is the coefficient of $T^6$01 in

$T^6$02

A subtlety arises when one or both centers are “small” charges with square $T^6$03. In that case one imposes the identification rule called bound-state metamorphosis. When one center $T^6$04 satisfies $T^6$05 and $T^6$06, the split $T^6$07 must be identified with the $T^6$08-dual split obtained by shifting $T^6$09, $T^6$10. Explicitly, if $T^6$11, $T^6$12, and $T^6$13, then

$T^6$14

When both centers have square $T^6$15, one must continue to identify an infinite $T^6$16-duality orbit of such splits. After subtracting all two-center contributions subject to these identifications, one finds

$T^6$17

All negative-discriminant states are therefore accounted for by two-centered black holes (Sen, 2011).

6. Generating functions, polar parts, and modular structure

For fixed magnetic invariant $T^6$18, the Fourier–Jacobi coefficient of the inverse Igusa cusp form is

$T^6$19

Physically this generating function decomposes as

$T^6$20

where $T^6$21 captures the single-centered black holes and is holomorphic in $T^6$22, while $T^6$23 captures two-centered bound states and is meromorphic in $T^6$24. The theorem of Dabholkar–Murthy–Zagier gives a unique non-holomorphic completion $T^6$25 such that

$T^6$26

and both completed pieces transform as true Jacobi forms of weight $T^6$27 and index $T^6$28. In particular,

$T^6$29

so the completion of the two-center part is controlled by $T^6$30 (Murthy et al., 2018).

In the heterotic $T^6$31 setting, the two-centered generating function $T^6$32 is defined as the sum of all terms subtracted from $T^6$33 in order to remove the two-center contributions. Near $T^6$34,

$T^6$35

which is the generating function of two-center bound states formed from two half-BPS cores. Subtracting the leading double pole and its $T^6$36-images yields $T^6$37, and one defines

$T^6$38

The Fourier coefficients of $T^6$39 in the attractor chamber are moduli-independent, vanish whenever $T^6$40 or $T^6$41 or $T^6$42, and are free of wall-crossing (Bhand et al., 6 Oct 2025).

In $T^6$43 CHL models, the two-centered generating function is

$T^6$44

and the single-centered generating function is

$T^6$45

The explicit decomposition $T^6$46 subtracts the $T^6$47 poles of $T^6$48, while the last term implements the quotient by the metamorphosis subgroup $T^6$49. For $T^6$50, each of the four sums $T^6$51 is absolutely convergent and uniformly so on compact subsets of the attractor chamber, and $T^6$52 is holomorphic in the domain $T^6$53 with only $T^6$54 poles (Singh, 17 Jun 2026).

7. Finite-temperature supergravity saddles and the disappearance of the bound state

A direct gravitational derivation of the two-centered contribution to the index uses the Euclidean path integral with periodic fermions,

$T^6$55

Because of supersymmetry, the index localizes semiclassically onto finite-temperature BPS saddles in which the geometry caps off smoothly in the bulk at each center. For a two-center split $T^6$56, one introduces doubled poles

$T^6$57

with $T^6$58 fixed by the new attractor equations, and the harmonic function

$T^6$59

The metric and one-form satisfy

$T^6$60

where

$T^6$61

Smoothness imposes the regularity conditions

$T^6$62

These are the finite-temperature analogues of the Denef integrability conditions (Boruch et al., 9 Jul 2025).

For $T^6$63, the regularity conditions can be solved in closed form, leaving a single real modulus $T^6$64 in a finite interval $T^6$65. In the extremal limit,

$T^6$66

which reproduces the Denef formula, and the full four-dimensional moduli space after quotienting translations has topology

$T^6$67

The on-shell action is

$T^6$68

so

$T^6$69

As the asymptotic moduli vary, the interval $T^6$70 shrinks and collapses when

$T^6$71

At that point the two-center index saddle ceases to exist. The net jump in the index is then

$T^6$72

which matches the standard primitive wall-crossing formula. A plausible implication is that the two-centered index admits a unified interpretation across microscopic partition functions, supersymmetric quantum mechanics, and finite-temperature supergravity saddles: in each description, the disappearance of the bound state is controlled by the same wall of marginal stability (Boruch et al., 9 Jul 2025).