Dynamical kickback and noncommuting impurities in a spin chain

Dynamical kickback and non commuting impurities in a spin chain


In an interacting continuous time quantum walk, while the walker (the cursor) is moving on a graph, computational primitives (unitary operators associated with the edges) are applied to ancillary qubits (the register). The model with one walker was originally proposed by R. Feynman, who thus anticipated many features of the Continuous Time Quantum Walk (CTWQ) computing paradigm. In this note we examine the behaviour of an interacting CTQW with two walkers and examine the interaction of the walkers with noncommuting primitives. We endow such a walk with a notion of trajectory, in the sense of sample path of an associated Markov process, in order to use such notions as sojourn time and first passage time as heuristic tools for gaining intuition about its behaviour.
Keywords: Continuous time quantum walks; birth and death processes; quantum annealing.

1 Introduction

We consider a collection of spin systems , , coupled with an ancilla qubit through a Hamiltonian of the form:


where . The integers are supposed to satisfy ; the unitary operators act on the state space of the ancilla qubit. We will take, in this note,


and will suppose that all the remaining are the identity operator. We will consider an initial condition in the eigenspace belonging to the eigenvalue 2 of the conserved number operator


and we will refer the system to the orthonormal basis , where and , formed by the simultaneous eigenstates of , , and , belonging, respectively, to the eigenvalue of and , to the eigenvalue of the remaining and to the eigenvalue of .
We will look at the above system from two points of view:

  1. as an Anderson model [1] with noise on the hopping parameters relative to the links and , with the peculiarity that the “random values” of these parameters are determined by the non commuting observables and ;

  2. as a version of Feynman’s model of a quantum computer [2], where the motion of spin-up excitations of the field (the clock) administers the primitives and to the ancilla (the register).

For an extensive analysis of related models in the subspace belonging to the eigenvalue of we refer the reader to [3]. In the subspace (because of Peres’ conservation laws [4]) the presence of an ancilla cannot affect the motion of the clocking excitation. In the subspace we will, on the contrary, give evidence of a peculiar three-body effect involving the two clocking excitations and the ancilla qubit, related to the fact that and do not commute.

2 Dynamical kickback

The study of the system introduced in the previous section is made easy by the fact that the two following projectors are constants of motion:


where if then 1, else 0. If, as we will always do in this note, we consider the evolution of the system from the initial state


we will be interested only in the matrix elements of the Hamiltonian between states belonging to the range of :


Looking at the model from the point of view (ii), it performs quite a trivial computational task: starting with the two clocking excitations in positions 1 and 2 and with the register “up”, it returns, if the two clocking excitations are found in the terminal positions and , the register again “up”.
In the case the only track of the fact that, in applying the identity to the register, the machine has temporarily flipped it (by applying ), can be seen by comparing the corresponding probability amplitude with the one for the free case in which all the are the identity operator: the inversion of phase (“dynamical kickback”) shown in figure figure 1 is easily understood by thinking that, while the two clocking excitations move to the right, they restore the register into the “up” state by applying minus the identity operator to the register in the successive steps (from right to left) .

Figure 1: Frame (a): . Frame (b): . Solid lines: the probability amplitude as a function of time, under the initial condition (5). The dashed lines refer to the free case.
Figure 2: Graphical representation of the Hamiltonian (2) on the weighted graph having the set of vertices with edges between nearest neighbour sites. Edges to which (2) attributes a positive weight are represented by thick lines. Frame (a): . Frame (b): .

The above simplistic description of the evolution of the “computation” with two clocking excitations holds only in the particular case considered up to now, as made clear by the example with shown in figure 1.
By direct inspection of the Hamiltonian (2) and, in particular, of the weighted graph on which our quantum walk takes place (figure 2), the role of the positions and of the two noncommuting impurities and is easily understood in the context of an interference phenomenon. The term

in shows that in the situation of figure 1 all the amplitudes , are phase inverted with respect to the free case evolving according to the finite difference Laplacian

In the situation of figure 1, on the contrary, only the amplitudes , are phase inverted and deviations with respect to the free case take place because of interference with the uninverted signal : that this interference can be destructive, suppressing the probability of the two excitations ever going beyond the noncommuting impurities, is shown by comparison of the two frames of figure 1. A notational remark: for a wave function in the range of we are suppressing explicit indication of the argument .

3 Sample paths

It is fairly intuitive to attribute the effect shown in figure 1 to the fact that, in the case , not only the computational path is available, but also the path , corresponding to the fact that the rightmost cursor can wait in for the leftmost cursor to jump in ( being thus applied to the register) and then both of them can jump to the right of ( being thus applied).
The above intuition (involving correlations between positions of the cursors at different times) can be made more precise in terms of the stochastic process associated, according to the prescription of [5] (as specialized to the present context in [6]) to the time evolution , in , of the initial condition (5). The transition probability per unit time from site to site is given by


We will indicate by the analogously defined quantities in the absence of interaction. Simulation of processes with the above transition fields are performed, in what follows, according to the first order algorithm outlined in [6].
The free process (some sample paths of which are shown in the insets of figure 3) can, in the region , be described in much the same way as the paradigmatic example of [6]: each component starts as a pure birth process (only steps to the right or upwards are allowed in an initial time interval); at each instant each link (edge between nearest neighbour lattice sites) can be traversed only in one direction; the allowed direction along a link is inverted each time the probability mass at one vertex of the link vanishes.
Mastering the slalom at in the interacting case, as shown by figure 3, requires subtle time correlations. There is an instant of time at which most of the trajectories that hit are there simultaneously. At a later time they radiate from in many different directions. As figure 4 shows, given that the process hits , it stays there, on the average, for a longer time in the interacting case than in the free case.

Figure 3: . Frames (a) and (b): sample paths of the interacting processes and respectively. Only the trajectories, in our sample of size , that hit the site are shown. In the insets, the corresponding trajectories of the free process are shown for comparison purposes.
Figure 4: . Conditional cumulative distribution function of the first passage time (Frame (a)) and sojourn time (Frame (b)) at site , given that in the time interval the process visits . Solid line: the interacting case. Dashed line: the free case. Both distribution functions are estimated from the subsample of those trajectories that do hit .

4 Conclusions and outlook

We have tried, in this note, to contribute to the effort of looking at quantum mechanics as a source of metaphors suggesting Markov processes with interesting dynamical behaviour, interesting from the point of view of, say, efficiently crossing (in the sense of [7]) a graph or a decision tree, or sampling a function to be minimized [8, 9, 10].
In [6] we have, in this spirit, shown how to mimic, by a Markov process of the class proposed in [5], the diffraction effect due to a sharp initial position. Here we have tried to formulate in the same stochastic language the interference effects due to different localizations, as shown in figure 2, of discontinuities of the Hamiltonian.
In the process of doing so we have explored the notions of first passage and sojourn times for a quantum walk, which might prove useful from the point of view of suggesting heuristics of quantum algorithms, in a context, such as Feynman’s [2], in which timing and synchronization issues play a major role.
From the point of view of physics, our analysis raises the question of finding, in the quantum mechanics of Anderson localization by non commuting impurities, an analog of the time dependent phenomenon shown by the stochastic process in figures 3 and 4: the sudden formation (see figure 4), in the situation of figure 2, of a probability bubble at and its delayed (see figure 4) bursting in random directions.


  1. P. Lee and T. Ramakrishnan. Disordered electronic systems. Rev.Mod.Phys., 57:287–337, 1985.
  2. R.P. Feynman. Quantum mechanical computers. Found. Phys., 16(6):507–31, 1986.
  3. D. de Falco and D. Tamascelli. Entropy generation in a model of reversible computation. RAIRO: Inf. Theor. Appl., 40:93–105, 2006.
  4. A. Peres. Reversible logic and quantum computers. Phys. Rev. A, 32(6):3266–3276, 1985.
  5. F. Guerra and R. Marra. Discrete stochastic variational principles and quantum mechanics. Phys. Rev. D, 29(8):1647–1655, 1984.
  6. D. de Falco and D. Tamascelli. Quantum walks: a Markovian perspective. In V. Geffert and et al., editors, SOFSEM 08, LNCS 4910, pages 519–530. Springer, 2008.
  7. E. Farhi and S. Gutmann. Quantum computation and decision trees. Phys. Rev. A, 58(2):915–928, 1998.
  8. B. Apolloni, C. Carvalho, and D. de Falco. Quantum stochastic optimization. Stoc. Proc. and Appl., 33:223–244, 1989.
  9. B. Apolloni, N. Cesa-Bianchi, and D. de Falco. A numerical implementation of Quantum Annealing. In Albeverio et al. (Eds.), Stochastic Processes, Physics and Geometry, Proceedings of the Ascona/Locarno Conference, 4-9 July 1988, pages 97–111. World Scientific, 1990.
  10. G.E. Santoro and E. Tosatti. Optimization using quantum mechanics: quantum annealing through adiabatic evolution. J. Phys. A: Math. Gen., 39:R393–R431, 2006.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description