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:2605.08042

The two-flavor Schwinger model at 50: Solving Coleman's puzzles

Gabriel Cuomo, Ross Dempsey, Andrei Katsevich, Igor R. Klebanov, Ilia V. Kochergin, Silviu S. Pufu, Benjamin T. Søgaard

hep-thhep-latquant-phSubmitted May 8, 2026• Announced May 11, 2026
arXivHTMLPDF

Abstract

In his 1976 paper "More about the massive Schwinger model", Coleman introduced $1+1$-dimensional Quantum Electrodynamics coupled to two charged massive fermions. By applying Abelian bosonization, he elucidated much of the physics of this two-flavor Schwinger model, but he listed three puzzles at the end of his paper. We present new analytical and numerical calculations to solve Coleman's three puzzles and thereby deepen our understanding of this model. These puzzles pertain to the theory with equal fermion masses at $θ= 0$ and at $θ= π$, as well as the size of isospin-breaking effects when the fermion masses are unequal. For the puzzle at $θ= π$, the solution is related to the structure of the zero-temperature phase diagram arXiv:2305.04437: for equal fermion masses $m$, the model exhibits spontaneous breaking of charge conjugation symmetry and absence of confinement for any value of the gauge coupling $g$, so that there is a smooth interpolation from weak to strong coupling. Using two-loop Renormalization Group and integrability methods, we show that the mass gap behaves as $\sim m e^{-0.111 g^2/m^2}$ in the strong coupling regime $m\ll g$. Our numerical results using the lattice Hamiltonian are in good agreement with this behavior. For the puzzle at $θ= 0$, the solution is related to a level crossing between two isosinglet particles with different discrete quantum numbers; we demonstrate the necessity of such a crossing by comparing integrability and weak coupling calculations, and we also exhibit the crossing numerically. Finally, we provide a new estimate for the size of isospin-breaking effects caused by different fermion masses at strong coupling.

Comments
68 pages, 14 figures

Related papers

Fermion discretization effects in the two-flavor lattice Schwinger model: A study with matrix product states

Tim Schwägerl, Karl Jansen, Stefan Kühn

We present a comprehensive tensor network study of staggered, Wilson, and twisted mass fermions in the Hamiltonian formulation, using the massive two-flavor Schwinger model as a benchmark. Particular emphasis is placed on twisted mass fermions, whose properties in this context have not been systematically explored before. We confirm the expected O ( a ) improvement in the free theory and observe that this improvement persists in the interacting case. By leveraging an electric-field-based method for mass renormalization, we reliably tune to maximal twist and establish the method’s applicability in the two-flavor model. Once mass renormalization is included, the pion mass exhibits rapid convergence to the continuum limit. Finite-volume effects are addressed using two complementary approaches: dispersion relation fits and finite-volume scaling. Our results show excellent agreement with semiclassical predictions and reveal a milder volume dependence for twisted mass fermions compared to staggered and Wilson discretizations. In addition, we observe clear isospin-breaking effects, suggesting intriguing parallels with lattice QCD. These findings highlight the advantages of twisted mass fermions for Hamiltonian simulations and motivate their further exploration—particularly in view of future applications to higher-dimensional lattice gauge theories.

arXiv:2509.02329•Physical review. D/Physical review. D.•Oct 16, 2025•1 citation

Phase transitions in Zn gauge models: Towards quantum simulations of the Schwinger-Weyl QED

Elisa Ercolessi, Paolo Facchi A, Giuseppe Magnifico +2 more

We study the ground-state properties of a class of ${\mathbb{Z}}_{n}$ lattice gauge theories in $1+1$ dimensions, in which the gauge fields are coupled to spinless fermionic matter. These models, stemming from discrete representations of the Weyl commutator for the U(1) group, preserve the unitary character of the minimal coupling and have, therefore, the property of formally approximating lattice quantum electrodynamics in one spatial dimension in the large-$n$ limit. The numerical study of such approximated theories is important to determine their effectiveness in reproducing the main features and phenomenology of the target theory, in view of implementations of cold-atom quantum simulators of QED. In this paper, we study the cases $n=2\textdiv{}8$ by means of a DMRG code that exactly implements Gauss's law. We perform a careful scaling analysis and show that, in absence of a background field, all ${\mathbb{Z}}_{n}$ models exhibit a phase transition which falls in the Ising universality class, with spontaneous symmetry breaking of the $CP$ symmetry. We then perform the large-$n$ limit and find that the asymptotic values of the critical parameters approach the ones obtained for the known phase transition of the zero-charge sector of the massive Schwinger model, which occurs at negative mass.

arXiv:1705.11047•Physical review. D/Physical review. D.•Oct 5, 2018•81 citations

Calculating composite-particle spectra in Hamiltonian formalism and demonstration in 2-flavor QED$_{1+1\text{d}}$

Etsuko Itou, Akira Matsumoto, Yuya Tanizaki

We consider three distinct methods to compute the mass spectrum of gauge theories in the Hamiltonian formalism: (1) correlation-function scheme, (2) one-point-function scheme, and (3) dispersion-relation scheme. The first one examines spatial correlation functions as we do in the conventional Euclidean Monte Carlo simulations. The second one uses the boundary effect to efficiently compute the mass spectrum. The third one constructs the excited states and fits their energy using the dispersion relation with selecting quantum numbers. Each method has its pros and cons, and we clarify such properties in their applications to the mass spectrum for the 2-flavor massive Schwinger model at $m/g=0.1$ and $θ=0$ using the density-matrix renormalization group (DMRG). We note that the multi-flavor Schwinger model at small mass $m$ is a strongly coupled field theory even after the bosonizations, and thus it deserves to perform the first-principles numerical calculations. All these methods mostly agree and identify the stable particles, pions $π_a$ ($J^{PG}=1^{-+}$), sigma meson $σ$ ($J^{PG}=0^{++}$), and eta meson $η$ ($J^{PG}=0^{--}$). In particular, we find that the mass of $σ$ meson is lighter than twice the pion mass, and thus $σ$ is stable against the decay process, $σ\to ππ$. This is consistent with the analytic prediction using the WKB approximation, and, remarkably, our numerical results are so close to the WKB-based formula between the pion and sigma-meson masses, $M_σ/M_π=\sqrt{3}$.

arXiv:2307.16655•Jul 31, 2023•0 citations

Berezinskii-Kosterlitz-Thouless transition in lattice Schwinger model with one flavor of Wilson fermion

Yuya Shimizu, Y. Kuramashi

We have made a detailed study of the phase structure for the lattice Schwinger model with one flavor of Wilson fermion on the $(m,g)$ plane. For numerical investigation, we develop a decorated tensor renormalization method for lattice gauge theories with fermions incorporating the Grassmann tensor renormalization. Our algorithm manifestly preserves rotation and reflection symmetries. We find not only a parity-broken phase but also a Berezinskii-Kosterlitz-Thouless (BKT) transition by evaluating the central charge and an expectation value of a projection operator into the parity-odd subspace. The BKT phase boundaries converge into the degenerated doubler pole $(m,g)=(\ensuremath{-}2,0)$, while the parity-breaking transition line ends at the physical pole $(m,g)=(0,0)$. In addition, our analysis of scaling dimensions indicates that a conformal field theory with SU(2) symmetry arises on the line of $m=\ensuremath{-}2$.

arXiv:1712.07808•Physical review. D/Physical review. D.•Feb 6, 2018•65 citations

Calculating composite-particle spectra in Hamiltonian formalism and demonstration in 2-flavor QED1+1d

Etsuko Itou, Akira Matsumoto, Yuya Tanizaki

Abstract We consider three distinct methods to compute the mass spectrum of gauge theories in the Hamiltonian formalism: (1) correlation-function scheme, (2) one-point-function scheme, and (3) dispersion-relation scheme. The first one examines spatial correlation functions as we do in the conventional Euclidean Monte Carlo simulations. The second one uses the boundary effect to efficiently compute the mass spectrum. The third one constructs the excited states and fits their energy using the dispersion relation with selecting quantum numbers. Each method has its pros and cons, and we clarify such properties in their applications to the mass spectrum for the 2-flavor massive Schwinger model at m/g = 0.1 and θ = 0 using the density-matrix renormalization group (DMRG). We note that the multi-flavor Schwinger model at small mass m is a strongly coupled field theory even after the bosonizations, and thus it deserves to perform the first-principles numerical calculations. All these methods mostly agree and identify the stable particles, pions πa (JPG = 1−+), sigma meson σ (JPG = 0++), and eta meson η (JPG = 0−−). In particular, we find that the mass of σ meson is lighter than twice the pion mass, and thus σ is stable against the decay process, σ → ππ. This is consistent with the analytic prediction using the WKB approximation, and, remarkably, our numerical results are so close to the WKB-based formula between the pion and sigma-meson masses, Mσ/Mπ = $$ \sqrt{3} $$ 3 .

Journal of High Energy Physics•Nov 30, 2023•14 citations

Flavoured Lattice Schwinger Model with Chiral Anomaly

Dogukan Bakircioglu

We introduce the \emph{flavoured lattice Schwinger model}, a $(1{+}1)$-dimensional $U(1)$ lattice gauge theory in which the fermion doubling problem is resolved by staggering a $\mathbb{Z}_{2}$ flavour degree of freedom rather than staggering chirality. Unlike the standard approaches, this construction preserves an exact axial $U(1)$ symmetry at finite lattice spacing. We derive the continuum limit, showing that the model reduces to the \emph{two-flavour} massless Schwinger model, with flavours $α\in\{0,1\}$ sharing one dynamical $U(1)$ gauge field. The central result is a well-defined, regularised, gauge-invariant lattice axial charge $Q_{G}^{A}$ whose continuum non-conservation $\langle dQ_{G}^{A}/dt\rangle = -(2g/π)\!\int\! dx\,\langle E(x)\rangle$ arises as a direct dynamical consequence of minimal gauge coupling. A particle-hole transformation on the $χ$ flavour exposes a hidden $U_{L}(2)\times U_{R}(2)$ chiral symmetry; non-Abelian bosonisation then identifies the model with a massive abelian Schwinger sector tensored with the level-$1$ $SU(2)$ Wess--Zumino--Witten model. Finally, we show that embedding the flavoured fermions in a ribbon-shaped $(2{+}1)$D Bernevig--Hughes--Zhang topological insulator and gauging the bulk in a constant background field factorises the boundary theory into \emph{two decoupled} single-flavour Schwinger models, one on each edge, identifying the lattice factor of $2$ as one quantum of Schwinger anomaly per edge.

hep-latCross-listedarXiv:2604.13146•Apr 15, 2026
In his 1976 paper "More about the massive Schwinger model", Coleman introduced 1+11+11+1-dimensional Quantum Electrodynamics coupled to two charged massive fermions. By applying Abelian bosonization, he elucidated much of the physics of this two-flavor Schwinger model, but he listed three puzzles at the end of his paper. We present new analytical and numerical calculations to solve Coleman's three puzzles and thereby deepen our understanding of this model. These puzzles pertain to the theory with equal fermion masses at θ=0θ= 0θ=0 and at θ=πθ= πθ=π, as well as the size of isospin-breaking effects when the fermion masses are unequal. For the puzzle at θ=πθ= πθ=π, the solution is related to the structure of the zero-temperature phase diagram arXiv:2305.04437: for equal fermion masses mmm, the model exhibits spontaneous breaking of charge conjugation symmetry and absence of confinement for any value of the gauge coupling ggg, so that there is a smooth interpolation from weak to strong coupling. Using two-loop Renormalization Group and integrability methods, we show that the mass gap behaves as ∼me−0.111g2/m2\sim m e^{-0.111 g^2/m^2}∼me−0.111g2/m2 in the strong coupling regime m≪gm\ll gm≪g. Our numerical results using the lattice Hamiltonian are in good agreement with this behavior. For the puzzle at θ=0θ= 0θ=0, the solution is related to a level crossing between two isosinglet particles with different discrete quantum numbers; we demonstrate the necessity of such a crossing by comparing integrability and weak coupling calculations, and we also exhibit the crossing numerically. Finally, we provide a new estimate for the size of isospin-breaking effects caused by different fermion masses at strong coupling.
We present a comprehensive tensor network study of staggered, Wilson, and twisted mass fermions in the Hamiltonian formulation, using the massive two-flavor Schwinger model as a benchmark. Particular emphasis is placed on twisted mass fermions, whose properties in this context have not been systematically explored before. We confirm the expected O ( a ) improvement in the free theory and observe that this improvement persists in the interacting case. By leveraging an electric-field-based method for mass renormalization, we reliably tune to maximal twist and establish the method’s applicability in the two-flavor model. Once mass renormalization is included, the pion mass exhibits rapid convergence to the continuum limit. Finite-volume effects are addressed using two complementary approaches: dispersion relation fits and finite-volume scaling. Our results show excellent agreement with semiclassical predictions and reveal a milder volume dependence for twisted mass fermions compared to staggered and Wilson discretizations. In addition, we observe clear isospin-breaking effects, suggesting intriguing parallels with lattice QCD. These findings highlight the advantages of twisted mass fermions for Hamiltonian simulations and motivate their further exploration—particularly in view of future applications to higher-dimensional lattice gauge theories.
We study the ground-state properties of a class of Zn{\mathbb{Z}}_{n}Zn​ lattice gauge theories in 1+11+11+1 dimensions, in which the gauge fields are coupled to spinless fermionic matter. These models, stemming from discrete representations of the Weyl commutator for the U(1) group, preserve the unitary character of the minimal coupling and have, therefore, the property of formally approximating lattice quantum electrodynamics in one spatial dimension in the large-nnn limit. The numerical study of such approximated theories is important to determine their effectiveness in reproducing the main features and phenomenology of the target theory, in view of implementations of cold-atom quantum simulators of QED. In this paper, we study the cases $n=2\textdiv{}8$ by means of a DMRG code that exactly implements Gauss's law. We perform a careful scaling analysis and show that, in absence of a background field, all Zn{\mathbb{Z}}_{n}Zn​ models exhibit a phase transition which falls in the Ising universality class, with spontaneous symmetry breaking of the CPCPCP symmetry. We then perform the large-nnn limit and find that the asymptotic values of the critical parameters approach the ones obtained for the known phase transition of the zero-charge sector of the massive Schwinger model, which occurs at negative mass.
Calculating composite-particle spectra in Hamiltonian formalism and demonstration in 2-flavor QED1+1d_{1+1\text{d}}1+1d​
We consider three distinct methods to compute the mass spectrum of gauge theories in the Hamiltonian formalism: (1) correlation-function scheme, (2) one-point-function scheme, and (3) dispersion-relation scheme. The first one examines spatial correlation functions as we do in the conventional Euclidean Monte Carlo simulations. The second one uses the boundary effect to efficiently compute the mass spectrum. The third one constructs the excited states and fits their energy using the dispersion relation with selecting quantum numbers. Each method has its pros and cons, and we clarify such properties in their applications to the mass spectrum for the 2-flavor massive Schwinger model at m/g=0.1m/g=0.1m/g=0.1 and θ=0θ=0θ=0 using the density-matrix renormalization group (DMRG). We note that the multi-flavor Schwinger model at small mass mmm is a strongly coupled field theory even after the bosonizations, and thus it deserves to perform the first-principles numerical calculations. All these methods mostly agree and identify the stable particles, pions πaπ_aπa​ (JPG=1−+J^{PG}=1^{-+}JPG=1−+), sigma meson σσσ (JPG=0++J^{PG}=0^{++}JPG=0++), and eta meson ηηη (JPG=0−−J^{PG}=0^{--}JPG=0−−). In particular, we find that the mass of σσσ meson is lighter than twice the pion mass, and thus σσσ is stable against the decay process, σ→ππσ\to ππσ→ππ. This is consistent with the analytic prediction using the WKB approximation, and, remarkably, our numerical results are so close to the WKB-based formula between the pion and sigma-meson masses, Mσ/Mπ=3M_σ/M_π=\sqrt{3}Mσ​/Mπ​=3​.
We have made a detailed study of the phase structure for the lattice Schwinger model with one flavor of Wilson fermion on the (m,g)(m,g)(m,g) plane. For numerical investigation, we develop a decorated tensor renormalization method for lattice gauge theories with fermions incorporating the Grassmann tensor renormalization. Our algorithm manifestly preserves rotation and reflection symmetries. We find not only a parity-broken phase but also a Berezinskii-Kosterlitz-Thouless (BKT) transition by evaluating the central charge and an expectation value of a projection operator into the parity-odd subspace. The BKT phase boundaries converge into the degenerated doubler pole $(m,g)=(\ensuremath{-}2,0)$, while the parity-breaking transition line ends at the physical pole (m,g)=(0,0)(m,g)=(0,0)(m,g)=(0,0). In addition, our analysis of scaling dimensions indicates that a conformal field theory with SU(2) symmetry arises on the line of $m=\ensuremath{-}2$.
Abstract We consider three distinct methods to compute the mass spectrum of gauge theories in the Hamiltonian formalism: (1) correlation-function scheme, (2) one-point-function scheme, and (3) dispersion-relation scheme. The first one examines spatial correlation functions as we do in the conventional Euclidean Monte Carlo simulations. The second one uses the boundary effect to efficiently compute the mass spectrum. The third one constructs the excited states and fits their energy using the dispersion relation with selecting quantum numbers. Each method has its pros and cons, and we clarify such properties in their applications to the mass spectrum for the 2-flavor massive Schwinger model at m/g = 0.1 and θ = 0 using the density-matrix renormalization group (DMRG). We note that the multi-flavor Schwinger model at small mass m is a strongly coupled field theory even after the bosonizations, and thus it deserves to perform the first-principles numerical calculations. All these methods mostly agree and identify the stable particles, pions πa (JPG = 1−+), sigma meson σ (JPG = 0++), and eta meson η (JPG = 0−−). In particular, we find that the mass of σ meson is lighter than twice the pion mass, and thus σ is stable against the decay process, σ → ππ. This is consistent with the analytic prediction using the WKB approximation, and, remarkably, our numerical results are so close to the WKB-based formula between the pion and sigma-meson masses, Mσ/Mπ = 3 \sqrt{3} 3​ 3 .
We introduce the \emph{flavoured lattice Schwinger model}, a (1+1)(1{+}1)(1+1)-dimensional U(1)U(1)U(1) lattice gauge theory in which the fermion doubling problem is resolved by staggering a Z2\mathbb{Z}_{2}Z2​ flavour degree of freedom rather than staggering chirality. Unlike the standard approaches, this construction preserves an exact axial U(1)U(1)U(1) symmetry at finite lattice spacing. We derive the continuum limit, showing that the model reduces to the \emph{two-flavour} massless Schwinger model, with flavours α∈{0,1}α\in\{0,1\}α∈{0,1} sharing one dynamical U(1)U(1)U(1) gauge field. The central result is a well-defined, regularised, gauge-invariant lattice axial charge QGAQ_{G}^{A}QGA​ whose continuum non-conservation ⟨dQGA/dt⟩=−(2g/π) ⁣∫ ⁣dx ⟨E(x)⟩\langle dQ_{G}^{A}/dt\rangle = -(2g/π)\!\int\! dx\,\langle E(x)\rangle⟨dQGA​/dt⟩=−(2g/π)∫dx⟨E(x)⟩ arises as a direct dynamical consequence of minimal gauge coupling. A particle-hole transformation on the χχχ flavour exposes a hidden UL(2)×UR(2)U_{L}(2)\times U_{R}(2)UL​(2)×UR​(2) chiral symmetry; non-Abelian bosonisation then identifies the model with a massive abelian Schwinger sector tensored with the level-111 SU(2)SU(2)SU(2) Wess--Zumino--Witten model. Finally, we show that embedding the flavoured fermions in a ribbon-shaped (2+1)(2{+}1)(2+1)D Bernevig--Hughes--Zhang topological insulator and gauging the bulk in a constant background field factorises the boundary theory into \emph{two decoupled} single-flavour Schwinger models, one on each edge, identifying the lattice factor of 222 as one quantum of Schwinger anomaly per edge.