Lattice Quantum Chromodynamics for Quantum Simulations

Technical Summary: Lattice Quantum Chromodynamics for Quantum Simulations

Problem Statement
Quantum Chromodynamics (QCD) is the fundamental theory of the strong interaction, yet its non-perturbative regime—essential for understanding hadron structure and high-density matter—remains computationally challenging. Standard lattice QCD approaches rely on Euclidean path integrals and importance sampling, which fail in the presence of a "sign problem," such as in real-time dynamics or when a theta ($\theta$) angle is included. While Hamiltonian formulations offer a pathway to real-time simulation, the exponential growth of the Hilbert space with lattice volume renders classical simulation intractable for large systems. Consequently, there is a critical need for frameworks that map lattice gauge theories onto quantum computers, specifically addressing the encoding of gauge and fermionic degrees of freedom, the construction of physical (gauge-invariant) states, and the development of efficient time-evolution algorithms.

Methodology
The authors develop a comprehensive framework for quantum simulations of $SU(N_c)$ lattice gauge theory with quarks, focusing on $N_c=3$ (QCD) in two and three spatial dimensions. The methodology is built upon the following pillars:

1. Hamiltonian Formulation and Discretization:
   The work utilizes the Kogut-Susskind Hamiltonian formulation. The authors consider both staggered and Wilson fermions to address the fermion doubling problem. The theory is discretized on a cubic spatial lattice with continuous time. The Hamiltonian is decomposed into electric ($H_E$), magnetic ($H_B$), theta ($H_\theta$), kinetic ($H_{kin}$), and mass ($H_{mass}$) terms. The $\theta$ term is included as a coupling to the topological charge density, introducing complex coefficients in the Hamiltonian.

2. Representation Basis and Physical Hilbert Space:
   A central contribution is the formulation of the physical Hilbert space in a representation basis. Instead of mapping group elements directly to qubits, the gauge degrees of freedom on links and the fermionic degrees of freedom on sites are encoded using irreducible representations (irreps) of $SU(N_c)$.

   • Link Hilbert Space: Links are described by states $|(R, r_R, r_L)\rangle$, where $R$ is an irrep and $r_{R,L}$ are basis vectors for the right and left half-links. This basis diagonalizes the quadratic Casimir operators.
   • Site Hilbert Space: Fermions are treated via an occupation number basis, which is then transformed into a representation basis using Clebsch-Gordan coefficients (CGCs). This allows particle and antiparticle states to be classified by their transformation properties under $SU(N_c)$.
   • Gauge Invariance: Physical states are constructed as singlets (trivial representations) formed by the tensor product of all irreps meeting at a lattice site (links and fermions). The authors explicitly define the "site order" convention required to unambiguously construct these singlet states, ensuring the Gauss law constraint is satisfied.

3. Truncation and Encoding:
   Since the Hilbert space of $SU(N_c)$ is infinite-dimensional, the authors implement truncation schemes. These include cutoffs on the highest weight of the irreps ($n_1$) or the quadratic Casimir eigenvalue ($C_2$) for gauge links, and limits on color occupation numbers for fermions. The physical states are encoded into qubits using a local encoding scheme, where qubits represent specific irreps and singlet multiplicity indices. This approach is noted to be more qubit-efficient than direct encodings for large Hilbert spaces.

4. Quantum Simulation Algorithms:
   The authors propose a Trotterization algorithm for time evolution. The time evolution operator $e^{-iHt}$ is approximated as a product of unitaries corresponding to the decomposed Hamiltonian terms.

   • Gate Construction: The evolution of plaquette and link operators is mapped to Givens rotations (two-level unitaries) and phase shifts. The authors detail the construction of quantum circuits for these operations, including the handling of Jordan-Wigner strings (phase flips) required for fermionic anticommutation relations.
   • Matrix Elements: The paper provides explicit formulae for the matrix elements of all Hamiltonian operators in the physical basis, expressed as products of "site factors" involving CGCs. These formulae are essential for calculating the rotation angles and phases required for the quantum gates.

Key Contributions

• Framework for $SU(3)$ with $\theta$ Angle: The paper presents the first formulation and simulation of lattice QCD with a $\theta$ angle in three spatial dimensions using a Hamiltonian approach on a quantum computer.
• Representation Basis Formalism: It provides a rigorous construction of the physical Hilbert space using irreps and singlet CGCs, offering a scalable alternative to group-element bases.
• Algorithmic Details: The work details the specific quantum circuits for Trotterized evolution, including the decomposition of complex Hamiltonian terms (like the $\theta$ term and Wilson kinetic terms) into elementary gates.
• Matrix Element Calculations: The authors derive and tabulate the matrix elements for plaquette, theta, and fermion kinetic operators, providing the necessary data for implementing the simulation circuits.

Results
The authors performed noiseless simulations on small lattices (up to a $2\times2$ periodic lattice and a $2\times2\times2$ open boundary cube) using up to 32 qubits on classical simulators (Qiskit Aer and CUDA-Q). The simulations demonstrated:

• Hilbert Space Truncation Effects: Analysis of how different gauge and fermion truncations affect observables like the average electric field ($\langle \Pi^2 \rangle$) and color occupation number ($\langle N_q \rangle$). Results indicated that strong coupling suppresses high-population states, making modest truncations effective for staggered fermions, while Wilson fermions showed greater sensitivity to truncation.
• Parameter Dependence: The evolution of observables was studied as a function of fermion mass and the Wilson parameter $r$. Heavier masses were found to suppress oscillations into fermion-rich states.
• Theta Angle Dynamics: Simulations varying the $\theta$ angle revealed aperiodic behavior in the strong-coupling regime, with distinct dynamics for positive and negative $\theta$.
• Hadronic States and String Dynamics: The authors demonstrated the preparation of "charged pion" states and the evolution of string states (both loop and fermi strings). They observed string breaking dynamics, where flux tubes decay into fermion-antifermion pairs, and tracked the probability of string breaking versus unbroken configurations.
• Baryon Chemical Potential: Using exact diagonalization on a single cube, the authors computed the energy cost of adding a baryon, serving as a proof-of-principle for studying high-density matter.

Significance and Claims
The paper positions itself as a foundational contribution to the development of practical quantum simulations for fundamental physics. The authors claim that their work provides the necessary "ingredients" for designing quantum simulations of lattice QCD in the Hamiltonian formulation.

• Proof of Concept: The primary significance lies in demonstrating that the complex machinery of $SU(3)$ lattice gauge theory, including the $\theta$ angle and fermionic matter, can be encoded and evolved on quantum hardware (simulated here).
• Scalability Path: The local encoding and representation basis are presented as a viable path toward larger simulations, potentially requiring fewer qubits than global encodings.
• Future Outlook: The authors modestly state that the current results are exploratory and do not yet make serious physical statements about real-time QCD due to the small lattice sizes, strong coupling regime, and aggressive truncations used. They identify the need for improved Hamiltonians (counterterms), better state preparation algorithms, and fault-tolerant hardware to reach the continuum limit and simulate physically relevant regimes. The work is framed as a stepping stone toward simulating dense matter (e.g., neutron stars) and axion-QCD dynamics in the future.