# A discrete relativistic spacetime formalism for 1 + 1-QED with continuum limits

Kevissen Sellapillay, Pablo Arrighi, Di Molfetta Giuseppe

► To cite this version:

Kevissen Sellapillay, Pablo Arrighi, Di Molfetta Giuseppe. A discrete relativistic spacetime formalism for 1 + 1-QED with continuum limits. Scientific Reports, 2022, 12 (1), pp.2198. 10.1038/s41598-022-06241-4. hal-03594721

HAL Id: hal-03594721

https://amu.hal.science/hal-03594721

Submitted on 28 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# A discrete relativistic spacetime formalism for 1 + 1-QED with continuum limits

Kevissen Sellapillay<fnref n="1" />,⁵✉, Pablo Arrighi²,³,⁵ & Giuseppe Di Molfetta⁴,⁵✉

We build a quantum cellular automaton (QCA) which coincides with 1 + 1 QED on its known continuum limits. It consists in a circuit of unitary gates driving the evolution of particles on a one dimensional lattice, and having them interact with the gauge field on the links. The particles are massive fermions, and the evolution is exactly U(1) gauge-invariant. We show that, in the continuous-time discrete-space limit, the QCA converges to the Kogut-Susskind staggered version of 1 + 1 QED. We also show that, in the continuous spacetime limit and in the free one particle sector, it converges to the Dirac equation—a strong indication that the model remains accurate in the relativistic regime.

Quantum physical phenomena can always be modelled classically by means of matrices and vectors. But, as far as we know, the dimension of these vectors grows exponentially with the number of particles, making these models intractable for classical computers. To simulate quantum physical phenomena efficiently, it seems we have no choice but to harness the laws of quantum mechanics themselves, as Feynman first suggested<fnref n="1" />. Quantum simulation could be applied to better understand condensed matter problems², simulate molecules, find ground states of Hamiltonians, or even simulate the dynamics of quantum field theories (QFT)³⁻⁵. It is the latter application that motivates this paper.

Amongst QFT, gauge theories are of fundamental importance to Physics, as they capture the fundamental interactions. Some of them have been recast in discrete space. Lattice QCD⁶ is the most famous example as it is extensively used to obtain theoretical numerical values, to be compared against experimental values coming out of particle accelerators: this procedure is partly how physicists are searching for new physics. Simulation has therefore taken a central role in the scientific method of particle physics. But these techniques are computationally heavy: finding a way to simulate lattice gauge theories efficiently and accurately with a quantum simulation device would be a game changer. Lattice gauge theories are also key for condensed matter through their application in spin liquids, and for quantum error correction e.g. via Kitaev's toric code⁷,⁸.

The 1 + 1 QED, also known as the Schwinger model⁹, is a good candidate for a first step towards the quantum simulation of the dynamics of a gauge theory. Indeed, it is based on the U(1) gauge group just like 3 + 1 QED. It captures many non trivial physical properties such as a mass gap, fermion confinement and chiral symmetry breaking¹⁰. It is exactly solvable in the massless limit⁹. These features explain why it is often used as a testbed for new techniques and ideas.

The standard ways to quantum simulate QFT are fundamentally non-relativistic, as they all begin by expressing the theory in continuous-time discrete-space Hamiltonian form, using Kogut-Susskind methods¹¹,¹². They then map the matter (fermions) and the gauge field (bosons) onto quantum systems on a lattice, whose interactions will mimic those of the target Hamiltonian³. Sometimes these interactions are implemented as discrete-time products of quantum gates, but even then these are obtained by approximating the target Hamiltonian via the Trotter formula, an approximation which remains valid only in the non-relativistic regime $\Delta_t \ll \Delta_x$. This approach was recently realized experimentally on an ion trap architecture¹³. Numerical techniques exist that come to complement the standard approach, based on tensors networks. Those use compact, approximate descriptions of quantum states¹⁴ such as the Density Matrix Renormalization Group (DMRG)¹⁵,¹⁶, discarding unwanted information about the states as they evolve, so that their description remain of manageable size—whilst attempting to keep track of the interesting physical ingredients.

In order to achieve quantum simulation in the relativistic regime $\Delta_x \approx \Delta_t$, one must keep space and time on an equal footing, discretizing both at the same time. This is sometimes referred to as digital quantum simulation.

<footnote n="1">¹Aix-Marseille Université, CPT, Campus de Luminy, case 907, 13288 Marseille, France. ²Université Paris-Saclay, CNRS, Laboratoire de recherche en informatique, 91405 Orsay, France. ³IXXI, Lyon, France. ⁴Aix-Marseille Université, Université de Toulon, CNRS, LIS, 13288 Marseille, France. ⁵These authors contributed equally: Kevissen Sellapillay, Pablo Arrighi and Giuseppe Di Molfetta. ✉email: kevissen.sellapillay@univ-amu.fr; giuseppe.dimolfetta@lis-lab.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Kogut–Susskind version of the Schwinger model. The gauge (boson) field is represented by the links in red. These links are states $|l\rangle$ where $l$ takes value in $\mathbb{Z}$. The operators acting on them are $L_p$ and $e^{\pm i\theta_p}$. The matter (fermion) field $\phi_p, \phi_p^\dagger$ are on the nodes in black. Even (odd) sites correspond to upper (lower) component of a spinor field. $a$ is the lattice spacing.</figcaption></figure>

Digital simulation has indeed been very successful at describing relativistic particles in different fields¹⁷, but so far it has not been able to produce simulation scheme for interacting QFT in the sector of more than 2 particles.

A proposal was made in Ref.¹⁸. From a QCA simulating relativistic Dirac equation a $U(1)$ gauge-invariant model was obtained. The convergence to the Schwinger model was not shown however and in fact, using the method presented in this paper, it can be shown that the Hamiltonian in the discrete space continuous time limit is not the correct one. We present here the circuit that gives the correct Hamiltonian in this limit.

In Refs.¹⁹,²⁰ a quantum walk (i.e. the one particle sector of a QCA) was proposed which unifies non-relativistic analog quantum simulation with relativistic digital quantum simulation. Just by imposing $\Delta_x = \Delta_t^{1-\alpha}$ and tuning the $\alpha$, the operator is found to have well defined limits lattice fermions both in continuous-time discrete-space, and the relativistic Dirac equation in the continuous spacetime limit—a property referred to as plasticity.

The QCA presented in this paper is closely related to these last two models. It is again based upon a QCA that recovers the relativistic Dirac equation, extended to become natively discrete gauge-invariant as in Ref.¹⁸. But this time, the QCA is plastic, allowing us to prove its continuum limit towards 1 + 1 QED, in the regime where 1 + 1 QED does have a limit, i.e. the non-relativistic regime. In other words, we recover the Hamiltonian of the Kogut–Susskind Schwinger model in the continuous-time discrete-space limit. In the continuous spacetime limit we show that the QCA yields the Dirac equation in the free one particle sector, allowing to make the bridge between the non-relativistic and the relativistic regimes. Altogether, the QCA coincides with 1 + 1 QED on its mathematical continuum limits, whenever these are defined.

The natively discrete digital circuit for staggered Schwinger model we propose is not seen in the literature. The QCA is staggered which is not usual in the QCA formalism²¹. One may wonder whether the approach, beyond the quantum simulation application, could be used to reframe QFT. Indeed, the fact that the QCA is gauge-invariant by construction, contains explicit relativistic and non-relativistic limits, is expressible by means of path integrals²²,²³, suggests that Quantum Computing point of view upon QFT may bring both rigour and pedagogy to the table—reviving the line of thought initiated by Feynman with his checkerboard propagator for 1 + 1 Dirac equation²⁴.

The paper is organized as follows. We first define the QCA model, that is the spacetime structure and the gates. Second we show the continuous-time and discrete-space limit towards the Kogut–Susskind version of the Schwinger model, by means of the Jordan–Wigner transformation from qubits to fermions. Third we show that in a continuous spacetime limit, we recover the Dirac equation for the free one particle sector. Finally, we prove that the model is gauge-invariant and conclude by giving some perspectives.

## Model

**The Kogut–Susskind staggered version of the Schwinger model.** The Schwinger model⁹ is a (1 + 1)D model invariant under the $U(1)$ gauge group. It models spinless electrons and their antiparticles, positrons, propagating on a 1D lattice and interacting with a $U(1)$ gauge field. We briefly summarize it by giving its Hamiltonian, which can be written using a temporal gauge ($A_0(x) = 0$, and $A(x) = A_1(x)$) as :

$$H = \int dx \left( \psi^{\dagger}(x) \left[ (i\partial_x + igA(x))\sigma_z + m\sigma_x \right] \psi(x) + \frac{1}{2}E^2(x) \right), \quad (1)$$

where $E(x)$ is the electric field observable at $x$, and $A(x)$ is its conjugate momentum, meaning

$$[A(x), E(y)] = i\delta(x - y). \qquad (2)$$

Here $\psi(x) = (\psi_1(x), \psi_2(x))^T$ is a two components fermion field satisfying

$$\{\psi_{\alpha}(x), \psi_{\beta}^{\dagger}(y)\} = \delta_{\alpha\beta}\delta(x - y). \qquad (3)$$

We now describe the staggered Kogut–Susskind version of the Schwinger model. This Kogut–Susskind procedure consists in putting fermion fields on the nodes of an infinite 1D lattice and bosonic gauge fields on the links between them, as depicted in Fig. 1. We can interpret occupied odd sites as electrons and unoccupied even sites as positrons, therefore particles and antiparticles are described by a single fermion field¹³ (staggered picture). This procedure allows to give a continuous time, discrete space formulation of a continuous spacetime model and is also a way to partly resolve the fermion doubling problem.

The Hamiltonian of this staggered version reads¹¹,¹²:
<figure><img src="image_4.png" /><figcaption>**Figure 2.** QCA structure. One application of G, the evolution operator, corresponds to two rows, first a row of W and then of W* gates. Black wires represent the fermionic fields, green wires represent the gauge field. The qubits are separated in space by Δx and the gates are separated in time by Δt.</figcaption></figure>

$$H_S = \frac{i}{2a} \sum_p (\phi_{p+1}^\dagger e^{-i\theta_p} \phi_p - h.c.) + m \sum_p (-1)^p \phi_p^\dagger \phi_p + \frac{ag^2}{2} \sum_p L_p^2, \quad (4)$$

where a is the lattice spacing, g the strength of the interaction (the charge of the particles), m the mass and φ the fermion field. To construct the spinor Ψ$_p$ from the fermion fields, we group fields in pairs (p, p + 1), where p is even. In between the fermion fields at p and p + 1, there is a gauge field link which takes values in $\mathbb{Z}$. The operators e$^{\pm i\theta_p}$ raises or lowers the value of this link [p, p + 1] such that :

$$e^{\pm i\theta_p} |l\rangle_p = |l \pm 1\rangle_p \quad (5)$$

and the electric field is given by $E_p = gL_p$ where

$$L_p |l\rangle_p = l |l\rangle_p. \quad (6)$$

**A quantum cellular automaton for 1 + 1 QED.** The model we propose consists in having one qubit per site p separated by Δx and gauge fields located on the links between each sites, at half step, modeled by states taking values in $\mathbb{Z}^{18}$. This gauge field could be experimentally represented by qudits or harmonic oscillators. We choose the evolution operator of the QCA to be :

$$G = \bigotimes_{p' \text{ even}} W_p^* e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2} \otimes_{p \text{ odd}} W_p e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2}, \quad (7)$$

where

$$W_p = \begin{pmatrix} I & 0 & 0 & 0 \\ 0 & e^{-i\zeta} \sin \theta I & \cos \theta V_{p+\frac{1}{2}} & 0 \\ 0 & -\cos \theta V_{p+\frac{1}{2}}^{\dagger} & e^{i\zeta} \sin \theta I & 0 \\ 0 & 0 & 0 & I \end{pmatrix} \quad (8)$$

corresponds to the fermions dynamics and the exponential applied before each $W_p$ codes the interaction with the bosons field at the same position. In particular, as depicted in Fig. 2, the gate $W_p$ and the gate $W_p^*$ are located in the space-time grid in between the qubit at position (p, p + 1) respectively for odd and even p. Each gate W and its conjugate acts on the local gauge field. Moreover, in order to avoid the fermion doubling problem we choose to work with a staggered QCA, in which occupied odd sites are interpreted as electrons and unoccupied even sites as positrons. This is consistent with the same interpretation that Kogut and Susskind gave of the Schwinger model. From a dynamical point of view, the matrix $W_p$ can be interpreted as follows: diagonal terms correspond to staying on the same site, eventually picking up a phase related to the mass. The non diagonal part correspond to hopping terms: a right-moving |1⟩ will decrease the gauge field it passes through, a left-moving one will increase it. The gauge field operators read:

$$V_{p+\frac{1}{2}} |l\rangle_{p+\frac{1}{2}} = |l-1\rangle_{p+\frac{1}{2}}, \quad (9)$$
and

$$V_{p+\frac{1}{2}}^{\dagger}|l\rangle_{p+\frac{1}{2}} = |l+1\rangle_{p+\frac{1}{2}}. \qquad (10)$$

The identities $\mathbb{I}$ means that the gauge link is left invariant. The operator $L$ in the exponentials of Eq. (7) acts on the states of the gauge field $|l\rangle_{p+\frac{1}{2}}$ in the following way:

$$L|l\rangle_{p+\frac{1}{2}} = l|l\rangle_{p+\frac{1}{2}}, \qquad (11)$$

for $l \in \mathbb{Z}$.

Now, in order to map the above automaton with the Schwinger model we move to the second quantization formalism. Let us start with the gate $W_p$ (8), which acts on a pair of qubits $(p, p+1)$ and a gauge link $|p, p+1\rangle$. It is useful to introduce $\tilde{W}_p = I_{<p} \otimes W_p \otimes I_{>p}$. Each $W_p$ can be written in terms of single qubit operators $E_{ij} = |i\rangle\langle j|$. Using the order $(0_p 0_{p+1}, 0_p 1_{p+1}, 1_p 0_{p+1}, 1_p 1_{p+1})$ and omitting identities on the gauge link for clarity, we have:

$$W_p = (E_{00})_p \otimes (E_{00})_{p+1} + e^{-i\zeta} \sin \theta (E_{00})_p \otimes (E_{11})_{p+1} \\ + e^{i\zeta} \sin \theta (E_{11})_p \otimes (E_{00})_{p+1} + \cos \theta V_{p+\frac{1}{2}} (E_{01})_p \otimes (E_{10})_{p+1} \\ - \cos \theta V_{p+\frac{1}{2}}^{\dagger} (E_{10})_p \otimes (E_{01})_{p+1} + (E_{11})_p \otimes (E_{11})_{p+1}. \qquad (12)$$

Thus, we ought to transform qubit operators to fermionic operators, via the standard Jordan-Wigner transformations, and finally introduce the following annihilation and creation operators:

$$\phi_p = I_{<p} \otimes (E_{01})_p \otimes (\sigma_z)_{p+1} \otimes (\sigma_z)_{p+2} \dots \\ \phi_p^{\dagger} = I_{<p} \otimes (E_{10})_p \otimes (\sigma_z)_{p+1} \otimes (\sigma_z)_{p+2} \dots, \qquad (13)$$

with $\{\phi_p, \phi_p^\dagger\} = \delta_{p,p'}I$ and $\{\phi_p, \phi_{p'}^\dagger\} = 0$. Notice that the $\sigma_z$ to the right are specifically there to ensure the correct anti-commutation relation (Supplementary Information 1), together they form the string operator of the Jordan-Wigner transformation. Looking at these operators, we see that creating or annihilating a fermion on a qubit lattice is very much a non-local operation and we could be worried that the resulting dynamics using such operators would be non-local and unphysical. However, the terms involving these operators in the Hamiltonian are quadratic, therefore they become perfectly local and physical. This can already be seen at the level of the gates written in this formalism. In fact, for $\tilde{W}_p$ we have:

$$\tilde{W}_p = \phi_p \phi_p^\dagger \phi_{p+1} \phi_{p+1}^\dagger + e^{-i\zeta} \sin \theta \phi_p \phi_p^\dagger \phi_{p+1}^\dagger \phi_{p+1} + e^{i\zeta} \sin \theta \phi_p^\dagger \phi_p \phi_{p+1} \phi_{p+1}^\dagger \\ - \cos \theta V_{p+\frac{1}{2}} \phi_p \phi_{p+1}^\dagger - \cos \theta V_{p+\frac{1}{2}}^\dagger \phi_p^\dagger \phi_p \phi_{p+1} + \phi_p^\dagger \phi_p \phi_{p+1}^\dagger \phi_{p+1}. \qquad (14)$$

Putting together these transformed gates, the global evolution reads:

$$G = \prod_{p' \text{ even}} \tilde{W}_{p'}^* e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2} \prod_{p \text{ odd}} \tilde{W}_p e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2}. \qquad (15)$$

## Methods

**Continuous limits.** In order to prove that the above staggered QCA, reformulated in terms of fermionic operators, converges to the Schwinger Hamiltonian, we introduce the following parametrization:

$$\begin{aligned} \Delta t &= \epsilon \\ \Delta x &= \epsilon^{1-\alpha} \\ \kappa &= \epsilon^{\alpha} \\ \theta &= \arccos(c\kappa) \\ \zeta &= m \frac{(-1)^{\kappa} \epsilon}{\sin(\theta)}, \end{aligned} \qquad (16)$$

where the case $\alpha = 1$ and $\alpha = 0$ correspond respectively to the continuous time and discrete space limit and the continuous spacetime limit.

In the non relativistic limit, the evolution is continuous in time and discrete in space. The dynamics is driven by a Hamiltonian, which is recovered looking at the first order of the global evolution operator of the QCA, as follows:

$$G = e^{-2iH_{QCA}\Delta t} \approx 1 - 2i\Delta t H_{QCA}. \qquad (17)$$

We then take the limit $\epsilon \to 0$ on $G$, using the parametrisation (16), for $\alpha = 1$. Finally we get:
$$ 
\begin{aligned}
\tilde{W}_p & \approx \phi_p \phi_p^\dagger \phi_{p+1} \phi_{p+1}^\dagger + (1 - im\epsilon) \phi_p \phi_p^\dagger \phi_{p+1}^\dagger \phi_{p+1} + (1 + im\epsilon) \phi_p^\dagger \phi_p \phi_{p+1} \phi_{p+1}^\dagger \\
& - \epsilon V_{p+\frac{1}{2}} \phi_p \phi_{p+1}^\dagger - \epsilon V_{p+\frac{1}{2}}^\dagger \phi_p^\dagger \phi_{p+1} + \phi_p^\dagger \phi_p \phi_{p+1}^\dagger \phi_{p+1} \\
& = 1 + \epsilon \left[ \phi_{p+1}^\dagger \phi_p V_{p+\frac{1}{2}} - V_{p+\frac{1}{2}}^\dagger \phi_p^\dagger \phi_{p+1} + im(\phi_p^\dagger \phi_p \phi_{p+1} \phi_{p+1}^\dagger - \phi_p \phi_p^\dagger \phi_{p+1}^\dagger \phi_{p+1}) \right].
\end{aligned}
\quad (18) $$

We can rewrite the mass term:

$$ 
\begin{aligned}
\phi_p^\dagger \phi_p \phi_{p+1} \phi_{p+1}^\dagger - \phi_p \phi_p^\dagger \phi_{p+1}^\dagger \phi_{p+1} = \phi_p^\dagger \phi_p \phi_{p+1} \phi_{p+1}^\dagger - (1 - \phi_p^\dagger \phi_p) \phi_{p+1}^\dagger \phi_{p+1} \\
\phi_p^\dagger \phi_p (\phi_{p+1} \phi_{p+1}^\dagger + \phi_{p+1}^\dagger \phi_{p+1}) - \phi_{p+1}^\dagger \phi_{p+1} = \phi_p^\dagger \phi_p - \phi_{p+1}^\dagger \phi_{p+1}.
\end{aligned}
$$

Finally we have

$$ \tilde{W}_p \approx 1 + \epsilon \left[ \phi_{p+1}^{\dagger} \phi_p V_{p+\frac{1}{2}} - V_{p+\frac{1}{2}}^{\dagger} \phi_p^{\dagger} \phi_{p+1} + im(\phi_p^{\dagger} \phi_p - \phi_{p+1}^{\dagger} \phi_{p+1}) \right], \quad (19) $$

and the interaction with the gauge field can be developed as

$$ e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2} \approx 1 - \frac{i}{2} \epsilon g^2 \sum_i L_i^2. \quad (20) $$

Combining everything, we get for G:

$$ 
\begin{aligned}
G &= \prod_{p' \text{ even}} \tilde{W}_p^* e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2} \prod_{p \text{ odd}} \tilde{W}_p e^{-\frac{i}{2} \Delta x \Delta t g^2 L^2} \\
&= \prod_{p'} (1 + \epsilon [\phi_{2p'}^\dagger \phi_{2p'} V_{2p'+\frac{1}{2}} - \phi_{2p'}^\dagger \phi_{2p'+1} V_{2p'+\frac{1}{2}} - im(\phi_{2p'}^\dagger \phi_{2p'} - \phi_{2p'+1}^\dagger \phi_{2p'+1})]) \\
&\quad (1 - \frac{i}{2} \epsilon g^2 \sum_i L_i^2) \prod_p (1 + \epsilon [\phi_{2p+2}^\dagger \phi_{2p+1} V_{2p+\frac{3}{2}} - \phi_{2p+1}^\dagger \phi_{2p+2} V_{2p+\frac{3}{2}} \\
&\quad + im(\phi_{2p+1}^\dagger \phi_{2p+1} - \phi_{2p+2}^\dagger \phi_{2p+2})]) (1 - \frac{i}{2} \epsilon g^2 \sum_i L_i^2).
\end{aligned}
\quad (21) $$

A straightforward calculation leads us to the leading order of the series:

$$ 
\begin{aligned}
G \approx 1 + \epsilon \sum_p [\phi_{2p+1}^\dagger \phi_{2p} V_{2p+\frac{1}{2}} - \phi_{2p}^\dagger \phi_{2p+1} V_{2p+\frac{1}{2}} + \phi_{2p+2}^\dagger \phi_{2p+1} V_{2p+\frac{3}{2}} - \phi_{2p+1}^\dagger \phi_{2p+2} V_{2p+\frac{3}{2}} \\
- ig^2 L_p^2 - im(\phi_{2p}^\dagger \phi_{2p} - 2\phi_{2p+1}^\dagger \phi_{2p+1} + \phi_{2p+2}^\dagger \phi_{2p+2})].
\end{aligned}
\quad (22) $$

Notice that the mass term can be rewritten as: $-2im(\phi_{2p}^\dagger \phi_{2p} - \phi_{2p+1}^\dagger \phi_{2p+1})$ Moreover the four hopping terms simplify because each pair of two terms is separated by one step:

$$ \sum_p \phi_{p+1}^\dagger \phi_p V_{p+\frac{1}{2}} - \text{h.c.} \quad (23) $$

Finally, we identify $V_{p+\frac{1}{2}}$ to $e^{-i\theta_p}$ and $V_{p+\frac{1}{2}}^\dagger$ to $e^{i\theta_p}$. Identifying the Hamiltonian in (22) using $G \approx 1 - 2i\epsilon H_{QCA}$, we find the Hamiltonian of the QCA to be:

$$ H_{QCA} = \sum_p \left[ \frac{i}{2} (\phi_{p+1}^\dagger \phi_p e^{-i\theta_p} - \text{h.c.}) + m(-1)^p \phi_p^\dagger \phi_p + \frac{g^2}{2} L_p^2 \right]. \quad (24) $$

The above one coincides with the Kogut-Susskind Hamiltonian of the Schwinger model $H_S$ with $a = 1$. We have thus identified a QCA-based quantum simulator for a QED toy model, namely a theory for both spinless electrons and positrons and their interaction with a dynamical gauge field.

We ought to be sure that in the relativistic limit our simulator reproduces the right dynamics. Here, we give a proof of that in the simplest scenario, the non-interacting case. Starting from the the one particle sector of the staggered QCA, we take the $\alpha = 0$ relativistic limit, namely for $\Delta t = \Delta x = \epsilon \to 0$ (16). The local gate which drives the automaton simplifies as follow

$$ W'_p = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & e^{-i\zeta} \sin\theta & \cos\theta & 0 \\ 0 & -\cos\theta & e^{i\zeta} \sin\theta & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix}, \quad (25) $$

and the global evolution operator reads:
$$G = \bigotimes_{p' \text{ even}} W_{p'}^* \otimes_{p \text{ odd}} W_p. \qquad (26)$$

Let's start with a general 1 particle state: $|\Psi(t)\rangle = \sum_x \psi(t, x)|1\rangle_x$ and separate the one particle state into pairs of even and odd sites:

$$|\Psi(t)\rangle = \sum_{x \text{ even}} \psi^l(t, x)|1\rangle_x + \psi^r(t, x + \epsilon)|1\rangle_{x+\epsilon}. \qquad (27)$$

After one time step evolution of the automaton, the recurrence relations on the amplitudes $s \equiv (\psi^l, \psi^r)$ which governs the single fermion, reads:

$$\mathbf{s}(t + 2\epsilon, x) = \begin{pmatrix} 0 & -\cos\theta e^{-i\zeta} \sin\theta \\ 0 & \cos^2\theta \end{pmatrix} \mathbf{s}(t, x - 2\epsilon) + \begin{pmatrix} e^{-2i\zeta} \sin^2\theta & \cos\theta e^{i\zeta} \sin\theta \\ -\cos\theta e^{-i\zeta} \sin\theta & e^{2i\zeta} \sin^2\theta \end{pmatrix} \mathbf{s}(t, x) \qquad (28)$$

$$+ \begin{pmatrix} \cos^2\theta & 0 \\ \cos\theta e^{i\zeta} \sin\theta & 0 \end{pmatrix} \mathbf{s}(t, x + 2\epsilon). \qquad (29)$$

Taking the limit $\epsilon \to 0$, we find the following differential equation:

$$\partial_t \mathbf{s}(t, x) = P \partial_x \mathbf{s}(t, x) + Q \mathbf{s}(t, x), \qquad (30)$$

where the operators $P$ and $Q$ are represented in the computational basis:

$$P = \begin{pmatrix} c^2 & c\sqrt{1-c^2} \\ c\sqrt{1-c^2} & -c^2 \end{pmatrix} \qquad (31)$$

and

$$Q = \begin{pmatrix} im\sqrt{1-c^2} & -cim \\ -cim & -im\sqrt{1-c^2} \end{pmatrix}. \qquad (32)$$

The operator $P$ is self-adjoint and its eigenvalues are $\pm c$. Two eigenvectors associated to these eigenvalues are:

$$b_{-} = -\sqrt{\frac{1-c}{2}} b_{0} + \sqrt{\frac{1+c}{2}} b_{1}, \quad b_{+} = \sqrt{\frac{1+c}{2}} b_{0} + \sqrt{\frac{1-c}{2}} b_{1}. \qquad (33)$$

The family $(b_-, b_+)$ forms an orthonormal basis of the two dimensional spin Hilbert space. Let us now rewrite Eq. (30) in this new orthonormal basis. A straightforward computation leads to:

$$i\gamma^0 \partial_0 \bar{s}(t, x) + i\gamma^1 \partial_1 \bar{s}(t, x) - m\bar{s}(t, x) = 0, \qquad (34)$$

where $\gamma^0 = \sigma_x$, $\partial_0 = \partial_t$, $\gamma^1 = \sigma_x \sigma_z$ and $\partial_0 = c\partial_x$.

**Remark.** Let us shortly discuss the last term of operator $W(p)$. In general, when two creation operators get exchanged, a minus sign is produced so as to respect the anti-commutation of their creation operators. The abrupt minus appearance of minus signs makes it harder to compute quantities about fermions—an issue which is sometimes referred to as the ‘sign problem’. Here, when $\Delta t \approx \Delta x$, we must indeed put a $-1$ in the bottom right coefficient of the gate, since two creation operators are crossing during the lapse of one W gate, as was shown in detail¹⁸. However, when $\Delta t \ll \Delta x$ we should not. Physically, this is because the two fermions now hardly have time to cross in the lapse of one gate. Mathematically, this shows through the fact that placing a $-1$ at this position forbids the development of the gate around identity, ruining any effort to obtain a non-relativistic continuous-time discrete-space limit towards the Kogut–Susskind Hamiltonian, or any other Hamiltonian for that matter. To get the best of both worlds, we use scaling factor $e^{i(\frac{\Delta t}{\Delta x})^2\pi} = e^{i\epsilon^{2\alpha\pi}}$, making the coefficient go to 1 in the non-relativistic parametrization ($\alpha = 1$) and to $-1$ in the relativistic parametrization ($\alpha = 0$). This is compatible with unitarity, plasticity, and fermionic computation.

$$W_p'' = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & \sin\theta & -\cos\theta V_{p+\frac{1}{2}} & 0 \\ 0 & \cos\theta V_{p+\frac{1}{2}}^{\dagger} & \sin\theta & 0 \\ 0 & 0 & 0 & e^{i\epsilon^{2\alpha\pi}} \end{pmatrix}. \qquad (35)$$

It is the choice that yields the Kogut–Susskind Hamiltonian in the non-relativistic regime. In the relativistic regime, however, $\alpha = 0$ and so at order 0, we see that:

$$\begin{aligned} \bar{W}_p & \approx \phi_p \phi_p^\dagger \phi_{p+1} \phi_{p+1}^\dagger + \phi_p \phi_p^\dagger \phi_{p+1}^\dagger \phi_{p+1} + \phi_p^\dagger \phi_p \phi_{p+1} \phi_{p+1}^\dagger - \phi_p^\dagger \phi_p \phi_{p+1}^\dagger \phi_{p+1} \\ &= 1 - 2\phi_p^\dagger \phi_p \phi_{p+1}^\dagger \phi_{p+1}. \end{aligned} \qquad (36)$$

Since this is a control-Z on the 2 qubits, the order of $G$ cannot be the identity. As expected, we cannot recover a many-body interacting Hamiltonian in this regime.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** The gauge transformation. Black dots are the fermions sites. Green dots in between sites represent the gauge field values. At a single site p, the gauge transformation will apply a gate $R_{\phi(p)}$ that gives a phase when the qubit there is in state $|1\rangle$ and apply gates $T_{\phi(p)}$ on the left and right gauge field points that produce phases according to the values of the gauge field. This transformation is applied on every sites, therefore a single gauge field is acted upon twice.</figcaption></figure>

## The staggered QCA is gauge invariant

In this section, we show that our QCA-based quantum simulator is gauge invariant at finite scale. We define the gauge transformation in Fig. 3, where the gates are defined as follows:

$$P_{\phi} = \otimes_{p} P_{\phi(p)} = \otimes_{p} (T_{\phi(p)} \otimes R_{\phi(p)} \otimes T_{-\phi(p)}), \qquad (37)$$

where $R_{\phi(p)}$ acts on the qubit at site p such that:

$$R_{\phi(p)} : |0\rangle \rightarrow |0\rangle \qquad (38)$$

$$|1\rangle \rightarrow e^{i\phi(p)}|1\rangle, \qquad (39)$$

and $T_{\phi(p)}$ acts on the gauge field states such that:

$$T_{\phi(p)}|l\rangle = e^{il\phi(p)}|l\rangle. \qquad (40)$$

From definition (37), we see that a gauge field at a given site $p + \frac{1}{2}$ will be acted upon twice, once by $-\phi(p)$ and once by $+\phi(p + 1)$. Is our QCA invariant under the above gauge transformation, or namely:

$$P_{\phi} G \stackrel{?}{=} G P_{\phi}. \qquad (41)$$

Equation (41) means that whatever the field $\phi(p)$ is, the evolution operator should give the same result if we apply it after the gauge transformation or before the gauge transformation, in other words, the dynamics should not be changed by this gauge field, the latter should have no physical consequences. However if we impose this $U(1)$ local phase, invariance is only possible if we introduce some interaction between qubits and the gauge fields, which is done through the $V$ and $V^\dagger$ operators, that will compensate the missing pieces of phase appearing during the particle movement, as done in$^{18}$. To show the gauge invariance (41), we start from a ket $|1\rangle$ at a given site (odd or even), with only gauge field values given in the rest of the space. We compute explicitly the results after applying first $P_{\phi}$ then $G$ and the other way around, and show that the result are exactly the same, therefore the gauge transformation $P_{\phi}$ has no physical consequences on the dynamics $G$. The model is thus gauge invariant (Supplementary Information 2).

## Conclusion

We have described a quantum cellular automaton (QCA) that simulates $1 + 1$ QED. It consists in a lattice of qubits encoding whether a fermion is present at a given site. These interact with the gauge field that lives on the links between those sites. The QCA was shown to coincide with the Kogut-Susskind version of the Schwinger model in the continuous-time discrete-space limit, and with the Dirac equation in the continuous spacetime limit in the one-particle sector. We go from one limit to the other just by imposing $\Delta_x = \Delta_t^{1-\alpha}$ and tuning the $\alpha$.

We still cannot ascertain the QCA recovers the Schwinger model in the interactive regime of the continuous spacetime limit, but then again it is not even clear that the lattice QFT has such a limit in the first place. The QCA coincides with the Schwinger model wherever it has a known, mathematically defined limit.

It coincides in the story it tells; of fermions propagating relativistically updating the gauge field, which in turn simply triggers a phase—thereby turning on the interaction.

Finally, it coincides in terms of its construction: the QCA retains the fundamental $U(1)$ gauge-invariance of QED, even in discrete spacetime. This gauge-invariance construction was originally proposed in Ref.$^{18}$ for QCA, and in$^{25}$ for Reversible CA. Both drew inspirations from gauge-invariant Quantum Walks$^{17,26,27}$. The question of extending these constructions to non-abelian gauge theories has been treated for Quantum Walks$^{28,29}$ and for Reversible CA$^{30}$. Lifting this to obtain $U(N)$ gauge-invariant QCA is no doubt one of the next steps lying ahead towards digital quantum simulation schemes for Yang Mills theories, such as QCD. An observation was
made in³¹ that the path integral of such a theory can be written as a real-time transfer operator which is itself a finite-depth local quantum circuit.

Another obvious next step is the extension to 2 + 1 and 3 + 1 dimensions. Of particular concern is the fact that we relied upon the Jordan Wigner mapping to encode fermions into qubits. This transformation is known to suffer locality issues³²⁻³⁴, which did not affect us because our terms were quadratic, cancelling out all non-local effects. In further dimensions we may not be so lucky however, and will have to rely on alternative transformations³⁵⁻³⁷. A possible route of investigation could be to use formulations of QED in 2 + 1D or 3 + 1 dimensions where the Hilbert space dimension is reduced by either approximating the gauge group³⁸ or by going in a rotating frame³⁹ which decouples the matter to the gauge field and keeps only local constraints on the latter.

The problem of preparing the ground state of such QCA is puzzling as was pointed out in Ref.¹⁸. We do not solve the problem but by recovering a Hamiltonian in the continuous-time discrete-space limit, we make the problem well-defined.

QCA are closely related to path integrals²². Lately a formalism was developed for dealing with interactions in a perturbative manner²³, within the QCA framework—which would be an interesting application here. Another vast topic for exploration is trying to understand the link between renormalization theory, and the way the parameters of the unitary gates must be made to vary with the lattice as we take our limits.

The plastic Quantum Walk¹⁹ upon which this plastic QCA is built is deformable to the point that a curved spacetime limit can been obtained just by putting a spacetime dependence in the parameter c. A fair question to ask is whether allowing for the same spacetime dependence here, would yield the Schwinger model on a curved background⁴⁰.

Received: 22 September 2021; Accepted: 20 January 2022

Published online: 09 February 2022

## References

1. Feynman, R. P. Simulating physics with computers. *Int. J. Theor. Phys.* **21**, 467–488. https://doi.org/10.1007/BF02650179 (1982).

2. Lloyd, S. Universal quantum simulators. *Science* **273**, 1073–1078. https://doi.org/10.1126/science.273.5278.1073 (1996).

3. Bañuls, M. C. et al. Simulating lattice gauge theories within quantum technologies. *Eur. Phys. J. D* **74**, 165. https://doi.org/10.1140/epjd/e2020-100571-8 (2020).

4. Preskill, J. Simulating quantum field theory with a quantum computer. In *Proceedings of The 36th Annual International Symposium on Lattice Field Theory—PoS (LATTICE 2018)*, 024. https://doi.org/10.22323/1.334.0024 (Sissa Medialab, Michigan State University, 2019).

5. Jordan, S. P., Lee, K. S. & Preskill, J. Quantum algorithms for quantum field theories. *Science* **336**, 1130–1133 (2012).

6. Knechtli, F., Günther, M. & Peardon, M. *Lattice Quantum Chromodynamics: Practical Essentials. SpringerBriefs in Physics* (Springer Netherlands, 2017).

7. Savary, L. & Balents, L. Quantum spin liquids: A review. *Rep. Prog. Phys.* **80**, 016502. https://doi.org/10.1088/0034-4885/80/1/016502 (2016) (**Publisher: IOP Publishing**).

8. Kitaev, A. Y. Fault-tolerant quantum computation by anyons. *Ann. Phys.* **303**, 2–30. https://doi.org/10.1016/S0003-4916(02)00018-0 (2003).

9. Schwinger, J. Gauge invariance and mass. II. *Phys. Rev.* **128**, 2425–2429. https://doi.org/10.1103/PhysRev.128.2425 (1962).

10. Melnikov, K. & Weinstein, M. Lattice Schwinger model: Connement, anomalies, chiral fermions, and all that. *Phys. Rev. D* **14**, 094504 (2000).

11. Kogut, J. & Susskind, L. Hamiltonian formulation of Wilson's lattice gauge theories. *Phys. Rev. D* **11**, 395–408. https://doi.org/10.1103/PhysRevD.11.395 (1975).

12. Banks, T., Susskind, L. & Kogut, J. Strong-coupling calculations of lattice gauge theories: (1 + 1)-dimensional exercises. *Phys. Rev. D* **13**, 1043–1053. https://doi.org/10.1103/PhysRevD.13.1043 (1976).

13. Martinez, E. A. et al. Real-time dynamics of lattice gauge theories with a few-qubit quantum computer. *Nature* **534**, 516–519. https://doi.org/10.1038/nature18318 (2016) (**Number: 7608 Publisher: Nature Publishing Group**).

14. Orús, R. A practical introduction to tensor networks: Matrix product states and projected entangled pair states. *Ann. Phys.* **349**, 117–158. https://doi.org/10.1016/j.aop.2014.06.013 (2014).

15. Byrnes, T. M. R., Sriganesh, P., Bursill, R. J. & Hamer, C. J. Density matrix renormalization group approach to the massive Schwinger model. *Phys. Rev. D* **66**, 013002. https://doi.org/10.1103/PhysRevD.66.013002 (2002).

16. Zapp, K. & Orús, R. Tensor network simulation of QED on infinite lattices: Learning from (1+1)d, and prospects for (2+1)d. *Phys. Rev. D* **95**, 114508. https://doi.org/10.1103/PhysRevD.95.114508 (2017) (**Publisher: American Physical Society**).

17. Di Molfetta, G., Brachet, M. & Debbasch, F. Quantum walks in artificial electric and gravitational fields. *Phys. A* **397**, 157–168 (2014).

18. Arrighi, P., Bény, C. & Farrelly, T. A quantum cellular automaton for one-dimensional QED. *Quantum Inf. Process.* **19**, 88. https://doi.org/10.1007/s11128-019-2555-4 (2020).

19. Di Molfetta, G. & Arrighi, P. A quantum walk with both a continuous-time limit and a continuous-spacetime limit. *Quantum Inf. Process.* **19**, 47. https://doi.org/10.1007/s11128-019-2549-2 (2019).

20. Manighalam, M. & Di Molfetta, G. Continuous time limit of the dtqw in 2d+ 1 and plasticity. *Quantum Inf. Process.* **20**, 1–24 (2021).

21. Arrighi, P. An overview of quantum cellular automata. *Nat. Comput.* **18**, 885–899. https://doi.org/10.1007/s11047-019-09762-6 (2019).

22. D'Ariano, G. M., Mosco, N., Perinotti, P. & Tosini, A. Path-sum solution of the Weyl Quantum Walk in 3 + 1 dimensions. *Philos. Trans. R. Soc. A Math. Phys. Eng. Sci.* **375**, 20160394. https://doi.org/10.1098/rsta.2016.0394 (2017). ArXiv:1705.08552.

23. Bisio, A., Mosco, N. & Perinotti, P. Scattering and perturbation theory for discrete-time dynamics. arXiv:1912.09768 [quant-ph] (2019).

24. Feynman, R. P., Hibbs, A. R. & Styer, D. F. *Quantum Mechanics and Path Integrals* (Courier Corporation, 2010).

25. Arrighi, P., Di Molfetta, G. & Eon, N. A gauge-invariant reversible cellular automaton. In *International Workshop on Cellular Automata and Discrete Complex Systems*, 1–12 (Springer, 2018).

26. Di Molfetta, G. & Debbasch, F. Discrete-time quantum walks: Continuous limit and symmetries. *J. Math. Phys.* **53**, 123302 (2012).

27. Márquez-Martín, I., Arnault, P., Di Molfetta, G. & Pérez, A. Electromagnetic lattice gauge invariance in two-dimensional discrete-time quantum walks. *Phys. Rev. A* **98**, 032333 (2018).

28. Arnault, P., Di Molfetta, G., Brachet, M. & Debbasch, F. Quantum walks and non-Abelian discrete gauge theory. *Phys. Rev. A* **94**, 012335. https://doi.org/10.1103/PhysRevA.94.012335 (2016).
29. Di Molfetta, G. & Pérez, A. Quantum walks as simulators of neutrino oscillations in a vacuum and matter. *New J. Phys.* **18**, 103038 (2016).

30. Arrighi, P., Molfetta, G. D. & Eon, N. Non-abelian gauge-invariant cellular automata. In *Theory and Practice of Natural Computing—8th International Conference, TPNC 2019, Kingston, ON, Canada, December 9-11, 2019, Proceedings*, 211–221. https://doi.org/10.1007/978-3-030-34500-6_15 (2019).

31. Farrelly, T. & Streich, J. Discretizing quantum field theories for quantum simulation. arXiv:2002.02643 [hep-lat, physics:hep-th, physics:quant-ph] (2020). ArXiv:2002.02643.

32. Huerta, L. & Zanelli, J. Bose-Fermi transformation in three-dimensional space. *Phys. Rev. Lett.* **71**, 3622–3624. https://doi.org/10.1103/PhysRevLett.71.3622 (1993).

33. D'Ariano, G. M., Manessi, F., Perinotti, P. & Tosini, A. Fermionic computation is non-local tomographic and violates monogamy of entanglement. *EPL (Europhys. Lett.)* **107**, 20009, https://doi.org/10.1209/0295-5075/107/20009 (2014). ArXiv:1307.7902.

34. D'Ariano, G. M., Manessi, F., Perinotti, P. & Tosini, A. The Feynman problem and Fermionic entanglement: Fermionic theory versus qubit theory. *Int. J. Mod. Phys. A* **29**, 1430025. https://doi.org/10.1142/S0217751X14300257 (2014). ArXiv:1403.2674.

35. Farrelly, T. Insights from Quantum Information into Fundamental Physics. arXiv:1708.08897 [quant-ph] (2017).

36. Brun, T. A. & Mlodinow, L. Quantum cellular automata and quantum field theory in two spatial dimensions. *Phys. Rev. A* **102**, 062222. https://doi.org/10.1103/PhysRevA.102.062222 (2020). ArXiv:2010.09104.

37. Mlodinow, L. & Brun, T. A. Fermionic and bosonic quantum field theories from quantum cellular automata in three spatial dimensions. arXiv:2011.05597 [quant-ph] (2020).

38. Haase, J. F. et al. A resource efficient approach for quantum and classical simulations of gauge theories in particle physics. *Quantum* **5**, 393. https://doi.org/10.22331/q-2021-02-04-393 (2021). ArXiv:2006.14160.

39. Bender, J. & Zohar, E. Gauge redundancy-free formulation of compact QED with dynamical matter for quantum and classical computations. *Phys. Rev. D* **102**, 114517. https://doi.org/10.1103/PhysRevD.102.114517 (2020).

40. Parker, L. & Toms, D. *Quantum Field Theory in Curved Spacetime: Quantized Fields and Gravity* (Cambridge University Press, 2009).

## Acknowledgements

We benefited from discussions with Cédric Bény, Nathanaël Éon, Terry Farrelly, Alberto Verga. This publication was made possible through the support of the ID# 61466 grant from the John Templeton Foundation, as part of the “The Quantum Information Structure of Spacetime (QISS)” Project (qiss.fr). The opinions expressed in this publication are those of the author(s) and do not necessarily reflect the views of the John Templeton Foundation.

## Author contributions

All authors contributed equally to this work.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-06241-4.

**Correspondence** and requests for materials should be addressed to K.S. or G.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022