Relativistic Wigner Function for Quantum Walks

TL;DR

This paper introduces a relativistic Wigner function for 2D quantum walks on a lattice, deriving its transport equation and analyzing lattice-induced corrections, bridging discrete quantum models with continuous relativistic fermion descriptions.

Contribution

It defines a relativistic Wigner function for discrete quantum walks and derives its transport equation, including lattice correction terms, connecting discrete models with continuous relativistic quantum mechanics.

Findings

Derived the transport equation for the relativistic Wigner function.

Identified the continuous limit matching 2D Dirac fermions.

Computed first-order lattice corrections to the transport equation.

Abstract

A relativistic Wigner function for free Discrete Time Quantum Walks (DTQWs) on the square $2D$ space-time lattice is defined. Useful concepts such as discrete derivatives and discrete distributions are also introduced. The transport equation obeyed by the relativistic Wigner function is obtained and degenerates at the continuous limit into the transport equation obeyed by the Wigner function of $2D$ Dirac fermions. The first corrections to the continuous equation induced by the discreteness of the lattice are also computed.

Equations

$$(U_j{\Psi }_j{)}_{j,p}=U(\theta ){\psi }_{j,p}$$

$$U(\theta)=\left(\begin{array}[]{cc}\cos\theta&-i\sin\theta\\ -i\sin\theta&\cos\theta\end{array}\right),$$

$$(D_{j}f{)}_{j,p}$$

$$(D_{jj}f{)}_{j,p}$$

$$(D_{p}f{)}_{j,p}$$

$$(D_{pp}f{)}_{j,p}$$

$$f_{j+1,p}$$

$$f_{j-1,p}$$

$$f_{j,p+1}$$

$$f_{j,p-1}$$

$$(D_{j}\Psi^{A})_{j,p}=(U\sigma_{3})^{A}_{B}(D_{p}\Psi^{B})_{j,p}+(U-\mathbbm{1})^{A}_{B}\Psi^{B}_{j,p}\color[rgb]{0,0,0}-\delta^{A}_{B}(D_{jj}\Psi^{B})_{j,p}+U^{A}_{B}(D_{pp}\Psi^{B})_{j,p}\ ,$$

$$D_j(fg)=(D_{j}f)g+f(D_{j}g)+(D_{j}f)(D_{jj}g)+(D_{jj}f)(D_{j}g),$$

$$D_p(fg)=(D_{p}f)g+f(D_{p}g)+(D_{p}f)(D_{pp}g)+(D_{pp}f)(D_{p}g).$$

$$\hat{h}(k)=\sum _n\exp (ikn)h_n.$$

$$h_n=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }\exp (-ikn)\hat{h}(k)dk.$$

$$f(h)=<f,h>_{n}\color[rgb]{0,0,0}=\sum_{n}f_{n}h_{n}.$$

$$\hat{f}(H)$$

$$<u,v{>}_k=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }dku(k)v(k),$$

$$<\hat{F},H{>}_k$$

$$\hat{F}(k)=-i(\sin k)\hat{f}(k).$$

$$ϵ\widehat{f^{\prime }}(K)=-i\sin (ϵK)\hat{f}(K)$$

$$\widehat{f^{\prime }}(K)=-iK\hat{f}(K)$$

$$\displaystyle{\widehat{f^{\prime}}}(K)=-iK\left(1-\frac{\epsilon^{2}K^{2}}{3}\color[rgb]{0,0,0}\right){\hat{f}}(K).$$

$$W^{AB}_{j,p,k_{j},k_{p}}\color[rgb]{0,0,0}=\sum_{(n_{j},n_{p})\in{{\mathbb{Z}}^{2}}}\exp{\left(+ik_{j}n_{j}+ik_{p}n_{p}\right)}\times(\Psi^{A})^{*}_{j-n_{j},p-n_{p}}(\Psi^{B})_{j+n_{j},p+n_{p}}.\color[rgb]{0,0,0}$$

$$(D_j{\Omega }^{AB}{)}_{j,p,n_j,n_p}=(D_j({\Psi }^A{)}^{\ast }{)}_{j-n_j,p-n_p}({\Psi }^B{)}_{j+n_j,p+n_p}+({\Psi }^A{)}_{j-n_j,p-n_p}^{\ast }(D_j{\Psi }^B{)}_{j+n_j,p+n_p}+({\Delta }_j^{AB}{)}_{j,p,n_j,n_p},$$

$$(D_{n_{j}}\Omega^{AB})_{j,p,n_{j},n_{p}}\color[rgb]{0,0,0}=-\color[rgb]{0,0,0}(D_{j}(\Psi^{A})^{*})_{j-n_{j},p-n_{p}}(\Psi^{B})_{j+n_{j},p+n_{p}}+(\Psi^{A})^{*}_{j-n_{j},p-n_{p}}(D_{j}\Psi^{B})_{j+n_{j},p+n_{p}}+(\Delta_{n_{j}}^{AB})_{j,p,n_{j},n_{p}}\ ,$$

$$({\Delta }_j^{AB}{)}_{j,p,n_j,n_p}=(D_j({\Psi }^A{)}^{\ast }{)}_{j-n_j,p-n_p}(D_{jj}{\Psi }^B{)}_{j+n_j,p+n_p}+(D_{jj}({\Psi }^A{)}^{\ast }{)}_{j-n_j,p-n_p}(D_j{\Psi }^B{)}_{j+n_j,p+n_p},$$

$$(\Delta_{n_{j}}^{AB})_{j,p,k_{\alpha},k_{\beta}}=-(D_{j}(\Psi^{A})^{*})_{j-n_{j},p-n_{p}}(D_{jj}\Psi^{B})_{j+n_{j},p+n_{p}}+\color[rgb]{0,0,0}(D_{jj}(\Psi^{A})^{*})_{j-n_{j},p-n_{p}}(D_{j}\Psi^{B})_{j+n_{j},p+n_{p}}.$$

$${\left((D_j+D_{n_j}){\Omega }^{AB}\right)}_{j,p,n_j,n_p}$$

$${\left((D_p+D_{n_p}){\Omega }^{AB}\right)}_{j,p,n_j,n_p}$$

$$\displaystyle\left((D_{j}+D_{n_{j}})\Omega^{AB}\right)_{j,p,n_{j},n_{p}}-(U\sigma_{3})^{B}_{C}\left((D_{p}+D_{n_{p}\color[rgb]{0,0,0}})\Omega^{AC}\right)_{j,p,n_{j},n_{p}}=$$

$$\left({\delta }_C^B{D}_j-(U{\sigma }_3{)}_C^B{D}_p\right)W^{AC}$$

$${\left(K^{AB}\left[\Omega \right]\right)}_{j,p,k_j,k_p}=\sum _{(n_j,n_p)\in {\mathbb{Z}}^2}\exp \left(+i{k}_j{n}_j+i{k}_p{n}_p\right)\times \left({\delta }_C^B{D}_{n_j}-(U{\sigma }_3{)}_C^B{D}_{n_p}\right){\Omega }_{j,p,n_j,n_p}^{AC}$$

Full text

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

I Introduction

Discrete Time Quantum Walks (DTQWs) are unitary quantum automata that can be viewed as formal generalisations of classical random walks. Following the seminal work of Feynman Feynman (1965)Grössing and Zeilinger Grössing (1988) and Aharonov Aharonov (1993) they were considered in a systematic way by Meyer Meyers (1996). DTQWs have been realized experimentally with a wide range of physical objects and setups Schmitz (2009); Zähringer (2010); Schreiber (2010); Genske (2013); Karski (2009); Sansoni (2012); Sanders (2003); Perets (2008), and are studied in a large variety of contexts, ranging from quantum optics Perets (2008) to quantum algorithmics Ambainis (2007); Magniez (2007), condensed matter physics Aslangul (2005); Bose (2003); Burgarth (2006); Bose (2007); DiMolfetta (2015), hydrodynamics Hatifi (2019) and biophysics Collini (2010); Engel (2007).

It is well known that several DTQWs can be viewed as discrete versions of the continuous Dirac and Weyl fermion dynamics Bialynicki-Birula (1994); Yepez (2005); Strauch (2006, 2007); Arrighi (2014); Bisio (2015, 2015); D’Ariano (2017).

These DTQWs have a continuous limit which coincides with the Dirac equation and they even display exact discrete gauge invariance properties Cedzich (2013); DiMolfetta (2013, 2014); Arnault (2016, 2016). What however remains unknown is the phase-space behaviour of these DTQWs. If one follows the standard procedure adopted for Dirac fermions Elze (1986); Vasak (1987); Gao (2017), one should first build a relativistic Wigner function for DTQWs and then describe the phase-space dynamics by the transport equation obeyed by that function. Until now, the only Wigner functions that have been considered for DTQWs Hinarejos (2012); Hinajeros (2013, 2015); Alberti (2014) are non relativistic Curtright (2014). Thus, these function and the equation they obey do not coincide, at the continuous limit, with the usual Wigner function and phase-space transport equation for Dirac fermions.

The aim of this article is to fill this gap for free DTQWs in flat $2D$ space-time. We first generalize basic concepts of continuous mathematics such as derivation and distribution theory to analysis on discrete lattice. We then define a discrete relativistic Wigner function for DTQWs and derive the corresponding discrete relativistic transport equation. At the continuous limit, the discrete Wigner function tends towards the Wigner function of Dirac particles and the transport equation tends towards the equation obeyed by the Wigner function of Dirac fermions in continuous space-time. We also compute the first correction to the transport equation induced by the discreteness of the lattice on which the DTQWs propagate. We finally discuss all results in the last section of the article.

II A simple Dirac QW

We work with two-component wave-functions $\Psi$ defined in $2D$ discrete space-time where instants are labeled by $j\in\mathbb{Z}$, spatial positions are labeled by $p\in\mathbb{Z}$ and $\Psi_{j}=(\psi_{j,p})_{p\in\mathbb{Z}}$. We introduce a basis $(b_{A})=(b_{L},b_{R})$ in Hilbert-space space and the components $\Psi^{A}=(\Psi^{L},\Psi^{R})$ of the arbitrary wave-function $\Psi$ in this basis. The Hilbert product is defined by $<\psi,\phi>=\sum_{A,j,p}(\psi^{A})^{*}_{j,p}(\phi^{A})_{j,p}$, which makes the basis $(b_{A})$ orthonormal. Consider now the quantum walk $\Psi_{j+1}=U_{j}T\Psi_{j}$ where $T$ is the spatial-translation operator defined by $(T\Psi_{j})_{j,p}=(\psi^{L}_{j,p+1},\psi^{R}_{j,p-1})^{T}$ and $U_{j}$ is an $SU(2)$ operator defined by

[TABLE]

where

[TABLE]

with constant $\theta$. It was shown in DiMolfetta (2014) that this quantum walk, at the continuous limit defined by $t_{j}=j\epsilon$, $x_{p}=p\epsilon$, $\theta=\epsilon m$, tends to the Dirac equation for a $2$D spinor of mass $m$ in flat Minkovski space-time with coordinates $(t,x)$ as $\epsilon$ tends to zero.

III Basic tools

III.1 Discrete derivatives

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

[TABLE]

where $f$ is an arbitrary $j$- and $p$-dependent quantity. These are discrete versions of the usual partial derivatives. Inverting the above equations delivers:

[TABLE]

The equation of motion of the QW can then be rewritten as:

[TABLE]

where $\sigma_{3}$ is the operator represented by the third Pauli matrix in the basis $(b_{A})$ i.e. $\sigma_{3}$ is represented by the matrix $\mbox{diag}(1,-1)$ in the basis $(b_{A})$ and $(\mathbbm{1})^{A}_{B}=\delta^{A}_{B}=1$ if $A=B$ and [math] otherwise. At the continuous 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 because they are of higher order.

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

[TABLE]

[TABLE]

Note that the terms which vanish at the continuous limit are of order $3$, and not $2$ in $\epsilon$.

III.2 Discrete distributions

In what follows, we will consider discrete Fourier transforms of quantities which do vanish at infinity. To give meaning to these Fourier transforms, one has to extend the theory of distributions to the discrete case. For simplicity sakes, all ideas are now presented for functions and distributions of one single discrete variable $n\in{\mathbb{Z}}$, which plays to role of $j$ or $p$. The extension to 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 ${\hat{h}}$ defined by

[TABLE]

All these functions vanish at infinity. The conjugate momentum $k$ takes value in the first Brillouin zone of the lattice i.e. $(-\pi,\pi)$. The inverse Fourier transform is thus defined by:

[TABLE]

Any function $f$ defined on $\mathbb{Z}$ can now be considered a distribution acting on this space of test functions and we define

[TABLE]

We now define the discrete Fourier transform $\hat{f}$ of the distribution $f$ by its action on functions on the variable $k$:

[TABLE]

Introducing the natural product in Fourier-space

[TABLE]

the definition of ${\hat{f}}$ can be rewritten as $<{\hat{f}},H>_{k}=<f,{\hat{H}}>_{n}$.

Consider the Fourier transform of the derivative $F=D_{n}f$.

[TABLE]

where $G(k)=-i(\sin k){\hat{f}}(k)$. Thus,

[TABLE]

At the continuous limit, the continuous position $N$ and wave-vector $K$ are related to $n$ and $k$ by $N=\epsilon n$ and $K=k/\epsilon$, so (14) becomes

[TABLE]

where $f^{\prime}=df/dN$. At orders [math] and $1$, (15) reads simply

[TABLE]

but one finds, for example at second order in $\epsilon$:

[TABLE]

IV Discrete Wigner function

We now introduce the ‘density’ in space-time \Omega^{AB}_{j,p,n_{j},n_{p}}=(\Psi^{A})^{*}_{j-n_{j},p-n_{p}}(\Psi^{B})_{j+n_{j},p+n_{p}}\, $(j,n_{j},p,n_{p})\in{\mathbb{Z}}^{4}$, $(A,B)\in\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^{AB}$ as the Fourier transform of $\Omega^{AB}$ with respect to $(n_{j},n_{p})$:

[TABLE]

Let us now 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 $\Omega^{AB}$ with respect to $j$ ($D_{j}$), $p$ ($D_{p}$), $n_{j}$ ($D_{n_{j}}$) and $n_{p}$ ($D_{n_{p}}$), where the derivatives $D_{n_{j}}$ and $D_{n_{p}}$ are defined as $D_{j}$ and $D_{p}$ above. The derivatives of $\Omega$ are best computed using the identities (6,7). One obtains

[TABLE]

[TABLE]

where

[TABLE]

[TABLE]

Similar relations can be written for $D_{p}$ and $D_{n_{p}}$.

Putting all this together delivers

[TABLE]

and

[TABLE]

Using the equation of motion (5) leads to

[TABLE]

Taking the Fourier transform delivers:

[TABLE]

where

[TABLE]

and

[TABLE]

with

[TABLE]

[TABLE]

[TABLE]

The quantity $K^{AB}\left[\Omega\right]$ 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:

[TABLE]

The final form of the discrete evolution equation obeyed by $W$ is thus:

[TABLE]

In contrast with the transport equation obeyed by the Wigner function of the continuous Dirac field, this equation does not involve only $W$ but also a functional of $\Psi$, $\Psi^{*}$. In the discrete case, this equation should therefore be viewed primarily as a discrete integro-differential equation obeyed by $\Psi$ and $\Psi^{*}$, not $W$. The lowest order correction to the continuous limit is presented in the next section as an example.

IV.1 Corrections to the continuous transport equation

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

At second order in $\epsilon$, the operator $U(\theta)$ reads

[TABLE]

where $\sim$ designates equality at second order. Since $D_{p}=\epsilon\partial_{x}$, one finds that

[TABLE]

leading to

[TABLE]

where $\sigma_{2}$ is the second Pauli matrix.

Since $k_{j}=\epsilon k_{t}$ and $k_{p}=\epsilon k_{x}$,

[TABLE]

The first part of the mass term reads

[TABLE]

The other part of the mass term depends linearly on the second derivatives $D_{jj}$ and $D_{pp}$. To compute these at second order in $\epsilon$, one can use the first order expressions of $D_{j}$ and $D_{p}$, which are $D_{j}\sim\epsilon\partial_{j}$ and $D_{p}\sim\epsilon\partial_{p}$, combined with the zeroth order expression for $U$, which is $U\sim{\mathbbm{1}}$. Using then the Dirac equation to eliminate the second derivatives leads to:

[TABLE]

Including the lowest order corrections to the continuous case, the transport equation obeyed by $W$ thus reads:

[TABLE]

V Conclusion

We have defined a relativistic Wigner function for free DTQWs on the square $2D$ space-time lattice. This definition uses the concepts of discrete derivatives and distributions which we have also introduced. We have established the transport equation obeyed by the relativistic Wigner function and proved that this equation degenerates at the continuous limit into the transport equation 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.

These results 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 Arrighi (2018); Jay (2019). This could be done by using Fourier series defined on general spectral sets (see for example Xu (2010) and references therein). An extension to DTQWs defined on graphs should also be envisaged Bru (2016); D’Ariano (2016).

Since several DTQWs with non constant mixing operators can be interpreted as fermions coupled to discrete gauge fields Cedzich (2013); DiMolfetta (2013, 2014); Arrighi (2016); Bru (2016); Arnault (2016); Arrighi (2016); Arnault (2016, 2017, 2016); Cedzich (2019), 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 Vasak (1987) and one wonders how the problem translates to DTQWs. Finally, DTQWs defined through unitaries which present a time randomness decohere and behave asymptotically like non quantum diffusions (see for example DiMolfetta (2016) 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 Debbasch (1997); Chevalier (2008). This should be confirmed and the asymptotic fully analyzed.