Maximal Torus Topological Entanglement Entropy in WZW Theories at Bounded Ground-State Degeneracy

Abstract: For an untwisted Wess--Zumino--Witten modular tensor category $\cC(\mathfrak g,k)$, let $r(\mathfrak g,k)$ be the number of simple objects and $\cD(\mathfrak g,k)$ its total quantum dimension. The vacuum-flux state on a torus divided into two cylinders has positive universal entropy contribution $Γ<em>{T<sup>2}=2\log\cD$. We maximize this quantity over all simple Lie algebras and positive integral levels subject to $r(\mathfrak g,k)\leq R$. If $q_R=\log_2R$, then $Γ</em>{\max}(R)\sim[7ζ(3)/(4π<sup>2)]q_R<sup>2$ as $R\to\infty$. The sequence $Sp(2n)<em>n$ attains the asymptotic coefficient. At equal rank and level, the Fourier expansion of the type-$C$ root product loses its even modes, leaving $$2π^{-2}\sum</em>{m\,\mathrm{odd}}m<sup>{-3}=7ζ(3)/(4π<sup>2)$$. An entropy--spectral inequality proves optimality among the classical families, while rank--level duality and fixed-rank estimates control unbalanced and exceptional sequences.

• The paper finds the maximal topological entanglement entropy in specified WZW the-ories to grow quadratically with binary capacity.
• The optimal scenario is achieved through balanced rank–level sequences such as Sp(2n)n under a given GSD cutoff.
• A key analytical discovery links the quadratic scaling factor to specific Fourier modes within root interactions.

Problem and setting

The paper addresses a quantitative question in the theory of $(2+1)$-dimensional topological order: given a bound $R$ on the number of anyon superselection sectors — equivalently, on the torus ground-state degeneracy (GSD) — how large can the universal topological entanglement entropy (TEE) contribution be? The optimization is restricted to untwisted Wess–Zumino–Witten (WZW) modular tensor categories $(\mathfrak g,k)$, i.e., simply connected Chern–Simons theories based on simple Lie algebras at positive integral level. The relevant observable is the vacuum-flux TEE for a torus cut into two cylinders,

$\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$

where $\mathcal D$ is the total quantum dimension; the factor of two reflects the two entanglement circles of the bipartition. The paper defines $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$ and proves

$$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$

with coefficient $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$, measured in nats per squared topological qubit. The balanced symplectic sequence $Sp(2n)_n$ attains this asymptotic coefficient.

State-dependent bounds and reduction to $\mathcal D$

For a general torus state $R$0 with flux probabilities $R$1, the universal contribution is $R$2, where $R$3. Relative entropy with respect to the quantum-dimension distribution $R$4 yields $R$5, so $R$6. The upper bound is saturated by the Kirby color state; among definite-flux states, all Abelian fluxes maximize $R$7 while any non-Abelian flux lowers it by $R$8. The vacuum flux state is therefore taken as canonical representative, reducing the optimization to maximizing $R$9 subject to $(\mathfrak g,k)$0. For Abelian categories $(\mathfrak g,k)$1 grows only linearly in binary capacity $(\mathfrak g,k)$2; quadratic growth requires non-Abelian quantum dimensions.

Half filling and the odd-mode constant

The mechanism is most transparent for type $(\mathfrak g,k)$3. Integrable highest weights are Young diagrams in an $(\mathfrak g,k)$4 rectangle, bijective to $(\mathfrak g,k)$5-element subsets of $(\mathfrak g,k)$6 ordered sites, so $(\mathfrak g,k)$7 and the filling fraction is $(\mathfrak g,k)$8; rank–level exchange acts as particle–hole exchange. Meanwhile the Kac–Peterson formula expresses $(\mathfrak g,k)$9 as a product over positive roots. With $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$0, Stirling's formula and a uniform root-height Riemann sum give

$\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$1

so $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$2 while $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$3: the quadratic scaling of TEE versus capacity originates from the $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$4 pair-root factors ($\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$5 and $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$6 channels), whereas the $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$7 long roots are subleading. A Fourier expansion converts $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$8 into a spectral form involving $\Gamma_{T^2}^{\rm vac}=2\log\mathcal D,$9, the same trilogarithm function appearing in the large-$\mathcal D$0 Chern–Simons free energy on $\mathcal D$1.

At balance $\mathcal D$2, i.e., $\mathcal D$3, the even Fourier modes vanish identically and the odd modes double:

$\mathcal D$4

Combining these estimates gives $\mathcal D$5 along $\mathcal D$6, with $\mathcal D$7. The power $\mathcal D$8 traces to the logarithmic Weyl kernel ($\mathcal D$9) plus two integrations ($$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$0).

Entropy–spectral inequality and global optimality

The proportional-limit efficiency reduces to maximizing $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$1, where $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$2. The key analytic input is the sharp inequality

$$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$3

with equality only at $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$4. Particle–hole symmetry makes $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$5 stationary; the inequality establishes uniqueness. The proof writes both sides as power series in $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$6, shows all coefficients of $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$7 beyond degree ten are strictly positive via explicit lower bounds on convolution coefficients and upper bounds on tail sums, and verifies the first ten coefficients by interval arithmetic using rational bounds on $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$8, $$\Gamma_{\max}(R)=\sup_{r(\mathfrak g,k)\le R}2\log\mathcal D$$9, and $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$0. This yields family coefficients $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$1 and $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$2: type $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$3 reaches exactly half the optimum because its root system contains only one pair channel, while types $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$4 contain both.

Boundary regimes are controlled separately. At fixed rank, $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$5 while $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$6, so the normalized ratio decays as $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$7; the exceptional families $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$8 have zero quadratic coefficient (for example $$\Gamma_{\max}(R)\sim \frac{7\zeta(3)}{4\pi^2}\,(\log_2 R)^2,\qquad R\to\infty,$$9 against $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$0). Unbalanced sequences with effective ratio tending to zero or infinity also give vanishing coefficient, using rank–level duality identities proved directly from vacuum products via Jacobi's complementary-minor theorem applied to discrete sine/cosine transforms, together with low-level orthogonal evaluations. Finiteness of the supremum at each fixed $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$1 follows from binomial sector bounds, and achievability holds for every sufficiently large cutoff since $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$2 along the balanced sequence. The same leading coefficient persists if the admissible class is enlarged to finite tensor products (stacks) of modular categories: stacking cannot improve the leading coefficient.

Limitations and open questions

The result is a statement about effective topological field theories; the paper explicitly does not address microscopic realizability of arbitrary $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$3 phases in lattice or condensed-matter systems. The optimization excludes fermionic phases described by super-modular categories, which require a choice of modular extension, as well as coset and twisted constructions and modular categories outside the WZW class. Whether comparable capacity bounds hold in those settings remains open. Subleading terms, which could distinguish the asymptotically tied $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$4 families, are not determined.

Conclusion

Within untwisted WZW modular tensor categories, the maximal vacuum torus TEE under a GSD cutoff $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$5 grows quadratically in binary capacity, $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$6, attained by balanced rank–level sequences such as $7\zeta(3)/(4\pi^2)=0.2131391994\ldots$7. The constant arises from the odd Fourier modes of the half-filled root interaction, and global optimality rests on an entropy–spectral inequality whose unique maximizer is half filling. The work provides an asymptotic benchmark for numerical studies of WZW trial states and isolates rank–level balance as the structural feature controlling universal entanglement capacity.