We numerically simulate a non-Abelian lattice gauge theory in two spatial dimensions, with tensor networks (TN), up to intermediate sizes ( > 30 matter sites) well beyond exact diagonalization. We focus on the SU(2) Yang-Mills model in Hamiltonian formulation, with dynamical matter and minimally truncated gauge field (hardcore gluon). Thanks to the TN sign-problem-free approach, we characterize the phase diagram of the model at zero and finite baryon number as a function of the quark bare mass and color charge. At intermediate system sizes, we detect a liquid phase of quark-pair bound-state quasiparticles (baryons), whose mass is finite towards the continuum limit. Interesting phenomena arise at the transition boundary where color-electric and color-magnetic terms are maximally frustrated: For low quark masses, we see traces of potential deconfinement, while for high masses, signatures of a possible topological order. Published by the American Physical Society 2024
We perform a zero temperature analysis of a non-Abelian lattice gauge model corresponding to an SU(3) Yang-Mills theory in
1+1D at low energies. Specifically, we characterize the model ground states via gauge-invariant matrix product states, identifying its phase diagram at finite density as a function of the matter-gauge interaction coupling, the quark filling, and their bare mass. Overall, we observe an extreme robustness of baryons: For positive free-field energy couplings, all detected phases exhibit colorless quasiparticles, a strong numerical hint that QCD does not deconfine in 1D. Additionally, we show that, having access to finite-density properties, it is possible to study the stability of composite particles, including multibaryon bound states such as the deuteron.
We investigate the finite-density phase diagram of a non-abelian SU(2) lattice gauge theory in (1+1) -dimensions using tensor network methods. We numerically characterise the phase diagram as a function of the matter filling and of the matter-field coupling, identifying different phases, some of them appearing only at finite densities. For weak matter-field coupling we find a meson BCS liquid phase, which is confirmed by second-order analytical perturbation theory. At unit filling and for strong coupling, the system undergoes a phase transition to a charge density wave of single-site (spin-0) mesons via spontaneous chiral symmetry breaking. At finite densities, the chiral symmetry is restored almost everywhere, and the meson BCS liquid becomes a simple liquid at strong couplings, with the exception of filling two-thirds, where a charge density wave of mesons spreading over neighbouring sites appears. Finally, we identify two tri-critical points between the chiral and the two liquid phases which are compatible with a SU(2)2 Wess-Zumino-Novikov-Witten model. Here we do not perform the continuum limit but we explicitly address the global U(1) charge conservation symmetry.
This thesis develops advanced Tensor Network (TN) methods to address Hamiltonian Lattice Gauge Theories (LGTs), overcoming limitations in real-time dynamics and finite-density regimes. A novel dressed-site formalism is introduced, enabling efficient truncation of gauge fields while preserving gauge invariance for both Abelian and non-Abelian theories. This formalism is successfully applied to SU(2) Yang-Mills LGTs in two dimensions, providing the first TN simulations of this system and revealing critical aspects of its phase diagram and non-equilibrium behavior, such as a Quantum Many-Body (QMB) scarring dynamics. A generalization of the dressed-site formalism is proposed through a new fermion-to-qubit mapping for general lattice fermion theories, revealing powerful for classical and quantum simulations. Numerical innovations, including the use of optimal space-filling curves such as the Hilbert curve to preserve locality in high-dimensional simulations, further enhance the efficiency of these methods. Together with high-performance computing techniques, these advances open current and future development pathways toward optimized, efficient, and faster simulations on scales comparable to Monte Carlo state-of-the-art.
We present a quantum simulation strategy for a (1+1)-dimensional SU(2) non-Abelian lattice gauge theory with dynamical matter, a hardcore-gluon Hamiltonian Yang-Mills, tailored to a six-level trapped-ion-qudit quantum processor, as recently experimentally realized [Nat. Phys. 18, 1053 (2022)]. We employ a qudit encoding fulfilling gauge invariance, an SU(2) Gauss’s law. We discuss the experimental feasibility of generalized Mølmer-Sørensen gates used to efficiently simulate the dynamics. We illustrate how a shallow circuit with these resources is sufficient to implement scalable digital quantum simulation of the model. We also numerically show that this model, albeit simple, can dynamically manifest physically relevant properties specific to non-Abelian field theories, such as baryon excitations. Published by the American Physical Society 2024
Gauge theories are of paramount importance in our understanding of fundamental constituents of matter and their interactions. However, the complete characterization of their phase diagrams and the full understanding of non-perturbative effects are still debated, especially at finite charge density, mostly due to the sign-problem affecting Monte Carlo numerical simulations. Here, we report the Tensor Network simulation of a three dimensional lattice gauge theory in the Hamiltonian formulation including dynamical matter: Using this sign-problem-free method, we simulate the ground states of a compact Quantum Electrodynamics at zero and finite charge densities, and address fundamental questions such as the characterization of collective phases of the model, the presence of a confining phase at large gauge coupling, and the study of charge-screening effects.
We study non-stabilizerness (magic) in the ground state of
(2+1)-dimensional
SU(2) Hamiltonian lattice gauge theory with matter, formulated in the dressed-site basis in the hardcore-gluon truncation and restricted to the zero baryon-number sector. Using matrix product states, we compute three facets of magic: the second-order stabilizer Rényi entropy (SRE)
M2, its non-local component
M2NL, estimated by the SRE
M2(ψ~) of the computational-basis Schmidt representative
ψ~, and a lower bound in terms of the anti-flatness
F of the entanglement spectrum. We also prove a stronger form of the sandwich relation:
−log2(1−4F)≤M2NL≤M2; the lower bound rests on a stronger inequality that we obtain for arbitrary Schmidt bases and rank, with an optimal constant. Combined with recent results, this yields the unconditional bracket
max{−log2(1−4F),M2(ψ~)−4}≤M2NL≤M2(ψ~). We emphasize a structural distinction that makes the non-local quantities the physically preferred diagnostics: whereas the full SRE depends on the (non-unique) encoding of the gauge-invariant local Hilbert space into qubits, the non-local magic, its estimate
M2(ψ~), and the anti-flatness are invariant under site-local re-encodings and are therefore intrinsic to the state and bipartition. By varying the gauge coupling for two values of the fermion mass on lattices up to
8×8 (
∼288 qubits) with bond dimension up to
128, we find that the non-local magic, as estimated by
M2(ψ~), furnishes a sharper and more bond-dimension-friendly probe of the gauge-matter delocalization crossover compared to the full SRE or the gauge-invariant entanglement entropy, retaining a clear signal at bond dimensions well below those needed to converge the ground state itself.