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:2103.15615 · SciPost Physics

Generalized eigenproblem without fermion doubling for Dirac fermions on a lattice

M. J. Pacholski, G. Lemut, J. Tworzydło, C. W. J. Beenakker

Published Dec 14, 2021• 16 citations
arXivHTMLPDF

Abstract

The spatial discretization of the single-cone Dirac Hamiltonian on the surface of a topological insulator or superconductor needs a special ``staggered’’ grid, to avoid the appearance of a spurious second cone in the Brillouin zone. We adapt the Stacey discretization from lattice gauge theory to produce a generalized eigenvalue problem, of the form \bm{\mathcal H}\bm{\psi}=\bm{E}\bm{\mathcal P}\bm{\psi} ℋ 𝛙 = 𝐄 𝒫 𝛙 , with Hermitian tight-binding operators \bm{\mathcal H} ℋ , \bm{\mathcal P} 𝒫 , a locally conserved particle current, and preserved chiral and symplectic symmetries. This permits the study of the spectral statistics of Dirac fermions in each of the four symmetry classes A, AII, AIII, and D.

DOI
10.21468/scipostphys.11.6.105

Related papers

Quantum circuit decomposition of the tangent-fermion Dirac operator

C. W. J. Beenakker

The Dirac operator on a lattice cannot be both local and free of fermion doubling, at least not without breaking fundamental symmetries. Non-local, symmetry-preserving discretizations that avoid doubling have a quantum circuit representation as a linear-combination-of-unitaries (LCU) in which both the number of terms and their norm (the subnormalization factor) grow with the lattice size, compromising the efficiency of a quantum algorithm. We show that the tangent-fermion discretization escapes this obstruction when the Dirac equation is written as a generalized eigenvalue problem with a local operator pencil: Each member of the pencil has an exact LCU, with term count that is independent of lattice size and with subnormalization factor of order unity, on a par with elliptic operators. This provides an efficient block-encoding primitive for Dirac spectra and Green functions without fermion doubling.

quant-pharXiv:2606.19020•Jun 18, 2026•0 citations

Fermion-doubling problem in Chiral discretizations of Quantum field theory: Definitive proof, Fixing, and Computation of two-point correlation function

Dogukan Bakircioglu, Pablo Arnault

We give the definitive proof that the Dirac Quantum Cellular Automaton (QCA) used for both quantum simulation and algorithmic foundations of Quantum Field Theory (QFT), and especially of Quantum Electrodynamics (QED), as put forward in References https://doi.org/10.1007/s11128-019-2555-4 and https://doi.org/10.22331/q-2023-11-08-1179, does exhibit Fermion Doubling (FD), albeit thrice as less severe as discrete-time standard Lattice Gauge Theories (LGTs) [arXiv:2505.0790], which are naive regarding the spacetime discretization of differential operators acting on fermionic fields. The proof is done for the (1 + 1)D Dirac-QCA model. We show that the (one-time-step) two-point correlation function, also called Green's function (GF), of the Dirac QCA, is of astonishing simplicity, which is in contrast with the GF of the Dirac equation. We also compare, both qualitatively and quantitatively, this Dirac QCA to the continuous-time-LGT spatial discretization of Dirac fermions regarding how well these two lattice models approximate their naive continuum limit$\unicode{x2014}$which is nothing but the Dirac equation$\unicode{x2014}$even when far away from that limit, a situation which must be considered because of experimental limitations in quantum simulation: the Dirac QCA is better for ultrarelativistic regimes, whereas continuous-time LGT is better for non-relativistic regimes. Then, we compute the GF of the FD-fixed model put forward in the last cited reference, called Flavored Dirac QCA (FDQCA)$\unicode{x2014}$which staggers an extra, artificial flavor $only$, on a diamond spacetime lattice, and does not stagger chirality as staggered fermions in usual LGT. The structure of this FDQCA two-point correlation function is of extreme simplicity, and can be expressed in a very simple manner in terms of the four chiral components of the FD-suffering, original-model GF.

hep-latCross-listedarXiv:2607.14874•Jul 17, 2026•0 citations

Artificial staggered magnetic field for ultracold atoms in optical lattices

Lih-King Lim, Andreas Hemmerich, C. Morais Smith

A time-dependent optical lattice with staggered particle current in the tight-binding regime was considered that can be described by a time-independent effective lattice model with an artificial staggered magnetic field. The low-energy description of a single-component fermion in this lattice at half-filling is provided by two copies of ideal two-dimensional massless Dirac fermions. The Dirac cones are generally anisotropic and can be tuned by the external staggered flux $\ensuremath{\phi}$. For bosons, the staggered flux modifies the single-particle spectrum such that in the weak coupling limit, depending on the flux $\ensuremath{\phi}$, distinct superfluid phases are realized. Their properties are discussed, the nature of the phase transitions between them is established, and Bogoliubov theory is used to determine their excitation spectra. Then the generalized superfluid--Mott-insulator transition is studied in the presence of the staggered flux and the complete phase diagram is established. Finally, the momentum distribution of the distinct superfluid phases is obtained, which provides a clear experimental signature of each phase in ballistic expansion experiments.

arXiv:0910.2610•Physical Review A•Feb 4, 2010•61 citations

Single-cone real-space finite difference schemes for the Dirac von Neumann equation

Magdalena Schreilechner, W. Pötz

Two finite difference schemes for the numerical treatment of the von Neumann equation for the (2+1)D Dirac Hamiltonian are presented. Both utilize a single-cone staggered space-time grid which ensures a single-cone energy dispersion to formulate a numerical treatment of the mixed-state dynamics within the von Neumann equation. The first scheme executes the time-derivative according to the product rule for "bra" and "ket" indices of the density operator. It therefore directly inherits all the favorable properties of the difference scheme for the pure-state Dirac equation and conserves positivity. The second scheme proposed here performs the time-derivative in one sweep. This direct scheme is investigated regarding stability and convergence. Both schemes are tested numerically for elementary simulations using parameters which pertain to topological insulator surface states. Application of the schemes to a Dirac Lindblad equation and real-space-time Green's function formulations are discussed.

arXiv:1503.02685•Mar 6, 2015•3 citations

Fragile single-cone Dirac quantum walks in two dimensions

C. W. J. Beenakker, J. Sánchez Férnan, J. Tworzydło

It is known that a one-dimensional (1D) quantum walk gives a local space-time discretization of the massless Dirac equation with a single quasi-energy cone (no fermion doubling at low energies), keeping the fundamental symmetries (chiral and time-reversal) of the continuum theory. We show that the analogous 2D construction is fundamentally more fragile. Local two-band quantum walks can have an unpaired Dirac cone, but the protecting symmetries then cease to be ordinary on-site symmetries: They become non-symmorphic, involving half-lattice translations, and are broken by generic spatial inhomogeneities. In particular, we demonstrate that the 2D Dirac quantum walk based on the Ho-Chalker network model can be gapped by potential scattering, as is evidenced by the breakdown of reflectionless Klein tunneling. The symmorphic symmetry obstruction is sharp: If we allow for couplings that are no longer strictly finite-range but decay exponentially with distance, we can construct a 2D quantum walk with a single Dirac cone and on-site chiral and time-reversal symmetries.

quant-pharXiv:2607.05112•Jul 7, 2026

Interacting chiral fermions on the lattice with matrix product operator norms

Jutho Haegeman, Laurens Lootens, Quinten Mortier +3 more

We develop a Hamiltonian formalism for simulating interacting chiral fermions on the lattice while preserving unitarity and locality and without breaking the chiral symmetry. The fermion doubling problem is circumvented by constructing a Fock space endowed with a semi-definite norm. When projecting our theory on the the single-particle sector, we recover the framework of Stacey fermions, and we demonstrate that the scaling limit of the free model recovers the chiral fermion field. Technically, we make use of a matrix product operator norm to mimick the boundary of a higher dimensional topological theory. As a proof of principle, we consider a single Weyl fermion on a periodic ring with Hubbard-type nearest-neighbor interactions and construct a variational generalized DMRG code to demonstrate that the ground state for large system sizes can be determined efficiently. As our tensor network approach does not exhibit any sign problem, we can add a chemical potential and study real-time evolution.

arXiv:2405.10285•May 16, 2024•1 citation
The spatial discretization of the single-cone Dirac Hamiltonian on the surface of a topological insulator or superconductor needs a special ``staggered’’ grid, to avoid the appearance of a spurious second cone in the Brillouin zone. We adapt the Stacey discretization from lattice gauge theory to produce a generalized eigenvalue problem, of the form \bm{\mathcal H}\bm{\psi}=\bm{E}\bm{\mathcal P}\bm{\psi} ℋ 𝛙 = 𝐄 𝒫 𝛙 , with Hermitian tight-binding operators \bm{\mathcal H} ℋ , \bm{\mathcal P} 𝒫 , a locally conserved particle current, and preserved chiral and symplectic symmetries. This permits the study of the spectral statistics of Dirac fermions in each of the four symmetry classes A, AII, AIII, and D.
The Dirac operator on a lattice cannot be both local and free of fermion doubling, at least not without breaking fundamental symmetries. Non-local, symmetry-preserving discretizations that avoid doubling have a quantum circuit representation as a linear-combination-of-unitaries (LCU) in which both the number of terms and their norm (the subnormalization factor) grow with the lattice size, compromising the efficiency of a quantum algorithm. We show that the tangent-fermion discretization escapes this obstruction when the Dirac equation is written as a generalized eigenvalue problem with a local operator pencil: Each member of the pencil has an exact LCU, with term count that is independent of lattice size and with subnormalization factor of order unity, on a par with elliptic operators. This provides an efficient block-encoding primitive for Dirac spectra and Green functions without fermion doubling.
We give the definitive proof that the Dirac Quantum Cellular Automaton (QCA) used for both quantum simulation and algorithmic foundations of Quantum Field Theory (QFT), and especially of Quantum Electrodynamics (QED), as put forward in References https://doi.org/10.1007/s11128-019-2555-4 and https://doi.org/10.22331/q-2023-11-08-1179, does exhibit Fermion Doubling (FD), albeit thrice as less severe as discrete-time standard Lattice Gauge Theories (LGTs) [arXiv:2505.0790], which are naive regarding the spacetime discretization of differential operators acting on fermionic fields. The proof is done for the (1 + 1)D Dirac-QCA model. We show that the (one-time-step) two-point correlation function, also called Green's function (GF), of the Dirac QCA, is of astonishing simplicity, which is in contrast with the GF of the Dirac equation. We also compare, both qualitatively and quantitatively, this Dirac QCA to the continuous-time-LGT spatial discretization of Dirac fermions regarding how well these two lattice models approximate their naive continuum limit$\unicode{x2014}$which is nothing but the Dirac equation$\unicode{x2014}$even when far away from that limit, a situation which must be considered because of experimental limitations in quantum simulation: the Dirac QCA is better for ultrarelativistic regimes, whereas continuous-time LGT is better for non-relativistic regimes. Then, we compute the GF of the FD-fixed model put forward in the last cited reference, called Flavored Dirac QCA (FDQCA)$\unicode{x2014}$which staggers an extra, artificial flavor onlyonlyonly, on a diamond spacetime lattice, and does not stagger chirality as staggered fermions in usual LGT. The structure of this FDQCA two-point correlation function is of extreme simplicity, and can be expressed in a very simple manner in terms of the four chiral components of the FD-suffering, original-model GF.
A time-dependent optical lattice with staggered particle current in the tight-binding regime was considered that can be described by a time-independent effective lattice model with an artificial staggered magnetic field. The low-energy description of a single-component fermion in this lattice at half-filling is provided by two copies of ideal two-dimensional massless Dirac fermions. The Dirac cones are generally anisotropic and can be tuned by the external staggered flux $\ensuremath{\phi}$. For bosons, the staggered flux modifies the single-particle spectrum such that in the weak coupling limit, depending on the flux $\ensuremath{\phi}$, distinct superfluid phases are realized. Their properties are discussed, the nature of the phase transitions between them is established, and Bogoliubov theory is used to determine their excitation spectra. Then the generalized superfluid--Mott-insulator transition is studied in the presence of the staggered flux and the complete phase diagram is established. Finally, the momentum distribution of the distinct superfluid phases is obtained, which provides a clear experimental signature of each phase in ballistic expansion experiments.
Two finite difference schemes for the numerical treatment of the von Neumann equation for the (2+1)D Dirac Hamiltonian are presented. Both utilize a single-cone staggered space-time grid which ensures a single-cone energy dispersion to formulate a numerical treatment of the mixed-state dynamics within the von Neumann equation. The first scheme executes the time-derivative according to the product rule for "bra" and "ket" indices of the density operator. It therefore directly inherits all the favorable properties of the difference scheme for the pure-state Dirac equation and conserves positivity. The second scheme proposed here performs the time-derivative in one sweep. This direct scheme is investigated regarding stability and convergence. Both schemes are tested numerically for elementary simulations using parameters which pertain to topological insulator surface states. Application of the schemes to a Dirac Lindblad equation and real-space-time Green's function formulations are discussed.
It is known that a one-dimensional (1D) quantum walk gives a local space-time discretization of the massless Dirac equation with a single quasi-energy cone (no fermion doubling at low energies), keeping the fundamental symmetries (chiral and time-reversal) of the continuum theory. We show that the analogous 2D construction is fundamentally more fragile. Local two-band quantum walks can have an unpaired Dirac cone, but the protecting symmetries then cease to be ordinary on-site symmetries: They become non-symmorphic, involving half-lattice translations, and are broken by generic spatial inhomogeneities. In particular, we demonstrate that the 2D Dirac quantum walk based on the Ho-Chalker network model can be gapped by potential scattering, as is evidenced by the breakdown of reflectionless Klein tunneling. The symmorphic symmetry obstruction is sharp: If we allow for couplings that are no longer strictly finite-range but decay exponentially with distance, we can construct a 2D quantum walk with a single Dirac cone and on-site chiral and time-reversal symmetries.
We develop a Hamiltonian formalism for simulating interacting chiral fermions on the lattice while preserving unitarity and locality and without breaking the chiral symmetry. The fermion doubling problem is circumvented by constructing a Fock space endowed with a semi-definite norm. When projecting our theory on the the single-particle sector, we recover the framework of Stacey fermions, and we demonstrate that the scaling limit of the free model recovers the chiral fermion field. Technically, we make use of a matrix product operator norm to mimick the boundary of a higher dimensional topological theory. As a proof of principle, we consider a single Weyl fermion on a periodic ring with Hubbard-type nearest-neighbor interactions and construct a variational generalized DMRG code to demonstrate that the ground state for large system sizes can be determined efficiently. As our tensor network approach does not exhibit any sign problem, we can add a chemical potential and study real-time evolution.