Exponential speedup in quantum simulation of Kogut-Susskind Hamiltonian via orbifold lattice

Technical Summary: Exponential Speedup in Quantum Simulation of Kogut-Susskind Hamiltonian via Orbifold Lattice

Problem Statement
The Kogut-Susskind (KS) Hamiltonian has long served as the canonical framework for simulating Yang-Mills theory and Quantum Chromodynamics (QCD) on quantum computers. However, the authors identify a fundamental incongruity: the KS formalism, developed prior to the advent of quantum computation, relies on compact unitary link variables. This reliance creates a Hilbert space with structural complexity that becomes intractable for $SU(N \geq 2)$ theories in more than one spatial dimension. Existing protocols face two primary constraints:

1. Exponential Scaling: Classical preprocessing requirements and quantum circuit depths scale exponentially with the number of qubits per link variable.
2. Implementation Barriers: There is no known quantum simulation protocol that allows for the continuum limit of the KS Hamiltonian without exponential resource growth. Even asymptotic gate counting is difficult, and explicit circuit constructions for $SU(3)$ in $3+1$ dimensions remain elusive or resource-prohibitive.

The core issue is the use of compact variables. The authors posit that theories based on non-compact variables (such as scalar field theories) are straightforward to simulate, whereas the compact nature of the KS Hamiltonian introduces unnecessary complexity.

Methodology
The paper proposes a methodological synthesis: deriving the Kogut-Susskind Hamiltonian as a specific limiting case of the orbifold lattice Hamiltonian. The orbifold lattice formulation utilizes non-compact complex link variables ($Z_{j,\vec{n}}$) rather than compact unitary ones.

The methodology proceeds as follows:

1. Theoretical Framework: The orbifold lattice Hamiltonian is defined using complex matrices $Z_{j,\vec{n}}$, which decompose into a Hermitian matrix $W_{j,\vec{n}}$ and a unitary link variable $U_{j,\vec{n}}$. The Hamiltonian includes mass terms ($m^2$) and a linear term ($\gamma$) designed to force $W_{j,\vec{n}}$ toward the identity matrix and the determinant of $U_{j,\vec{n}}$ toward 1.
2. The Infinite Mass Limit: The authors demonstrate that in the limit where the scalar mass parameters ($m^2$ and $m^2_{U(1)}$) approach infinity, the scalar fields and the $U(1)$ part of the gauge field decouple. In this limit, the orbifold lattice action reduces exactly to the standard Wilson action (and consequently the KS Hamiltonian) at finite lattice spacing.
3. Numerical Verification: To validate this equivalence at finite lattice spacing and strong coupling (beyond the continuum limit), the authors employ Euclidean path integral methods and Hybrid Monte Carlo (HMC) simulations. They simulate $(2+1)$-dimensional pure Yang-Mills theory for $SU(2)$ and $SU(3)$ gauge groups.
4. Extrapolation Strategy: Instead of simulating the KS Hamiltonian directly, the protocol involves:
   • Simulating the orbifold lattice Hamiltonian at various finite mass values.
   • Systematically extrapolating the results to the infinite-mass limit.
   • Comparing observables (spatial/temporal plaquettes, Polyakov loops) against standard Wilson action results to confirm convergence.

Key Contributions

• Analytical Proof of Equivalence: The paper provides a rigorous analytical treatment showing that the orbifold lattice action converges to the Wilson action (and thus the KS Hamiltonian) as the scalar mass goes to infinity, even at finite lattice spacing. This establishes the KS Hamiltonian as a controlled limit of the orbifold formulation.
• Resolution of Technical Obstacles: By utilizing the orbifold lattice Hamiltonian, the authors bypass the exponential scaling associated with compact variables. The orbifold Hamiltonian belongs to a class of bosonic systems with polynomial potential terms, allowing for explicit quantum circuit construction using only CNOT gates, single-qubit rotations, and Quantum Fourier Transforms. This leads to a polynomial scaling of gate complexity with respect to the number of qubits per degree of freedom.
• Numerical Evidence: The authors present the first numerical evidence confirming that the KS limit can be smoothly approached from the orbifold lattice. They demonstrate convergence of plaquette values and Polyakov loops for $SU(2)$ and $SU(3)$ in $2+1$ dimensions, showing that the orbifold results match the Wilson action results as $m^2 \to \infty$.
• Handling of Linear Counterterms: The study investigates the role of the linear term $\hat{H}_{linear}$ (parameterized by $\gamma$). They show that tuning $\gamma$ to a critical value $\gamma_c$ (where $\langle \text{Tr}(W-1) \rangle = 0$) significantly improves convergence, effectively removing lattice spacing shifts and allowing the KS limit to be evaluated using mass scales comparable to the inverse lattice spacing ($m \sim a^{-1}$).

Results

• Convergence of Observables: Numerical simulations for $SU(2)$ and $SU(3)$ on $8^3$ and $4 \times 16^2$ lattices show that spatial and temporal plaquettes, as well as the Polyakov loop, converge to the values predicted by the Wilson action as $1/m^2 \to 0$.
• Decoupling of Scalars: The data confirms that as the mass increases, the scalar fields ($W$) and the $U(1)$ components decouple from the $SU(N)$ gauge dynamics, leaving only the pure gauge theory.
• Extrapolation Validity: The paper demonstrates that the KS limit can be accurately recovered by extrapolating data from simulations where the dimensionless product $a \times m$ is of order $O(1)$. This implies that no new energy scales beyond the cutoff are required to estimate the KS limit.
• Comparison of Limits: The results confirm that the orbifold lattice formulation yields the same continuum limit as the KS formulation, validating the orbifold approach as a viable alternative for quantum simulation.

Significance and Claims
The paper claims to provide the first viable method to make the Kogut-Susskind Hamiltonian programmable on digital quantum computers for arbitrary gauge groups and dimensions. The significance lies in the "methodological synthesis":

• Exponential Speedup: By simulating the orbifold lattice Hamiltonian (which admits efficient quantum circuits) and extrapolating to the infinite-mass limit, one can effectively simulate the KS Hamiltonian with an exponential speedup compared to direct KS simulation methods.
• Preservation of Theoretical Benefits: This approach retains the direct correspondence with Wilson's Euclidean path-integral formulation (a standard framework in lattice QCD) while resolving the technical implementation disadvantages of the KS Hamiltonian on quantum hardware.
• Universality: The framework is applicable to arbitrary $SU(N)$ groups and spatial dimensions, offering a path forward for simulating QCD and Yang-Mills theories that was previously blocked by resource constraints.

The authors emphasize that this is not a new physical theory but a reformulation that leverages the simplicity of non-compact variables to achieve computational efficiency, effectively transferring the advantages of the orbifold lattice to the standard KS framework.