# Bound Molecules in an Interacting Quantum Walk

###### Abstract

We investigate a system of two atoms in an optical lattice, performing a quantum walk by state-dependent shift operations and a coin operation acting on the internal states. The atoms interact, e.g., by cold collisions, whenever they are in the same potential well of the lattice. Under such conditions they typically develop a bound state, so that the two atoms effectively perform a quantum walk together, rarely moving further from each other than a few lattice sites. The theoretical analysis is based on a theory of quantum walks with a point defect, applied to the difference variable. We also discuss the feasibility of an experimental realization in existing quantum walk experiments.

###### pacs:

42.50.Ex, 31.15.aq, 03.65.GeQuantum walks are the discrete-time and discrete-space analog of free quantum particle motion. They exhibit many of the well-known features of quantum dynamics: In the translation invariant case their motion is governed by a dispersion relation, from which the position distribution at large times is obtained as the distribution of group velocity. By modifying the coin operation of the quantum walk on a small number of sites, one can create the analog of potentials, with non-trivial scattering cross sections and bound states. In this paper we consider the possibility of an interaction between two walking particles, which becomes effective whenever the particles occupy the same site. In the simplest case, the Hadamard quantum walk, the interaction consists of just an additional phase. We show that with such interaction one can expect the formation of molecules, with all the features known from ordinary quantum dynamics of interacting two-particle systems: The two atoms stay together forever, performing a quantum walk in its own right. The distance between the particles remains bounded in the sense that they occupy nearby lattice sites with high probability, with estimates valid for all times.

Quantum walks have recently been demonstrated experimentally in setups based on controlled neutral atoms Karski et al. (2009), ions Schmitz et al. (2009); Zähringer et al. (2010), individual photons Broome et al. (2010), wave guides Peruzzo et al. (2010); Bromberg et al. (2009) or light pulses Schreiber et al. (2010). While all these experiments show the expected one-particle dynamics, their ability to demonstrate the type of interaction we describe differs very much. In some realizations interaction is very natural if not unavoidable: when two atoms sit in the same potential well of an optical lattice they are expected to pick up a collision phase depending on the duration of the contact. In spite of the fact that the atom-atom interaction would be much too weak for chemical binding in free space, the lattice creates a complex interference effect, which nevertheless leads to binding. The relevant conditions and parameters for an implementation of molecules of two atoms walking in a 1D optical lattice are discussed below. In a suitable system molecules are easy to prepare by starting two particles on the same site. The initial wave function then typically has a good overlap with the bound state, and the group velocity of unbound pairs is higher than that of molecules. Hence one can just wait for the unbound atom pairs to move away. Our conclusion is that an experiment should be possible in the near future, perhaps even by extending an existing experiment.

Our present work was partly inspired by analogous work in the continuous-time case Winkler et al. (2006), in which the motion between neighboring sites is given by tunneling. The spectral structure of continuous-time lattice dynamics is indeed analogous to that of discrete-time walks, except for the fact that in discrete time the spectrum of the unitary one-step operator is on the unit circle, and in the continuous-time case the spectrum of the Hamiltonian is on the real line, so it makes sense to speak of high energies and low energies. In contrast, for discrete-time quantum dynamics there is no distinction between repulsive and attractive interaction potentials. Hence, in contrast to Winkler et al. (2006), we cannot speak of “repulsively bound” molecules. In the discrete time case, a recent study of correlations in two-particle quantum walks Štefaňák et al. (2011), numerically found enhanced correlations related to the molecule states studied in this paper.

### Free Quantum Walks.—

Let us begin by describing quantum walks of a single particle. A walk runs on an -dimensional lattice, which we take as the cubic lattice . The walking particles generally have an additional quantum degree of freedom, described in a -dimensional Hilbert space . We can think of this as the internal states, or a “coin” degree of freedom, which can be “flipped” by a suitable unitary operation. The basis states of the system Hilbert space are thus with position and internal state label . The dynamical time step is given by a unitary operator , which commutes with lattice translations. A crucial condition is locality: has non-zero matrix elements only between states on lattice sites closer together than some fixed finite distance, called the neighborhood size of the walk. Since commutes with translation it is convenient to partially diagonalize this operator by Fourier transform, i.e., to look at the walk in momentum space. In the momentum representation wave functions depend on , i.e., the Brillouin zone of the lattice, with . The walk in this representation acts by multiplying the wave function at with a unitary -matrix . The locality condition then becomes the statement that each entry must be a polynomial in the variables . The standard example is the so-called Hadamard walk in one dimension (), for which

(1) |

Here the first factor is usually called the (internal state-dependent) shift and the second the coin. In one lattice dimension a decomposition into (usually several) such factors is always possible, but in higher dimension this is an open question. There is no intrinsic connection between the space dimension and the internal state space dimension . In particular, it is not necessary to choose a different coin for each dimension (which would give Inui et al. (2004)).

For a non-interacting walk of two distinguishable particles, we just take where and are the momentum coordinates of the two particles. Apart from a doubled lattice dimension () and a squared dimension of the internal state space (), this is a walk exactly as described above. Obviously, since commutes with the particle permutation, we may restrict the two particle walk to either the Bose or Fermi subspace. The influence of entangled initial states on the position distribution of the non-interacting walk has been studied theoretically in Omar et al. (2006), and experimentally with photons in Peruzzo et al. (2010).

### Interacting Quantum Walks .—

To make a two-particle walk interacting we need to introduce a dependence on the particle coordinates. The general way to introduce such modifications is to define a new transition operator as , where is the unperturbed, fully translation invariant two particle walk, and is a space dependent coin operation, i.e., a unitary which acts (in a possibly different way) on the internal degree of freedom at each site. We choose an operator , which differs from the identity operator only when the two particles are on the same lattice site (), in which case it acts on the joint coin space as some fixed operator . We denote by the projection operator onto the set of collision points, i.e., the diagonal in the -diagram. The overall interacting walk operator is

(2) |

Note that only acts on the translation degrees of freedom and acts on the coin space, so these operators commute and the second factor is unitary. Since is taken to be the same for every collision point, we thus break the separate translation invariance, but retain the invariance of the walk under joint translations. Therefore, the total momentum is conserved by .

We will develop the basic theory for such walks quite generally, for any in any coin dimension and also any lattice dimension (see Appendix A for the salient conditions). However, for illustration we will focus especially on the simplest case, which we call a walk with singlet collisions. This is defined by having the coin space at the collision points one-dimensional, so the coin acts just by a phase . In the Fermi case, and when the single particle coin space is two dimensional as for the Hadamard walk, this condition is automatically fulfilled even for the free walk. In the Bose case of the Hadamard walk, particles can only end up on the same site if they come from opposite directions, i.e., in a “symmetric singlet” . If the combination of free coin and collision coin leave this state invariant, we are again in the singlet collision case.

In order to define the interaction for more than two particles, or for a quantum lattice gas of such particles, more work is required. We describe this in the supplement Appendix B, mostly for the case of singlet collisions.

### Numerical Phenomenology.—

We start the walk (2) from the initial state in which both particles are at the origin, and their joint internal state is antisymmetric. Since the Fermi subspace is invariant under , all subsequent collisions will also be of singlet type. Then by multiplying the initial vector times with (2), and taking the modulus square with respect to the internal states, we get the joint position distribution of the two particles at time . The result is shown in Fig. 1, once for the non-interacting case (), and once for the interacting case with value .

Note that due to Fermi symmetry both diagrams are symmetric with respect to reflection at the diagonal. The peaks near are expected from the theory of the free Hadamard walk Ahlbrecht et al. (2011), and are near . The striking concentration of probability near the diagonal is the hallmark of the bound state. The width of this distribution in the off-diagonal direction remains constant in time, whereas along the diagonal it shows the characteristic behavior of a 1D walk. This will be proved below, and the peak velocity of this walk will be determined via Eq. (9) to be . For comparison, squares with edges at and are drawn in Fig. 1.

Turning from space and time diagrams to momentum and energy, let us consider the spectral properties of the operator (2). Since it commutes with joint translations, we can diagonalize it together with total momentum . For each value of we get an operator , whose eigenvalues we can compute. To treat this numerically, we close the system to a ring, which also discretizes . The eigenvalues of , as a function of , are shown in Fig. 2. Again, one part of this diagram is expected from the non-interacting case. Indeed, in that case also commutes with translations in the difference variable , which implies absolutely continuous spectrum for . These bands are seen in discretized form in Fig. 2, and the indication of continuous bands is that the spectral density increases with the size of the ring. In addition, however, we see a single line in the gap between the bands. This is the bound state.

### Analysis: quantum walks w ith a point defect.—

We first rewrite the problem in center of mass and relative coordinates. In momentum space this corresponds to the passage from the variables to and . Since is conserved we can fix and study the operator as an operator acting on the coordinate and the internal states. In this problem (for which is just an external parameter) the interaction appears as a modification of a free walk acting only at the origin. The analysis of point defects in a general quantum walk is of some interest in its own right. We therefore drop the parameter for the moment, and consider this general problem. So let

(3) |

denote the free walk matrix in relative momentum space, which corresponds, for each value of , to a single particle 1D walk with a four-dimensional coin space. We thus have to analyze the family of 1D quantum walks , where is the projection onto the subspace at the origin, and is the coin modification at the origin.

Since is a finite rank projection, the bands, i.e., the continuous spectrum of , are the same as those of (Kato, 1976, Thm.IV.5.35). To find the eigenvalues, suppose that is the momentum space representation of an eigenvector with eigenvalue . Then . Note that is just the value of at the origin in position representation, so is a function independent of . This determines the whole function by

(4) |

The condition for an eigenvalue is thus that is normalizable and . Since is the projection onto the zero Fourier component, i.e., integration with respect to , we introduce the operator

(5) |

and get the eigenvalue condition for in (4) to determine the eigenvector for eigenvalue in the form . This can be rewritten as

(6) |

Now it turns out that for not in the spectrum of , the operator on the right is unitary (see Appendix C). Therefore, the eigenvalue condition (6) can be satisfied for any in the gap in Fig. 2, for a suitable collision operator .

Applying this to our standard example, the Hadamard walk with singlet collisions, the space is one-dimensional, and the integral (5) is readily evaluated by the residue theorem, turning (6) into a formula for the collision parameter in , given the total momentum and the phase . It turns out that this relation can be solved for , giving

(7) | |||||

Here the constraint results from picking the correct pole for the residue evaluation in (5). It implies that we do not get two bound states for every momentum. Surprisingly, the branches forbidden by this constraint do not run inside the bands or off the unit circle but also in the gap. The results are shown in Fig. 3. The overlap of the initial wave function with the bound states, i.e. the probability to observe a molecule, can be computed from (5) and (7), see Appendix D. It achieves its maximum of at .

### Coin states of the molecule.—

If we restrict the Hilbert space of the relative motion of the two particles just to the bound states, we have a picture resembling very much the spectrum of a single walking particle with two internal states. Is it actually the same? That is, is there a quantum walk with two internal states and only nearest neighbor shifts which reproduces exactly the walking dynamics of the molecule? For this it is only necessary to compare the spectra , , because if these coincide, we can find a unitary commuting with translations and mapping one walk to the other. Since we have an explicit form of the spectrum, we can indeed identify an appropriate coin. The walk of the molecules with interaction phase is unitarily isomorphic to a walk with the second factor in (1) replaced by

(8) |

The exact form of the unitary isomorphism is determined by the bound states which are computed in Appendix E. The walk uses both branches in (7), though only “virtually” if one or both of them are forbidden by the constraint. It is not clear whether the bound states of general interacting quantum walks always allow such an interpretation.

### Velocities.—

It is apparent from the right panel of Fig. 1 that the molecules are slower than the free particles. The square in that figure corresponds to the propagation of one site per step. The free Hadamard walk has maximal group velocity , corresponding to the off-diagonal peaks. The group velocity of the molecules (with ) is

(9) |

with the same constraint as in (7). The maximal speed (always at ) is for the walk with , as shown in Fig. 1. For smaller the maximal speed according to (9) approaches one, but becomes forbidden by the constraint. Nevertheless it is possible to design molecules, which are faster than the free atoms, see Appendix F.

### Prospects for experimental realization.—

We plan to realize walking molecules by extending an existing experiment Karski et al. (2009), in which the walk of single Cs atoms has been demonstrated. These happen to be Bosons, so the singlet collision case discussed above is not automatically realized. The internal coin states are taken as the two hyperfine states , , split by the atomic clock transition at 9.2 GHz. Spin-dependent shift operations work, because the up and down state see complementary circular polarization components of the optical lattice formed by counterpropagating linearly polarized waves. Internal coin operations are Rabi rotations driven by microwaves.

On-site interactions via s-wave cold collisions can be accurately described using pseudopotential methods Stock and Deutsch (2006). The use of Cs atoms is especially favorable, because its large triplet scattering length allows fast interactions Leo et al. (2000). Roughly, given a 3D optical lattice with isotropic harmonic confinement kHz, the state acquires a collisional phase shift in about 10 s. This is just below the current step times of our experiment, so could well be incorporated. A further effect of the interaction is a level shift for , which may render the coin operations at 9.2 GHz ineffective at collision points, thus suppressing the transitions to , . In this way a singlet collision case could also be realized for Bosons. For the collision phase to be well controlled and coherent, it is necessary to cool the motional states to the ground state not only along the axis but also in the lateral direction, in which at the moment the confinement is much weaker. We plan to enhance the lateral confinement by an additional blue-detuned doughnut-mode laser. Then cooling to average vibrational quantum numbers in all three directions seems feasible, which should allow sufficiently many coherent steps to clearly demonstrate the molecule binding effect.

### Acknowledgments.—

Both groups acknowledge the support by the DFG (grant Forschergruppe For 635), and the EU projects CORNER, COQUIT (Hannover) and AQUTE (Bonn). A. Alberti acknowledges support by the A. v. Humboldt Foundation.

## Appendix A: Molecules in higher dimension

Generally speaking, a point perturbation has less and less effect in higher dimension. One therefore expects the formation of molecules to become more difficult in higher dimension. This intuition is basically correct, but since points on a lattice have finite thickness, there is no critical dimension above which molecule formation disappears.

From the analysis given in the paper a sufficient condition is clear: Assume that, for some given total momentum , the walk in the difference coordinate has a band gap, i.e., there is a on the unit circle so that is invertible for all . Choosing (see (6)) we then get a walk with an eigenvalue of maximal degeneracy precisely at . By perturbation theory this will typically split into a collection of bound state in a neighborhood of . This gives a (possibly small) range of and suitable initial preparations, which will show the molecule behavior.

The band gap condition can be met by constructing the free walk from a product (not tensor product) of “slow” walks in the coordinate directions. Here a slow walk is defined as one with nearly constant eigenvalues . For example, consider a walk of the form

By choosing appropriate unitary matrices , a large variety of -dimensional walks can be implemented. However, when we take all , we get a rather trivial walk which only has steps along the main diagonal and, for all , only the eigenvalues . Then if we take all close to the identity, we get a non-trivial walk in dimensions, whose eigenvalues are and . The same is true for the non-interacting pair walk . Clearly, this now satisfies the band gap condition. Intuitively this construction slows down the group velocity to such an extent that the perturbation at a single point can take effect.

It should be noted, however, that the band gap condition is by no means necessary. For example, if we always take time steps together, we are looking at the walk , whose spectrum, on the one hand, will still show molecule states if does. On the other hand, any bit of continuous spectrum will be spread out by a factor , and eventually wraps around the circle, maybe with multiplicity . The molecule states will thus be embedded in the continuum.

In the continuous-time version of the problem Winkler et al. (2006) the band gap condition is always satisfied , since the lattice kinetic energy is bounded, and a sufficiently strong point interaction can always create a bound state. Hence, as already remarked in Winkler et al. (2006), suitably constructed point interactions lead to molecule formation even in high dimension.

## Appendix B: Interaction in Quantum Cellular automata

Unlike the continuous time case, where one can just add Hamiltonian interaction terms, the interaction unitaries in the discrete time case have to be multiplied in some order. This leads to ambiguities for triple collisions. The way to treat these is to consider the many-particle system as a “lattice gas”, a special kind of quantum cellular automaton(QCA) Schumacher and Werner (). This is analogous to a lattice spin system with the spin up/down states interpreted as occupied/unoccupied. The special feature of a lattice gas (as opposed to a general QCA) is that there is a notion of particle number for each cell, and that the local dynamical rules preserve particle number. The characteristic change of perspective is to focus not so much on the different particles, but on what can happen in a single cell.

One standard transition from single particle walk to a non-interacting lattice gas is Bose/Fermi second quantization. Since for every basis state of the one-particle space we have a number operator , we can define the number of particles in the cell as . In the Fermi case this is clearly bounded by , but is unbounded in the Bose case. This is not quite consistent with the idea of finite volume atoms in a potential well, nor with the usual assumption for QCAs that there should be only finitely many states per cell. Indeed the QCA’s cell finiteness can be seen as a hard core repulsion condition, so a QCA tends to impose some interaction. There is a procedure for systematically introducing an interaction for a walk which is given in shift&coin factorization Vogts (2009). We illustrate it for the factorization (1) of the Hadamard walk.

The QCA system consists of a double chain of qubits, which we may call the right moving chain and the left moving chain. This describes the action of the shift operation: In the shift step the entire right (resp. left) moving chain is shifted on step to the right (resp. left). The coin operation acts at each site. We denote the basis states for all qubits by for empty and for occupied, so at each site we have the basis vectors . In order to guarantee particle number conservation, the coin has to act in block diagonal form on the number subspaces. The empty space will just be left invariant. There is no loss of generality in this choice, because we are free to choose a global phase. On the basis vectors and of the one particle space we just use the coin operation of the free walk. Then it is easy to see that the subspace of total particle number is invariant, and isomorphic to the free one-particle walk. The only free parameter in this construction is the phase for the two-particle state . Again, it is easy to check that on the -particle subspace we exactly get the singlet collision case discussed in the paper. However, the construction also extends to arbitrarily many particles.

This construction can be extended to other walks with an explicit decomposition into shift and coin operations (allowing also for many factors). Since such decompositions are not unique there are some trade-offs possibleVogts (2009): For example, one can use a decomposition in which only one if the states is ever shifted. The multiple chain picture then only has a qubit chain for transport, and a stationary chain, with an empty state and, for particle number , the remaining states of the one-particle walk. This leads to a QCA with states per cell. However, the interaction neighborhood size of this QCA may be much larger than that of the one-particle walk. In contrast, taking a moving chain for every one of the states suggests a picture with states per cell, and in this case (at least in 1D) the neighborhood can be arranged to be the same as for the walk. Further options are introduced by allowing more than single occupation of the one-particle coin states, as in the free Bose walk.

Which of these many possibilities is realized is, of course, a question that has to be answered by analyzing the physics of the given system. The point to be made here is only that the language of QCAs makes clear what physical information is required to specify an “interaction” completely.

## Appendix C: A Unitarity Lemma

In this appendix we prove the claim that the operator appearing in (6) is unitary. We will do this by proving the following, more general claim:

Consider a unitary operator on a Hilbert space and a point on the unit circle (, ) not in the spectrum of , which just means that the inverse exists and is a bounded operator. Now let be an arbitrary projection and set , considered as an operator on . We claim that this operator is invertible and that is unitary.

The crucial step for the proof is to show the following relation

(11) |

We can show this for the special case , and then multiply the above equation with from both sides. Since is unitary it has a spectral decomposition , where is the spectral measure on the unit circle. We can then express the operator (11) directly in the functional calculus (we use for hermitian adjoint and an overbar for complex conjugate):

This determines the hermitian part of , and we can write for some hermitian operator on .

For the next we use the functional calculus of , with the spectral measure of , which is now supported on the real axis. We note that due to the projection , which need not commute with , this measure cannot be easily obtained from the spectral measure . Then

Now, since the integrand has absolute value , and so the integral represents a unitary operator.

## Appendix D: Overlap with bound states

The time evolution of a two-particle quantum state under the considered quantum walk is essentially determined from its overlap with the bound states in Fourier space. This is also important for experimental implementations since it describes how efficient some preparation, e.g., “both particles in singlet state at the origin”, is for generating molecules. In other words, we want to estimate the splitting in Fig. 1 between the total probability near the diagonal and the rest, i.e., the probability for ballistically moving away from the diagonal. Since the total momentum is a conserved quantity we omit the -dependence of in the following discussion.

First, let us calculate the normalization of the eigenstate of the quantum walk corresponding to eigenvalue as determined from (4). Abbreviating we get

(12) |

where we used the definition of from equation (5) and the unitarity of . This yields now the probability to capture the walker initially prepared in an internal state located at the origin in the bound state

(13) |

since respects (6). We want to evaluate (13) for our standard example, the Hadamard Walk with singlet collisions and the singlet state as the initial state, where is given by . By (Appendix D: Overlap with bound states) we have to calculate the reciprocal of

where we used (6) and the convention in the last step. Defining via one can deduce from (7) the relation

(14) |

Hence, we get for the capture probability in the bound state corresponding to the eigenvalue the expression

(15) |

Note that this is a function of the total momentum since as well as depend on it. Analyzing this dependency, we see that as approaches a band edge, tends to zero and therefore such values are suppressed as Fourier components of the bound state. That behavior can also be seen in Fig. 4.

For fixed interaction phase the ratio between the occurrence of molecules and unbound two particle states in an experiment is determined from the integral of with respect to the total momentum . As can be seen from Fig. 5 there is a fairly high probability to observe molecules with a maximum of at .

Expression (15) also allows us to study the long time behavior of the molecule’s position distribution scaled by , where is the discrete time parameter. As outlined in corollary 7 in Ahlbrecht et al. (2011), this asymptotic distribution can be evaluated with knowledge of the group velocity (9) and the capture probability (15). The result is shown in Fig. 6.

## Appendix E: Analytic form of bound states

The goal of this section is to derive an explicit formula for the bound states of a quantum walk with a point defect at the origin. We then apply this procedure to our standard example, the Hadamard walk with singlet collisions. This allows us to derive an explicit form of the dependent transformation needed to construct the virtual walk (8) of the molecule as explained in the main part of the paper.

The components of the unnormalized eigenvector corresponding to lattice site can be determined from (4), with the substitution the formula reads

(16) | |||||

Here, we have to substitute the correct form (7) of the eigenvalue . The initial state in our example is the singlet state which effectively reduces the problem to the Fermi sector. Of course, the Fermi symmetry is reflected in the eigenstate . If we denote the operator which exchanges tensor factors by , that is, for all vectors , then and hence for all .

In order to determine for negative we use residual calculus which requires knowledge of the poles of the integral kernel in (16). The operator may have singularities at the zeros of the quadratic polynomial and possibly at . Conjugation of by yields the equality , which implies that the singularities of are either of first order and inverses with respect to each other or is constant, in which case the operator has a pole of first order at . Analyzing the case where is not constant we get the following expression for the singularity with modulus smaller than one, denoted by ,

(17) | |||||

with as in (14). Clearly, and if only contributes to the integral (16). Since the integral kernel is analytic at it remains to calculate the component wise residual operator

at . The coefficients of the normalized eigenstate in terms of this operator read

Note that since the modulus of is smaller than one the state decays exponentially in . This implies that the distance of the two particles in the molecule state is exponentially concentrated around zero.

In the case that the singularities of are at the eigenstate is strictly localized on a finite set of lattice sites. This is because for the integral kernel in (16) is analytic if , hence for these cases. For our example this happens exactly at . This property leads to an interesting feature of the considered quantum walk: By engineering the initial state of the quantum walk on few lattice sites it is possible to generate states which are close to a momentum eigenstate at points . Here, the spread of the initial state in the relative coordinate is bounded by the spread of the eigenstate and in the center of mass direction by the desired accuracy of the momentum preparation.

With the help of the eigenstates corresponding to the two branches we can construct the one dimensional quantum walk mimicking the time evolution of the molecule explicitly. This is done by identifying the states with the dependent eigenvectors of some one dimensional quantum walk with the same dispersion relation .

## Appendix F: Fast Molecules

So far only the two particle Hadamard Walk with singlet collisions was considered as explicit exemplification of the ideas of this paper. From Fig. 1 and Fig. 6 one could think that the bound states are slower or at most as fast as the free walkers.

In this section we will give two numerical examples to show that this is not true in general and that it is possible to design interactions that allow for molecules spreading faster then the free walk. Since in the case of singlet collisions the only freedom lies in the selection of the interaction phase we turn to the Bose case where we can choose an arbitrary three by three unitary matrix in symmetric subspace.

A trivial possibility to generate such fast molecules is choose an interaction matrix that counteracts the coin operation. In that case the walk operator for the two particle Hadamard walk is given by

which acts on the walkers only with the shift at the collision point. If we now prepare the particles at the same lattice side for example in the state they will be just shifted to the left with maximal velocity.

Two nontrivial examples of fast molecules are given in Fig. 7, were we have choosen two unitary interactions for which the bound state in the symmetric subspace spreads faster than the individual particles of the free quantum walk. Moreover it is even possible to design the interaction in such a way that the resulting molecule travels nearly at the maximal possible velocity as in the trivial case.

## References

- Karski et al. (2009) M. Karski, L. Förster, J.-M. Choi, A. Steffen, W. Alt, D. Meschede, and A. Widera, Science 325, 174 (2009).
- Schmitz et al. (2009) H. Schmitz, R. Matjeschk, C. Schneider, J. Glueckert, M. Enderlein, T. Huber, and T. Schaetz, Phys. Rev. Lett. 103, 090504 (2009).
- Zähringer et al. (2010) F. Zähringer, G. Kirchmair, R. Gerritsma, E. Solano, R. Blatt, and C. F. Roos, Phys. Rev. Lett. 104, 100503 (2010).
- Broome et al. (2010) M. A. Broome, A. Fedrizzi, B. P. Lanyon, I. Kassal, A. Aspuru-Guzik, and A. G. White, Phys. Rev. Lett. 104, 153602 (2010).
- Peruzzo et al. (2010) A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, K. Wörhoff, Y. Bromberg, Y. Silberberg, M. G. Thompson, and J. L. OBrien, Science 329, 1500 (2010).
- Bromberg et al. (2009) Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silberberg, Phys. Rev. Lett. 102, 253904 (2009).
- Schreiber et al. (2010) A. Schreiber, K. N. Cassemiro, V. Potoček, A. Gábris, P. J. Mosley, E. Andersson, I. Jex, and C. Silberhorn, Phys. Rev. Lett. 104, 050502 (2010).
- Winkler et al. (2006) K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. H. Denschlag, A. J. Daley, A. Kantian, H. P. Büchler, and P. Zoller, Nature 441, 853 (2006).
- Štefaňák et al. (2011) M. Štefaňák, S. M. Barnett, B. Kollár, T. Kiss, and I. Jex, New J. Phys. 13, 033029 (2011).
- Inui et al. (2004) N. Inui, Y. Konishi, and N. Konno, Phys. Rev. A 69, 052323 (2004).
- Omar et al. (2006) Y. Omar, N. Paunković, L. Sheridan, and S. Bose, Phys. Rev. A 74, 042304 (2006).
- Ahlbrecht et al. (2011) A. Ahlbrecht, H. Vogts, A. Werner, and R. Werner, J. Math. Phys. 52, 042201 (2011).
- Kato (1976) T. Kato, Perturbation theory for linear operators, 2nd ed. (Springer, 1976).
- Stock and Deutsch (2006) R. Stock and I. H. Deutsch, Phys. Rev. A 73, 032701 (2006).
- Leo et al. (2000) P. J. Leo, C. J. Williams, and P. S. Julienne, Phys. Rev. Lett. 85, 2721 (2000).
- (16) B. Schumacher and R. F. Werner, “Reversible quantum cellular automata,” arXiv:quant-ph/0405174.
- Vogts (2009) H. Vogts, “Discrete time quantum lattice systems,” Doctoral Thesis, Technische Universität Braunschweig (2009).