Hamiltonian Large-Nc Expansion

• Hamiltonian large-Nc expansion is a systematic method that expands gauge theory Hamiltonians in inverse powers of the number of colors to simplify complex Hilbert spaces.
• It employs a leading order truncation by retaining only single-plaquette contributions while subleading terms introduce multi-plaquette correlations with controlled corrections.
• This framework enhances quantum simulations and effective field theories by reducing computational complexity and improving predictions in baryon spectroscopy and low-energy QCD.

The Hamiltonian large-$N_c$ expansion refers to the systematic expansion of Hamiltonians for quantum chromodynamics (QCD) and related gauge theories in inverse powers of the number of colors, $N_c$. This method leverages $N_c$ as a formal expansion parameter to organize contributions in the Hamiltonian and to drastically simplify the structure of gauge-invariant Hilbert spaces and interaction terms, especially in lattice regularizations. The expansion yields tractable models, often with substantial reduction in the dimensionality and complexity of the Hilbert space, while maintaining accuracy through a power-counting hierarchy in $1/N_c$. It has become a critical framework for both non-perturbative analytic studies and for efficient quantum simulation of lattice gauge theories.

1. Hamiltonian Formulation and Large-$N_c$ Power Counting

The canonical Hamiltonian for an $SU(N_c)$ lattice gauge theory in the Kogut–Susskind formulation is

$$\hat H \;=\; \sum_{l}\frac{g^2}{2}\,\hat E^2_l \;-\; \frac{1}{2g^2}\sum_{p} \left(\hat\Box_p+\hat\Box_p^\dagger\right),$$

where $g$ is the gauge coupling, $\hat E^2_l$ is the quadratic Casimir (electric energy) on each link $l$, and $\hat\Box_p$ is the oriented product of link operators around each elementary plaquette $p$ (Ciavarella et al., 2024).

Power counting in $1/N_c$ is governed by the group structure of SU($N_c$):

• Each Wilson loop wrapping $m$ plaquettes contributes a factor $N_c^{1-m}$.
• The electric term has typical matrix elements $\sim N_c$, whereas the magnetic (plaquette) term is $\mathcal{O}(1/N_c)$.
• The full Hamiltonian can thus be organized as

$$\hat H = \hat H^{(0)} + \frac{1}{N_c}\,\hat H^{(1)} + \mathcal{O}(1/N_c^2).$$

The leading-order structure (LO) isolates all processes at $N_c\rightarrow\infty$, discarding multi-plaquette loops and non-local correlations (Ciavarella et al., 14 Mar 2025, Ciavarella et al., 2024). Subleading corrections couple to longer loops and more intricate correlations.

2. Large-$N_c$ Expansion of the Hilbert Space

Physical gauge-invariant states are constructed by the action of closed loops (Wilson loops) on the vacuum: $\bigl|\{P_p,\bar P_p\}\bigr\> = \prod_{p}(\hat\Box_p)^{P_p}(\hat\Box_p^\dagger)^{\bar P_p}|0\>,$ with $P_p$, $\bar P_p\in\mathbb N$. At leading order, only single-plaquette loops (length one) are retained; multiple-plaquette loops are suppressed by $1/N_c$ per additional enclosed plaquette.

This hierarchy creates a drastic reduction in Hilbert space dimensionality. For $P$ plaquettes:

• Unconstrained, the local Hilbert space per plaquette is a qutrit (no loop, fundamental, or anti-fundamental).
• Imposing the “no-two-adjacent-excited” constraint further reduces the Hilbert space to a restricted dimension $\lesssim 2^P$ (Ciavarella et al., 2024, Ciavarella et al., 2024).
• The reduction is $N_c$-independent at leading order.

Subleading $1/N_c$ corrections systematically introduce states with two-plaquette (or longer) loops and the associated dynamical couplings.

3. Leading-Order and Subleading Hamiltonians

At $N_c\rightarrow\infty$, the Hamiltonian reduces to a form involving only single-plaquette objects: $$\hat H_\mathrm{LO} = \sum_p\left(\tfrac{8}{3}g^2 - \frac{1}{2g^2}\right)\,\hat P_{1,p} - \frac{1}{g^2\sqrt{2}}\sum_p\left(\hat P_{0,\partial p}\right)\,\hat X_p,$$ where $\hat P_{1,p}$ projects onto the charge-conjugation-even ("qubit") excitation at $p$, $\hat X_p$ is the Pauli $X$ operator, and $\hat P_{0,\partial p}$ is a projector enforcing the local constraint that surrounding plaquettes are unexcited (Ciavarella et al., 2024).

At subleading order ($O(1/N_c)$), the Hamiltonian acquires new terms allowing neighbor-coupled excitations: $$\hat H_{1/N_c} = \sum_p \left[\frac{8}{3}g^2\,\hat P_{1,p} - \frac{4}{3}g^2\,\hat P_{1,p}\sum_{\hat e}\hat P_{1,p+\hat e}\right] -\frac{1}{g^2\sqrt{2}}\sum_p \left(\hat P_{0,\partial p}\right)\hat X_p + \ldots$$ Further “conditional-flip” and multi-plaquette terms arise, suppressed by $1/g^2$ and $1/N_c$ (Ciavarella et al., 2024, Ciavarella et al., 14 Mar 2025).

4. Loop Representation, Krylov Subspaces, and Truncation Strategies

A loop or arrow representation parametrizes physical states in terms of collections of non-crossing, oriented single-plaquette loops plus internal routing data. At leading order, each plaquette is labeled by a qutrit: $$|0\rangle,\; |\circlearrowleft\rangle,\; |\circlearrowright\rangle,$$ where $|\circlearrowleft\rangle$ and $|\circlearrowright\rangle$ correspond to fundamental and anti-fundamental fluxes, respectively (Ciavarella et al., 2024).

Hilbert space truncation can be implemented via local Krylov subspaces, generated by repeated action of the plaquette operator on the vacuum, with maximal occupations and link representations constrained by the order in $1/N_c$. Truncation to $(1,1,1)$ allows only single-plaquette loops for each plaquette and link, giving qubit-per-plaquette models in the charge-conjugation-even sector. Systematic improvement is possible by incorporating higher representations and multiple occupancy, with all truncation errors entering at $\mathcal{O}(1/N_c^2)$ (Ciavarella et al., 14 Mar 2025).

5. Physical Implications: Baryons, Operators, and Model Reduction

The Hamiltonian large-$N_c$ expansion also arises in baryon spectroscopy and in the effective Hamiltonians for nucleon-nucleon and hyperon-nucleon systems. Negative-parity baryon multiplets in QCD exhibit simplified spectra at leading order, with only three degenerate “tower” eigenvalues, each highly degenerate and entirely determined by the coefficients of three dominant operators: $$O_1 = N_c\,\mathbb{1},\qquad O_2 = \boldsymbol{\ell}\cdot\boldsymbol{s},\qquad O_3 = \frac{15}{N_c} L^{(2)ij} G^{i,a} G^{j,a}$$ This linear structure reflects the compatibility of quark-shell, resonance-scattering, and core-excited frameworks in the large-$N_c$ limit (Matagne et al., 2012, Matagne et al., 2011).

For effective field theories, the $1/N_c$ expansion of the Hartree Hamiltonian reorganizes the operator basis, implying sum rules among low-energy constants and drastically reducing free parameters in chiral effective theories of baryon interactions (Liu et al., 2017). This model reduction increases the predictivity and robustness of low-energy QCD models.

6. Applications to Quantum Simulation and Computational Gains

The Hamiltonian large-$N_c$ expansion is essential to enable quantum simulation of non-Abelian gauge theories:

• At leading order, the mapping requires only qubits (or qutrits) per plaquette, with constrained local dynamics, rather than the exponentially larger Hilbert space of the full gauge theory.
• Superconducting and trapped-ion quantum devices have implemented $SU(3)$ "PXP-like" models for real-time evolution on $5\times5$ and $8\times8$ lattices, with circuit depths and qubit counts orders of magnitude smaller than brute-force encodings (Ciavarella et al., 2024).
• New Krylov-inspired truncations following this expansion have enabled 17–19 orders of magnitude reduction in quantum hardware complexity compared to previous approaches while retaining accuracy to $\mathcal O(1/N_c^2)$ (Ciavarella et al., 14 Mar 2025).

7. Extensions, Corrections, and Outlook

Corrections at $O(1/N_c)$ and beyond allow systematic restoration of multi-plaquette correlations, higher representation content, and full group-theoretic structure. For high-energy QCD (e.g., Reggeon field theory), subleading $1/N_c^2$ corrections correspond to merging vertices (e.g., $2\to1$ Reggeon fusion), which become critical in dense limits for the unitarization of scattering amplitudes (Altinoluk et al., 2014). In nucleon-nucleon potentials, the hierarchy of operators imposed by the large-$N_c$ expansion establishes the dominance of specific interactions and motivates minimal models for time-reversal and parity-violating observables (Samart et al., 2016).

The Hamiltonian large-$N_c$ expansion thus provides a systematic, physically motivated, and computationally essential framework for both analytic and quantum simulation treatments of non-Abelian gauge theories, underpinning both a deeper theoretical understanding and the realization of tractable quantum computations of strongly coupled systems.