---
title: 1/16-BPS Superconformal Index in N=4 SYM
url: https://www.emergentmind.com/topics/1-16-bps-superconformal-index
type: topic
---

# 1/16-BPS Superconformal Index in N=4 SYM

Searching arXiv for recent and foundational papers on the 1/16-BPS superconformal index in 4d \(\mathcal N=4\) SYM.
The \(1/16\)-BPS superconformal index is a protected trace that isolates states or local operators annihilated by a chosen supercharge and its conjugate in radially quantized superconformal field theory. In four-dimensional \(\mathcal N=4\) super-Yang–Mills theory, it is the standard protected quantity associated with the most general single-supercharge BPS sector, and it has been studied both as a cohomological counting problem for actual BPS operators and as an index with \((-1)^F\) cancellations. Its modern significance lies in the conjunction of three themes: weak-coupling \(Q\)-cohomology, finite-\(N\) matrix-integral and representation-theoretic formulas, and the emergence of supersymmetric AdS\(_5\) black-hole entropy in Cardy-like limits [1305.6314].

## 1. Definition and protected-sector structure

A standard weak-coupling definition fixes the supercharge pair
\[
Q=Q_4,\qquad S=S^4,
\]
with radial quantization identifying \(S=Q^\dagger\). The relevant BPS bound is
\[
2\{Q,Q^\dagger\}=A=E-2J_1-H_1-H_2-H_3.
\]
A \(1/16\)-BPS operator is then annihilated by \(Q\) and \(Q^\dagger\), equivalently a cohomology class of \(Q\) on gauge-invariant operators saturating this bound [1305.6314].

In this formulation, the index is obtained from a refined partition function \(Z(x,a,b,u,v,w)\) by setting the bonus-charge fugacity to \(x=-1\),
\[
I(a,b,u,v,w)=Z(-1,a,b,-u,-v,-w),
\]
so that the protected trace becomes an index with the usual \((-1)^F\) grading. The same paper relates this quantity to the standard \(\mathcal N=4\) superconformal index of Kinney et al. by an explicit change of variables [1305.6314].

In later finite-\(N\) and Cardy-limit treatments, the same protected sector is often written in matrix-integral form. One convenient refined notation is
\[
I_N(t,u,v;p,q),
\qquad pq=tuv,
\]
with single-letter index
\[
f(t,u,v;p,q)=1-\frac{(1-t)(1-u)(1-v)}{(1-p)(1-q)},
\]
and unitary-matrix integral
\[
I_N(t,u,v;p,q)=\int_{U(N)} dU \,
\exp\!\left(\sum_{n=1}^{\infty}\frac{f(t^n,u^n,v^n;p^n,q^n)}{n}\,
\mathrm{Tr}(U^n)\mathrm{Tr}[(U^\dagger)^n]\right).
\]
This presentation makes the fugacity-refined \(1/16\)-BPS counting amenable to representation theory, asymptotics, and integrable-system methods [2604.00924].

## 2. Cohomological formulation at weak coupling

The weak-coupling problem was formulated not merely as an index computation with signs, but as an attempt to determine the actual number and structure of \(1/16\)-BPS operators. To organize the cohomology, one restricts to the BPS “letters” that saturate the free BPS bound and packages them into a \((2|3)\) superfield \(\Upsilon(Z)\), depending on auxiliary commuting variables \(z_\alpha\) and anticommuting variables \(\theta_a\), with
\[
\{Q,\Upsilon(Z)\}=\Upsilon(Z)^2.
\]
This algebraic relation is the key simplification behind the cohomological analysis [1305.6314].

The physical \(1/16\)-BPS states are the \(Q\)-cohomology classes on gauge-invariant polynomials in the coefficients of \(\Upsilon(Z)\). The same problem can be rephrased as relative Lie algebra cohomology,
\[
H^\ast(\mathfrak G_N,\mathfrak{sl}_N;\mathbb C),
\qquad
\mathfrak G_N=\mathbb C[z_1,z_2]\otimes \wedge(\theta_1,\theta_2,\theta_3)\otimes \mathfrak{su}_N,
\]
which simultaneously captures the \(Q\)-action and gauge invariance [1305.6314].

This cohomological perspective is conceptually distinct from a bare index calculation. The index is protected under continuous deformations and is insensitive to many bose-fermi cancellations, whereas the cohomology problem attempts to enumerate the actual protected operators. The distinction became important in later discussions of black-hole entropy, because the protected trace and the absolute cohomology count need not show the same growth [1305.6314].

## 3. Infinite-\(N\) spectrum and the multi-graviton correspondence

At infinite \(N\), the weak-coupling analysis produced a complete classification of \(1/16\)-BPS operators. The basic single-trace cohomology representatives are of the schematic form
\[
\mathrm{Tr}\Big[\partial^{k_1}_z\Upsilon\;\partial^{k_2}_z\Upsilon\;\cdots\Big]_{z=\theta=0},
\]
with fermionic derivative insertions \(k_i\in\{0,1\}\) and the \(z\)-derivatives symmetrized or antisymmetrized appropriately. Since trace relations can be ignored at infinite \(N\), these single-trace classes generate the full cohomology by taking products [1305.6314].

The corresponding single-trace generating function \(Z_{\rm st}(x,a,b,u,v,w)\) is an explicit rational expression with subtraction terms that remove double counting due to derivative identities and equations of motion. Its significance is that it gives the exact single-particle counting function for \(1/16\)-BPS operators built from the cohomology representatives [1305.6314].

A central result is the equality between this gauge-theory single-trace counting and the single-particle BPS graviton spectrum in \(AdS_5\times S^5\):
\[
Z_{\rm sp} = Z_{\rm st}(x\mid z, x^2yz, x^2z/y, xz/w, xzw/v, xzv).
\]
In the “blind” unrefined limit, the resulting single-trace counting again agrees with the graviton partition function. The infinite-\(N\) conclusion is therefore sharp: the complete set of \(1/16\)-BPS operators is exactly the free multi-graviton spectrum in the bulk [1305.6314].

The finite-\(N\) extrapolation proposed in the same work was that all \(1/16\)-BPS operators remain of multi-graviton form, with trace relations reducing the number of independent states but never generating genuinely new cohomology classes. Numerical tests in low-dimensional sectors of \(SU(2)\), \(SU(3)\), and some \(SU(4)\) examples found no counterexample. The same analysis emphasized a puzzle: if no new cohomology appears at large charge, then the growth remains too small to account for large supersymmetric AdS\(_5\) black holes, since the free-graviton estimate gives \(S\sim E^{5/6}\), and even the symmetric-polynomial upper bound gives \(S\sim N^{1/3}E^{2/3}\), both parametrically below \(N^2\) black-hole entropy [1305.6314].

## 4. Finite-\(N\) index, exact small-charge theorem, and interpolation

A later finite-\(N\) analysis shifted attention from absolute cohomology to the Hamiltonian index
\[
\mathcal I_N(x)=\mathrm{Tr}_{\mathcal H_{\rm phys}}\!\left[(-1)^F\,x^n\,e^{-\beta\{\mathcal Q,\bar{\mathcal Q}\}}\right]
=\sum_n d_N(n)\,x^n,
\]
where in \(\mathcal N=1\) language
\[
n=3(2J+Q).
\]
The corresponding matrix integral is
\[
\mathcal I_N(x)=\int DU\, \exp\!\left( \sum_{j=1}^\infty \frac1j\, s(x^j)\,\mathrm{Tr}\,U^j\,\mathrm{Tr}\,(U^\dagger)^j \right),
\]
with
\[
s(x)=1-\frac{(1-x^2)^3}{(1-x^3)^2}.
\]
This provides an exact finite-\(N\) definition of the indexed degeneracies \(d_N(n)\) [2005.10843].

The multi-graviton comparison is explicit. The single-graviton index is
\[
\mathrm{grav}(x)=\frac{3x^2}{1-x^2}-\frac{2x^3}{1-x^3},
\]
and the multi-graviton index is its plethystic exponential,
\[
\mathcal I_{\text{multi-grav}}(x)=
\exp\!\left(\sum_{k=1}^\infty \frac1k\,\mathrm{grav}(x^k)\right)
=
\prod_{n=1}^\infty \frac{(1-x^{3n})^2}{(1-x^{2n})^3}.
\]
Its asymptotic growth is
\[
\log d_{\text{grav}}(n)\sim \frac{\pi}{3}\sqrt{5n}.
\]
This is the graviton-gas regime [2005.10843].

One of the main rigorous results is the threshold theorem
\[
n<2(N+1)\quad\Longrightarrow\quad d_N(n)=d_{\text{grav}}(n).
\]
The proof uses the matrix-integral expansion together with the Frobenius character formula and orthogonality of \(U(N)\) and symmetric-group characters. The first possible deviation occurs at
\[
n\ge 2(N+1),
\]
which is precisely the onset of finite-\(N\) trace relations in representation-theoretic language [2005.10843].

Numerically, for \(N=2,\dots,10\), the finite-\(N\) index displays three regimes: exact agreement with the multi-graviton index at small charge, a short crossover region, and large-\(n\) growth that approaches the black-hole entropy curve. The microscopic data also show regular “bumps” of size \(\sim N\) around the smooth black-hole curve. The overall conclusion is that the \(U(N)\) SYM index interpolates between a multi-graviton regime at small charge and a black-hole regime at large charge [2005.10843].

## 5. Cardy-like limits, black-hole entropy, and the second sheet

The large-charge regime is naturally described by Cardy-like asymptotics. In the finite-\(N\) index analysis, the fixed-\(N\) limit \(n\to\infty\) is governed by
\[
\log \mathcal I_N(\tau\to 0) =
-\frac{\pi N^2}{27}\left(\frac1{\tau^2}+\frac6\tau\right)+O(1),
\]
and the large-\(N\) saddle reproduces the supersymmetric AdS\(_5\) black-hole entropy exactly, not merely asymptotically in \(N\) [2005.10843].

A complementary viewpoint emphasizes that the index is multivalued and admits several Cardy-like limits. On the “second sheet,” shifting one of the chemical potentials by \(2\pi i\) changes the effective bose-fermi cancellations, and the index becomes the “R-charge index”
\[
\mathcal I={\rm Tr}\, e^{-\pi i R}
\,e^{-\{ \mathcal Q,{\mathcal Q}^\dagger \}
-\omega_1\left(J_1 + \frac{1}{2}R \right)- \omega_2 \left(J_2 + \frac{1}{2}R \right)}.
\]
In the resulting three-dimensional effective field theory, chiral multiplets generically do not have zero modes, the vector multiplet does, and the Kaluza–Klein tower induces a dynamical Chern–Simons level \(k_{\rm dynamical}=N\), so the effective theory is \(SU(N)_N\) and flows to a trivially gapped, confined vacuum. This yields a robust derivation of the \(\beta^{-2}\), \(\beta^{-1}\), and \(\beta^{0}\) terms in the second-sheet asymptotics [2104.01464].

A more general modular approach uses \(SL(3,\mathbb Z)\) covariance of the elliptic gamma function to generate a three-integer family of generalized Cardy limits labeled by \(\mathbf m=(m,n_1,n_2)\). The leading free energy becomes
\[
\log I_N = -i\pi\frac{N^2}{m}
\frac{[m\phi_1][m\phi_2]\big([m\phi_1]+[m\phi_2]+1\big)}
{(m\tau+1-mn_1)(m\sigma+1-mn_2)},
\]
and its Legendre transform gives
\[
S=\frac{2\pi}{m}\sqrt{
Q_1Q_2+Q_2Q_3+Q_1Q_3-\frac12 N^2(J_1+J_2)}.
\]
For \(m=1\) one recovers the standard AdS\(_5\) BPS black-hole result; for general \(m\) the entropy is reduced by \(1/m\), with an interpretation in terms of an \(m\)-divisible BPS subsector and a \(\mathbb Z_m\) quotient of the Hopf fiber [2104.07030].

On the gravity side, quantum corrections around nearly \(1/16\)-BPS AdS\(_5\) black holes lead to a partition function governed by an \(\mathcal N=2\) super-Schwarzian with \(SU(1,1|1)\) symmetry. In the supersymmetric specialization \(\alpha=1/2\), the large-\(N\) result is
\[
\mathcal I(\beta,\mathfrak j_1,\mathfrak j_2)=(-1)^{R^\star}e^{S^\star(\mathfrak j_1,\mathfrak j_2)},
\]
together with a gap
\[
\Delta_{\rm gap}=\frac{\tilde\Delta(\mathcal J_1,\mathcal J_2)}{N^2}
\]
between the \(1/16\)-BPS black holes and the lightest near-BPS states in the same charge sector. This formulation gives a gravitational explanation of why the index can reproduce the extremal entropy while the non-BPS continuum cancels in supermultiplets [2203.01331].

## 6. Reformulations and refinements

Recent work has recast the finite-\(N\) \(\mathcal N=4\) index in terms of the elliptic Ruijsenaars–Schneider integrable system and elliptic Macdonald polynomials. After rewriting the matrix integral as an elliptic hypergeometric integral with a Macdonald-type weight and kernel, the index becomes the discrete sum
\[
I_N=\chi'_N\sum_{\lambda\in\Lambda^N}u^{|\lambda|}
B_\lambda(p,q,t)\,\mathcal N_\lambda(p,q,t),
\]
where
\[
\Lambda^N=\{\lambda\in\mathbb Z^N\mid \lambda_1\ge\lambda_2\ge\cdots\ge\lambda_N\}
\]
is the set of generalized partitions, \(B_\lambda\) are structure constants from the elliptic Cauchy expansion, and \(\mathcal N_\lambda\) are elliptic norms. The elliptic deformation parameter \(p\) organizes a perturbative expansion, while the \(p\to0\) limit reproduces ordinary Macdonald data and the deformed Schur index. In the special “deformed \(1/2\)-BPS” limit \(q=t\), the zeroth-order index collapses to
\[
I_{N,0}(u)=\frac{1}{(u;u)_N}=\sum_{\ell(\lambda)\le N}u^{|\lambda|},
\]
which is exactly the \(1/2\)-BPS counting [2604.00924].

A different refinement inserts charge conjugation:
\[
\mathcal I_C=\mathrm{Tr}\!\left[(-1)^F\, C\, e^{\mu_i q_i}\right].
\]
For \(\mathcal N=4\) \(SU(N)\) SYM, this “charged superconformal index” tracks \(C\)-parity in the same protected \(1/16\)-BPS sector. Its construction relies on a “charged character” \(K_C\), conjecturally identified with an ordinary character of \(SO(N+1)\) for even \(N\) or \(SP(N-1)\) for odd \(N\), which yields a matrix-integral formula over the corresponding classical group. Like the ordinary index, the charged index is independent of \(N\) in the large-\(N\) limit, but it refines the protected spectrum by separating \(C\)-even and \(C\)-odd states [1111.1773].

Taken together, these developments show that the \(1/16\)-BPS superconformal index is not a single computational object but a family of closely related protected observables and formulations. Weak-coupling cohomology, finite-\(N\) representation theory, Cardy-like asymptotics, modular properties, integrable-system expansions, and charge-conjugation refinements all probe the same protected sector from different directions. The persistent interpretive issue is the relation between indexed growth, actual cohomology, and black-hole microstates: early weak-coupling counting emphasized a sparsity puzzle, whereas later finite-\(N\), second-sheet, and gravitational analyses exhibited black-hole entropy directly in the index. This suggests that the \(1/16\)-BPS sector is controlled not by a single approximation but by several distinct protected regimes whose relation remains structurally rich [1305.6314].

Source: https://www.emergentmind.com/topics/1-16-bps-superconformal-index