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:2312.16865 · International Symposium on Lattice Field Theory (LATTICE)

Cluster-Algorithm-Amenable Models of Gauge Fields and Matter

Emilie Huffman, Debasish Banerjee

Published May 3, 2024• 1 citation
arXivHTMLPDF

Abstract

Typical fermion algorithms require the computation (or sampling) of the fermion determinant. We focus instead on cluster algorithms which do not involve the determinant and involve a more physically relevant sampling of the configuration space. We develop new cluster algorithms and design classes of models for fermions coupled to $\mathbb{Z}_2$ and $U(1)$ gauge fields that are amenable to being simulated by these cluster algorithms in a sign-problem free way. Such simulations should contain rich phase diagrams and are particularly relevant for quantum simulator experiments.

DOI
10.22323/1.453.0016

Related papers

Quantum Monte Carlo for gauge fields and matter without the fermion determinant

Debasish Banerjee, Emilie Huffman

Monte Carlo simulations of strongly interacting fermionic systems are plagued by the fermion sign problem, making the nonperturbative study of many interesting regimes of dense quantum matter, or of theories of odd numbers of fermion flavors, challenging. Moreover, typical fermion algorithms require the computation (or sampling) of the fermion determinant. We focus instead on the meron cluster algorithm, which can solve the fermion sign problem in a class of models without involving the determinant. We develop and benchmark new meron algorithms to simulate fermions coupled to Z2 and U(1) gauge fields in the presence of appropriate four-fermi interactions. Such algorithms can be used to uncover potential exotic properties of matter, particularly relevant for quantum simulator experiments. We demonstrate the emergence of the Gauss’ law at low temperatures for a U(1) model in (1+1)D . Published by the American Physical Society 2024

arXiv:2305.08917•Physical review. D/Physical review. D.•Feb 27, 2024•5 citations

Meron-Cluster Algorithms for Quantum Link Models

João C. Pinto Barros, Thea Budde, Marina Krstić Marinković

State-of-the-art algorithms for simulating fermions coupled to gauge fields often rely on integrating fermion degrees of freedom. While successful in simulating QCD at zero chemical potential, at finite density these approaches are hindered by the sign problem. We discuss the simulation of lattice gauge theories in the Hamiltonian formalism and present a generalized meron-cluster algorithm for the simulation of the $U\left(1\right)$ Quantum Link Model for spin $1/2$. This enables the study of models directly relevant to current quantum simulators and is a promising first step toward constructing new efficient algorithms for more complicated gauge theories.

arXiv:2402.01039•International Symposium on Lattice Field Theory (LATTICE)•May 4, 2024•4 citations

Variational Monte Carlo algorithm for lattice gauge theories with continuous gauge groups: A study of ( 2 + 1 ) -dimensional compact QED with dynamical fermions at finite density

Julian Bender, Patrick Emonts, J. I. Cirac

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 $\text{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.

arXiv:2304.05916•Physical Review Research•Nov 7, 2023•13 citations

Field and Gauge Theories with Ultracold Gauge Potentials and Fields

João C. Pinto Barros

In the last decade there has been an intense activity aimed at the quantum simulation of interacting many-body systems using cold atoms [1, 2]. The idea of quantum simulations traces back to Feynman [3], who argued that the ideal setting to study quantum systems would be a quantum experimental setup rather then a classical one - the latter one being fundamentally limited due to its hardware classical structure. This is a particularly important problem given the intrinsic complexity of interacting many-body problems, and the difficulties that arise when tackling them with numerical simulations - two paradigmatic examples being the sign(s) problem affecting Monte Carlo simulations of fermionic systems, and the real time dynamics in more than one spatial dimension. Ultracold atoms offer a very powerful setting for quantum simulations. Atoms can be trapped in tailored optical and magnetic potentials, also controlling their dimensionality. The inter-atomic interactions can be tuned by external knobs, such as Feshbach resonances. This gives a large freedom on model building and, with suitable mappings, they allow the implementation of desired target models. This allowed an impressive exploitation of quantum simulators on the context of condensed matter physics. The simulation of high-energy physics is an important line of research in this field and it is less direct. In particular it requires the implementation of symmetries like Lorentz and gauge invariance which are not immediately available in a cold atomic setting. Gauge fields are ubiquitous in physics ranging from condensed matter [4–6] and quantum computation [7, 8] to particle physics [9], an archetypical example being Quantum Chromodynamics (QCD) [10,11], the theory of strong nuclear forces. Currently open problems in QCD, providing a long-term goal of cold atomic simulations, include confinement/deconfinement and the structure of color superconducting phases at finite chemical potential [12]. Even though QCD is a very complicated theory (to simulate or study), it is possible to envision a path through implementation of simpler models. Furthermore, it is also expected that interesting physics is found on such “intermediate models” which may deserve attention irrespectively of the QCD study. A very relevant model in this regard is the Schwinger model (Quantum Electrodynamics in 1+1 dimensions) [13]. This theory exhibits features of QCD, such as confinement [14], and is at the very same time amenable to both theoretical studies and simpler experimental schemes. This model was the target of the first experimental realization of a gauge theory with a quantum simulator [15]. The work on this Thesis is, in part, motivated by the study of toy models which put in evidence certain aspects that can be found in QCD. Such toy models provide also intermediate steps in the path towards more complex simulations. The two main aspects of QCD which are addressed here are symmetry-locking and confinement. The other main motivation for this study is to develop a systematic framework, through dimensional mismatch, for theoretical understanding and quantum simulations of long-range theories using gauge theories. The model used to study symmetry-locking consists of a four-fermion mixture [16]. It has the basic ingredients to exhibit a non-Abelian symmetry-locked phase: the full Hamiltonian has an SU (2) × SU (2) (global) symmetry which can break to a smaller SU (2) group. Such phase is found in a extensive region of the phase diagram by using a mean-field approach and a strong coupling expansion. A possible realization of such system is provided by an Ytterbium mixture. Even without tuning interactions, it is shown that such mixture falls inside the the locked-symmetry phase pointing towards a possible realization in current day experiments. The models with dimensional mismatch investigated here have fermions in a lower dimensionality d + 1 and gauge fields in higher dimensionality D + 1. They serve two purposes: establish mappings to non-local theories by integration of fields [17] and the study of confinement [18]. In the particular case of d = 1 and D = 2 it is found that some general non-local terms can be obtained on the Lagrangian [17]. This is found in the form of power-law expansions of the Laplacian mediating either kinetic terms (for bosons) or interactions (for fermions). The fact that such expansions are not completely general is not surprising since constraints do exist, preventing unphysical features like breaking of unitarity. The non-local terms obtained are physically acceptable, in this regard, since they are derived from unitary theories. The above mapping is done exactly. In certain cases it is shown that it is possible to construct an effective long-range Hamiltonian in a perturbative expansion. In particular it is shown how this is done for non-relativistic fermions (in d = 1) and 3 + 1 gauge fields. These results are relevant in the context of state of the art experiments which implement models with long-range interactions and where theoretical results are less abundant than for the case of local theories. The above mappings establish a direct relation with local theories which allow theoretical insight onto these systems. Examples of this would consist on the application of Mermin- Wagner-Hohenberg theorem [19, 20] and Lieb-Robinson bounds [21], on the propagation of quantum correlations, to non-local models. In addition they can also provide a path towards implementation of tunable long-range interactions with cold atoms. Furthermore, in terms of quantum simulations, they are in between the full higher dimensional system and the full lower dimensional one. Such property is attractive from the point of view of a gradual increase of complexity for quantum simulations of gauge theories. Another interesting property of these models is that they allow the study of confinement beyond the simpler case of the Schwinger model. The extra dimensions are enough to attribute dynamics to the gauge field, which are no longer completely fixed by the Gauss law. The phases of the Schwinger model are shown to be robust under variation of the dimension of the gauge fields [18]. Both the screened phase, of the massless case, and the confined phase, of the massive one, are found for gauge fields in 2+1 and 3+1 dimensions. Such results are also obtained in the Schwinger- Thirring model. This shows that these phases are very robust and raises interesting questions about the nature of confinement. Robustness under Thirring interactions are relevant because it shows that errors on the experimental implementation will not spoil the phase. Even more interesting is the case of gauge fields in higher dimensions since confinement in the Schwinger model is intuitively atributed to the dimensionality of the gauge fields (creating linear potentials between particles). This Thesis is organized as follows. In Chapter 1, some essential background regarding quantum simulations of gauge theories is provided. It gives both a brief introduction to cold atomic physics and lattice gauge theories. In Chapter 2, it is presented an overview over proposals of quantum simulators of gauge potentials and gauge fields. At the end of this Chapter, in Section 2.3 it is briefly presented ongoing work on a realization of the Schwinger model that we term Half Link Schwinger model. There, it is argued, some of the generators of the gauge symmetry on the lattice can be neglected without comprimising gauge invariance. In Chapter 3, the results regarding the phase diagram of the fourfermion mixture, exhibiting symmetry-locking, is presented. In Chapter 4, the path towards controlling non-local kinetic terms and interactions is provided, after a general introduction to the formalism of dimensional mismatch. At the end the construction of effective Hamiltonians is described. Finally, Chapter 5 concerns the study of confinement and the robustness of it for 1+1 fermions. The first part regards the Schwinger-Thirring with the presence of a -term while, in the second part, models with dimensional mismatch are considered. The thesis ends with conclusions and perspectives of future work based on the results presented here.

•Oct 4, 2017•0 citations

Quantum Simulation of Gauge Theories for Particle and Nuclear Physics

Zohreh Davoudi

Lattice field theory, along with its algorithmic and hardware ecosystems, has been at the forefront of computational particle and nuclear physics. It continues to deliver impressive results on the hadronic spectrum, structure, decays, and reactions. Yet, this vigorous campaign has fallen short in addressing a range of problems involving dense matter and general dynamical phenomena. The reason is that such problems require an exponential scaling of computing time and space in system size. Quantum simulation, enabled by quantum-computing algorithms and hardware technology, promises a way forward by offering several polynomially efficient algorithms compared with their inefficient classical counterparts. Lattice gauge theorists have engaged in a multi-pronged program to leverage such new possibilities, and have steadily advanced the state of theory, algorithm, and hardware implementations and co-design. In this talk, I motivate the quantum-computational lattice-field-theory program; introduce the questions such a program is expected to address and the strategies it involves; report on recent progress; and end with a note on challenges and opportunities ahead.

hep-latCross-listedarXiv:2605.20417•May 21, 2026•0 citations

Meron-Cluster Solution of Fermion Sign Problems

Shailesh Chandrasekharan A, U.‐J. Wiese

We present a general strategy to solve the notorious fermion sign problem using cluster algorithms. The method applies to various systems in the Hubbard model family as well as to relativistic fermions. Here it is illustrated for nonrelativistic lattice fermions. A configuration of fermion world lines is decomposed into clusters that contribute independently to the fermion permutation sign. A cluster whose flip changes the sign is referred to as a meron. Configurations containing meron clusters contribute $0$ to the path integral, while all other configurations contribute $1$. The cluster representation describes the partition function as a gas of clusters in the zero-meron sector.

arXiv:cond-mat/9902128•Physical Review Letters•Oct 18, 1999•223 citations
Typical fermion algorithms require the computation (or sampling) of the fermion determinant. We focus instead on cluster algorithms which do not involve the determinant and involve a more physically relevant sampling of the configuration space. We develop new cluster algorithms and design classes of models for fermions coupled to Z2\mathbb{Z}_2Z2​ and U(1)U(1)U(1) gauge fields that are amenable to being simulated by these cluster algorithms in a sign-problem free way. Such simulations should contain rich phase diagrams and are particularly relevant for quantum simulator experiments.
Monte Carlo simulations of strongly interacting fermionic systems are plagued by the fermion sign problem, making the nonperturbative study of many interesting regimes of dense quantum matter, or of theories of odd numbers of fermion flavors, challenging. Moreover, typical fermion algorithms require the computation (or sampling) of the fermion determinant. We focus instead on the meron cluster algorithm, which can solve the fermion sign problem in a class of models without involving the determinant. We develop and benchmark new meron algorithms to simulate fermions coupled to Z2 and U(1) gauge fields in the presence of appropriate four-fermi interactions. Such algorithms can be used to uncover potential exotic properties of matter, particularly relevant for quantum simulator experiments. We demonstrate the emergence of the Gauss’ law at low temperatures for a U(1) model in (1+1)D . Published by the American Physical Society 2024
State-of-the-art algorithms for simulating fermions coupled to gauge fields often rely on integrating fermion degrees of freedom. While successful in simulating QCD at zero chemical potential, at finite density these approaches are hindered by the sign problem. We discuss the simulation of lattice gauge theories in the Hamiltonian formalism and present a generalized meron-cluster algorithm for the simulation of the U(1)U\left(1\right)U(1) Quantum Link Model for spin 1/21/21/2. This enables the study of models directly relevant to current quantum simulators and is a promising first step toward constructing new efficient algorithms for more complicated gauge theories.
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)(2+1)(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)\text{U}(1)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.
In the last decade there has been an intense activity aimed at the quantum simulation of interacting many-body systems using cold atoms [1, 2]. The idea of quantum simulations traces back to Feynman [3], who argued that the ideal setting to study quantum systems would be a quantum experimental setup rather then a classical one - the latter one being fundamentally limited due to its hardware classical structure. This is a particularly important problem given the intrinsic complexity of interacting many-body problems, and the difficulties that arise when tackling them with numerical simulations - two paradigmatic examples being the sign(s) problem affecting Monte Carlo simulations of fermionic systems, and the real time dynamics in more than one spatial dimension. Ultracold atoms offer a very powerful setting for quantum simulations. Atoms can be trapped in tailored optical and magnetic potentials, also controlling their dimensionality. The inter-atomic interactions can be tuned by external knobs, such as Feshbach resonances. This gives a large freedom on model building and, with suitable mappings, they allow the implementation of desired target models. This allowed an impressive exploitation of quantum simulators on the context of condensed matter physics. The simulation of high-energy physics is an important line of research in this field and it is less direct. In particular it requires the implementation of symmetries like Lorentz and gauge invariance which are not immediately available in a cold atomic setting. Gauge fields are ubiquitous in physics ranging from condensed matter [4–6] and quantum computation [7, 8] to particle physics [9], an archetypical example being Quantum Chromodynamics (QCD) [10,11], the theory of strong nuclear forces. Currently open problems in QCD, providing a long-term goal of cold atomic simulations, include confinement/deconfinement and the structure of color superconducting phases at finite chemical potential [12]. Even though QCD is a very complicated theory (to simulate or study), it is possible to envision a path through implementation of simpler models. Furthermore, it is also expected that interesting physics is found on such “intermediate models” which may deserve attention irrespectively of the QCD study. A very relevant model in this regard is the Schwinger model (Quantum Electrodynamics in 1+1 dimensions) [13]. This theory exhibits features of QCD, such as confinement [14], and is at the very same time amenable to both theoretical studies and simpler experimental schemes. This model was the target of the first experimental realization of a gauge theory with a quantum simulator [15]. The work on this Thesis is, in part, motivated by the study of toy models which put in evidence certain aspects that can be found in QCD. Such toy models provide also intermediate steps in the path towards more complex simulations. The two main aspects of QCD which are addressed here are symmetry-locking and confinement. The other main motivation for this study is to develop a systematic framework, through dimensional mismatch, for theoretical understanding and quantum simulations of long-range theories using gauge theories. The model used to study symmetry-locking consists of a four-fermion mixture [16]. It has the basic ingredients to exhibit a non-Abelian symmetry-locked phase: the full Hamiltonian has an SU (2) × SU (2) (global) symmetry which can break to a smaller SU (2) group. Such phase is found in a extensive region of the phase diagram by using a mean-field approach and a strong coupling expansion. A possible realization of such system is provided by an Ytterbium mixture. Even without tuning interactions, it is shown that such mixture falls inside the the locked-symmetry phase pointing towards a possible realization in current day experiments. The models with dimensional mismatch investigated here have fermions in a lower dimensionality d + 1 and gauge fields in higher dimensionality D + 1. They serve two purposes: establish mappings to non-local theories by integration of fields [17] and the study of confinement [18]. In the particular case of d = 1 and D = 2 it is found that some general non-local terms can be obtained on the Lagrangian [17]. This is found in the form of power-law expansions of the Laplacian mediating either kinetic terms (for bosons) or interactions (for fermions). The fact that such expansions are not completely general is not surprising since constraints do exist, preventing unphysical features like breaking of unitarity. The non-local terms obtained are physically acceptable, in this regard, since they are derived from unitary theories. The above mapping is done exactly. In certain cases it is shown that it is possible to construct an effective long-range Hamiltonian in a perturbative expansion. In particular it is shown how this is done for non-relativistic fermions (in d = 1) and 3 + 1 gauge fields. These results are relevant in the context of state of the art experiments which implement models with long-range interactions and where theoretical results are less abundant than for the case of local theories. The above mappings establish a direct relation with local theories which allow theoretical insight onto these systems. Examples of this would consist on the application of Mermin- Wagner-Hohenberg theorem [19, 20] and Lieb-Robinson bounds [21], on the propagation of quantum correlations, to non-local models. In addition they can also provide a path towards implementation of tunable long-range interactions with cold atoms. Furthermore, in terms of quantum simulations, they are in between the full higher dimensional system and the full lower dimensional one. Such property is attractive from the point of view of a gradual increase of complexity for quantum simulations of gauge theories. Another interesting property of these models is that they allow the study of confinement beyond the simpler case of the Schwinger model. The extra dimensions are enough to attribute dynamics to the gauge field, which are no longer completely fixed by the Gauss law. The phases of the Schwinger model are shown to be robust under variation of the dimension of the gauge fields [18]. Both the screened phase, of the massless case, and the confined phase, of the massive one, are found for gauge fields in 2+1 and 3+1 dimensions. Such results are also obtained in the Schwinger- Thirring model. This shows that these phases are very robust and raises interesting questions about the nature of confinement. Robustness under Thirring interactions are relevant because it shows that errors on the experimental implementation will not spoil the phase. Even more interesting is the case of gauge fields in higher dimensions since confinement in the Schwinger model is intuitively atributed to the dimensionality of the gauge fields (creating linear potentials between particles). This Thesis is organized as follows. In Chapter 1, some essential background regarding quantum simulations of gauge theories is provided. It gives both a brief introduction to cold atomic physics and lattice gauge theories. In Chapter 2, it is presented an overview over proposals of quantum simulators of gauge potentials and gauge fields. At the end of this Chapter, in Section 2.3 it is briefly presented ongoing work on a realization of the Schwinger model that we term Half Link Schwinger model. There, it is argued, some of the generators of the gauge symmetry on the lattice can be neglected without comprimising gauge invariance. In Chapter 3, the results regarding the phase diagram of the fourfermion mixture, exhibiting symmetry-locking, is presented. In Chapter 4, the path towards controlling non-local kinetic terms and interactions is provided, after a general introduction to the formalism of dimensional mismatch. At the end the construction of effective Hamiltonians is described. Finally, Chapter 5 concerns the study of confinement and the robustness of it for 1+1 fermions. The first part regards the Schwinger-Thirring with the presence of a -term while, in the second part, models with dimensional mismatch are considered. The thesis ends with conclusions and perspectives of future work based on the results presented here.
Lattice field theory, along with its algorithmic and hardware ecosystems, has been at the forefront of computational particle and nuclear physics. It continues to deliver impressive results on the hadronic spectrum, structure, decays, and reactions. Yet, this vigorous campaign has fallen short in addressing a range of problems involving dense matter and general dynamical phenomena. The reason is that such problems require an exponential scaling of computing time and space in system size. Quantum simulation, enabled by quantum-computing algorithms and hardware technology, promises a way forward by offering several polynomially efficient algorithms compared with their inefficient classical counterparts. Lattice gauge theorists have engaged in a multi-pronged program to leverage such new possibilities, and have steadily advanced the state of theory, algorithm, and hardware implementations and co-design. In this talk, I motivate the quantum-computational lattice-field-theory program; introduce the questions such a program is expected to address and the strategies it involves; report on recent progress; and end with a note on challenges and opportunities ahead.
We present a general strategy to solve the notorious fermion sign problem using cluster algorithms. The method applies to various systems in the Hubbard model family as well as to relativistic fermions. Here it is illustrated for nonrelativistic lattice fermions. A configuration of fermion world lines is decomposed into clusters that contribute independently to the fermion permutation sign. A cluster whose flip changes the sign is referred to as a meron. Configurations containing meron clusters contribute 000 to the path integral, while all other configurations contribute 111. The cluster representation describes the partition function as a gas of clusters in the zero-meron sector.