# Quantum walks and quantum simulations with Bloch oscillating spinor atoms

###### Abstract

We propose a scheme for the realization of a quantum walker and a quantum simulator for the Dirac equation with ultracold spinor atoms in driven optical lattices. A precise control of the dynamics of the atomic matter wave can be realized using time-dependent external forces. If the force depends on the spin state of the atoms, the dynamics will entangle the inner and outer degrees of freedom which offers unique opportunities for quantum information and quantum simulation. Here, we introduce a method to realize a quantum walker based on the state-dependent transport of spinor atoms and a coherent driving of the internal state. In the limit of weak driving the dynamics is equivalent to that of a Dirac particle in 1+1 dimensions. Thus it becomes possible to simulate relativistic effects such as Zitterbewegung and Klein tunneling.

###### pacs:

03.75.Mn,03.67.Ac,03.65.Pm## I Introduction

Quantum simulators aim at the simulation of complex quantum systems in well controllable laboratory experiments Bulu09 (). Such a simulation is especially useful when the original quantum system is experimentally not accessible and classical simulations are impossible due to the exponential size of the Hilbert space. Furthermore, quantum simulators offer the possibility to tune the experimental parameters to explore novel physical phenomena. Important examples include the simulation of solid state systems with ultracold atoms Bloc08 (), Dirac dynamics with graphene Novo05 (); Kats06 () or trapped ions Lama07 (); Gerr10 () and sonic black holes in Bose-Einstein condensates Gara00 ().

Ultracold atoms in optical lattices are especially suited for such a task, since their dynamics can be controlled with an astonishing precision and their dynamics can be measured in situ. In the present paper we propose to use spinor atoms in tilted or driven optical lattices to realize a quantum walker – a paradigmatic system in quantum information science Ahar93 (); Kars09 (); Schm09 (); Schr10 (); Zahr10 () – and a quantum simulator for relativistic Dirac dynamics. In contrast to previous proposals which were based on the realization of artificial gauge fields using atoms with a tripod internal structure Zhu07 (); Otte09 (), we focus on simple spinor atoms with two internal states, which are used routinely in ongoing experiments. The necessary correlations between the internal and the external dynamics can be realized with state dependent external forces, in particular by magnetic gradient fields Gorl03 (); Gust08 (); Weld09 (); Bohi09 (); Hall10 (). Using a suitable driving of this external field and microwave transitions between the internal states Fors09 (); Misc10 (), one can obtain control over the full dynamics of the atoms.

## Ii Bloch oscillations and transport in optical lattices

A Bloch oscillation is the counter-intuitive dynamics of a quantum particle in a periodic potential subject to a static external force. The force accelerates the particle until it reaches the edge of the Brillouin zone, where it is Bragg reflected. If Landau-Zener tunneling to higher bands can be neglected, this leads to a fully periodic motion. This peculiar kind of dynamics was postulated by Bloch already in 1928 in the context of electrons in crystals Bloc28 () but never observed because of the strong scattering of the electrons. Thus it took another seven decades before Bloch oscillations could be demonstrated for electrons in semiconductor superlattices Feld92 (), photons in waveguide arrays Pert99 () and ultracold atoms in optical lattices Daha96 (). Especially the latter realization shows an astonishing level of coherence and offers various possibilities to precisely control the atomic dynamics. Recent experiments demonstrated Bloch oscillations over tenth of thousands of periods and macroscopic distances Gust08 (). An example of Bloch oscillations of an atomic matter wave is shown in Fig. 1 (a).

Extensive possibilities to control the atomic motion can be realized in driven optical lattices, i.e. in lattices with a time-dependent external field Hart04 (); Klum07 (); Brei07 (); Brei08 (); Sias08 (); Zene09 (); Hall10 (). A particular interesting case is a periodic driving, when the direction of the field is reversed before the matter wave reaches the edge of the Brillouin zone, decelerating the wave packet back to a standstill. Thus the atoms always have a positive momentum such that directed transport is realized, as illustrated in Fig. 1 (b) for a sinusoidal driving.

To be precise, we have simulated the dynamics of ultracold Cs atoms in an optical lattice with a wavelength of and a depth of . The strength of the external field was assumed to be , which can be easily realized by accelerating the optical lattice Daha96 () or by a magnetic gradient field Gust08 (); Hall10 (). The dynamics of the atomic matter wave is then given by the Wannier-Stark Hamiltonian

(1) |

where and is the atomic recoil energy. The initial state has been chosen as a Bloch state with zero momentum weighted by a Gaussian with a width . For a static field one then finds the celebrated Bloch oscillations with a period given by the Bloch time . The maximum displacement of the wave packet is given by , being the width of the ground Bloch band Hart04 ().

A directed transport of the atoms can be realized by means of a time-periodic driving. During one period of the driving the atoms are first accelerated and the decelerated back to a standstill, such that the atomic wavefunction is displaced in space but otherwise unaffected. An example for directed transport in a driven optical lattice is shown in Fig. 1 (b) for a sinusoidal driving

(2) |

assuming a driving strength of with a period of . The transport properties of a driven optical lattice can be calculated analytically within the tight-binding approximation Hart04 (); Klum07 (). It is found that transport is possible only for a pure ac-driving, while it is generally forbidden in the case of a combined dc- and ac-field except for the case of a resonant driving . The actual transport velocity depends on the initial state, but an upper bound can easily be found:

(3) |

where is an ordinary Bessel function. The case of a pure ac-driving corresponds to . The dispersion of a wave packet is generally negligible for an initially broad wave packet as it vanishes with Klum07 (); Brei07 (); Hall10 (). This approach to controlling transport in driven optical lattices was experimentally demonstrated in Sias08 (); Zene09 ().

## Iii Quantum Walks

Now consider the dynamics of a spinor atom, where the atoms experience a different field strength depending on their initial state. Using Bloch oscillations and the directed transport described above, it is possible to entangle the position and the internal degrees of freedom. This provides the basis for the implementation of a quantum walk or a Dirac quantum simulator.

One possibility to induce a state dependent transport of spinor atoms is the use of a magnetic gradient field. In particular we consider the dynamics of ultracold cesium atoms with the internal states and Gorl03 (). The magnetic moments of these two internal states are opposite, such that they move into different directions in a magnetic gradient field. The effective potential is given by

(4) |

assuming a linear magnetic gradient field . The effective Landé factor of the two internal states is given by and Arim77 () and denotes the Bohr magneton. The magnetic gradient field can also be varied in time to control the atomic motion Hall10 (). Furthermore, we assume that the internal state of the atoms can be manipulated coherently by resonant microwave pulses with Rabi frequency Fors09 (). The dynamics of the atomic state is thus given by the Hamiltonian

(5) | |||||

A coherent quantum walker can be implemented by a sinusoidal driving of the gradient field . During one period, the atom are first accelerated and then decelerated back to a standstill as illustrated in Fig. 1 (b). However, the transport direction depends on the direction of the external field, which is opposite for both internal states, thus realizing an effective quantum walk. The coin operation is then realized by a pulse coupling the two hyperfine levels after integer multiples of the driving time . An example of a quantum walk is shown in Fig. 2 (c), where we have plotted the dynamics of the complete atomic density . The time-dependence of the magnetic gradient field and the microwave pulses is illustrated in the upper panels of the figure. Initially, all atoms are assumed to be in the internal state , where the spatial wavefunction is a Bloch state with , weighted with a Gaussian envelope of width . The depth of the optical lattice is given by .

It has to be noted that the magnitude of the force is different for the two spin states leading to a trivial overall displacement of the atomic density pattern. Furthermore this introduces a phase shift between the two spinor components, which has more severe consequences. This phase shift can be avoided by choosing internal states with the same , i.e instead of , or it can be compensated actively by a static magnetic field.

## Iv Dirac Dynamics

For weak driving, the dynamics of a quantum walker is equivalent to that of a Dirac particle in 1+1 dimension Stra06 (), which can be seen as follows. During one period , the atoms are first displaced depending on their internal state and then the internal state is rotated by an angle . Assuming that the internal rotation is fast compared to , the dynamics is given by the evolution operator with

(6) |

where denotes the translation operator over a distance . In the limit of weak driving, i.e small values of and , the dynamics can be described by an effective Hamiltonian, , which is given by

(7) |

This is just the Dirac Hamiltonian with an effective mass and the momentum operator .

According to these results, ultracold spinor atoms show an effective relativistic dynamics if the displacement and the rotation angle are reduced. This is simply realized by increasing the lattice depth and decreasing the strength of the microwave driving field. An example is shown in Fig. 4 for atoms initially prepared in the internal state with a width . The optical lattice is now deeper as in the previous examples, . Transitions to the other spin state are induced at integer multiples of the period with a variable rotation angle . For , the internal state does not change such that the atoms are steadily transported by the time-periodic external force. Increasing and thus the effective mass of the Dirac particle leads to a Zitterbewegung of the atoms, which is illustrated in Fig. 4 (a). As expected, the period of the oscillations decreases with the effective mass . Microscopically, the oscillations of the position expectation value result from the interference of the two components of the atomic state, which are shown in Fig. 4 (b,c)

One of the most surprising predictions of Dirac theory is the Klein paradoxon: Relativistic particles are not repelled by a strong repulsive potential, but perfectly transmitted if the height of the potential step exceeds Thal92 (). Signatures of this effect can be readily explored with ultracold atoms in driven optical lattices. We consider a situation where an additional blue detuned tophat laser beam is focused aside of the atom cloud. We model the induced optical dipole potential by a tanh profile such that the total optical potential is given by

If the additional potential is weak enough, the directed transport mechanism remains mostly unaffected. However, the atomic wave packets accumulate a dynamical phase depending on their position which is equivalent to a potential in the effective Dirac hamiltonian (7). In this regime, the dynamics is described by the discrete evolution operator

(9) |

and given in Eqn. (6). The dynamics in the presence of an additional potential is shown in Fig. 5, comparing a full simulation to the Dirac approximation (9). A good agreement is observed for weak potential steps, .

A pronounced signature of Klein tunneling observed in Fig. 5 is that the repulsive potential step enhances the transmission of the atomic matter wave. Without the potential, the atoms exibit a Zitterbewegung as discussed above. If a weak additional potential step is included, e.g. , the oscillating motion is gone and the atoms move steadily above the barrier. An even stronger potential, however, suppresses the directed transport underlying the Dirac quantum simulator. The effective description by the Dirac hamiltonian (7) is no longer appropriate and the atoms are reflected just like ordinary Schrödinger particles. In this regime the Dirac approximation clearly deviates from the exact simulation results.

## V Decoherence

Real experiments with spinor atoms are limited by a loss of coherence due a coupling to the environment. The schemes discussed above make use on the opposite magnetic moment of the two hyperfine states and to transport the atoms in different directions. However, this makes the atomic state vulnerable to dephasing caused by magnetic field fluctuations Kars09 (), while spatial coherence can be conserved over seconds Gust08 (). To explore the effects of decoherence, we simulate the dynamics in the presence of pure dephasing described by the Master equation

(10) |

For the numerical calculations we use the quantum jump method Dali92 (), averaging over 100 trajectories.

Decoherence of the atomic states essentially turns a quantum random walk into a classical one, which is confirmed by the simulations shown in Fig. 6. After a short initial period the standard deviation grows much slower in the the presence of dephasing. The linear increase gradually changes into the classical diffusion law if is increased. The different scaling is analyzed quantitatively in part (b), where is plotted on a loglog scale together with a linear fit, omitting the initial expansion of the atomic matter wave. The different slope of the two curves directly reveals the different diffusion exponents. The fit yields an exponent of and for and , respectively, which is in good agreement with the theoretical expectation for a quantum and for a classical random walk. A quantitative analysis of the decoherence process is provided in Fig. 6 (c), where the standard deviation after a fixed propagation time is plotted as function of . One observes that decreases with until the random walk becomes completely classical at .

Also the Dirac quantum simulator discussed in the previous section is vulnerable to decoherence of the atomic hyperfine states. Figure 7 shows an example of the atomic Zitterbewegung for the same parameters as in Fig. 4 for different values of the dephasing rate . The oscillations of the mean position of the atoms are much less pronounced for . We find that the dephasing rate must not exceed to observe a clear signature of a Zitterbewegung. For stronger dephasing no interference effects and thus no oscillations are visible.

## Vi Conclusion

In the present paper we have explored the rich dynamics of spinor atoms in tilted and driven optical lattices. These systems are nowadays routinely realized experimentally and can be controlled with a high accuracy. We have discussed how a state-dependent transport can be realized by a periodic driving and how this can be used to implement a quantum simulator. A quantum random walk is realized when transitions between the internal states are driven by microwave or Raman pulses. In the continuum limit of weak driving, the dynamics is given by the Dirac equation in 1+1 dimension. Thus it is possible to investigate relativistic effects such as Zitterbewegung or Klein tunneling in a well controllable laboratory experiment.

All elements needed for the implementation of the proposed quantum simulator have been demonstrated quite recently such that an experimental realization is possible with current technology. A major obstacle can be decoherence of the internal atomic state. On the other hand, decoherence processes offer the possibility to experimentally study the quantum classical transition for a Dirac particle.

###### Acknowledgements.

This work has been supported by the German Research Foundation (DFG) through the research fellowship program (grant no WI 3415/1). I thank F. Trimborn and H. J. Korsch for valuable comments.## References

- (1) I. Bulutaand and F. Nori, Science 326, 106 (2009).
- (2) I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
- (3) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005).
- (4) M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. Physics 2, 620 (2006).
- (5) L. Lamata, J. León, T. Schätz, and E. Solano, Phys. Rev. Lett. 98, 253005 (2007).
- (6) R. Gerritsma, F. Zähringer, E. Solano, R. Blatt, and C. F. Roos, Nature 463, 68 (2010).
- (7) L. J. Garay, J. R. Anglin, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 85, 4643 (2000).
- (8) Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. A 48, 1687 (1993).
- (9) M. Karski, L. Förster, J.-M. Choi, A. Steffen, W. Alt, D. Meschede, and A. Widera, Science 325, 174 (2009).
- (10) H. Schmitz, R. Matjeschk, Ch. Schneider, J. Glueckert, M. Enderlein, T. Huber, and T. Schaetz, Phys. Rev. Lett. 103, 090504 (2009) .
- (11) A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, P. J. Mosley, E. Andersson, I. Jex, and Ch. Silberhorn, Phys. Rev. Lett 104, 050502 (2010).
- (12) F. Zähringer, G. Kirchmair, R. Gerritsma, E. Solano, R. Blatt, and C. F. Roos, Phys. Rev. Lett 104, 100503 (2010).
- (13) S.-L. Zhu, B. Wang, and L.-M. Duan, Phys. Rev. Lett. 98, 260402 (2007).
- (14) J. Otterbach, R. G. Unanyan, and M. Fleischhauer, Phys. Rev. Lett. 102, 063602 (2009).
- (15) E. Haller, R. Hart, M. J. Mark, J. G. Danzl, L. Reichsöllner, and H.-C. Nägerl, Phys. Rev. Lett. 104, 200403 (2010).
- (16) A. Görlitz, T. L. Gustavson, A. E. Leanhardt, R. Löw, A. P. Chikkatur, S. Gupta, S. Inouye, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 90, 090401 (2003).
- (17) M. Gustavsson, E. Haller, M. J. Mark, J. G. Danzl, G. Rojas-Kopeinig, and H.-C. Nägerl, Phys. Rev. Lett. 100, 080404 (2008).
- (18) D. M. Weld, P. Medley, H. Miyake, D. Hucul, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 103, 245301 (2009).
- (19) P. Böhi, M. F. Riedl, J. Hoffrogge, J. Reichl, T. W. Hänsch, and P. Treutlein, Nature Phys. 5, 592 (2009).
- (20) L. Förster, M. Karski, J.-M. Choi, A. Steffen, W. Alt, D. Meschede, A. Widera, E. Montano, J. H. Lee, W. Rakreungdet, and P. S. Jessen, Phys. Rev. Lett. 103, 233001 (2009).
- (21) B. Mischuck, I. H. Deutsch, and P. S. Jessen, Phys. Rev. A 81, 023403 (2010).
- (22) F. Bloch, Z. Phys. 52, 555 (1928).
- (23) J. Feldmann, K. Leo, J. Shah, D. A. B. Miller, J. E. Cunningham, T. Meier, G. von Plessen, A. Schulze, P. Thomas, and S. Schmitt-Rink, Phys. Rev. B 46, R7252 (1992).
- (24) T. Pertsch, P. Dannberg, W. Elflein, A. Bräuer, and F. Lederer, Phys. Rev. Lett. 83, 4752 (1999).
- (25) M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, Phys. Rev. Lett. 76, 4508 (1996).
- (26) T. Hartmann, F. Keck, H. J. Korsch, and S. Mossmann, New J. Phys. 6, 2 (2004).
- (27) A. Klumpp, D. Witthaut, and H. J. Korsch, J. Phys. A: Math. Theor. 40, 2299 (2007).
- (28) B. M. Breid, D. Witthaut, and H. J. Korsch, New J. Phys. 8, 110 (2006).
- (29) B. M. Breid, D. Witthaut, and H. J. Korsch, New J. Phys. 9, 62 (2007).
- (30) C. Sias, H. Lignier, Y. P. Singh, A. Zenesini, D. Ciampini, O. Morsch, and E. Arimondo, Phys. Rev. Lett. 100, 040404 (2008).
- (31) A. Zenesini, H. Lignier, D. Ciampini, O. Morsch, and E. Arimondo, Phys. Rev. Lett. 102, 100403 (2009).
- (32) E. Arimondo, M. Inguscio, and P. Violino, Rev. Mod. Phys. 49, 31 (1977).
- (33) F. W. Strauch, Phys. Rev. A 73, 054302 (2006).
- (34) B. Thaller, The Dirac eqation, Springer, Berlin Heidelberg New York (1992).
- (35) J. Dalibard, Y. Castin, and K. Mølmer, Phys. Rev. Lett. 68, 580 (1992).