Quantum Simulation of Generalized Parton Distributions in the Schwinger Model

Technical Summary: Quantum Simulation of Generalized Parton Distributions in the Schwinger Model

Problem Statement
Generalized Parton Distributions (GPDs) are fundamental to understanding the internal structure of hadrons, encoding correlations between transverse impact parameters and longitudinal momentum. However, calculating GPDs from first principles in Quantum Chromodynamics (QCD) is hindered by the sign problem in Euclidean lattice QCD, which prevents the direct computation of real-time, light-cone correlation functions. While classical approaches like Large Momentum Effective Theory (LMET) offer workarounds via quasi-PDFs, they require extreme momentum boosting and suffer from power corrections. Quantum computing offers a paradigm shift by enabling direct simulation of Hamiltonian dynamics and real-time evolution, potentially bypassing the sign problem. Despite this potential, simulating GPDs remains a formidable challenge, particularly regarding the preservation of discrete symmetries (Charge Conjugation $C$, Parity $P$, Time Reversal $T$) and the efficient preparation of hadronic states with non-zero momentum.

Methodology
The authors propose a comprehensive quantum algorithm to simulate GPDs in the (1+1)-dimensional Schwinger model (quantum electrodynamics), serving as a proxy for full QCD. The methodology is distinguished by several key technical choices:

1. Wilson Fermions: Unlike current quantum simulations that predominantly utilize staggered fermions, this work employs Wilson fermions for lattice discretization. This choice is critical for the strict preservation of exact charge conjugation ($C$) symmetry on the lattice. The authors argue that $C$-symmetry is vital for studying charge-neutral mesons, as it eliminates unphysical real-part contaminations in correlation functions and ensures GPDs remain odd functions of the momentum fraction $x$.
2. Hamiltonian Formulation: The Schwinger model is formulated in the axial gauge. The continuum Hamiltonian is discretized using Wilson fermions, augmented by a Wilson term to circumvent the fermion doubling problem while preserving $C$, $P$, and $T$ symmetries. The temporal gauge field is eliminated via Gauss's law, resulting in a non-local interaction term in the Hamiltonian.
3. State Preparation: The algorithm utilizes a Quantum-Number-Resolving Variational Quantum Eigensolver (VQE). To prepare hadronic states with specific non-zero momenta ($|h(k)\rangle$) and the vacuum ($|\Omega\rangle$), the authors construct a set of orthogonal reference states. A momentum-injection circuit, $D(k)$, maps a reference state to a state with definite momentum. The VQE ansatz, $U(\theta)$, is designed to be symmetry-preserving, optimizing a multi-state loss function to simultaneously extract the vacuum and excited hadronic states within the uncharged sector ($Q=0$).
4. Correlation Function Measurement: The core of the algorithm evaluates the gauge-invariant light-cone correlation function, which includes Wilson lines. The authors implement a controlled version of the momentum translation operator and controlled gauge links ($W_0$) using an auxiliary qubit. By measuring the off-diagonal expectation values ($\sigma_x$ and $\sigma_y$) of the auxiliary qubit after real-time evolution, the spatial correlation functions $\tilde{H}_{q/h}(N, \xi)$ are extracted.
5. Complexity Analysis: The resource requirements are analyzed using the second-order Suzuki-Trotter product formula for time evolution. The total complexity scales polynomially with the system size $N$ and the inverse precision $\epsilon$, specifically as $O(N^p + N^{5/2}/\epsilon^{1/2})$.

Key Contributions

• Unified Framework: The paper establishes a unified quantum algorithmic framework capable of computing both Parton Distribution Functions (PDFs) and GPDs, rather than treating them as disjoint evaluations.
• Symmetry Preservation: It demonstrates that implementing Wilson fermions in a quantum computing framework does not incur higher computational complexity compared to staggered fermions but offers superior symmetry properties. Specifically, it rigorously preserves exact $C$-symmetry, which is shown to be essential for obtaining physically meaningful results for neutral mesons.
• Algorithmic Construction: The work details the construction of specific quantum circuits for momentum injection ($D(k)$), gauge link implementation, and the symmetry-preserving VQE ansatz required for the Schwinger model.

Results
The authors benchmarked their approach using Exact Diagonalization (ED) on a spatial lattice of size $N=18$ (and up to $N=22$ for specific checks):

• Mass Spectrum: The calculated mass spectrum of the lightest $q\bar{q}$ bound state reproduced the known analytical continuum limit ($m_h/g \approx 0.564$) via linear extrapolation, validating the Hamiltonian formulation.
• Correlation Functions: The spatial light-cone correlation functions $\tilde{H}_{q/h}(N, \xi)$ were extracted. The results confirmed that the real part of the correlation function vanishes asymptotically as the lattice volume increases, consistent with $C$-symmetry.
• Comparison with Staggered Fermions: A direct comparison showed that the ratio of the real to imaginary parts of the correlation function is significantly suppressed for Wilson fermions compared to staggered fermions. This confirms that the Wilson formulation effectively eliminates unphysical contaminations arising from finite-volume effects and symmetry violations.
• GPDs and PDFs: Fourier transforms of the correlation functions yielded momentum-space GPDs. The resulting PDFs (limit $\xi=0$) peaked at $x=0.5$, correctly reflecting the two-valence-quark nature of the bound state in 1+1 dimensions. The GPDs were found to be odd functions of $x$, satisfying the theoretical requirement that the elastic charge form factor for a charge-neutral state must vanish.

Significance and Claims
The paper claims that this work provides a vital precursor to full (3+1)-dimensional QCD simulations. By demonstrating that Wilson fermions can be efficiently implemented on quantum hardware without increasing complexity, the authors argue that this approach offers a robust path toward simulating non-perturbative phenomena while strictly adhering to fundamental symmetries. The polynomial scaling of the algorithm with system size and precision highlights the theoretical advantage of quantum computers over classical methods for directly simulating light-cone correlations. The authors position this work as a foundational step toward generalizing these algorithms to non-Abelian gauge theories (e.g., SU(2) or SU(3) Yang-Mills) and, ultimately, to full QCD on future fault-tolerant quantum hardware.