Sorbonne Universit´e, Observatoire de Paris, Universit´e PSL, CNRS, LERMA, F-75005, Paris, France

* **Correspondence: ** Fabrice Debbasch

**Academic Editors: **Hossein Hosseinkhani and Sotirios Baskoutas

**Special Issue:** Quantum Confinement Effects in Nano Material

**Received:** October 15, 2019 | **Accepted:** December 08, 2019 | **Published: **January 06, 2020

Recent Progress in Materials **2020**, Volume 2, Issue 1, doi:10.21926/rpm.2001002

**Recommended citation: **Debbasch F. Relativistic Wigner Function for Quantum Walks. *Recent Progress in Materials* **2020**;2(1):13; doi:10.21926/rpm.2001002.

© 2020 by the authors. This is an open access article distributed under the conditions of the Creative Commons by Attribution License, which permits unrestricted use, distribution, and reproduction in any medium or format, provided the original work is correctly cited.

**Abstract**

A relativistic Wigner function for free Discrete Time Quantum Walks (DTQWs) on the square 2*D *space-time lattice is defined. The transport equation obeyed by the relativistic Wigner function is derived in terms of discrete derivatives and degenerates in the continuum limit into the transport equation obeyed by the Wigner function of 2*D *Dirac fermions. The first corrections to the continuous equation induced by the discreteness of the lattice are also determined and the similarities with non-quantum relativistic transport in phase space are discussed further.

**Keywords **

Discrete time quantum walks; relativistic phase-space; relativistic transport equation

**1. Introduction**

Discrete Time Quantum Walks (DTQWs) are unitary quantum automata that can be viewed as the formal generalizations of classical random walks. Following the seminal work of Feynman [**1**], GrÖssing and Zeilinger [**2**] and Aharonov [**3**], they were considered in a systematic way by Meyer [**4**]. DTQWs have been realized experimentally with a wide range of physical objects and setups [**5**,**6**,**7**,**8**,**9**,**10**,**11**] which include quantum dots [**12**] and silicon nano-photonic chips [**13**]. DTQWs are studied in several diverse of contexts, ranging from quantum optics [**14**] to quantum algorithmics [**15**,**16**], condensed matter physics [**17**,**18**,**19**,**20**,**21**], including transport in graphene [**22**,**23**,**24**], hydrodynamics [**25**] and biophysics [**26**,**27**].

It is well known that several DTQWs can be viewed as discrete versions of the continuous Dirac and Weyl fermion dynamics [**28**,**29**,**30**,**31**,**32**,**33**,**34**,**35**]. These DTQWs have a continuum limit which coincides with the Dirac equation and they even display exact discrete gauge invariance properties [**36**,**37**,**38**,**39**,**40**].

However, the phase-space behavior of these DTQWs still remains unknown, which could be investigated well with the help of Wigner functions. Due to several fundamental and practical reasons, building Wigner functions for DTQWs is indeed an interesting problem. Quite generally speaking, there are essentially three possible logically complete formulations of quantum mechanics. Formulation 1 is the canonical quantization. It starts by considering position and momentum as operators which do not commute. The commutation relation is then realized by choosing position and momentum to be operators acting on time-dependent wave-functions defined over physical space and the time-evolution on the wave-function is fixed by the Hamiltonian. Formulation 2 is based on Feynman’s path integral, which attributes to each possible path between an initial state and a final state with a certain probability amplitude defined, at least formally, in terms of the classical *i.e. *non-quantum action of the system. Formulation 3 is based on the Wigner function [**41**], which is a quantum generalization of classical phase-space densities and whose time-evolution is fixed by its Moyal bracket [**41**,**42**] with the Hamiltonian.

DTQWs are historically connected to Formulation 2 because they originate in Feynman’s attempts to find a suitable discretization to compute path-integrals for fermions. DTQWs also connect with Formulation 1, because they can be viewed as wave-functions defined over discrete space-time and evolving unitarily under the action of their so-called pseudo-Hamiltonian. However, the connection with Formulation 3 has not been established. Indeed, even if Wigner functions have already been considered for Dirac quantum walks [**43**,**44**,**45**,**46**], the proposed definitions are not manifestly covariant and do not coincide in the continuum limit, with the standard, manifestly covariant Wigner functions for Dirac fermions [**47**,**48**,**49**]. Nor does their evolution equation coincides, in the continuum limit, with the usual phase-space transport equation for Dirac fermions.

Moreover, defining Wigner functions for DTQWs is not a purely academic topic. For example, interesting results on trapped non-relativistic fermions obtained recently in terms of the non-relativistic Wigner function [**50**] suggest that a relativistic Wigner function is the proper tool to study confinement of Dirac fermions and Dirac quantum walks.

More generally, the Wigner function of a quantum system can be obtained by Fourier transformation of the space-time representation of the density operator of this system with respect to one of its variables. Wigner functions can thus be used to study how density operators evolve in time and, in particular, to follow their time-evolutions in the presence of noise. Working with Wigner functions is actually often easier than working with density operators as using phase-space is rather intuitive. For example, Wigner functions are usually the best way to write down and study the Master equations describing decoherence, if the classical limit of the transport equation coincides with standard Kolmogorov equations describing phase-space transport in stochastic non-quantum systems. Owing to the relevance of both DTQWs and decoherence to quantum information, constructing Wigner functions for DTQWs and determining their evolution equations are thus important issues in the field.

The aim of this article is to define a relativistic Wigner function for free DTQWs in flat 2*D *space-time. We first extend continuous distribution theory to analyze a discrete lattice. We then define a discrete relativistic Wigner function for DTQWs and derive the corresponding discrete relativistic transport equation, which we express in terms of discrete analogs to usual continuous derivatives. In the continuum limit, the discrete Wigner function tends towards the Wigner function of Dirac particles and the transport equation tends toward the equation obeyed by the Wigner function of Dirac fermions in continuous space-time. We also suggest the first correction to the transport equation induced by the discreteness of the lattice on which the DTQWs propagate and finally discuss the results in the last section.

**2. A Simple Dirac QW**

We work with two-component wave-functions Ψ defined in 2*D *discrete space-time where instants are labeled by *j *∈Z, spatial positions are labeled by *p *∈Z and Ψ* _{j }*= (

*ψ*)

_{j,p}

_{p}_{∈}

_{Z}. We introduce a basis (

*b*) = (

_{A}*b*) in Hilbertspace space and the components Ψ

_{L}, b_{R}*= (Ψ*

^{A }*Ψ*

^{L},*) of the arbitrary wave-function Ψ in this basis. The Hilbert product is defined by*

^{R}*< ψ, φ >*=$∑_{A,j,p}(ψ^A )_{j,p }^* (φ^A)_{j,p}$, which makes the basis (

*b*) orthonormal. Consider now the quantum walk Ψ

_{A}

_{j}_{+1 }=

*U*Ψ

_{j}T*where*

_{j }*T*is the spatial-translation operator defined by (

*T*Ψ

*)*

_{j}*=$(Ψ_{j,p+1}^L,Ψ_{j,p-1}^R)^T$and*

_{j,p }*U*is an

_{j }*SU*(2) operator defined by (

*U*Ψ

_{j}*)*

_{j}*=*

_{j,p}*U*(

*θ*)Ψ

*(1) where*

_{j,p}*U*(

*θ*) =$\begin{pmatrix} {\rm cos}\theta &-i {\rm sin}\theta \\ -i {\rm sin}\theta & {\rm cos}\theta \end{pmatrix}$

*,*(2) with constant

*θ*. It was shown earlier [

**38**] that this quantum walk, in the continuum limit defined by

*t*=

_{j}*jϵ*,

*x*=

_{p}*pϵ*,

*θ*=

*ϵm*, tends to the Dirac equation for a 2D spinor of mass

*m*in the flat Minkovski space-time with coordinates (

*t, x*) as

*ϵ*tends to zero.

**3.**** Basic Tools **

*3.1 Discrete Derivatives*

*3.1 Discrete Derivatives*

Let us now rewrite the above equations with the help of discrete covariant derivative. We define:

(*D _{j}f* )

*=$\frac{1}{2}$(*

_{j,p}*f*

_{j}_{+1,p}−

*f*

_{j}_{−1,p})

(*D _{jj}f*)

*= $\frac{1}{2}$(*

_{j,p}*f*

_{j}_{+1,p}+

*f*

_{j}_{−1,p}− 2

*f*

*)*

_{j,p}(*D _{p}f*)

*= $\frac{1}{2}$(*

_{j,p}*f*

_{j,p}_{+1}−

*f*

_{j,p}_{−1})

(*D _{pp}f*)

*= $\frac{1}{2}$(*

_{j,p}*f*

_{j,p}_{+1}+

*f*

_{j,p}_{−1}− 2

*f*

*) (3), where*

_{j,p}*f*is an arbitrary

*j*- and

*p*-dependent quantity. These are discrete versions of the usual partial derivatives.

Inverting the above equations delivers:

*f*_{j}_{+1,p} = *f** _{j,p}* + (

*D*

_{j}*f*)

*+ (*

_{j,p}*D*

_{jj}*f*)

_{j,p}*f*_{j}_{−1,p} = *f** _{j,p}* − (

*D*

_{j}*f*)

*+ (*

_{j,p}*D*

_{jj}*f*)

_{j,p}*f*_{j,p}_{+1} = *f** _{j,p} *+ (

*D*

*)*

_{p}f*+(*

_{j,p}*D*

*)*

_{pp}f

_{j,p}*f _{j,p}*

_{−1}=

*f*− (

_{j,p}*D*)

_{p}f*+ (*

_{j,p}*D*)

_{pp}f*(4).*

_{j,p}Let us now rewrite the equation of motion of the walk, Ψ_{j}_{+1 }= *U _{j}T*Ψ

*, in terms of the above discrete derivatives. One starts by subtracting Ψ*

_{j}*from both sides of the equation, getting Ψ*

_{j }

_{j}_{+1 }− Ψ

*= (*

_{j }*U*−1)Ψ

_{j}T*(5) which can be rewritten as (*

_{j},*D*Ψ)

_{j}*+ (*

_{j,p}*D*Ψ)

_{jj}*= ((*

_{j,p}*U*−1)Ψ

_{j}T*)*

_{j}*(6), or, equivalently, (*

_{j,p}*D*Ψ)

_{j}*= ((*

_{j,p}*U*−1)Ψ

_{j}T*)*

_{j}*− (*

_{j,p}*D*Ψ)

_{jj}*(7).*

_{j,p.}Let us now express *T*Ψ* _{j }*in terms of the discrete derivatives. Since $(TΨ_j)_{j,p}^L$=$(Ψ_j)_{j,p+1}^L$ , one can write $(TΨ_j)_{j,p}^L$=$(Ψ_j)^L$ $+(D_p Ψ_j^L)_{j,p} +(D_{pp} Ψ_j^L)_{j,p} $(8). Similarly,$(TΨ_j)_{j,p}^R$ =$(Ψ_j)^R$ -$(D_p Ψ_j^R)_{j,p} $ +$(D_{pp} Ψ_j^R)_{j,p} $ (9). It follows that ((

*U*−1)Ψ

_{j}T*)*

_{j}*= ((*

_{j,p}*U*−1)Ψ

_{j}*)*

_{j}*+ (*

_{j,p}*U*σ

_{j}_{3}

*D*Ψ

_{p}*)*

_{j}*+ (*

_{j,p }*U*Ψ

_{j}D_{pp}*)*

_{j}*(10), where*

_{j,p}*σ*

_{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*). The equation of motion of the QW can thus be rewritten as:

_{A}$(D_J Ψ^A)_{j,p} $=$(Uσ_3)_B^A (D_p Ψ^B)_{j,p} $ + (U-1)$\frac{A}{B} Ψ_{j,p}^B$ -$δ_B^A (D_{jj} Ψ^B)_{j,p}$+$U_B^A (D_{PP} Ψ^B)_{j,p} $ (11). Where (1)$\frac{A}{B}=δ_B^A=$1 if A=B and 0 otherwise. In the continuum limit, the left-hand side and the first term on the right-hand side deliver the differential terms in the Dirac equation, the second contribution to the right-hand side delivers the mass term while the other two terms vanish as they are of a higher order.

Let us finally mention the following identities, which will be used in the next sections:

*D _{j}*(

*fg*) = (

*D*)

_{j}f*g*+

*f*(

*D*) + (

_{j}g*D*)(

_{j}f*D*) + (

_{jj}g*D*)(

_{jj}f*D*) (12),

_{j}g*D*(

_{p}*fg*) = (

*D*)

_{p}f*g*+

*f*(

*D*) + (

_{p}g*D*)(

_{p}f*D*) + (

_{pp}g*D*)(

_{pp}f*D*) (13)

_{p}g*.*Note that the terms which disappear in the continuum limit are of order of 3, and not 2 in

*ϵ*.

*3**.2 Discrete **D**istributions*

*3*

*.2 Discrete*

*D*

*istributions*

We are interested in defining a relativistic Wigner function for DTQWs. This will involve Fourier analysis of the density operator of the walk with respect to one of its time- and one of its space-variables. The space-time grid is infinite in both time and space and there is no reason to suppose that the density operator of the walk vanishes when either of these variables, and especially the time variable, tends to be infinite. It follows that the standard discrete Fourier transform, which operates on a finite number of points, cannot *a priori *be used for the Fourier analysis of the density operator. A similar problem arises in continuous analysis, when one wants to perform the Fourier analysis of functions of a real variable which does not decrease sufficiently at infinity. In the continuous analysis, the standard remedy is to work with distributions. Therefore extend to the discrete setting to the usual continuous theory of distributions. In a case where the density operator *does *has compact support, the theory will give compatible with the usual harmonic analysis defined by discrete Fourier transforms. For the sake of simplicity, all ideas are now presented for functions and distributions of one single discrete variable *n *∈Z, which plays the role of *j *or *p*. The extension of these two discrete variables used in the following sections is straightforward.

We introduce as test functions the space of all functions *h *of *n *which admit a discrete Fourier transform *ĥ*^{ }defined by $ĥ(k) =∑_nexp(ikn)h_n$. (14). All these functions have compact support. The conjugate momentum *k *takes value in the first Brillouin zone of the lattice *i.e. *(−*π, π*). The inverse Fourier transform is thus defined by:

Any function *f *defined on Z can now be considered a distribution acting on this space of test functions and we define *f(h)=*$∑_nf_n h_n $ (16). We now define the discrete Fourier transform $\hat{f}$ of the distribution

_{n}=

*f*by its action on functions on the variable

*k*. Let

*H*be such a function. We define $\hat{f}$(

*H*) by

$$\hat{f}(H)=∑_n \frac{1}{2π} ∫_{-π}^πdkexp(ikn) f_n H(k)$$

$=< f, Ĥ >_n $(17). Introducing the natural product in Fourier-space *< u, v > _{k}*

_{ }=$\frac{1}{2π} ∫_{-π}^πdku(k) v(k)$ (18), the definition of $\hat{f}$ can be rewritten as

*<*$\hat{f}$

*, H>*=

_{k}*< f, Ĥ*

^{ }

*>*.

_{n}Consider the Fourier transform $\hat{F}$ of the derivative *F *= *D _{n}f*.

*<$\hat{F}$ * *, H> _{k}*$=< F,Ĥ> _n$

$=∑_n\frac{1}{2π} ∫_{-π}^πdkexp(ikn) F_n H(k)$

$=∑_n\frac{1}{2π }∫_{-π}^πdkexp(ikn) \frac{1}{2} (f_{n+1}-f_{n-1})H(k)$

$=∑_n\frac{1}{2π} ∫_{-π}^πdkexp(ikn) \frac{ 1}{2} (exp(-ik)-exp(ik))f_n H(k)$

$=∑_n\frac{1}{2π} ∫_{-π}^πdkexp(ikn) f_n (-i {\rm sin} k)H(k)$

$=< G,H> _k$(19), where *G*(*k*) = −*i*( sin*k*)$\hat{f}$(*k*). Thus,$\hat{F}$(*k*) = −*i*(sin*k*)$\hat{f}$(*k*) (20). In the continuum limit, the continuous position *N *and wave-vector *K *are related to *n *and* k* by *N* = *ϵn *and *K* = *k/ϵ*, so (20) becomes *ϵ * $\hat{{f}'}$(*K*)= - *i*sin(*ϵK*) $\hat{f}$(*K*) (21), where *f’*^{ }= *df/dN*. At orders 0 and 1, (21) reads simply ^{ }$\hat{{f}'}$(*K*) = −*iK*$\hat{f}$(*K*) (22), but one finds, for example at second order in *ϵ*:

**4****. Discrete Wigner Function**

We now introduce the ‘density’ in space-time$Ω_{j,p,n_j,n_p}^{AB}=(Ψ^A)_{j-n_j,p-n_p}^* (Ψ^B)_{j+n_j,p+n_p} ,(j,n_j,p,n_p )∈ Z^4,(A,B) ∈ \left \{ L,R \right \}^2$. We consider that this object is, at fixed (*j, p*), a discrete distribution acting on functions of (*n _{j},n_{p}*) which admit a discrete Fourier transform with respect to these variables and we define the discrete Wigner ‘function’

*W*as the Fourier transform of Ω

^{AB }*with respect to (*

^{AB }*n*):

_{j},n_{p}$W_{j,p,k_j,k_p}^{AB}=∑_{(n_j,n_p)∈Z^2}exp(+ik_j n_j + ik_p n_p) × (Ψ^A)_{j-n_j, p-n_p}^* (Ψ^B)_{j+n_j, p+n_p} $.(24)

Now we obtain from the equations of motion of the QW, an equation of motion for *W*. Following the computation carried out in the continuous case *i.e. *for the Dirac equation, we first compute the discrete derivatives of Ω* ^{AB }*with respect to

*j*(

*D*),

_{j}*p*(

*D*),

_{p}*n*(

_{j }*D*) and

_{nj}*n*(

_{p }*D*), where the derivatives

_{np}*D*and

_{nj }*D*are defined as

_{np }*D*and

_{j }*D*above. The derivatives of Ω are best computed using the identities (12, 13). One obtains

_{p }$(D_j Ω^{AB})_{j,p,n_j,n_p}=(D_j (Ψ^A)^*)_{j-n_j,p-n_p} (Ψ^B)_{j+n_j,p+n_p } + (Ψ^A)_{j-n_j,p-n_p}^* (D_j Ψ^B) )_{j+n_j, p+n_p}+(∆_j^{AB})_{j,p,n_j,n_p}$(25),

$(D_{nj} Ω^{AB})_{j,p,n_j,n_p}=-(D_j (Ψ^A)^*)_{j-n_j,p-n_p} (Ψ^B)_{j+n_j,_ p+n_p }+ (Ψ^A)_{j-n_j, p-n_p }^*(D_j Ψ^B))_{j+n_j,p+n_p}+\left (∆_{nj}^{AB}\right )_{j,p,n_j,n_p } $(26). Where $(∆_j^{AB})_{j,p,n_j,n_p }=(D_j (Ψ^A)^*)_{j-n_j,p-n_p} (D_{jj }Ψ^B)_{j-n_j,p-n_p}+ (D_{jj} (Ψ^A )^*)_{j-n_j,p-n_p} (D_j Ψ^B)_{j+n_j,p+n_p}$(27),$(∆_j^{AB})_{j,p,k_α,k_β}=-(D_j (Ψ^A)^*)_{j-n_j,p-n_p} (D_{jj} Ψ^B)_{j+n_j,p+n_p} + (D_{jj}(Ψ^A )^*)_{j-n_j,p-n_p} (D_j Ψ^B)_{j+n_j,p+n_p}$ (28). Similar relations can be written for *D _{p }*and

*D*. Putting all this together delivers $((D_j+D_{nj})Ω^{AB})_{j,p,n_j,n_p}={2(Ψ^A)^*}_{j-n_j,p-n_p } (D_j Ψ^B)_{j+n_j,p+n_p}+ 2(D_{jj }Ψ^A)_{j-n_j,p-n_p}^* (D_j Ψ^B)_{j+n_j,p+n_p}$(29). And $((D_p+D_{n_p })Ω^{AB})_{j,p,n_j,n_p}={2(Ψ^A)^*}_{j-n_j,p-n_p} (D_p Ψ^B)_{j+n_j,p+n_p} + 2(D_{pp} Ψ^A)_{j-n_j,p-n_p}^* (D_p Ψ^B)_{j+n_j,p+n_p}$(30). Using the equation of motion (11) leads to $((D_j+D_{nj})Ω^{AB})_{j,p,n_j,n_p}-(Uσ_3 )_C^B (D_p+D_{np})Ω^{AC})_{j,p,n_j,n_p }={2(Ψ^A )^*}_{j-n_j,p-n_p} (U-1)_C^B (Ψ^C )_{j+n_j,p+n_p} - 2(Ψ^A )_{j-n_j,p-n_p}^* ((δ_C^B D_{jj}-U_C^B D_{pp})(Ψ^C))_{j+n_j,p+n_p}+2(D_{jj} Ψ^A )_{j-n_j, p-n_p}^* D_j (Ψ^B))_{j+n_j,p+n_p}-2D_{pp} (Ψ^A)^*)_{j-n_j,p-n_p} (Uσ_3 )_C^B D_j (Ψ^C))_{j+n_j,p+n_p}$(31).

_{np}Taking the Fourier transform delivers:

$(δ_C^B D_j-(Uσ_3 )_C^B D_p)W^AC=K^{AB} [Ω]+M^{AB}$(32). Where $(K^{AB} [Ω])_{j,p,k_j,k_p} =∑_{(n_j,n_p)∈Z^2}exp(+ik_j n_j + ik_p n_p )× δ_C^B D_{nj}-(Uσ_3 )_C^B D_{np})Ω_{j,p,n_j,n_p}^{AC} $(33). And $(M^{AB})_{j,p,k_j,k_p} =(M_C^{AB} (W))_{j,p,k_j,k_p}+(M_S^{AB})_{j,p,k_j,k_p} $(34). With $(M_C^{AB} (W))_{j,p,k_j,k_p} =2(U-1)_C^B (W^{AC})_{j,p,k_j,k_p} $(35), $(M_S^{AB})_{j,p,k_j,k_p} =∑_{(n_j,n_p)∈Z^2}exp(+ik_j n_j + ik_p n_p ) μ_{j,p,n_j,n_p} $(36), $μ_{j,p,n_j,n_p} =-2(Ψ^A )_{j-n_j,p-n_p}^* ((δ_C^B D_{jj}-U_C^B D_{pp}))(Ψ^C))_{j+n_j,p+n_p}+2(D_{jj }(Ψ^A)^* )_{j-n_j,p-n_p}(D_j (Ψ^B))_{j+n_j,p+n_p}-2(D_{pp }(Ψ^A)^*)_{j-n_j,p-n_p} (Uσ_3 )_C^B (D_j (Ψ^C))_{j+n_j,p+n_p}$(37).

The quantity *K ^{AB }*[Ω] involves discrete Fourier transforms of discrete time- and space-derivatives. To rewrite these into a more appealing form, we use the computation carried out in the previous section. One finds:

$(K^{AB} [Ω])_{j,p,k_j,k_p} =∑_{(n_j,n_p)∈Z^2}exp(+ik_j n_j + ik_p n_p )× (-i( {\rm sin} k_j)δ_C^B+i( {\rm sin} k_p ) (Uσ_3 )_C^B ) Ω_{j,p,n_j,n_p}^{AC}=-i(( {\rm sin} k_j)δ_C^B-( {\rm sin} k_p)(Uσ_3 )_C^B)W_{j,p,k_j,k_p}^{AC}=(K^{AB} (W))_{j,p,k_j,k_p} $(38). The final form of the discrete evolution equation obeyed by *W *is thus:

$(δ_C^B D_j-(Uσ_3)_C^B D_p)W^{AC}-K^{AB }(W)-M_C^{AB} (W)=M_S^{AB}$(39).

*4**.1 Corrections to the **C**ontinuous **T**ransport **E**quation *

*4*

*.1 Corrections to the*

*C*

*ontinuous*

*T*

*ransport*

*E*

*quation*

Let us show how to perform an expansion around the continuum limit by computing the lowest order corrections to the transport equation induced by the discreteness of the space-time lattice.

At second order in, the operator *U*(*θ*) reads $U(θ)\sim \begin{pmatrix} 1-\frac{ϵ^2 m^2}{2}&-iϵm \\ -iϵm& 1-\frac{ϵ^2 m^2}{2} \end{pmatrix}$(40), where ∼ designates equality at second order. Since D_{p}=ϵ∂_{x}, one finds that $U(θ) σ_3 D_p\sim ϵ\begin{pmatrix} 1&+iϵm \\ -iϵm& -1 \end{pmatrix}∂_x $(41), leading to $δ_B^A D_j-(Uσ_3)_B^A D_p \sim ϵ(δ_B^A ∂_t-(σ_3 )_B^A ∂_x-ϵ^2 m(σ_2 )_B^A ∂_x $(42), where *σ*_{2 }is the second Pauli matrix.

Since $k_j=ϵk_t$ and $k_p=ϵk_x$ , $K^{AB} (W) \sim-iϵ(k_j δ_C^B-k_p (σ_3 )_C^B+ϵmk_p (σ_2 )_C^B )W^{AC}$ (43). The first part of the mass term reads $M_C^{AB} (W) \sim-iϵm(σ_1 )_C^B W^{AC}-\frac{ϵ^2 m^2}{2} W^{AB}$(44).

The other part of the mass term depends linearly on the second derivatives *D _{jj }*and

*D*. To compute these at second order in , one can use the first order expressions of

_{pp}*D*and

_{j }*D*, which are$D_j\simϵ∂_j$ and $D_p\simϵ∂_p$, combined with the zeroth order expression for

_{p}*U*, which is

*U*∼ 1. Using then the Dirac equation to eliminate the second derivatives leads to:

$M_S^{AB} \sim+2ϵ^2 m^2 W^{AB}$(45). Including the lowest order corrections to the continuous case, the transport equation obeyed by *W *thus reads:

$[(∂_t+ik_j ) δ_C^B-(∂_x+ik_p ) (σ_3 )_C^B+im(σ_1 )_C^B ] W^{AC}=ϵ[m(∂_x-ik_p ) (σ_2 )_C^B+\frac{3m^2}{2} δ_C^B ] W^{AC}$(46). This equation can be rewritten as $[K_j (γ^j )_C^B+K_p (γ^p )_C^B-mδ_C^B ] W^{AC}=-iϵ[mK¯_p (σ_3 )_C^B+\frac{3m^2}{2 }(σ_1 )_C^B ] W^{AC}$(47), where *K _{j }*=

*k*−

_{j }*i∂*,

_{t}*K*=

_{p }*k*−

_{p }*i∂*,

_{x}*K*

^{¯}

*=*

_{p }*k*+

_{p }*i∂*, and the Dirac matrices satisfying the 2

_{x}*D*Clifford algebra are

*γ*=

^{j }*σ*

_{1 }and

*γ*=

^{p }*iσ*

_{2}. Now, contract (47) on the left with $[K_j (γ^0 )_B^D+K_p (γ^1 )_B^D-mδ_B^D ]$.This delivers an equation of the form $(Μ+iK)_C^D W^{AC}=0$ where both

*M*and

*K*are hermitian. One thus deduces that both $M_C^D W^{AC}$and $K_C^D W^{AC}$vanish. The explicit computation delivers $M= M_0+ϵM_1$, $K=K_0+ϵK_1$(48). With $ M_0=k_j^2-k_{p-}^2 m^2-(∂_t^2-∂_x^2 )$(49), $M_1=m(k_j k_p+∂_t ∂_x ) σ_2+\frac{3m^2}{2} ∂_t 1+∂_x σ_3$(50),

*K*

_{0 }=−2

*k*+ 2

_{j}∂_{t }*k*(51), $K_1=-m(k_j ∂_x-k_p ∂_t ) σ_2-m(k_p^2+∂_x^2 ) σ_1+\frac{3m^2}{2 }k_j 1+k_p σ_3$(52).

_{p}∂_{x},It is clear that the equation *MW *= 0 is a mass-shell condition (hence the notation *M*) on the Wigner function *W*. Observe that, even for vanishing *i.e. *for the Dirac equation, the usual mass-shell condition is not exact, but only an approximation valid for Wigner functions that varies sufficiently slowly in both space and time. The *M _{1}*

_{ }term represents the correction to the mass-shell equation at the first order in the finite step of the lattice over which the DTQW is defined.

In order to understand the second equation, *KW *= 0, let us recall that, in flat 2*D *space-times, one particle relativistic non quantum transport can be described by introducing a scalar distribution *f *in 4*D *phase-space R^{4 }= {(*t,x,p ^{t},p^{x}*)} which verifies the equations in the form $p^μ \frac{∂f}{∂x^μ }+\frac{∂}{∂p^μ}(F^μ f)+\frac{∂^2}{∂p^μ ∂p^v } (D^{μv} f)=C(f,f)$(53). For example, in the relativistic Boltzmann or Boltzmann-Vlasov equations [51–53],

*F*is the force that acts on the particles,

*D*vanishes identically and

*C*(

*f,f*) is the relativistic Boltzmann collision kernel. For relativistic stochastic processes [

**54**], the transport equation is the so-called Kolmogorov equation (hence the notation K). For these processes,

*C*(

*f,f*) vanishes,

*F*represents the deterministic part of the force acting on the particles, and

*D*is the squared amplitude of a Gaussian noise.

Now we consider the free Dirac fermions, it is natural to set *F *= 0. Since there is no noise in the problem considered, *D *also vanishes, and so does *C *(*f,f*), because the one fermion under consideration does not collide with other fermions. The linear operator acting on *f *in the non-quantum transport equation then coincides formally with *K _{0}*

_{ }and describes the Dirac Hamiltonian flow in 4

*D*relativistic space-time. The operator

*K*

_{1}_{ }determines the correction induced by the finite lattice-step. If noise is added to the problem, extra terms involving at least the first and second order derivatives with respect to

*k*would appear in the transport equation, making it a quantum discrete Kolmogorov equation for fermions.

**5.**** Discussion**

We have defined a relativistic Wigner function for free DTQWs on the square 2*D *space-time lattice. This definition uses the concepts of discrete derivatives and distributions which have also been introduced. We have established the transport equation that is obeyed by the relativistic Wigner function and proved that this equation degenerates in the continuum limit into the transport equation that is obeyed by the Wigner function of Dirac fermions. We have finally computed the first corrections to this equation induced by the discreteness of the lattice and discussed how this equation connects with non-quantum relativistic transport in 4*D *phase-space.

It is intuitive to discuss how these results compare with others, such as non manifestly covariant definitions of relativistic Wigner functions. Let us focus for example on [**55**], which introduces non manifestly covariant Wigner functions for Dirac and Maxwell fields in the context of relativistic Quantum Field Theory (QFT). First, QFT is a quantum theory of fields, not of particles. In particular, the classical coordinates are classical fields, not positions in space-time, and the momenta are also fields. It may, therefore, appear unnatural to try and define Wigner functions for QFT which depend on coordinates in classical phase-space, and a more natural option would be to introduce Wigner ‘functions’ which are functionals of the fields and of their conjugate momenta. This route is certainly not trivial, but we believe it can produce relativistically covariant Wigner functionals. Leaving this to future publications, let us now suppose one nevertheless wants to construct, for relativistic QFT or for an approximate ‘relativistic quantum mechanics’, Wigner functions defined over a one-particle relativistic phase-space. Classical relativistic statistical physics offers several choices for the one particle phase-space. A first choice is the tangent bundle of the space-time manifold. A second choice is the cotangent bundle. In both cases, the phase-space has twice the dimension of the space-time manifold. Now, a classical particle lives on its mass-shell. One can therefore choose as classical relativistic phase-space the massshell bundle. At any given point in a *D*-dimensional space-time, the mass-shell is (one sheet of) a *D *− 1- dimensional hyperboloid in momentum space and the mass-shell bundle is thus of dimension 2*D *− 1.

All three choices are equivalent in classical statistical physics [**52**,**53**]. The transport equations obtained by choosing the tangent or the cotangent bundle as phase-space look simpler than those obtained by choosing the mass-shell bundle. But the dynamics encoded in these simpler-looking equations leave particles on their mass-shells, so restricting these equations to the mass-shells delivers the transport equations obtained by working directly with the mass-shell bundles. Because the associated transport equations look simpler, most authors working in classical relativistic statistical physics now choose to work in the tangent or cotangent bundle. Working with the mass-shell bundle may appear more natural in situations where the problem presents a preferred reference frame, as is the case for example with relativistic stochastic processes [**54**]. But, even in these cases, rewriting everything in a manifestly covariant manner by working in the tangent or cotangent bundle is interesting.

Things are more complex in quantum physics. The choice between working in a manifestly covariant manner in the tangent or cotangent bundle and working in a non manifestly covariant manner in a phase-space of lower dimension remains. However, both options now are not equivalent (in the mathematical sense). Consider for example the 2*D *Dirac equation. As shown by the equation M*W *= 0 derived and discussed in the previous section, the momentum *p *which appears as a variable of the Wigner function defined in a manifestly covariant treatment is *not *restricted to the classical mass-shell [**47**], nor to any hyper-surface in momentum space, even, which corresponds to 2*D *continuous Dirac fermions. The manifestly covariant treatment is thus not equivalent to a non manifestly covariant treatment developed in a phase-space of lower dimension. For the Dirac equation, both treatments exist. Indeed, references [47–49] use a manifestly covariant approach while reference [**55**] uses a non manifestly covariant one. This author believes that there are good reasons to stick to a manifestly covariant treatment, at least when no preferred reference frame exists, as is the case for the Dirac equation. For example, the fermionic transport equation obtained in [**55**] is much more complex than the transport equation obtained or used in [47–49] It, therefore, does not seem to admit a simple geometrical interpretation in terms of Hamiltonian flow in phase-space or mass-shell violation, making, in particular, the connexion to non-quantum transport hard to analyze. And changing reference frame in non-covariant transport equations is a really arduous task, as explicitly noted in [**55**].

Now QWs defined on a grid do exhibit a preferred reference frame, and this is why several authors have presented a non manifestly covariant construction of Wigner functions for QWs. The drawback of this approach is seen when one considers QWs which admit a continuous limit identical to the Dirac equation. The transport equation derived from a non manifestly covariant treatment does not admit as continuous limit the standard transport equation obeyed by the manifestly covariant Wigner function of the Dirac field, but rather the transport equation derived in [**55**]. As explained in the previous paragraph, this non manifestly covariant equation does not resemble classical relativistic transport equations and is, therefore, ill suited to study decoherence for these walks. We feel this is a serious drawback, especially because decoherence is a crucial issue in quantum information.

The results presented in this article can be extended in several directions. One should first address DTQWs defined, both on higher dimensional and on more general lattices, like for example planar triangular and hexagonal ones [**23**,**24**]. This could be done by using Fourier series defined on general spectral sets (see for example [**56**] and references therein). An extension to DTQWs defined on graphs should also be envisaged [57–59]. Since several DTQWs with non constant mixing operators can be interpreted as fermions coupled to discrete gauge fields [**36**,**37**,**38**,**39**,**40**,**60**,**61**,**62**,**63**,**64**,**65**], one should define a Wigner function which incorporates these gauge fields and, in particular the electromagnetic field and the gravitational field. For example, defining a gauge-invariant Wigner function for Dirac particles couples to electromagnetic fields is highly non trivial [**47**] and one wonders how the problem translates to DTQWs. Finally, DTQWs defined through unitaries that present a time randomness decohere and behave asymptotically like non-quantum diffusions (see for example [**66**] and references therein). One then expects the transport equation for the relativistic Wigner function to approach asymptotically relativistic transport equations similar to those obtained for relativistic stochastic processes [**54**,**67**]. This should be confirmed and the asymptotics fully analyzed.

**Author Contributions**

Fabrice Debbasch did all the work.

**Competing Interests**

The author has declared that no competing interests exist.

**References**

- Feynman RP, Hibbs AR, Styer DF. Quantum mechanics and path integrals. International Series in Pure and Applied Physics. McGraw-Hill Book Company. 1965.
- GrÖssing G, Zeilinger A. Quantum cellular automata. Compl Syst. 1988; 2: 197-208. [CrossRef]
- Aharonov Y, Davidovich L, Zagury N. Quantum random walks. Phys Rev A. 1993; 48: 1687. [CrossRef]
- Meyers DAT. From quantum cellular automata to quantum lattice gases. J Stat Phys. 1996; 85: 551-574. [CrossRef]
- Schmitz H, Matjeschk R, Schneider CH, Glueckert J, Enderlein M, Huber T, et al. Quantum walk of a trapped ion in phase space. Phys Rev Lett. 2009; 103: 090504. [CrossRef]
- Zahringer F, Kirchmair G, Gerritsma R, Solano E, Blatt R, Roos CF. Realization of a quantum walk with one and two trapped ions. Phys Rev Lett. 2010; 104: 100503. [CrossRef]
- Schreiber A, Cassemiro KN, Potocek V, Gabris A, Mosley PJ, Andersson E, et al. Photons walking the line: A quantum walk with adjustable coin operations. Phys Rev Lett. 2010; 104: 050502. [CrossRef]
- Karski M, Forster L, Cho JM, Steffen A, Alt W, Meschede D, et al. Quantum walk in position space with single optically trapped atoms. Science. 2009; 325: 174-177. [CrossRef]
- Genske M, Alt W, Steffen A, Werner AH, Werner RF, Meschede D, et al. Electric quantum walks with individual atoms. Phys Rev Lett. 2013; 110: 190601 [CrossRef]
- Sansoni L, Sciarrino F, Vallone G, Mataloni P, Crespi A, Ramponi R, et al. Two-particle bosonic-fermionic quantum walk via integrated photonics. Phys Rev Lett. 2012; 108: 010502. [CrossRef]
- Sanders BC, Bartlett SD, Tregenna B, Knight PL. Quantum quincunx in cavity quantum electrodynamics. Phys Rev A. 2003; 67: 042305. [CrossRef]
- Zhi-Hao B, Hao Q, Xiang Z, Jian L, Peng X. A quantum walk in phase space with resonator-assisted double quantum dots. Chinese Phys B. 2016; 25: 020307. [CrossRef]
- Wan LX, Yan LB, Huang JG, Zhang G, Kwek LC, Fitzsimons J, et al. A large-number and multilayer quantum walk using silicon nano-photonic chip. In Conference on Lasers and Electro-Optics, OSA Terchnical Digest (online); Optical Society of America, 2018. [CrossRef]
- Perets HB, Lahini Y, Pozzi F, Sorel M, Morandotti R, Silberberg Y. Realization of quantum walks with negligible decoherence in waveguide lattices. Phys Rev Lett. 2008; 100: 170506. [CrossRef]
- Ambainis A. Quantum walk algorithm for element distinctness. SIAM J Comput. 2007; 37: 210-239. [CrossRef]
- Magniez F, Nayak A, Roland J, Santha M. Search via quantum walk. In SIAM J. Comput - Proceedings of the thirty-ninth annual ACM symposium on Theory of computing; New York, 2007; 575-584. [CrossRef]
- Aslangul C. Quantum dynamics of a particle with a spin-dependent velocity. J Phys A. 2005; 38: 1-16. [CrossRef]
- Bose S. Quantum communication through an unmodulated spin chain. Phys Rev Lett. 2003; 91: 207901. [CrossRef]
- Burgarth D. Quantum state transfer with spin chains. PhD Thesis. 2006.
- Bose S. Quantum communication through spin chain dynamics: An Introductory Overview. Contemp Phys. 2007; 48: 13-30. [CrossRef]
- Di Molfetta G, Debbasch F, Brachet M. Non-linear optical Galton board: Thermalization and continuous limit. Phys Rev E. 2015; 92: 042923. [CrossRef]
- Bougroura H, Aissaoui H, Chancellor N, Kendon V. Quantum walk transport properties on graphene structures. Phys Rev A. 2016; 94: 062331. [CrossRef]
- Arrighi P, Di Molfetta G, M´arquez Mart´ın I, P´erez A. Dirac equation as a quantum walk over the honeycomb and triangular lattices. Phys Rev A. 2018; 97: 062111. [CrossRef]
- Jay G, Debbasch F, Wang JB. Dirac quantum walks on triangular and honeycomb lattices. Phys Rev A. 2019; 99: 032113. [CrossRef]
- Hatifi M, Di Molfetta G, Debbasch F, Brachet M. Quantum walk hydrodynamics. Scient Rep. 2019; 9: 2989. [CrossRef]
- Collini E, Wong CY, Wilk KE, Curmi PMG, Brumer P, Scholes GD. Coherently wired light-harvesting in photosynthetic marine algae at ambient temperature. Nature. 2010; 463: 644-647. [CrossRef]
- Engel GS, Calhoun TR, Read RL, Ahn TK, Manal T, Cheng YC, et al. Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature. 2007; 446: 782-786. [CrossRef]
- Strauch FW. Relativistic quantum walks. Phys Rev A. 2006; 73: 054302. [CrossRef]
- Strauch FW. Relativistic effects and rigorous limits for discrete- and continuous-time quantum walks. J Math Phys. 2007; 48: 082102. [CrossRef]
- Bialynicki-Birula I. Weyl, Dirac, and Maxwell equations on a lattice as unitary cellular automata. Phys Rev D. 1994; 49: 6920. [CrossRef]
- Yepez J. Relativistic path integral as a lattice-based quantum algorithm. Quantum Inf Proc. 2005; 4: 471-509. [CrossRef]
- Arrighi P., Nesme V. and Forets M. The Dirac equation as a quantum walk: Higher dimensions, observational convergence. J Phys A. 2014; 47: 465302. [CrossRef]
- Bisio A, D’Ariano G, Tosini A. Quantum field as a quantum cellular automaton: The Dirac free evolution in one dimension. Ann Phys. 2015; 354: 244-264. [CrossRef]
- Bisio A, D’Ariano G, Perinotti P, Tosini A. Weyl, Dirac and Maxwell quantum cellular automata. Found Phys. 2015; 45: 1203-1221. [CrossRef]
- D’Ariano G, Erba M, Perinotti P. Isotropic quantum walks on lattices and the Weyl equation. Phys Rev A. 2017; 96: 062101. [CrossRef]
- Cedzich C, Ryba´r T, Werner AH, Alberti A, Genske M, Werner RF. Propagation of quantum walks in electric fields. Phys Rev Lett. 2013; 111: 160601. [CrossRef]
- Di Molfetta G, Brachet M, Debbasch F. Quantum walks as massless Dirac fermions in curved space-time. Phys Rev A. 2013; 88: 042301. [CrossRef]
- Di Molfetta G, Brachet M, Debbasch F. Quantum walks in artificial electric and gravitational Fields. Physica A. 2008; 397: 157-168. [CrossRef]
- Arnault P, Debbasch F. Quantum walks and discrete gauge theories. Phys Rev A. 2016; 93: 052301. [CrossRef]
- Arnault P, Di Molfetta G, Brachet M, Debbasch F. Quantum walks and non-abelian discrete gauge theories. Phys Rev A. 2016; 94: 01233. [CrossRef]
- Curtright TL, Fairlie DB, Zachos ZK. A concise treatise on quantum mechanics in phase space. World Sci Pub Co Inc, 2014. [CrossRef]
- de Gosson M, Luef F. Moyal bracket and ehrenfest’s theorem in born-jordan quantization. Quantum Rep. 2019; 1: 71-81. [CrossRef]
- Hinarejos M, P´erez A, Banuls MC. Wigner function for a particle in an infinite lattice. New J Phys. 2012; 14: 103009. [CrossRef]
- Alberti A, Alt W, Werner R, Meschede D. Decoherence models for discrete-time quantum walks and their application to neutral atom experiments. New J Phys. 2014; 16: 123052. [CrossRef]
- Hinajeros M, Banuls MC, Perez A. Wigner formalism for a particle on an infinite lattice: Dynamics and spin. New J Phys. 2015; 17: 013037. [CrossRef]
- Hinajeros M, Banuls MC, Perez A. A study of wigner functions for discrete-time quantum walks. J Comput Theoret Nanosc. 2013; 10: 1626. [CrossRef]
- Vasak D, Gyulassy M, Helze HTh. Quantum transport theory for Abelian plasmas. Ann Phys. 1987; 173: 462-492. [CrossRef]
- Helze HTh, Gyulassy M, Vasak D. Transport equations for the QCD quark Wigner operator. Nucl Phys B. 1986; 276: 706-728. [CrossRef]
- Gao J, Pu S, Wang Q. Covariant chiral kinetic equation in Wigner function approach. Phys Rev D. 2017; 96: 016002. [CrossRef]
- Dean DS, Le Doussal P, Majumdar SN, Sher G. Wigner function of noninteracting trapped fermions. Phys Rev A. 2018; 97: 063614. [CrossRef]
- Israel W. Covariant fluid mechanics and thermodynamics: An introduction, in Relativistic Fluid Dynamics, A.M. Anile and Y. Choquet-Bruhat, eds, Lecture Notes in Mathematics 1385, 1987.
- Debbasch F, van Leeuwen W. General relativistic Boltzmann equation I: Covariant treatment. Physica A. 2009; 388: 1079-1104. [CrossRef]
- Debbasch F, van Leeuwen W. General relativistic Boltzmann equation II: Manifestly covariant treatment. Phys A. 2009; 388: 1818-1834. [CrossRef]
- Chevalier C, Debbasch F. Relativistic diffusions: A unifying approach. J Math Phys. 2008; 49: 043303. [CrossRef]
- Bialynicki–Birula I. Relativistic Wigner functions. EPJ Web of Conferences. 2014; 78. [CrossRef]
- Xu Y. Fourier series and approximation on hexagonal and triangular domains. Constr Approx. 2010; 31: 115. [CrossRef]
- Bru LA, de Valc´arcel GJ, Di Molfetta G, P´erez A, Rolda`n E, Silva F. Quantum walk on a cylinder. Phys Rev A. 2016; 94: 032328. [CrossRef]
- D’Ariano G, Erba M, Perinotti P, Tosini A. Virtually Abelian quantum walks. J Phys A. 2016; 50: 035301. [CrossRef]
- Bisio A, D’Ariano G, Erba M, Perinotti P, Tosini A. Quantum walks with a one–dimensional coin. Phys Rev A. 2016; 93: 062334. [CrossRef]
- Arrighi P, Facchini S, Forets M. Quantum walking in curved spacetime. Quantum Inf Proc. 2016; 15: 3467-3486. [CrossRef]
- Arnault P, Debbasch F. Landau Levels for discrete time quantum walks in artificial magnetic fields. Phys A. 2016; 443: 179-191. [CrossRef]
- Bru LA, Hinarejos M, Silva F, de Valca´rcel GJ, Rolda´n E. Electric quantum walks in two dimensions. Phys Rev A. 2016; 93: 032333. [CrossRef]
- Arrighi P, Patricot CH. Quantum walking in curved spacetime: (3+1) dimensions, and beyond. Quantum Inf Comput. 2016; 17.
- Arnault P, Debbasch F. Quantum walks and gravitational waves. Ann Phys. 2017; 383: 645-661. [CrossRef]
- Cedzich C, Geib T, Werner AH, Werner RF. Quantum walks in external gauge fields. J Mathemat Phys. 2019; 60: 012107. [CrossRef]
- Di Molfetta G, Debbasch F. Discrete-time Quantum Walks in random artificial Gauge Fields. Quant Stud Mathem Found. 2016; 3: 293-311. [CrossRef]
- Debbasch F, Mallick K, Rivet JP. Relativistic Ornstein-Uhlenbeck process. J Stat Phys. 1997; 88: 945. [CrossRef]