Generalized Parton Distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED2), within the Hamiltonian formulation of lattice field theory. Quasi-distributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasi-parton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.
Parton distribution functions (PDFs) describe the inner, nonperturbative structure of hadrons. Their computation involves matrix elements with a Wilson line along a direction on the light cone, posing significant challenges in Euclidean lattice calculations, where the time direction is not directly accessible. We propose implementing the light-front Wilson line within the Hamiltonian formalism using tensor network techniques. The approach is demonstrated in the massive Schwinger model (quantum electrodynamics in 1 + 1 dimensions), a toy model that shares key features with quantum chromodynamics. We present accurate continuum results for the fermion PDF of the vector meson at varying fermion masses, obtained from first-principle calculations directly in Minkowski space. Our strategy also provides a useful path for quantum simulations and quantum computing.
We present a quantum algorithm for simulating Generalized Parton Distributions (GPDs) in the Schwinger model. Unlike the staggered fermions widely utilized in current quantum simulations, we employ Wilson fermions for lattice discretization. This choice is critical for the quantum computation of GPDs due to their strict preservation of charge conjugation symmetry. We construct a comprehensive algorithmic framework that includes the preparation of hadronic states with non-zero momentum and the measurement of light-cone correlation functions incorporating Wilson lines. We provide a complexity analysis, demonstrating that the resources required for our algorithm scale polynomially with both the number of qubits and the desired precision
ε. Finally, we benchmark our approach using exact diagonalization, extracting mass spectra and GPDs (also parton distribution functions) that are consistent with theoretical expectations and fundamental physical constraints.
A bstract Calculations of the parton distribution function (PDF) and distribution amplitude (DA) are highly relevant to core experimental programs as they provide non-perturbative inputs to inclusive and exclusive processes, respectively. Direct computation of the PDFs and DAs remains challenging because they are non-perturbative quantities defined as light-cone correlators of quark and gluon fields, and are therefore inherently time-dependent. In this work, we use a uniform quantum strategy based on tensor network simulation techniques to directly extract these hadronic quantities from first principles using the matrix product state of the hadrons in the same setup. We present exemplary numerical calculations with the Nambu-Jona-Lasinio model in 1+1 dimensions and compare with available exact diagonalization and quantum circuit simulation results. Using tensor networks, we evaluate the PDF and DA at various strong couplings in the large-qubit limit, which is consistent with expectations at perturbative and non-relativistic limits.
Parton distribution functions (PDFs) describe universal properties of bound states and allow us to calculate scattering amplitudes in processes with large momentum transfer. Calculating PDFs involves the evaluation of matrix elements with a Wilson line in a light-cone direction. In contrast to Monte Carlo methods in Euclidean spacetime, these matrix elements can be directly calculated in Minkowski-space using the Hamiltonian formalism. The necessary spatial- and time-evolution can be efficiently applied using established tensor network methods. We present PDFs in the Schwinger model calculated with matrix product states.
Hadronic tensors encode the nonperturbative structure of hadrons probed in deep inelastic scattering (DIS), yet their direct evaluation requires real-time evolution that presents a challenge for traditional Euclidean lattice approaches. In this work, we present the first quantum simulation of real-time hadronic current-current correlators in a confining gauge theory, from which DIS-inspired structure functions are extracted as a proof-of-principle demonstration in the Schwinger model, i.e (1+1)-dimensional QED. Using two complementary quantum-simulation strategies -- quantum-circuit and tensor-network methods -- we compute the real-time current-current correlator directly on the lattice and validate our results against exact diagonalization where applicable. From this correlator, we compute the hadronic tensor and determine the longitudinal structure function, the sole nonvanishing DIS observable in two space-time dimensions. Our study demonstrates that quantum simulation offers a viable complementary pathway towards the evaluation of real-time observables relevant for hadronic structure. It also provides a foundation for extending the calculations from Schwinger model to other gauge theories.
The simulation of lattice gauge theories with tensor network (TN) methods is becoming increasingly fruitful. The vision is that such methods will, eventually, be used to simulate theories in (
3+1) dimensions in regimes difficult for other methods. So far, however, TN methods have mostly simulated lattice gauge theories in (
1+1) dimensions. The aim of this paper is to explore the simulation of quantum electrodynamics (QED) on infinite lattices with TNs, i.e., fermionic matter fields coupled to a
U(1) gauge field, directly in the thermodynamic limit. With this idea in mind we first consider a gauge-invariant infinite density matrix renormalization group simulation of the Schwinger model---i.e., QED in
(1+1) d. After giving a precise description of the numerical method, we benchmark our simulations by computing the subtracted chiral condensate in the continuum, in good agreement with other approaches. Our simulations of the Schwinger model allow us to build intuition about how a simulation should proceed in (
2+1) dimensions. Based on this, we propose a variational ansatz using infinite projected entangled pair states (PEPS) to describe the ground state of
(2+1) d QED. The ansatz includes
U(1) gauge symmetry at the level of the tensors, as well as fermionic (matter) and bosonic (gauge) degrees of freedom both at the physical and virtual levels. We argue that all the necessary ingredients for the simulation of
(2+1) d QED are, a priori, already in place, paving the way for future upcoming results.