We introduce a class of variational states to study ground-state properties and real-time dynamics in (2 + 1)-dimensional compact QED.These are based on complex Gaussian states which are made periodic to account for the compact nature of the U (1) gauge field.Since the evaluation of expectation values involves infinite sums, we present an approximation scheme for the whole variational manifold.We calculate the ground-state energy density for lattice sizes up to 20 × 20 and extrapolate to the thermodynamic limit for the whole coupling region.Additionally, we study the string tension both by fitting the potential between two static charges and by fitting the exponential decay of spatial Wilson loops.As the ansatz does not require a truncation in the local Hilbert spaces, we analyze truncation effects which are present in other approaches.The variational states are benchmarked against exact solutions known for the one plaquette case and exact diagonalization results for a Z 3 lattice gauge theory.Using the time-dependent variational principle, we study real-time dynamics after various global quenches, e.g., the time evolution of a strongly confined electric field between two charges after a quench to the weak-coupling regime.Up to the points where finite-size effects start to play a role, we observe equilibrating behavior.
We introduce a class of variational states to study ground-state properties and real-time dynamics in (2 + 1)-dimensional compact QED.These are based on complex Gaussian states which are made periodic to account for the compact nature of the U (1) gauge field.Since the evaluation of expectation values involves infinite sums, we present an approximation scheme for the whole variational manifold.We calculate the ground-state energy density for lattice sizes up to 20 × 20 and extrapolate to the thermodynamic limit for the whole coupling region.Additionally, we study the string tension both by fitting the potential between two static charges and by fitting the exponential decay of spatial Wilson loops.As the ansatz does not require a truncation in the local Hilbert spaces, we analyze truncation effects which are present in other approaches.The variational states are benchmarked against exact solutions known for the one plaquette case and exact diagonalization results for a Z 3 lattice gauge theory.Using the time-dependent variational principle, we study real-time dynamics after various global quenches, e.g., the time evolution of a strongly confined electric field between two charges after a quench to the weak-coupling regime.Up to the points where finite-size effects start to play a role, we observe equilibrating behavior.
Lattice gauge theories coupled to fermionic matter account for many interesting phenomena in both high-energy physics and condensed-matter physics. Certain regimes, e.g., at finite fermion density, are difficult to simulate with traditional Monte Carlo algorithms due to the so-called sign problem. We present a variational, sign-problem-free Monte Carlo method for lattice gauge theories with continuous gauge groups and apply it to (2+1)-dimensional compact QED with dynamical fermions at finite density. The variational ansatz is formulated in the full gauge-field basis, i.e., without having to resort to truncation schemes for the U(1) gauge-field Hilbert space. The ansatz consists of two parts: first, a pure gauge part based on Jastrow-type ansatz states (which can be connected to certain neural-network ansatz states) and, second, a fermionic part based on gauge-field-dependent fermionic Gaussian states. These are designed in such a way that the gauge-field integral over all fermionic Gaussian states is gauge-invariant and at the same time still efficiently tractable. To ensure the validity of the method we benchmark the pure gauge part of the ansatz against another variational method and the full ansatz against an existing Monte Carlo simulation where the sign problem is absent. Moreover, in limiting cases where the exact ground state is known we show that our ansatz is able to capture this behavior. Finally, we study a sign-problem affected regime by probing density-induced phase transitions.
This thesis investigates advances in quantum simulation methods for lattice gauge theories, with a focus on (2+1)-dimensional quantum electrodynamics (QED), which shares important features with (3+1)-dimensional quantum chromodynamics (QCD), such as confinement and asymptotic freedom. Two main projects are developed. The first project investigates the ongoing coupling in a compact pure U(1) gauge theory using a step-scaling function approach. A quantum variational method allows calculations over a wide range of couplings. Plaquette expectation value results are compared with classical Monte Carlo simulations to obtain, in principle, the physical value of the lattice spacing. This method can potentially be extended to the study of QED, which also includes fermionic matter. The second project investigates the static potential between two charges, analyzing the Coulomb, confining, and string-breaking regimes. It involves the development of a tailored variational approach that incorporates system symmetries and mutual information to represent fermionic degrees of freedom. Results from a trapped-ion quantum device agree well with classical simulations, validating the method. Furthermore, this work evaluates superconducting quantum devices, noting issues with noise but also demonstrating improvements in hardware and potential through error limitation. Overall, this work underscores the feasibility of variational quantum algorithms for lattice gauge theory simulations and demonstrates their promise for advancing theoretical physics with quantum computers. This thesis explores advancements in quantum simulation methods for lattice gauge theories, focusing on (2+1)-dimensional Quantum Electrodynamics (QED), which shares key features with (3+1)-dimensional Quantum Chromodynamics (QCD), such as confinement and asymptotic freedom. Two main projects are developed. The first investigates the running coupling in a compact U(1) pure gauge theory using a step scaling function approach. A quantum variational method enables computation across a wide range of couplings. The plaquette expectation value is matched with classical Monte Carlo simulations to obtain, in principle, the physical value of the lattice spacing. This method can be eventually extended to study QED, where fermionic matter is included. The second project examines the static potential between two charges, analyzing Coulomb, confining, and string-breaking regimes. It involves developing a custom variational Ansatz that incorporates system symmetries and mutual information to represent fermionic degrees of freedom. Results from a trapped-ion quantum device align well with exact classical simulations, validating the method. Additionally, the thesis evaluates superconducting devices, identifying noise challenges but noting improvements in hardware and the potential of error mitigation. Overall, this work highlights the viability of variational quantum algorithms for lattice gauge theory simulations and demonstrates their promise for advancing theoretical physics with quantum computing.
Lattice gauge theories (LGTs)1–4 can be used to understand a wide range of phenomena, from elementary particle scattering in high-energy physics to effective descriptions of many-body interactions in materials5–7. Studying dynamical properties of emergent phases can be challenging, as it requires solving many-body problems that are generally beyond perturbative limits8–10. Here we investigate the dynamics of local excitations in a Z2 LGT using a two-dimensional lattice of superconducting qubits. We first construct a simple variational circuit that prepares low-energy states that have a large overlap with the ground state; then we create charge excitations with local gates and simulate their quantum dynamics by means of a discretized time evolution. As the electric field coupling constant is increased, our measurements show signatures of transitioning from deconfined to confined dynamics. For confined excitations, the electric field induces a tension in the string connecting them. Our method allows us to experimentally image string dynamics in a (2+1)D LGT, from which we uncover two distinct regimes inside the confining phase: for weak confinement, the string fluctuates strongly in the transverse direction, whereas for strong confinement, transverse fluctuations are effectively frozen11,12. We also demonstrate a resonance condition at which dynamical string breaking is facilitated. Our LGT implementation on a quantum processor presents a new set of techniques for investigating emergent excitations and string dynamics. In a quantum simulation of a (2+1)D lattice gauge theory using a superconducting quantum processor, the dynamics of strings reveal the transition from deconfined to confined excitations as the effective electric field is increased.
Real-Time String Dynamics in 3+1D Lattice Quantum Electrodynamics
Understanding real-time string dynamics in three spatial dimensions is essential for connecting quantum simulations of lattice gauge theories (LGTs) to the physical dimensionality of QED and QCD, where transverse fluctuations and competing local processes proliferate. We present the first real-time simulations of string breaking in 3+1D lattice quantum electrodynamics. Using tree tensor networks, we simulate the quench dynamics of electric flux strings in a 3+1D U(1) LGT with dynamical matter. At strong coupling, the string breaks resonantly at a sharp resonance condition of mass and gauge couplings, converting electric energy into matter--antimatter pairs that screen the static charges. Off resonance, we classify all competing channels---pair production, string deformations and extensions, and flux loops---whose multiplicity, extensive for pair production and flux loops, depletes the string sector even far from resonance. A channel-resolved perturbation theory quantitatively reproduces these dynamics and their Fourier spectrum. Our results establish diagnostics and benchmarks for upcoming quantum simulators of higher-dimensional LGTs.
Compact lattice Quantum Electrodynamics is a complex quantum field theory with dynamical gauge and matter fields and it has similarities with Quantum Chromodynamics, in particular asymptotic freedom and confinement. We consider a (2+1)-dimensional lattice discretization of Quantum Electrodynamics with the inclusion of dynamical fermionic matter. We define a suitable quantum algorithm to measure the static potential as a function of the distance between two charges on the lattice and we use a variational quantum calculation to explore the Coulomb, confinement and string breaking regimes. A symmetry-preserving and resource-efficient variational quantum circuit is employed to prepare the ground state of the theory at various values of the coupling constant, corresponding to different physical distances, allowing the accurate extraction of the static potential from a quantum computer. We demonstrate that results from quantum experiments on the Quantinuum H1-1 trapped-ion device and emulator, with full connectivity between qubits, agree with classical noiseless simulations using circuits with 10 and 24 qubits. Moreover, we visualize the electric field flux configurations that mostly contribute in the wave-function of the quantum ground state in the different regimes of the potential, thus giving insights into the mechanisms of confinement and string breaking. These results are a promising step forward in the grand challenge of solving higher dimensional lattice gauge theory problems with quantum computing algorithms.
Gauge theories describe fundamental interactions in the standard model of particle physics and build effective theories in condensed matter physics. Being notoriously hard to simulate, these theories are commonly regularized as lattice gauge theories which can be evaluated with high precision through Monte Carlo algorithms or variational methods. However, Monte Carlo algorithms are not applicable in all regimes due to the sign problem; while variational methods, quantum and classical alike, only give upper-bounds on the ground-state energy. Their precision crucially depends on the chosen ansatz. Here, we present a framework based on a hierarchy of semidefinite programs to yield increasingly good lower-bounds on the ground-state energy and propagate these bounds to extended observables like Wilson loops and mesonic strings. We demonstrate the numerical capabilities by applying the algorithm to a (1+1)-dimensional Z2 theory with dynamic fermionic matter and its pure-gauge version in (2+1)-dimensions. In both cases, we obtain certified intervals for the ground-state energy. For a one-dimensional system with 129 sites, the certified interval has a relative spread of 0.04%. For a two-dimensional system with 97 sites (corresponding to a 7×8 lattice), the relative spread is 2.08%. Beyond energy estimates, the proposed method provides rigorous bounds on both short- and long-range observables such as Wilson loops over multiple plaquettes and the magnitude of the mesonic string. The method can be both used as a computational microscope into many-body physics and a way to certify results of quantum simulators.