Open access peer-reviewed chapter

Invariance in Quantum Walks

By Miquel Montero

Submitted: November 6th 2015Reviewed: March 7th 2016Published: August 24th 2016

DOI: 10.5772/62872

Downloaded: 1175


In this Chapter, we present some interesting properties of quantum walks on the line. We concentrate our attention in the emergence of invariance and provide some insights into the ultimate origin of the observed behavior. In the first part of the Chapter, we review the building blocks of the quantum-mechanical version of the standard random walk in one dimension. The most distinctive difference between random and quantum walks is the replacement of the random coin in the former by the action of a unitary operator upon some internal property of the later. We provide explicit expressions for the solution to the problem when the most general form for the homogeneous unitary operator is considered, and we analyze several key features of the system as the presence of symmetries or stationary limits. After that, we analyze the consequences of letting the properties of the coin operator change from site to site, and from time step to time step. In spite of this lack of homogeneity, the probabilistic properties of the motion of the walker can remain unaltered if the coin variability is chosen adequately. Finally, we show how this invariance can be connected to the gauge freedom of electromagnetism.


  • quantum walks
  • invariance
  • symmetry
  • Dirac equation
  • gauge transform

1. Introduction

In their origins [15], quantum walks (QWs) were thought as the quantum-mechanical generalization of the standard random walk in one dimension: the mathematical model describing the motion of a particle which follows a path that consists of a succession of jumps with fixed length whose direction depends on the random outcome of flipping a coin. In the quantum version, the coin toss is replaced by the action of a unitary operator upon some intrinsic degree of freedom of the system, a quantum observable with only two possible eigenvalues: for example, the spin of an electron, the polarization of a photon, or the chirality of a molecule.

After this preliminary analysis, it became clear that the similitude between these two processes was mainly formal and that random and QWs displayed divergent properties [6]. The most remarkable of these discrepancies is perhaps the ability of unbiased QWs to spread over the line, not as the square root of the elapsed time, the fingerprint of any diffusion process, but with constant speed [7]. This higher rate of percolation enables the formulation of quantum algorithms [8, 9] that can tackle some problems in a more efficient way than their classical analogs: For instance, QWs are very promising resources for optimal searching [1012]. Today, QWs have exceeded the boundaries of quantum computation and attracted the attention of researchers from other fields as, for example, information theory or game theory [1316].

As a consequence of this wide interest, diverse extensions of the discrete-time QW on the line have been considered in the past. Most of these variations are related with the properties of the unitary coin operator [17], backbone of the novel features of the process. Thus, one can find in the literature QWs whose evolution depends on more than one coin [1820], QWs that suffer from decoherence [21, 22], or QWs driven by inhomogeneous, site-dependent coins [2328]. There are also precedents where the temporal variability of the QW is explicit: in the form of a recursive rule for the coin selection, as in the so-called Fibonacci QWs [29, 30], through a given function that determines the value of the coin parameters [3133], or by means of an auxiliary random process that modifies properties of the coin [34].

The main goal in most of these seminal papers is to find out new and exciting features that the considered modifications introduce in the behavior of the system, like the emergence of quasiperiodic patterns or the induction of dynamic localization. Recent works [3537], however, have also regarded the issue from the opposite point of view, by exploring the conditions under which the evolution of the system results unchanged. In particular, Montero [37] considers the case of a discrete-time QW on the line with a time-dependent coin, a unitary operator with changing phase factors.

These phase factors are three parameters that appear in the definition of the coin operator whose relevance has been sometimes ignored in the past: When these phases are static magnitudes, they are superfluous [38], but if they are dynamic quantities, they can substantially modify the evolution of the system. This fact does not close the door to the possibility that a set of well-tuned variable phase factors can keep the process unchanged from a probabilistic perspective. This defines a control mechanism that can compensate externally induced decoherence and introduces a nontrivial invariance to be added to other well-known symmetries of QWs [3941].

In this Chapter, we will review the approach taken in [37] and consider a generalization of it. Now, the evolution of the discrete-time quantum walker on the line will be subjected to the introduction of a fully inhomogeneous coin operator: The properties of the unitary operator will depend both on the location and on the present time through the action of the aforementioned phase factors. This extra variability leads to additional constraints to be satisfied by these magnitudes if one wants to guarantee that the properties of the motion of the walker remain unaltered. Finally, we will connect our results with those appearing in the study of Di Molfetta et al. [36], where the authors considered how the inclusion of time- and site-dependent phase factors in the coin operator of a quantum walk on the line may induce some dynamics which, in the continuous limit, can be linked with the propagation of a Dirac spinor coupled to some external electromagnetic field. We will also explore the implications of this mapping here.


2. Fundamentals of QWs

We begin this Chapter with a survey of the fundamental concepts required in the designing of discrete QWs on the line. In its simplest version, the particle represented by the walker can occupy detached and numerable locations on a one-dimensional space. This space of positions may be just a topological space (a graph or a chain, for instance) or can be endowed with a metric. In such a case, it is usual to consider that the sites are separated by a fixed distance l, so that X=nl. Within this standard framework, time increases in discrete steps as well,

There is another kind of QW, called continuous quantum walk, in which the walker can modify its position at any time: this is the quantum counterpart of continuous-time random walk. The evolution of processes belonging to this category is ruled by a Hamiltonian and the corresponding Schrödinger equation. In spite they are different, discrete, and continuous QWs share common traits [42].

T=tτ, τbeing the sojourn time so that variable tbecomes a non-negative integer index, t{0,1,2,}, and the evolution of the system is just a sequence of states, |ψt.

Up to this point, there is no significant difference between random and quantum walks. The major distinction is found in the nature of the random event that determines the progress of the particle. While in a world governed by the laws of classical mechanics, randomness is the way in which we describe the uncertain effect of multiple (and usually uncontrollable) external agents acting upon a system, in the realms of quantum mechanics randomness is not an exogenous ingredient. This means that we can use some internal degree of freedom in the quantum system with two possible eigenvalues (the spin, the polarization, or the chirality) as a proxy for the coin and understand that any change in this inner property is the result of the act of tossing. Therefore, to represent the state of the walker, we need two different Hilbert spaces: P, the Hilbert space of particle positions spanned by the basis {|n:n}, and the Hilbert space of the coin states, C, which is spanned by the basis {|+,|}. The expression of |ψtin the resulting Hilbert space , CP, reads


where we have introduced the wave-function components ψ±(n,t), the two-dimensional projection of the state of the walker into the elements of the basis:


Now, we have to consider the mechanism that connects these two properties, position and quirality, which eventually leads to a model for the dynamics of ψ±(n,t). Evolution in the discrete-time, discrete-space QW can be regarded as the result of the action of operator

on the state of the system |ψt. As it can be observed, the practical implementation of operator
has two stages: In the first one, the unitary operator
modifies exclusively the internal degree of freedom of the quantum system, in what represents the throw of the coin as indicated earlier,


In a second step, the shift operator S^movesthe walker depending on the result obtained after the last toss:

With the present definition, the problem is spatially homogeneous and the system displays translational invariance. Therefore, alternative shift rules may be considered with equivalent results, as in the case of directed quantum walks [43, 44], where the particle can either remain still in the place or proceed in a fixed direction but never move backward.


Therefore, the state of the system at a later time |ψt+1is recovered by application of T^to the preset state:


and the complete evolution of the system is determined once |ψ0|ψt=0is selected. As in any quantum problem, one can consider for the initial state of the walker any combination of the elements in the basis of , a configuration that may lead to some degree of uncertainty in the position and/or the chirality of the system. However, the interest in establishing parallelisms between classical and quantum walkers encourages the choice in which, at the beginning, the particle position is known exactly, but its internal degree of freedom is aligned arbitrarily:


Needless to say that the linearity and the translational invariance of the problem ensure that the solution for a general initial state can be recovered by direct superposition of the evolution of Eq. (7), Eqs. (14) to (17) later.

The similarities and dissimilarities between classical and QWs must be grounded on the analysis of the probability mass function (PMF) of the process, ρ(n,t), the probability that the walker can be found in a particular position nat a given time t. The PMF for a random walk is


where pis the probability of obtaining a head as the result of flipping the coin. For the QW, ρ(n,t)is the sum of the squared modulus of the wave-function components,


On the basis of the values of the moduli of ψ±(n,t)we can also express the probability of obtaining a head value or a tail value when measuring the global coin state of the walker:


or the value of M(n,t),


another interesting magnitude that can be connected with the local magnetization of the system if the internal degree of freedom has its origin in the spin of the particle [45].

2.1. General solution

The evolution operator

induces the following set of recursive equations in the wave-function components,




whose general solution [38] can be written in a compact way by using ψ+(0,0)and ψ(0,0),


and the nonzero components of the wave function at time t=1,


since ψ+(1,1)=ψ(+1,1)=0, cf. Eqs. (12) and (13). In terms of the preceding quantities, and for n{t,t+2,,t2,t}, one has








It is noted that in this picture the evolution of each component depends only on their own initial values. In fact, it can be shown [38] that |ψ+(+1,1)|2can be understood as the “rightward initial velocity” of our quantum walker, whereas |ψ(1,1)|2would play the role of the “leftward initial velocity.”

Even though the expression for Λ(n,t)is completely explicit, Eq. (16), it may be instructive to show how the set of equations that cross correlate the evolution of the two components of the wave function, Eqs. (12) and (13), turns now into a single, two-step recursive formula that governs the whole dynamics:


Equation (16) is recovered from the above relationship once one considers the initial condition Λ(0,0)=1, together with the boundary conditions Λ(n,t)=Λ(n,t)=0, for nt1.

Observe how Λ(n,t)does not depend on χ, α, β, γ, or η. It is a function of θthrough the value of cosθ, a property that can be also observed in Eq. (18). One could infer from this feature that cos2θplays in QWs the same role of pin random walks and that the rest of parameters represent mathematical degrees of freedom without correspondence in the physical world. This impression can be strengthened by computing the value of the PMF in simple examples as, for instance, when ncoincides with t: in this case, ρclas.(t,t)=ptwhile ρ(t,t)=cos2tθ.

This conclusion is illusory, however. It is well known [19] that ρ(n,t)does not depend on χand that α, βand γappear in the PMF only in the following combination φ=α+βγ. But it is true as well that one needs to specify θ, φ, and ηto determine even the most basic aspects of the evolution of QWs. Figure 1 illustrates this fact. In the upper panel, we observe how the probability is distributed unevenly for positive and negative values of n, although θ=π/4. In the lower panel, we face the reversed situation, θ=π/8but ρ(n,t)shows no clear asymmetry.

2.2. Stationary PMF

Figure 1 also shows us that the disparity in the bias is not the most striking aspect that distinguishes QWs from their classical analogues. These differences can be appreciated more easily when one considers the stationary limit [34]. It can be shown [38] that for t1, the probability mass function ρ(n,t)is well described by ρ¯(n,t),


in the range tcosθ<n<tcosθ, 0<θ<π/2. As it can be seen in Figure 1, the agreement between ρ(n,t)and ρ¯(n,t)is greater for small values of n, whereas when |n| approaches to tcosθ, ρ(n,t)displays an oscillatory behavior around ρ¯(n,t). Regardless of this, Eq. (19) captures the essence of ρ(n,t): its U-shaped profile, with a central flat region and two local maxima in the vicinity of ±tcosθ. These traits are in clear contrast to the bell-shaped contour of the classical PMF, centered around (2p1)t, the mean value of the displacement of the random walker (Figure 1).

Figure 1.

Probability mass function aftert= 100 time steps. The dots correspond to the exact result for: (a)θ=π/4,η=π/16,φ=π; (b)θ=π/8,η=3π/16,φ=π; the boxes represent classical probabilities withp=cos2θ, whereas the black solid lines correspond toρ¯(n,t), cf. Eq.(19). We have only depicted probabilities for even values ofn, since in this case probabilities for odd values ofnare identically null.

Regarding the expectation value of the position of the quantum walker, Xt,


its magnitude does not stem from the location of the largest maximum of ρ(n,t), but has its origin in the skewness of the distribution. An elementary analysis of ρ¯(n,t)reveals that any bias in Xtis determined in the long run by the sign of the expression between parentheses in the right-hand side of Eq. (19). Therefore, as long as


the expectation value of the position of the walker will increase linearly with time:


as it can be checked in Figure 2. The converse is not true [40, 41]: in order to get quantum walkers that show an exact symmetry in the parity one has to demand that


Figure 2.

Expectation value of the position of the walker aftert =40 time steps. The dots correspond to the exact result for the QW withθ=π/6,η=π/6,φ=0, the boxes represent the classical mean position whenp=cos2θ=3/4, whereas the black solid line corresponds to the approximate law, Eq.(21), which in this case readsXtt/2whenl=1.

but also that

Eq. (23) implies |ψ+(+1,1)|2=|ψ(1,1)|2=1/2, see Eqs. (26) and (27) below. In other words, this is the condition that ensures the absence of bias in the “initial velocities.”


equations that have only three main families of solutions [38], being the most relevant of them the one corresponding to η=π/4, φ=π/2, for any choice of θ.


3. Inhomogeneous QWs

The fact that not only χbut even αand β(after a suitable choice of γ) can be completely ignored in the previous analysis can lead to the false conclusion that these phases can be disregarded in any other situation. We will devote the rest of this Chapter to the analysis of a framework where these magnitudes play a crucial role.

Consider a general inhomogeneous, time-dependent unitary operator U^t:


where αn,t, βn,t, χn,t, and θn,tare two-dimensional sets of real quantities. Now, we can define a new evolution operator

based on
and the standard shift operator
Eq. (5),
in such a way that the state of the particle at time t+1is the result of the application of T^tto |ψt,


In this case, the information supplied by the initial state of the system is not so important: Assume that |ψ0is of the form depicted in Eq. (7). Then, one has that


with φ0,0α0,0+β0,0γ. Note how this expression is invariant under the interchange


In practice, this means that we can modify θ0,0and α0,0+β0,0in order to obtain any desired value for |ψ+(+1,1)|and |ψ(1,1)|, irrespective of ηand γ. The complex arguments of ψ+(+1,1)and ψ(1,1)can be recovered with a suitable choice of χ0,0and α0,0β0,0.

The recursive equations of the wave-function components under the present dynamics induced by

are straightforward variations of Eqs. (12) and (13):




Since we have a specific interest in revealing a new kind of invariance, we will introduce ψ±(n,t), the solution to a certain inhomogeneous, time-dependent appealing problem




Therefore, our task is to find out nontrivial relationships connecting both set of parameters. Regarding this, note that θn,tis the same in both cases: as we have seen in Section 2, there are some features of the process, which are exclusively encoded in these magnitudes, and therefore, we will exclude them from the present analysis.


4. Invariance

The properties of the system enumerated up to this point are based on the moduli of the components of the wave function. This means, in particular, that if one has that ψ±(n,t)and ψ±(n,t)are linked through the following identities:




ρ(n,t)or M(n,t)will remain unchanged. The new magnitudes introduced in Eqs. (32) and (33), ξn,tand ζn,t, are two additional sets of arbitrary real constants, whose meaning will be discussed below.

If we assume the validity of Eqs. (32) and (33) and replace these expressions in Eqs. (28) and (29), the conditions to recover Eqs. (30) and (31) are


These equations lead to the following prescription to modify the phases leaving invariant the moduli of the components of the wave function:


4.1. Invariance of global observables

The first conclusion that can be drawn from Eqs. (34)(36) is that there is an infinite variety of choices for ξn,tand ζn,tthat does not modify the main properties of the quantum walker. The hard task is to identify those with a clear physical meaning or relevance. In a previous work [37], it has been considered one example that belongs to the following category:


This assumption simplifies enormously Eqs. (34)(36):


One particular choice that satisfies the above requirements is βn,t=β0, a constant value for all nand t, and the following functional forms for ξn,tand ζn,t:


a possible solution of Eqs. (37) and (38). The above expressions lead to the following homogeneous update rule for βn,t, t0,


where β1is an arbitrary constant, whose value cannot be assessed on the basis of the knowledge of ρ(t,n), P±(t), or M(n,t): it can only be inferred from the relative phase of the spinor components.

We illustrate in Figure 3 the invariance of ρ(t,n)in spite of the time- and site-inhomogeneous phase shifts that Eq. (44) introduces in the wave-function components, cf. Eqs. (42) and (43). Here, we have set θ=π/3, η=π/3, γ=0, χ=0, α=0, β0=0, and β1=1/10. With this choice, ψ±(n,t)are real functions that solve a stationary homogeneous problem, whereas ψ±(n,t)exhibit a complex, correlated behavior: for example, ψ(n,t)is a symmetric function around n=0, while neither the real part nor the imaginary part of ψ(n,t)shows this symmetry.

We can sketch a complementary picture that may help in the understanding the behavior of U^twhen βtfollows Eq. (44), through a geometrical analogy. Let us introduce ut, a time-dependent, unit-length vector in R3. Let us denote by θand βtits polar and azimuthal spherical coordinates, respectively. Then, we can recover the coin operator U^tthrough the scalar projection of the Pauli vector of operators, σ^, with Cartesian components


onto the utdirection, that is,


where I^Pis the identity operator defined in the position space P. The evolution of utis a step-like precession around the North Pole. Observe how, as in the example shown in Figure 2, when (β1β0)/πis an irrational number, the precession of utis not a periodic phenomenon at all. The absence of periodicity implies that vector utdefines an everywhere dense but enumerable subset of points in the ring associated with colatitude θon the sphere, and thus, the unconditional probability of choosing a particular value for βtis uniformly distributed in the stationary limit.

4.2. Exact invariance

Obviously, we can go further and demand exact invariance in the problem. This can be achieved by setting ζn,t=ξn,t. Eqs. (34)(36) read now [36]:


Figure 3.

Comparison of the wave function after t = 16 time steps. The red solid lines and dots correspond to a timehomogeneous QW. The blue-dotted lines show the real parts of the magnitudes associated with a time-dependent QW, while the imaginary parts are depicted by green-dashed lines.

As is shown below, these equations can be expressed in terms of finite differences which in turn lead to partial derivatives. In fact, in the expression of χn,t, it appears a time derivative, whereas the formulas for αn,tand βn,tcontain a spatial derivative. To illustrate these statements, consider the simple choice


Equations (46) to (48) read, as we have anticipated,


This means, in particular, that we can transform an inhomogeneous coin into a time-dependent one


4.3. Continuous limit

Let us express Eqs. (34) to (36) in a slightly different way. Consider the discrete difference operators Δnand Δtdefined as follows:


and similarly for Δnζn,tand Δtζn,t. In terms of these operators, Eqs. (34) to (36) now read:


Observe how the expression connecting βn,tand βn,tdepends explicitly on ξn,tand ζn,t, in the sense that it is not merely a function of the increments, cf. Eq. (41) above. In fact, we can rearrange the previous expressions in order to emphasize the distinct effects of ξn,tand ζn,t:


At this point, it is appropriate to note that we are not taking into account the issue of the parity of indexes nand t: since the instances of ξn,tand ζn,tthat appear in Eqs. (34) to (36) are those whose subscripts have the same parity, only one of the two terms in the right-hand side of Eqs. (53) and (54) is relevant or even well defined.

However, our interest in this Section is to analyze the continuous limit, τ0, l0. Up to the first order in τand l, one has that discrete difference operators Δnand Δtbecome partial derivatives:


We need to relate land τin order to obtain an unambiguous limit. We will assume that l=cτ, where cis the characteristic speed associated with the action of the shift operator S^upon the state of the walker. Therefore, depending on the physical nature of the system, crepresents the velocity at which the information is transferred, and it may coincide with the speed of light in vacuum. With this prescription, one has that Eqs. (55)(57) turn into


where ±are defined as follows,


and X±=cT±Xare the coordinates of the null geodesics in a flat (1+1) space and time. Observe how we have removed the subscripts: the dependency on Xand Tof all the magnitudes is implicitly assumed from now on.

The exact invariance, ζ=ξ, was analyzed in detail in the study of Di Molfetta et al [36]. There, it is shown how the recurrence equations of the wave-function components of the walker, Eqs. (28) and (29), can be mapped into equations describing the propagation of a Dirac spinor with charge eand masses m±coupled to a two-dimensional Maxwell potential A:


whose respective space-time components must change according to the formulas


Note that ζ=ξimplies that


and, when one introduces these relationships into Eqs. (67) and (68), one obtains the standard gauge transformations for the components of the potential A,


a transform that keeps invariant the electric field EXacting upon the system,


If we reconsider the example introduced at the end of Section 4.2,


we can conclude that it corresponds to a case in which the electric field EXis constant, where we are replacing the electric potential ϕ, ϕ=cAT, by a time-dependent magnetic potential AX,


In the most general case, when ζξ, the transformation rule for Ais


which departs from the gauge invariance of potential A. However, if we investigate the change in the electric field induced by Eqs. (75) and (76) we find


Clearly, ζ=ξis not the only solution to the constraint


that results in the invariance of the electric field. A possible choice is to demand that both ξand ζsatisfy the 2-dimensional wave equation by their own


Another alternative solution to Eq. (78) has appeared above, in Section 4.1. The equivalent expressions for Eqs. (37) and (38) in the continuous limit read:


what provides another solution to Eq. (78). Note that in this case Eqs. (65) and (66) show not merely covariance but perfect invariance in the mass-less case, m+=m=0, since



5. Conclusion

Along this Chapter, we have analyzed some interesting aspects of discrete-time QWs on the line, specifically those related with the emergence of invariance. In the first part, we have elaborated a succinct but comprehensive review covering the main features of the most elementary version of this process, when the unitary operator which assumes the function of the coin in the classical analog is kept fixed. We have described the dynamics that determines the evolution of the walker, supplied explicit formulas for assessing the precise state of the system at any time and approximate expressions that capture the main traits of the process in the stationary limit. These equations have been very useful to pinpoint the role played by the different parameters on the solution to the problem and put into context the generalization considered afterward.

The second part of the Chapter contemplates the situation in which the coin is time and site dependent. In particular, we have focused our interest on the phase parameters that define the unitary operator and determined the constraints that must be imposed in these changing phases if one wants to obtain invariance. This invariance can be demanded at two different levels: one can require that the invariance connects states belonging to the same ray of the Hilbert space or a milder condition, that the transformation modifies unevenly the two wave-function components. In this latter case, global properties (e.g., the probability that the particle is in a particular place or in a given spin state) remain unaltered but some other local quantum properties depending on the relative phase of these components can become modified.

The Chapter ends by analyzing the introduced invariance in the continuous limit. This approach unveils that the evolution of a time- and site-inhomogeneous quantum walk can be understood in terms of the dynamics of a particle coupled to an electromagnetic field and that the new symmetry shown by the walker can be interpreted as a manifestation of the well-known gauge invariance of electromagnetism.



The author acknowledges partial support from the Spanish Ministerio de Economa y Competitividad (MINECO) under Contract No. FIS2013-47532-C3-2-P and from Generalitat de Catalunya under Contract No. 2014SGR608.


  • There is another kind of QW, called continuous quantum walk, in which the walker can modify its position at any time: this is the quantum counterpart of continuous-time random walk. The evolution of processes belonging to this category is ruled by a Hamiltonian and the corresponding Schrödinger equation. In spite they are different, discrete, and continuous QWs share common traits [42].
  • With the present definition, the problem is spatially homogeneous and the system displays translational invariance. Therefore, alternative shift rules may be considered with equivalent results, as in the case of directed quantum walks [43, 44], where the particle can either remain still in the place or proceed in a fixed direction but never move backward.
  • Eq. (23) implies | ψ + ( + 1 , 1 ) | 2 = | ψ − ( − 1 , 1 ) | 2 = 1 / 2 , see Eqs. (26) and (27) below. In other words, this is the condition that ensures the absence of bias in the “initial velocities.”

© 2016 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Miquel Montero (August 24th 2016). Invariance in Quantum Walks, Research Advances in Quantum Dynamics, Paul Bracken, IntechOpen, DOI: 10.5772/62872. Available from:

chapter statistics

1175total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Quantum Walks

By Takuya Machida

Related Book

First chapter

Classical and Quantum Conjugate Dynamics – The Interplay Between Conjugate Variables

By Gabino Torres-Vega

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us