Discrete Geometry from Quantum Walks

A particular family of Discrete Time Quantum Walks (DTQWs) simulating fermion propagation in $2$D curved space-time is revisited. Usual continuous covariant derivatives and spin-connections are generalized into discrete covariant derivatives along the lattice coordinates and discrete connections. The concepts of metrics and $2$-beins are also extended to the discrete realm. Two slightly different Riemann curvatures are then defined on the space-time lattice as the curvatures of the discrete spin connection. These two curvatures are closely related and one of them tends at the continuous limit towards the usual, continuous Riemann curvature. A simple example is also worked out in full.

Quite remarkably, the DTQWs which converge towards fermions coupled to Yang-Mills fields also admit exact discrete gauge invariances and discrete field strength 'tensors'. These can be used e.g., to build new self-consistent discrete models of Dirac fermions interacting with Yang-Mills fields, where DTQWs are not only acted upon by gauge fields, but also act as sources to these fields. On the contrary, no exact discrete gauge invariance has been displayed for DTQWs which converge towards fermions coupled to gravitational fields, and no discrete field strength i.e., Riemann curvature has been defined either.
This article demonstrates how this gap can be filled. We focus on a certain family of DTQWs in discrete 2D space-time whose continuous limit coincides the dynamics of a (massless) Dirac fermion. We first show that these admit an exact discrete Lorentz gauge invariance and then present two alternate definitions of the Riemann curvature for these walks. More precisely, we define for each discrete walk discrete covariant derivatives in the direction of the grid coordinates. These derivatives generalise the usual covariant derivatives of differential geometry and allow the identification, not only of a discrete metric and a discrete 2-bein, but also of a discrete spin-connection defined on the lattice. The basic idea is then to define the Riemann curvature tensor of the space-time lattice (or of the DTQW) as the curvature of the spin connection using as gauge group the set of Lorentz transformations acting on spinors [32,33]. It turns out that there are actually two ways of implementing this idea and we therefore introduce two different discrete Riemann curvatures on the space-time lattice. The first Riemann curvature ρ * depends on a (nearly arbitrary) reference connection while the second one ρ s does not. It turns out that the curvature ρ * of the DTQW essentially represents the difference between the curvature ρ s of the DTQW and the curvature ρ s of the reference connection. We also show that, in the continuous limit, the Riemann curvature tends towards the usual, continuous Riemann tensor. We finally compute the curvature ρ s on a simple example before discussing all results.

Blueprint: The 2D Dirac Equation
The curved space-time Dirac equation is usually written in the form [34,35] iγ a e µ a D µ Ψ = mΨ, where Ψ is a spinor, and e µ a are the n-bein coefficients, which we suppose to be symmetrical, i.e., e µ a = e a µ (note that this relation makes sense because there are as many values of µ as there are values of a and because it equates coefficients i.e., tensor components in a given basis, and not tensors). The γ's are the so-called Dirac operators obeying the usual Clifford algebra, and In 2D space-time, Greek and Latin indices above only take two values, conveniently denoted by 0 and 1. The spin-connection ω has thus only two independent components ω 001 = −ω 010 and ω 101 = −ω 110 . The spinor Hilbert space is also two-dimensional and is equipped with the Hermitian product, where µ = (−det(g µν )) 1/2 where g µν are the metric components built from the 2-bein, which can be defined by g µν e µ a e ν b = η ab where (η ab ) = diag(1, −1) are the components of the 2D Minkovski metric in an orthonormal basis of the tangent space. We now choose an orthonormal basis (b − , b + ) in spinor space and represent an arbitrary spinor Ψ by its two components (Ψ − , Ψ + ). We also choose the Dirac operators γ 0 and γ 1 to ensure that their matrix representations in this basis coincide respectively with σ x and +iσ y where and are the first two Pauli matrices. The commutator then reads γ 0 , γ 1 = 2σ z = −2diag(1, −1). Expanding the compact notation above, the Dirac therefore equation reads: The 2-bein, the metric and the two non-vanishing connection coefficients can then be practically read off directly from the Dirac equation. Taking the continuous limit of the QWs considered in this article delivers this form of the Dirac equation [22]. In the next section, we will use discrete derivatives and put the QW equations in a form similar to (6) and thus identify in the discrete equations a 2-bein, a metric and a connection.
By definition, Lorentz transformations on spinors are generated by the commutator of the γ s. Thus, in 2D space-time, the Lorentz transform Ψ(Λ) of a spinor Ψ has components Ψ ± (Λ) = exp(±Λ)Ψ ± for an arbitrary, possibly space-and time-dependent Λ. In addition, the components of the spin connection transform according to ω µab (Λ) = ω µab + ∂ µ Λ. It follows from this that R µνab = ∂ µ ω νab − ∂ ν ω µab is invariant under Lorentz transformation. This quantity is the (µνab)-, so-called mixed component of the Riemann curvature tensor. The components R µναβ of the Riemann curvature tensor on the coordinate basis The Ricci tensor and the scalar curvature are defined from R µναβ in the standard manner. Note that the expression of R µνab is linear in the connection because the Lorentz group is abelian in 2D space-time. In what follows, a discrete Riemann curvature tensor will be computed by implementing Lorentz transformations on the discrete equations and identifying an invariant quantity.

A Simple Two-Step Quantum Walk
We work with two-component wave-functions Ψ defined in 2D discrete space-time where instants are labeled by j ∈ N and spatial positions are labeled by p ∈ Z and Ψ j = (ψ j,p ) p∈Z . We introduce a basis (b A ) = (b L , b R ) in Hilbert-space space and the components Ψ A = (Ψ L , Ψ R ) of the arbitrary wave-function Ψ in this basis. The Hilbert product is defined by < ψ, φ >= ∑ A,j,p (ψ A ) * j,p (φ A ) j,p , which makes the basis (b A ) orthonormal. Consider now the quantum walk Ψ j+1 = U j TΨ j where T is the spin-dependent spatial translation operator defined by (TΨ j ) j,p = (ψ L j,p+1 , ψ R j,p−1 ) T and U j is an SU(2) operator defined by where This article focuses on the two-step QW obtained by looking at the state of the original walk at only one in every two time steps, say the steps which correspond to even values of j (this is sometimes called the stroboscopic approach).
Written in full, the discrete equations of the two-step QW read: where c j,p = cos θ j,p and s j,p = sin θ j,p As shown in [22], this two-step QW admits a continuous limit if θ admits one and this limit coincides with the Dirac equation in a curved 2D space-time where the spinor connection and curvature depend on the derivatives of θ. The aim of this article is to show that the discrete equation can also be used to define a discrete metric, a discrete space-time connection and a discrete Riemann 'tensor' i.e., a full discrete geometry.

Covariant Discrete Derivatives
To define the geometry induced by this QW on the space-time lattice, it is necessary to change basis in the wave-function Hilbert space. The easiest way to do that is to write the equations of motion of the QW in an invariant, basis-independent manner by introducing covariant discrete derivatives in Hilbert space.
We start by defining the following simple, non covariant discrete derivatives: where f is an arbitrary jand p-dependent quantity. These are discrete versions of the usual partial derivatives. Inverting the above equations delivers: The equation of motion of the QW can then be rewritten as: where and σ 3 is the operator represented by the third Pauli matrix in the basis (b A ) i.e., σ 3 is represented by the matrix diag(1, −1) in the basis (b A ). Suppose now we change spin basis and rewrite (12) in a new, possibly jand p-dependent local basis b α = (b − , b + ). We need to introduce the operator r j,p which transforms the original basis b A into the basis b α , (b α ) j,p = (r j,p ) A α b A , and its inverse r −1 j,p . Thus, Let us now define covariant time-and space-derivatives, starting with derivation with respect to time. One has: This shows that (D j ψ A ) does not transforms as ψ A under a change of basis in Hilbert space, but this also suggests introducing a new, covariant time-derivative of the form where (A) = (A 0 , A 1 ) is an arbitrary jand p-dependent field. Using (15), one can write: where and and the time-and space-indices j and p have been omitted from the latest equations for readability purposes. Equation (17) proves that D j (A) is a covariant time-derivative. Space derivatives are slightly more complex. Using again (12), one can write: and As before, this suggests defining a spatial covariant derivative by: where (A) = (A 0 , A 1 , A 2 ) is an arbitrary jand p-dependent field and the transformation laws for A reads and The equation of motion (12) of the QW can be rewritten in terms of covariant derivatives. We introduce a time-connection A and a space-connection B, fixing only at this stage the values of their 1-and 2-components: We also introduce a 'mass'M and impose that thus ensuring that the equation of motion (12) can be written as: The 0-components of both connections and of the mass M will be specified in the next section. In addition, the status of M is discussed in the last section of this article.
Equation (30) is one step closer to the continuous Dirac equation that the original form of the equations of motion obeyed by the two-step walk. In particular, it shows that the time-connection A and the space-connection B are to be understood as two components of a single, space-time connection (A, B). this point of view will be adopted form here on.
The most important difference between (30) and (6) is that the operator Wσ 3 is not diagonal in the basis (b A ). Changing spin basis to make this operator diagonal is the goal of the next section.

Preliminary Gauge Change
Proceeding as in [22,23], we now change gauge i.e., spin basis by defining an operator r which puts W j,p σ 3 in diagonal form. The characteristic polynomial of W j,p σ 3 reads where δ j,p = c j,p−1 − c j,p+1 and π j,p = c j,p−1 c j,p+1 . Let (x α ) j,p , α = +, − be the two (possibly complex) roots of P j,p . From Equations (12) and (30), the eigenvalues (x α ) j,p actually determine two local transport velocities. More precisely, these eigenvalues actually define a set of local 2-bein coefficients (e µ 0 , e ν 1 ) jp (see above for details) on the space-time lattice. One finds e 0 0 = 1, e 0 This in turn defines the inverse metric 'components' on the space-time lattice In the usual differential, and thus continuous geometry, the Greek indices on n-bein coefficients, (inverse) metric components, etc. refer to components on the so-called coordinate basis (∂ µ ) = (∂ t , ∂ x ). In the discrete case, the equivalent of the basis (∂ µ ) is clearly the set (D j , D p ) and we therefore define accordingly the 2-bein 'vectors' e 0 = e j 0 D j + e p 0 D p and e 1 = e j 1 D j + e p 1 D p . The quantities g µν can be interpreted similarly as the components of the inverse metric g jj D Changes of space-time coordinates can then be implemented in the spirit of [36]. We finally define the discrete inverse 2-bein by the usual relations E a µ e µ b = δ a b where δ a b is the Kronecker symbol. We now recall that, in curved space-time, a spinor is normalized to unity, not with respect to the usual Lebesgue measure d 2 x , but with respect to the metric-induced measure (−detg)d 2 x where detg stands for the determinant of the metric components. This means that the usual Hilbert product < ψ, φ >= ∑ A,j,p (ψ A ) * j,p (φ A ) j,p , which makes the initial basis b A orthonormal, does not coincide with the natural Hilbertian product to be used in spinor space. We therefore define the new Hilbertian product by < ψ, φ > s = ∑ A,j,p µ j,p (ψ A ) * j,p (φ A ) j,p , a new basis (b α ) j,p made of two eigenvectors of Wσ 3 normalized with respect to < · > s and we define r j,p as the operator which transforms the original basis b A into the basis b α .

Choice of the Mass and Space-Time-Connection
Let us now specify the 0-components of the connections A and B as well as the mass M. Equations (29), (27) and (28) lead to: where We now define −iM α β as the non-diagonal part of N α β . This fully specifies M in any basis of the Hilbert space and it also leads to where O α β is the diagonal part of N α β . Since r was chosen to make (Wσ 3 ) α β diagonal, this last equation

Local Lorentz Transformations
Extending the definition of global Lorentz transformations for DTQWs proposed in [36], we now define the local Lorentz transform of the spinor Ψ by ψ − j,p → λ j,p ψ − j,p and ψ + j,p → λ −1 j,p ψ + j,p for an arbitrary, real and non-vanishing field λ defined on the 2D space-time lattice. Alternately, upon a Lorentz transformation, ψ → exp(Λσ z )ψ where λ = exp(Λ) and σ z is the operator represented by the third Pauli matrix in the basis (b − , b + ), and we use the practical notation which is not invariant under Lorentz transformation but becomes Note that the product M − + M + − , which can be interpreted as the squared mass of the walk, is invariant under Lorentz transformation.
The connection matrices also change under Lorentz transformation. Of particular interest are the diagonal parts of these connections because they obey a relatively simple transformation law. Indeed, The first two equations lead to where The following two equations lead to while the final two equations deliver Equating both expressions of exp(±2D pp Λ) delivers It is best to retain for tanh(+2D p Λ) an expression which does not favour a set of components over the other. We therefore choose as the final expression for tanh(+2D p Λ).

Riemann Curvature I
Assuming that (A 1 ) − − (A 1 ) + + does not vanish and inverting the functions sinh and tanh, Equations (43) and (52) can be rewritten under the form The identity [D j , D p ] = 0 then leads to Introduce now a reference connection (A * , B * ), with the sole constraint that L j (A * , A) and L p (B * , B) are both defined, and write Note that the identities L j (A, A) = L p (B, B) = 0 then imply We then define the discrete Riemann curvature ρ * jp (Λ) by By (56), which represents the discrete Riemann curvature ρ * of the connection (A, B) i.e., the curvature ρ * of the DTQW.

Riemann Curvature II
Suppose now that one is interested in a curvature which caracterizes only how the connection coefficients change under Lorentz transformations which vary slowly in time and space i.e., for which D j Λ, D p Λ and D pp Λ are all much smaller than unity. At the continuous limit, all Lorentz transformations are automatically slowly varying in both time and space because the time and space coordinates t and x are related to j and p by t j = j and x p = p, where is an infinitesimal [22,23], so that D p ∼ ∂ x and D pp ∼ 2 ∂ xx . However, slowly varying Lorentz transformations can also be considered outside the continuous limit (see the example in the next section).
The limit case of Lorentz transformations varying slowly in space is actually singular. Indeed, in the general case, Equations (37)-(42) relate the two independent variations B 0 (Λ) − B 0 and B 1 (Λ) − B 1 to the two independent discrete derivatives D p Λ and D pp Λ. Inverting these equations thus delivers D p Λ in terms of the two independent variables B 0 (Λ) − B 0 and B 1 (Λ) − B 1 . To study the limit case of slowly varying Lorentz transformations, suppose D p Λ = O( ) and D pp Λ = O( α ) with α > 1. Equations (37)-(42) then read: At the first order in , D pp Λ vanishes from the equations so both variationsB 0 (Λ) = B 0 (Λ) − B 0 andB 1 (Λ) = B 1 (Λ) − B 1 depend on the single variable D p Λ and they are therefore not independent. Indeed,B 1 (Λ) = 2B 2B0 (Λ). In this limit, the general problem, which depends on two variables, thus degenerates into a single variable problem, thus making the limit singular. To define curvature, one then needs only one of the two variations and it is natural to retainB 0 (Λ). The equation for D p Λ then reads and the equation for D j Λ becomes similarly From this choice and the identity [D j , D p ] = 0 follows where the s index stands for 'slow'. The 'slow' discrete Riemann curvature tensor ρ s jp (Λ) of a connection is then defined by: where the index 's' stands for slow, ensuring that ρ s jp (Λ) = ρ s jp (0). In addition, the Riemann of the DTQW is defined as ρ s jp (0).

Relation between the Two Riemann Curvatures
Let us now investigate how this second discrete Riemann tensor is related to the first one introduced in the previous section. To do so, suppose that both connections (A * , B * ) and (A(Λ), B(Λ)) are close to (A, B), in the sense that their coefficients in the basis (b α ) are close to those of (A, B). This implies in particular that D j Λ, D p Λ and D pp Λ are small (see Equations (37)-(42)) i.e., that Λ is slowly varying in time and space. To simplify the discussion, we also suppose that there exist a Λ * such that (A * , B * ) = (A(Λ * ), B(Λ * )), and Λ * is then also slowly varying in time and space. We now convert L j and L p into a functionL j of (Ā(Λ), A) and a functionL p of (B(Λ), B) and expand these two newly introduced functions in their first variable at first order around 0. This leads to: From this follows that which leads to whereρ Using again (67), this becomes Now, by definition,L j represents D j Λ andL p represents D p Λ. Since we are considering Λ's which vary slowly in time and space, Equations (63) and (64) are valid. Thus, In particular, ρ * jp (0) = ρ s jp (0) − ρ s jp (Λ * ), which links the two Riemann curvatures ρ * and ρ s of the space-time lattice.

Continuous Limit
Let us now discuss the continuous limit of ρ s . The continuous limit addresses situations where the operator U and the wave-function of the walk vary on time-and space-scale much larger than the grid cell. The physical time t and spatial coordinate x along the grid are related to j and p by t j = j and x p = p where is an infinitesimal. It has been shown in [22,23] that the continuous limit of the 2-step walk then coincides with the Dirac equation in a curved space-time with metric (g µν ) = diag(1, cos −2 θ). In particular, the matrices representing A 1 and B 1 in the basis (b α ) then tend towards unity while the matrix representing A 0 tends towards −(ω 001 /2) × 1 and B 0 tends towards −(ω 101 /2) × 1. The discrete Riemann curvature then tends towards 2 /2 × R µνab , where R µνab is the mixed component of the usual Riemann curvature tensor to µ = 0, ν = 1, a = 0, b = 1. This component contains all the information one needs about the Riemann tensor because this tensor, in 2D space-times, has only one independent component. The 1/2 in the multiplicative factor comes form the fact that the zeroth components of the discrete connection tend towards ω/2 (as opposed to ω). The 2 factor comes from the fact that curvatures are obtained by taking second discrete or continuous derivatives and that the above relation between (j, p) and (t, x) implies D j = ∂ t and D p = ∂ x . Finally, the components R µν01 of the continuous Riemann curvature tensor on the coordinate basis (∂ µ ) can be recovered by taking the continuous limit of E a µ E b ν ρ s where (E a µ ) is the discrete inverse 2-bein.

Example
The continuous limit of the walks studied in this article corresponds to the propagation of a Dirac spinor in a space-time metric of the form ds 2 = dt 2 − a 2 (t, x)dx 2 where t and x are the continuous coordinates corresponding to j and p and a(t) = 1/(cos θ). Fixing these coordinates i.e., retaining this form for the metric, the simplest space-times with non vanishing curvature are realized by choosing the function a independent of x. We now therefore choose an angle θ which depends only on j and proceed to compute, as an example, the first of the curvatures defined above. Since nothing depends on the spatial position, all quantities are now indexed by j only.
For such walks, the operators W and L take the simpler form and with ∆θ j = θ j+1 − θ j . A simple computation leads to (x ± ) j = ± | c j |. These values of (x ± ) j lead to (g jj ) j = 1, (g pp ) j = −c 2 j and (g jp ) j = 0. If θ j = π/2, the components of the discrete metric itself read (g jj ) j = 1, (g pp ) j = −c −2 j and (g jp ) j = 0. In addition, −µ 2 = −c 2 j . We now retain (assuming c j = 0) where κ j = cos(θ j+1 /2) and σ j = sin(θ j+1 /2). The matrix (r A α ) j can be read off these equations: and its inverse reads: The components of W + L are not modified by the change of basis i.e., and a direct computation delivers where the index j tracing the time-dependence of all quantities has been suppressed for readability purposes.
Since all angles depend only on j, only the connection B enters the curvature. One finds that This leads to ρ s j = D j (B 0 ) − − with (B 0 ) − − given by Equation (83).

Conclusions
We have revisited a particular family of DTQWs whose continuous limit coincides with the 2D curved space-time Dirac dynamics written in synchronous coordinates. We have defined discrete covariant derivatives of the spinor wave-function along the grid coordinates, thus introducing discrete spin-connections and also generalised the notions of metric and 2-bein to the discrete lattice. We have then defined two different discrete curvatures from the transformation properties of the discrete spin-connections under Lorentz transformations. Both curvatures are closely related and one of them coincides, in the continuous limit, with the usual Riemann curvature from differential geometry. We have finally computed this discrete Riemann curvature on a particularly simple example.
Let us now comment on these results. In an arbitrary space-time, the most complete characterization of curvature is given the Riemann tensor. This tensor is usually computed from the space-time connection, but it can also be obtained from spinor connection [32,33]. The definition and computation of discrete curvature presented in this article thus start with a definition of discrete spinor connections for DTQWs, which is itself based upon the definition of discrete first and second discrete partial derivatives with respect to the grid coordinates. In the discrete case, spinor connections have a richer structure than in the continuous case because they contain more coefficients. In 2D space-time, a continuous spinor connection is fully defined by two coefficients, whereas one needs five coefficients to fully define a discrete spinor connection. These five coefficients can be partitioned into two sets, one of two coefficients pertaining to discrete covariant derivatives with respect to the discrete time index j, and one of three coefficients pertaining to covariant derivatives with respect to the discrete space index p. Note that these two sets only mix if one performs discrete Lorentz transformations in space-time, and these have not been considered in this article, where only Lorentz transformations in spinor space are carried out. We have therefore chosen, for readability purposes, to use a different letter for each set of coefficients (A defines discrete covariant time-derivatives and B defines discrete covariant space-derivatives). In addition, the discrete space-time connection is thus represented by (A, B).
The computation of the Riemann tensor as the curvature of the spin connection coefficients using as gauge group the set of Lorentz transformations in spinor space does not deliver the usual space-time components of the tensor, but the so-called mixed components R µνab , from which the usual space-time components can be recovered through partial contraction with the inverse n-bein coefficients. This applies both to the continuous and the discrete case. In 2D, there is only one independent component to the usual continuous Riemann tensor, and the discrete one also has only one independent component.
We have proposed two different definitions of Riemann curvature for the DTQWs considered in this article. If one is mainly interested in quantum simulation of conventional continuous physics, the second definition, which makes use of slowly varying Lorentz transformations, is clearly the one of choice, if only because its continuous limit gives back the usual Riemann curvature of differential geometry. However, DTQWs are interesting in other contexts, for example in quantum computing and quantum algorithms, where the continuous limit is not necessarily of particular importance. It is therefore useful to develop, for these contexts, a very general notion of curvature which is not linked to what happens at the continuous limit or for slowly varying Lorentz transformations. This is why we have offered our first definition of Riemann curvature. For technicality reasons, the first definition makes it necessary to introduce a reference connection and the obtained Riemann curvature thus depends on this reference connection. In essence, the reference connection is the connection for which the first Riemann curvature vanishes. The easiest way to see this is to go back to Equation (72). This states that the first curvature of a connection is the difference between the second curvature of that connection and the second curvature of the reference connection. There is no canonical choice for the reference connection i.e., the reference connection must be chosen on case by case basis, according to the context and interests of the computation.
The whole approach developed in this article is close in spirit to work which has been done in the last fifteen years, where classical Markov chains are used to define Ricci curvatures of graphs [37][38][39]. Indeed, a Markov chain is essentially a discrete diffusion. It therefore defines a Laplace operator on the discrete structure where it lives and, thus, a Ricci curvature. Similarly, a DTQW is essentially a spin 1/2 wave propagating on the lattice. Since a spin 1/2 wave obeys the Dirac equation, a DTQW essentially defines discrete equivalents to all quantities appearing in the Dirac equation i.e., an n-bein, and thus a metric, and a spin-connection. Once one has a discrete equivalent of the spin-connection, one can compute its curvature (in the sense of gauge theories), which coincides with the Riemann curvature. It is remarkable that classical Markov chains thus provide only a generalization of the Ricci curvature while quantum walks deliver equivalents to all geometrical objects of usual interest, from the metric to the spin-connection and, thus to the full Riemann curvature tensor.
As already mentioned, curvature is often introduced in differential geometry without using spinors. A standard approach is to first define a notion of parallel transport for tensor (including vector and 1-form) fields and then introduce curvature as the natural object which measures how much parallel transport along an infinitesimal closed loop modifies a tensor field. The parallel transport generates a covariant derivative of tensor fields and, thus, a space-time connection, encoded for example in the Christoffel symbols, which can be used to compute the space-time curvature. However, the space-time connection can also be represented by the so-called rotation coefficients, which determine the covariant derivatives of the n-bein. Now, this space-time connection can be extended in a canonical way to spinor fields. One can thus define the covariant derivatives of spinor fields and these can be used to compute directly the space-time curvature. This approach presents the advantage of being very close in spirit to Yang-Mills gauge theories: at each point in space-time the fiber is the spinor Hilbert space, the theory is invariant under the local action of various gauge groups (U(1), the Lorentz group, etc.) in this fiber, and one computes for each group the curvature or field strength from the associated covariant spinor derivative.
Let us add a few comments about the mass M. If one focuses on usual physics, only the continuous limit counts and it has been shown in [22,23] that the mass M then vanishes. Thus, at the continuous limit, the DTQW under consideration describes the continuous dynamics of massless fermions. By definition, the DTQW does not describe known physics outside this limit i.e., in the discrete regime. In particular, we decided to call the matrix M the mass of the DTQW because the form (30) of the discrete equation resembles the form (6) of the continuous equation and the coefficient in front of Ψ in (6) is −i times the physical mass. However, this does not presuppose that the matrix M, outside the continuous limit, shares any property with physical masses. In particular, the two coefficients M − + and M + − are not necessarily identical outside the continuous limit, where both vanish. In this sense, the matrix M should be considered as a generalised mass, which becomes a physical mass only at the continuous limit. Note that it is possible to construct DTQWs close to those considered in this article and whose continuous limits describe fermions of non vanishing mass (see, for example, [29]). Let us also recall that some DTQWs not considered in this manuscript have, even in the continuous regime, complex mass terms [23] which thus do not share the properties of physical masses.
Let us now conclude by mentioning possible extensions of this work. One should first address more general DTQWs coupled to arbitrary Yang-Mills fields. The extension to both higher dimensional space-times and higher spins should also prove interesting, starting with walks defined on square lattices, then moving on to more general grids, the ultimate goal being DTQWs on graphs. For example: what are the necessary graph properties for a DTQW to define a curvature on the graph? Or, how can one use graph geometry to write more efficient quantum algorithms? One should finally extend all these computations to alternate, comparable discrete models such as Lattice Gauge Theories (LGTs) and compare the results with those obtained for DTQWs.
Funding: This research received no external funding.