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.
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.
The Dirac quantum walk (a 1+1 dimensional space-time discretization of the Dirac equation) has a 2μ mass gap both at the center and at the corner of the quasi-energy-momentum Brillouin zone. Gupta and Short recently noted [Quantum 9, 1845 (2025)] that the Dirac vacuum can create a particle-hole pair at the zone corner with the release of an energy 2μ to the environment. They removed this vacuum instability at the expense of fermion doubling, the appearance of a second low-energy Dirac cone. Here we show that an alternative discretization scheme, with a tangent rather than a sine dispersion relation, offers stability while retaining a single Dirac cone. The key step is the Cayley transformation from the unit circle of Floquet eigenvalues e−iε to the real line of unbounded energies E=2tan(ε/2). We compute the Schwinger effect (particle-hole pair creation in a uniform electric field) for tangent fermions and show that the pair-production rate agrees with the continuum result for a single Dirac cone. The zone corner is exactly decoupled if the scalar potential is coupled through the Hermitian generator of the quantum map. If it is coupled as a split operator, in order to preserve exact gauge invariance on the lattice, the corner does contribute - with a weight that vanishes quadratically with the lattice constants, in contrast to the Dirac quantum walk where the zone-corner instability survives the continuum limit.
We construct a bosonic analog of a two-dimensional topological Dirac semimetal (DSM). The low-energy description of the most basic 2D DSM model consists of two Dirac cones at positions $\ifmmode\pm\else\textpm\fi{}{\mathbf{k}}_{0}$ in momentum space. The local stability of the Dirac cones is guaranteed by a composite symmetry Z2TI, where T is time reversal and I is inversion. This model also exhibits interesting time-reversal and inversion symmetry breaking electromagnetic responses. In this work we construct a bosonic version by replacing each Dirac cone with a copy of the O(4) nonlinear sigma model (NLSM) with topological theta term and theta angle $\ensuremath{\theta}=\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi}$. One copy of this NLSM also describes the gapless surface termination of the 3D bosonic topological insulator (BTI). We compute the time-reversal and inversion symmetry breaking electromagnetic responses for our model and show that they are twice the value one gets in the DSM case matching what one might expect from, for example, a bosonic Chern insulator. We also investigate the stability of the BSM model and find that the composite Z2TI symmetry again plays an important role. Along the way we clarify many aspects of the surface theory of the BTI including the electromagnetic response, the charges and statistics of vortex excitations, and the stability to symmetry-allowed perturbations. We briefly comment on the relation between the various descriptions of the O(4) NLSM with $\ensuremath{\theta}=\ensuremath{\pi}$ used in this paper (a dual vortex description and a description in terms of four massless fermions) and the recently proposed dual description of the BTI surface in terms of 2+1-dimensional quantum electrodynamics with two flavors of fermion (N=2QED3). In a set of four Appendixes we review some of the tools used in the paper and also derive some of the more technical results.
We consider discrete spacetime models known as quantum walks, which can be used to simulate Dirac particles. In particular we look at fermion doubling in these models, in which high momentum states yield additional low energy solutions which behave like Dirac particles. The presence of doublers carries over to the `second quantised' version of the walks represented by quantum cellular automata, which may lead to spurious solutions when introducing interactions. Moreover, we also consider pseudo-doublers, which have high energy but behave like low energy Dirac particles, and cause potential problems regarding the stability of the vacuum. To address these issues, we propose a family of quantum walks, that are free of these doublers and pseudo-doublers, but still simulate the Dirac equation in the continuum limit. However, there remain a small number of additional low energy solutions which do not directly correspond to Dirac particles. While the conventional Dirac walk always has a zero probability for the walker staying at the same point after a single time step, we obtain the family of walks by allowing this probability to be non-zero.
The Luttinger model is a paradigm for the breakdown due to interactions of the Fermi liquid description of one-dimensional massless Dirac fermions. Attempts to discretize the model on a one-dimensional lattice have failed to reproduce the established bosonization results because of the fermion-doubling obstruction: a local and symmetry-preserving discretization of the Hamiltonian introduces a spurious second species of low-energy excitations, while a nonlocal discretization opens a single-particle gap at the Dirac point. Here, we show how to work around this obstruction by discretizing both space and time to obtain a local Lagrangian for a helical Luttinger liquid with Hubbard interaction. The approach enables quantum Monte Carlo simulations that preserve the topological protection of an unpaired Dirac cone.
Discrete-time quantum walks (QWs) are transportation models of single quantum particles over a lattice. Their evolution is driven through causal and local unitary operators. QWs are a powerful tool for quantum simulation of fundamental physics, as some of them have a continuum limit converging to well-known physics partial differential equations, such as the Dirac or the Schr\"odinger equation. In this paper, we show how to recover the Dirac equation in (3+1) dimensions with a QW evolving in a tetrahedral space. This paves the way to simulate the Dirac equation on a curved space-time. This also suggests an ordered scheme for propagating matter over a spin network, of interest in loop quantum gravity, where matter propagation has remained an open problem.
Abstract Two models are first presented, of a one-dimensional discrete-time quantum walk (DTQW) with temporal noise on the internal degree of freedom (i.e., the coin): (i) a model with both a coin-flip and a phase-flip channel, and (ii) a model with random coin unitaries. It is then shown that both these models admit a common limit in the spacetime continuum, namely, a Lindblad equation with Dirac-fermion Hamiltonian part and, as Lindblad jumps, a chirality flip and a chirality-dependent phase flip, which are two of the three standard error channels for a two-level quantum system. This, as one may call it, Dirac Lindblad equation, provides a model of quantum relativistic spatial diffusion, which is evidenced both analytically and numerically. This model of spatial diffusion has the intriguing specificity of making sense only with original unitary models which are relativistic in the sense that they have chirality, on which the noise is introduced: the diffusion arises via the by-construction (quantum) coupling of chirality to the position. For a particle with vanishing mass, the model of quantum relativistic diffusion introduced in the present work, reduces to the well-known telegraph equation, which yields propagation at short times, diffusion at long times, and exhibits no quantumness. Finally, the results are extended to temporal noises which depend smoothly on position.