Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2609.35722

Multilevel Plaquette-Space Sampling for Lattice Gauge Theories with Local Constraint Solves

Ankur Singha, Jacob Finkenrath, Karl Jansen, Vipul Arora, Shinichi Nakajima

hep-lathep-thmath-phquant-phSubmitted Sep 28, 2026• Announced Sep 29, 2026
arXivHTMLPDF

Abstract

Lattice gauge theories are an important class of physical models with high-dimensional structured distributions that underpin first-principles calculations in particle, nuclear, and condensed-matter physics. Recent advances in generative modeling have opened new avenues for sampling Boltzmann distributions in lattice gauge theories, offering the potential to alleviate limitations of traditional Monte Carlo methods, including critical slowing down and topological freezing. However, generative samplers for lattice gauge theories are typically constructed in link space, where correlations become increasingly long-ranged toward weak coupling, posing a challenge for learning. Plaquettes provide a more natural representation, as the action is local in these variables. However, exact Bianchi constraints restrict them to a lower-dimensional manifold, complicating direct generative modeling. We introduce a multilevel normalizing-flow construction that samples directly in plaquette space while satisfying these constraints exactly. The key idea is a coarse-to fine factorization that transforms a globally coupled constraint problem into a sequence of local solves: at each refinement, every determined plaquette depends on at most four newly generated variables, while the size of the constraint problem remains independent of the lattice size. We validate the construction for $U(1)$ in two and four dimensions and for $SU(2)$ in two dimensions. Our multilevel Plaquette-Space Sampler (PSS) substantially outperforms link-space baselines, with the advantage increasing toward weak coupling, where the gauge coupling becomes small. This regime is particularly challenging for generative sampling and, in asymptotically free gauge theories, is relevant to continuum studies.

Comments
9 pages main text, 38 pages total, 10 figures, 9 tables

Related papers

Matter-induced plaquette terms in a $\mathbb{Z}_2$ lattice gauge theory

Matjaž Kebrič, Fabian Döschl, Umberto Borla +4 more

Lattice gauge theories (LGTs) provide a powerful framework for studying confinement, topological order, and exotic quantum matter. In particular, the paradigmatic phenomenon of confinement, where dynamical matter is coupled to gauge fields and forms bound states, remains an open problem. In addition, LGTs can provide low-energy descriptions of quantum spin liquids, which is the focus of ongoing experimental research. However, the study of LGTs is often limited theoretically by their numerical complexity and experimentally in implementing challenging multi-body interactions, such as the plaquette terms crucial for the realization of many exotic phases of matter. Here we investigate a $(2+1)$D $\mathbb{Z}_2$ LGT coupled to hard-core bosonic matter featuring a global U(1) symmetry, and show that dynamical matter naturally induces sizable plaquette interactions even in the absence of explicit plaquette terms in the Hamiltonian. Using a combination of density matrix renormalization group simulations and neural quantum state calculations up to a system size of $20 \times 20$, we analyze the model across different fillings and electric field strengths. At small coupling strength, we find a large plaquette expectation value, independent of system size, for a wide range of fillings, which decreases in the presence of stronger electric fields. Furthermore, we observe signatures of a confinement-deconfinement transition at weak coupling strengths. Our results demonstrate that dynamical U(1) matter can induce complex multi-body interactions, suggesting a natural route to the realization of strong plaquette terms and paving the way for realizing a topological quantum spin liquid protected by a large gap.

cond-mat.quant-gasCross-listedarXiv:2602.13192•Feb 16, 2026•0 citations

Protecting Quantum Simulations of Lattice Gauge Theories through Engineered Emergent Hierarchical Symmetries

Zhanpeng Fu, Wei Zheng, Roderich Moessner +2 more

We present a strategy for the quantum simulation of many-body lattice models with constrained Hilbert spaces. We focus on lattice gauge theories (LGTs), which underlie a wide range of phenomena in particle physics, condensed matter, and quantum information. In present-day quantum computing platforms, perfect restrictions of the Hilbert space to the desired gauge sectors are beyond reach: for LGTs, violations of the local constraint are unavoidable, posing a formidable challenge for the emulation of the underlying physics. Here, we develop a Floquet-engineering framework that restructures departures from a target sector such that a series of emergent local symmetries occurs hierarchically in time and in a controllable way. This leads to a set of approximate dynamical selection rules that strongly restrict inter-sector couplings, resulting in a pronounced, symmetry-controlled hierarchy of lifetimes for the state population to spread among sectors. Concretely, this protects $U(1)$ LGTs against violations of the {defining} local symmetry. While some sectors remain very long-lived, others are destabilized on shorter timescales. We numerically verify our theory for the one-dimensional $U(1)$ quantum link model. In addition, we reveal that `defects', whose movement accounts for violations of the gauge constraint, are kinetically constrained, becoming mobile only through the assistance of intra-sector dynamics, which we describe using an effective quantum marble model. Our results can thus be used to extend the lifetime, in the spirit of passive error correction, of quantum simulations of complex many-body problems when emergent or desired local symmetries are only implemented approximately.

quant-pharXiv:2604.11085•Apr 14, 2026•0 citations

Deforming the Trail: Baseline Quantum Circuitry for $\text{SU(2)}_k$ Lattice Gauge Theory

Zoë Webb-Mack, Natalie Klco

Quantifying quantum resources for simulating the fundamental forces of Nature is sensitive to the mapping of gauge fields onto finite quantum computational architectures. When locally truncating lattice gauge theories in the irreducible representation basis, it has been proposed to further deform the theory via quantum groups. The purpose of this deformation is (1) to provide an infinite tower of finite-dimensional ($d = k+1$) groups systematically approximating the infinite-dimensional gauge links and (2) to restore the physical unitarity of a plaquette operator diagonalization procedure analytically derived from the field continuum by recontracting vertex pairs. For the SU(2)$_k$ Yang-Mills pure-gauge theory, we provide a constructive strategy of gauge-variant completions to extend this unitarity to the entire computational Hilbert space, leading to well-defined time evolution unitaries as targets for optimized circuit synthesis. Leveraging basic circuit decompositions and symmetries of the diagonalized plaquette operator, we report resource upper-bounds on the generalized-controlled-X two-qudit gates for arbitrary local truncation $d$, reducing estimates and scaling relative to the non-deformed theory by three polynomial powers from $O(d^8)$ to $O(d^5)$. Examining the stronger q-deformed gauge constraint, which softens the total flux at vertices, we show that the physical Hilbert space dimension of the deformed plaquette operator scales equivalently to its non-deformed counterpart with a constant factor $0.2563(5)$. Thus, despite affecting interactions at all scales as exemplified by the observed flux hierarchy inversion symmetry, q-deformation continues to pass scrutiny as a reliable truncation offering advantages in quantum circuit synthesis.

quant-pharXiv:2605.15076•May 15, 2026•0 citations

A resource efficient approach for quantum and classical simulations of gauge theories in particle physics

Jan F. Haase, Luca Dellantonio, Alessio Celi +4 more

Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use $2+1$ dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.

•Jan 1, 2021•116 citations

A resource efficient approach for quantum and classical simulations of gauge theories in particle physics

Jan F. Haase, Luca Dellantonio, Alessio Celi +4 more

Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2 + 1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.

arXiv:2006.14160•Quantum•Feb 4, 2021•129 citations

Flow-Based Generative Sampling for the Auxiliary-Field Hubbard Model

Dominic Schuh

Understanding strongly correlated electron systems remains one of the central challenges in condensed matter physics. Finding ab initio descriptions of non-perturbative systems is of great importance for modeling carbon-based materials such as graphene or fullerenes. In this regime, computational methods become indispensable tools for performing simulations, making predictions, and ultimately comparing to experiments. The Hubbard model offers a minimal description of the dynamics of electrons in such materials. As exact methods are limited to small system sizes by the exponentially growing Hilbert space, larger-scale first-principle studies rely on numerical approaches. Alongside powerful Hamiltonian-based methods, the auxiliary-field description of the model offers a natural setting for Lagrangian-based Monte Carlo sampling. However, substantial challenges, such as lack of ergodicity and the numerical sign problem, arise in these approaches, leading to long autocorrelation times, biased results, and large statistical uncertainties. In recent years, generative models, especially normalizing flows, have been applied extensively to various theories to address shortcomings of traditional Monte Carlo-based algorithms, such as critical slowing down or topological freezing. While these approaches have shown promise in alleviating these issues, a poorly understood and unfavorable scaling of these methods has stalled their progress, and whether they can reach state-of-the-art system sizes in lattice quantum chromodynamics remains to be seen. This motivates the application of these methods to condensed matter systems, where physically relevant system sizes are reachable at substantially smaller scales. In this work, we study for the first time the application of generative models to the finite-temperature auxiliary-field Hubbard model. We begin by employing an equivariant normalizing flow that respects multiple symmetries inherent to the Boltzmann distribution, demonstrating a substantial improvement in sampling efficiency over Hybrid Monte Carlo and standard non-equivariant flows. We further develop Symmetry-Enforcing Stochastic Modulation (SESaMo), a more scalable algorithm capable of modeling both exact and broken symmetries in normalizing flows, and show that it achieves significantly higher effective sample sizes than traditional equivariant techniques at large system sizes. Equipped with SESaMo, we compute physical observables at half-filling in excellent agreement with exact diagonalization at substantially lower computational cost than Hybrid Monte Carlo. We then turn to the doped Hubbard model, where the theory exhibits a numerical sign problem. We introduce an annealing scheme to simulate the Hubbard model at finite chemical potential in the spin basis. Making use of the less severe sign problem in this description, the annealing scheme enables ergodic sampling of the multi-modal distribution, yielding more accurate results with substantially reduced uncertainties compared to state-of-the-art optimized Hybrid Monte Carlo. Lastly, we employ a masked autoencoder for distribution estimation to learn the Boltzmann distribution of the discrete auxiliary-field Hubbard model. Preliminary results show a significantly higher effective sample size compared to normalizing flows with continuous degrees of freedom, while also exhibiting a milder sign problem than the continuous spin basis.

•Aug 31, 2026•0 citations
Lattice gauge theories are an important class of physical models with high-dimensional structured distributions that underpin first-principles calculations in particle, nuclear, and condensed-matter physics. Recent advances in generative modeling have opened new avenues for sampling Boltzmann distributions in lattice gauge theories, offering the potential to alleviate limitations of traditional Monte Carlo methods, including critical slowing down and topological freezing. However, generative samplers for lattice gauge theories are typically constructed in link space, where correlations become increasingly long-ranged toward weak coupling, posing a challenge for learning. Plaquettes provide a more natural representation, as the action is local in these variables. However, exact Bianchi constraints restrict them to a lower-dimensional manifold, complicating direct generative modeling. We introduce a multilevel normalizing-flow construction that samples directly in plaquette space while satisfying these constraints exactly. The key idea is a coarse-to fine factorization that transforms a globally coupled constraint problem into a sequence of local solves: at each refinement, every determined plaquette depends on at most four newly generated variables, while the size of the constraint problem remains independent of the lattice size. We validate the construction for U(1)U(1)U(1) in two and four dimensions and for SU(2)SU(2)SU(2) in two dimensions. Our multilevel Plaquette-Space Sampler (PSS) substantially outperforms link-space baselines, with the advantage increasing toward weak coupling, where the gauge coupling becomes small. This regime is particularly challenging for generative sampling and, in asymptotically free gauge theories, is relevant to continuum studies.
Matter-induced plaquette terms in a Z2\mathbb{Z}_2Z2​ lattice gauge theory
Lattice gauge theories (LGTs) provide a powerful framework for studying confinement, topological order, and exotic quantum matter. In particular, the paradigmatic phenomenon of confinement, where dynamical matter is coupled to gauge fields and forms bound states, remains an open problem. In addition, LGTs can provide low-energy descriptions of quantum spin liquids, which is the focus of ongoing experimental research. However, the study of LGTs is often limited theoretically by their numerical complexity and experimentally in implementing challenging multi-body interactions, such as the plaquette terms crucial for the realization of many exotic phases of matter. Here we investigate a (2+1)(2+1)(2+1)D Z2\mathbb{Z}_2Z2​ LGT coupled to hard-core bosonic matter featuring a global U(1) symmetry, and show that dynamical matter naturally induces sizable plaquette interactions even in the absence of explicit plaquette terms in the Hamiltonian. Using a combination of density matrix renormalization group simulations and neural quantum state calculations up to a system size of 20×2020 \times 2020×20, we analyze the model across different fillings and electric field strengths. At small coupling strength, we find a large plaquette expectation value, independent of system size, for a wide range of fillings, which decreases in the presence of stronger electric fields. Furthermore, we observe signatures of a confinement-deconfinement transition at weak coupling strengths. Our results demonstrate that dynamical U(1) matter can induce complex multi-body interactions, suggesting a natural route to the realization of strong plaquette terms and paving the way for realizing a topological quantum spin liquid protected by a large gap.
We present a strategy for the quantum simulation of many-body lattice models with constrained Hilbert spaces. We focus on lattice gauge theories (LGTs), which underlie a wide range of phenomena in particle physics, condensed matter, and quantum information. In present-day quantum computing platforms, perfect restrictions of the Hilbert space to the desired gauge sectors are beyond reach: for LGTs, violations of the local constraint are unavoidable, posing a formidable challenge for the emulation of the underlying physics. Here, we develop a Floquet-engineering framework that restructures departures from a target sector such that a series of emergent local symmetries occurs hierarchically in time and in a controllable way. This leads to a set of approximate dynamical selection rules that strongly restrict inter-sector couplings, resulting in a pronounced, symmetry-controlled hierarchy of lifetimes for the state population to spread among sectors. Concretely, this protects U(1)U(1)U(1) LGTs against violations of the {defining} local symmetry. While some sectors remain very long-lived, others are destabilized on shorter timescales. We numerically verify our theory for the one-dimensional U(1)U(1)U(1) quantum link model. In addition, we reveal that `defects', whose movement accounts for violations of the gauge constraint, are kinetically constrained, becoming mobile only through the assistance of intra-sector dynamics, which we describe using an effective quantum marble model. Our results can thus be used to extend the lifetime, in the spirit of passive error correction, of quantum simulations of complex many-body problems when emergent or desired local symmetries are only implemented approximately.
Deforming the Trail: Baseline Quantum Circuitry for SU(2)k\text{SU(2)}_kSU(2)k​ Lattice Gauge Theory
Quantifying quantum resources for simulating the fundamental forces of Nature is sensitive to the mapping of gauge fields onto finite quantum computational architectures. When locally truncating lattice gauge theories in the irreducible representation basis, it has been proposed to further deform the theory via quantum groups. The purpose of this deformation is (1) to provide an infinite tower of finite-dimensional (d=k+1d = k+1d=k+1) groups systematically approximating the infinite-dimensional gauge links and (2) to restore the physical unitarity of a plaquette operator diagonalization procedure analytically derived from the field continuum by recontracting vertex pairs. For the SU(2)k_kk​ Yang-Mills pure-gauge theory, we provide a constructive strategy of gauge-variant completions to extend this unitarity to the entire computational Hilbert space, leading to well-defined time evolution unitaries as targets for optimized circuit synthesis. Leveraging basic circuit decompositions and symmetries of the diagonalized plaquette operator, we report resource upper-bounds on the generalized-controlled-X two-qudit gates for arbitrary local truncation ddd, reducing estimates and scaling relative to the non-deformed theory by three polynomial powers from O(d8)O(d^8)O(d8) to O(d5)O(d^5)O(d5). Examining the stronger q-deformed gauge constraint, which softens the total flux at vertices, we show that the physical Hilbert space dimension of the deformed plaquette operator scales equivalently to its non-deformed counterpart with a constant factor 0.2563(5)0.2563(5)0.2563(5). Thus, despite affecting interactions at all scales as exemplified by the observed flux hierarchy inversion symmetry, q-deformation continues to pass scrutiny as a reliable truncation offering advantages in quantum circuit synthesis.
Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2+12+12+1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.
Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2 + 1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.
Understanding strongly correlated electron systems remains one of the central challenges in condensed matter physics. Finding ab initio descriptions of non-perturbative systems is of great importance for modeling carbon-based materials such as graphene or fullerenes. In this regime, computational methods become indispensable tools for performing simulations, making predictions, and ultimately comparing to experiments. The Hubbard model offers a minimal description of the dynamics of electrons in such materials. As exact methods are limited to small system sizes by the exponentially growing Hilbert space, larger-scale first-principle studies rely on numerical approaches. Alongside powerful Hamiltonian-based methods, the auxiliary-field description of the model offers a natural setting for Lagrangian-based Monte Carlo sampling. However, substantial challenges, such as lack of ergodicity and the numerical sign problem, arise in these approaches, leading to long autocorrelation times, biased results, and large statistical uncertainties. In recent years, generative models, especially normalizing flows, have been applied extensively to various theories to address shortcomings of traditional Monte Carlo-based algorithms, such as critical slowing down or topological freezing. While these approaches have shown promise in alleviating these issues, a poorly understood and unfavorable scaling of these methods has stalled their progress, and whether they can reach state-of-the-art system sizes in lattice quantum chromodynamics remains to be seen. This motivates the application of these methods to condensed matter systems, where physically relevant system sizes are reachable at substantially smaller scales. In this work, we study for the first time the application of generative models to the finite-temperature auxiliary-field Hubbard model. We begin by employing an equivariant normalizing flow that respects multiple symmetries inherent to the Boltzmann distribution, demonstrating a substantial improvement in sampling efficiency over Hybrid Monte Carlo and standard non-equivariant flows. We further develop Symmetry-Enforcing Stochastic Modulation (SESaMo), a more scalable algorithm capable of modeling both exact and broken symmetries in normalizing flows, and show that it achieves significantly higher effective sample sizes than traditional equivariant techniques at large system sizes. Equipped with SESaMo, we compute physical observables at half-filling in excellent agreement with exact diagonalization at substantially lower computational cost than Hybrid Monte Carlo. We then turn to the doped Hubbard model, where the theory exhibits a numerical sign problem. We introduce an annealing scheme to simulate the Hubbard model at finite chemical potential in the spin basis. Making use of the less severe sign problem in this description, the annealing scheme enables ergodic sampling of the multi-modal distribution, yielding more accurate results with substantially reduced uncertainties compared to state-of-the-art optimized Hybrid Monte Carlo. Lastly, we employ a masked autoencoder for distribution estimation to learn the Boltzmann distribution of the discrete auxiliary-field Hubbard model. Preliminary results show a significantly higher effective sample size compared to normalizing flows with continuous degrees of freedom, while also exhibiting a milder sign problem than the continuous spin basis.