# Bloch oscillations in lattice potentials with controlled aperiodicity

###### Abstract

We numerically investigate the damping of Bloch oscillations in a one-dimensional lattice potential whose translational symmetry is broken in a systematic manner, either by making the potential bichromatic or by introducing scatterers at distinct lattice sites. We find that the damping strongly depends on the ratio of lattice constants in the bichromatic potential and that even a small concentration of scatterers can lead to strong damping. Moreover, collisional interparticle interactions are able to counteract aperiodicity-induced damping of Bloch oscillations. The discussed effects should readily be observable for ultracold atoms in optical lattices.

###### pacs:

03.75.-b, 67.85.Hj, 72.10.Fk^{†}

^{†}thanks: present address: Dept. of Physics, Univ. Würzburg, Germany

^{†}

^{†}thanks: present address: Photon Research Associates, Port Jefferson, NY

## I Introduction

The oscillatory motion of particles in a periodic potential, when subject to an external force, was predicted by Felix Bloch in 1928 Bloch (1928). Bloch oscillations were first observed in the 1990s both in semiconductor superlattices Waschke et al. (1993) and in systems of laser-cooled atoms in optical lattices Wilkinson et al. (1996); Ben Dahan et al. (1996). Since then, they have also been studied in atomic quantum gases Morsch et al. (2001); Roati et al. (2004); Gustavsson et al. (2008); Fattori et al. (2008) and have found applications in ultracold atomic-physics-based precision measurements Battesti et al. (2004); Carusotto et al. (2005); Ferrari et al. (2006); Cladé et al. (2006).

Bloch oscillations are absent in solids due to fast damping from scattering by defects and phonons. Their observation in long-period semiconductor superlattices relies on oscillation periods that are shorter than the characteristic scattering lifetime, and even in those systems, damping of Bloch oscillations due to disorder is ubiquitous von Plessen et al. (1994); Diez et al. (1998). In contrast, such damping is absent in optical lattice systems which are inherently defect-free, allowing for the observation of a large number of oscillations Roati et al. (2004); Gustavsson et al. (2008); Fattori et al. (2008). Damping can however be induced by the introduction of disorder into the lattice potential, which can in principle be done with various techniques Horak et al. (1998); Diener et al. (2001); Gavish and Castin (2005); Ospelkaus et al. (2006); Wang et al. (2004). Also, in the case of quantum gases, a damping of Bloch oscillations arises from mean-field interactions between weakly interacting, Bose-condensed atoms Wu and Niu (2001); Holthaus (2000); Witthaut et al. (2005); Menotti et al. (2003); Modugno et al. (2004); Gustavsson et al. (2008); Fattori et al. (2008). In the presence of both interactions and disorder, a reduction of disorder-induced damping due to screening of disorder by the mean field has been predicted Schulte et al. (2008). Experiments with ultracold atoms thus not only constitute a versatile testbed for behavior expected in solid-state systems but may also display novel effects.

Recently, the damping of Bloch oscillations in a Bose-Einstein condensate Drenkelforth et al. (2008) has been observed for an optical lattice with a superimposed randomly corrugated optical field. A related theoretical investigation of disorder-induced damping of Bloch oscillations has considered the case of Gaussian spatial noise Schulte et al. (2008). In this paper, we numerically investigate the dynamics of an atomic wave packet in a potential with a systematically degraded translational symmetry, considering two scenarios: The first is based on the use of a weak bichromatic potential Diener et al. (2001) of a variable wavelength ratio, and the second considers scatterers (impurities) pinned at single sites of the potential Gavish and Castin (2005). We find that the damping strongly depends on the ratio of lattice constants in the bichromatic case and that even a small concentration of scatterers can lead to strong damping. We also include effects of the mean-field interaction and find that the rate at which damping of the Bloch oscillations occurs is reduced, similar to the case of Gaussian disorder Schulte et al. (2008). Both effects should be observable experimentally with existing ultracold-atom technology.

This paper is organized as follows: After a brief discussion of fundamental aspects of tilted lattices in Sec. II, Section III investigates the influence of aperiodicity on the damping dynamics of Bloch oscillations. Sec. IV addresses the interplay of aperiodicity and the interaction between atoms. Conclusions are given in Sec. V.

## Ii Tilted periodic potentials

The Hamiltonian for the motion of a particle in a one-dimensional periodic potential with lattice constant possesses a complete set of eigenfunctions that obey Bloch’s theorem . The corresponding eigenvalues are periodic in momentum space and form energy bands, (where is the band index, the quasimomentum, and the width of the first Brillouin zone). Under the influence of an externally applied constant force (“tilt”) , the quasimomentum evolves as . Due to the periodicity of the energy bands, this results in oscillations of the particle’s group velocity . These oscillations occur with a period and a maximum displacement in coordinate space, where is the width of the -th band.

In the following, we consider particles that are confined to the lowest Bloch band, =1, corresponding to sufficiently deep potentials and small enough tilts such that Zener tunneling Zener (1934) at the band edges is negligible. Using the split-operator method Feit et al. (1982), we perform numerical simulations of the dynamics of Bloch oscillations of a Gaussian wave packet

(1) |

that evolves according to the Hamiltonian

(2) |

for motion in a tilted periodic potential that is modified by a weak additional potential .

## Iii Effects of modified periodicity

This section investigates the influence of the additional potential , which either is a weak additional periodic potential with variable lattice constant, or arises from the local interaction with scatterers pinned at single sites of the tilted periodic potential. In the absence of , the energies of neighboring sites differ by a fixed amount , corresponding to a spatially homogeneous phase difference . In the presence of , generally becomes position-dependent, leading to global dephasing and thus to a broadening of the wave packet in momentum space, as illustrated in Fig. 1.

### iii.1 Bichromatic potentials

A tunable bichromatic potential is generated by the addition of

(3) |

with variable relative amplitude and lattice-constant ratio . If is a rational number , the total potential has a periodicity . If furthermore exceeds the spatial range probed by the wave packet, the potential can be considered disordered.

The evolution of the wave packet (Eq. 1) typically exhibits collapses and revivals of center-of-mass oscillations that are coupled to breathing-mode excitations at twice the Bloch frequency as is shown in Fig. 2 (visible after a transient phase immediately following the switch-on of ). The presence of a collapse and a revival of the Bloch oscillations is a consequence of the absence of dissipation in our model. To characterize the decay of Bloch oscillations, the numerical data for their initial decay, for a given and , are fitted with the function

(4) |

Results for the Gaussian-decay constant are shown in Fig. 3.

Clearly an increase in the perturbation amplitude leads to an overall increase in the damping of the oscillations. However, the dependence of the damping on the ratio is less trivial. For integer values of , the potential retains its original periodicity and no dephasing of Bloch oscillations occurs. It is also suppressed for half-integer values of (i.e. for ), which correspond to a doubling of the lattice constant, and correspondingly a halving of the Brillouin zone, and of the corresponding Bloch period. For this case, the dynamics of the wave packet in the modified band structure can easily be visualized. The shape of the band is essentially that of , but it is folded back into the new Brillouin zone, thus forming a closed loop, with a small splitting at the new zone boundaries (cf. Fig. 4). In one Bloch cycle, most of the wave function tunnels through the tiny gap from the lower to the upper portion at one boundary and then back to the lower portion at the other boundary (Bloch-Zener tunneling Breid et al. (2006, 2007)). The time needed for one such cycle is , that is the Bloch period for the unperturbed potential . This behavior can directly be seen in the time dependence of the overlap of the wave function with the initial wave packet , which has a periodicity of ; contributions at odd multiples of only grow very slowly over a large number of cycles.

The dynamics of Bloch oscillations in Figs. 2 and 3 are for a wave packet with width in a lattice with depth [where is the recoil energy], and a small external force . In coordinate space the amplitude of the Bloch oscillations is . The additional lattice is chosen to have a small relative amplitude , which is sufficient to cause noticeable damping already after a few Bloch cycles, depending on the lattice constant ratio .

Experimentally, such a bichromatic potential is straightforward to realize in the context of optical lattices, using two laser beams of different wavelength. With present-day tunable-laser technology, a large fraction of the range of shown in Fig. 3 can be accessed.

### iii.2 Scatterers on distinct sites

For scatterers pinned at a set of distinct sites of the lattice, we model the potential as a sum of Gaussians,

(5) |

For the simulation, the parameters of the optical lattice and the tilt are chosen such that the size of the wave packet and its oscillation amplitude in the unperturbed titled potential each cover a large number of lattice sites. A small number of scatterers are then randomly placed within the range of the wave packet’s motion. The amplitude and width of the scatterers are chosen such that the valleys of the unperturbed potential are effectively filled up where the scatterers are located. Their effect on the Bloch oscillations is shown in Fig. 5. The general trend with an increasing number of scatterers is a more rapid damping of the Bloch oscillations, with the details of the dynamics depending on their spatial arrangement. Clearly, already a small number of scatterers can lead to a rapid damping of Bloch oscillations.

The results shown in Fig. 5 are for a comparatively shallow lattice () and a tilting force , resulting in a large amplitude for undamped Bloch oscillations. This condition for the scatterers is well fulfilled by setting and .

In the context of ultracold atoms, the placement of scatterers pinned at single lattice sites can be achieved, for example by using atoms with two internal states in conjunction with a state-dependent lattice depth et al (2009), or by using two atomic species in a species-dependent optical lattice Catani et al. (2009).

## Iv Effect of Interactions

The discussion so far has been restricted to noninteracting particles. We now consider the case of a Bose-Einstein condensate of atoms with repulsive interparticle interaction . The evolution of the condensate wave function in the bichromatic potential is then determined by the mean-field Hamiltonian

(6) |

The simulation results (Fig. 6) show that an increase in the coupling strength leads to a reduction of the damping constant , up to a characteristic value , beyond which an increase in leads to an increase in . We interpret the reduction of as a partial screening of the potential corrugations by the mean field Schulte et al. (2008) as increases, which is eventually overcompensated by mean-field-induced dephasing Holthaus (2000); Witthaut et al. (2005).

The parameters for the simulation are the same as those in Sec. 3.1, leading to an associated coupling strength . This value should be compared with an effective one-dimensional interaction parameter for a trapped atomic Bose-Einstein condensate, where (with atomic s-wave scattering length and atomic mass ), and , where is the Thomas-Fermi radius. Already for a small condensate of Rb atoms with atoms and m in an isotropic 50 Hz trap Pertot et al. (2009), and an optical lattice with nm, we obtain , which is in the vicinity of . Hence, a significant modification of the damping rate in a bichromatic potential due to mean-field effects can be expected; the coupling depends, for example on the atom number in the condensate, which is variable. Alternatively, an investigation of the interplay is possible for species in which the mean-field interaction can be tuned via a Feshbach resonance. In this context we mention that Bloch oscillations with widely controllable mean-field interactions in a (monochromatic) optical lattice have recently been demonstrated with cesium condensates Gustavsson et al. (2008); Fattori et al. (2008).

## V Conclusions

We have numerically investigated the damping of Bloch oscillations resulting from a controlled breakdown of the periodicity of the lattice potential. The effects discussed here, including the effects of the mean-field interaction on disorder-induced dephasing, should be readily observable in experiments with ultracold atoms in optical lattices.

###### Acknowledgements.

We thank D. Pertot and B. Gadway for a critical reading of the manuscript. This work is supported by NSF Grant Nos DMR-0605919 (S.W. and A.C.D) and PHY-0855643 (D.S.) as well as a DAAD scholarship (S.W.).## References

- Bloch (1928) F. Bloch, Z. Phys. A 52, 555 (1929).
- Waschke et al. (1993) C. Waschke, H. G. Roskos, R. Schwedler, K. Leo, H. Kurz, and K. Köhler, Phys. Rev. Lett. 70, 3319 (1993).
- Wilkinson et al. (1996) S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Q. Niu, and M. G. Raizen, Phys. Rev. Lett. 76, 4512 (1996).
- Ben Dahan et al. (1996) M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, Phys. Rev. Lett. 76, 4508 (1996).
- Morsch et al. (2001) O. Morsch, J. H. Müller, M. Cristiani, D. Ciampini, and E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001).
- Roati et al. (2004) G. Roati, E. de Mirandes, F. Ferlaino, H. Ott, G. Modugno, and M. Inguscio, Phys. Rev. Lett. 92, 230402 (2004).
- Gustavsson et al. (2008) M. Gustavsson, E. Haller, M. J. Mark, J. G. Danzl, G. Rojas-Kopeinig, and H.-C. Nägerl, Phys. Rev. Lett. 100, 080404 (2008).
- Fattori et al. (2008) M. Fattori, C. D’Errico, G. Roati, M. Zaccanti, M. Jona-Lasinio, M. Modugno, M. Inguscio, and G. Modugno, Phys. Rev. Lett. 100, 080405 (2008).
- Battesti et al. (2004) R. Battesti, P. Cladé, S. Guellati-Khélifa, C. Schwob, B. Grémaud, F. Nez, L. Julien, and F. Biraben, Phys. Rev. Lett. 92, 253001 (2004).
- Carusotto et al. (2005) I. Carusotto, L. Pitaevskii, S. Stringari, G. Modugno, and M. Inguscio, Phys. Rev. Lett. 95, 093202 (2005).
- Ferrari et al. (2006) G. Ferrari, N. Poli, F. Sorrentino, and G. M. Tino, Phys. Rev. Lett. 97, 060402 (2006).
- Cladé et al. (2006) P. Cladé, E. de Mirandes, M. Cadoret, S. Guellati-Khélifa, C. Schwob, F. Nez, L. Julien, and F. Biraben, Phys. Rev. Lett. 96, 033001 (2006).
- von Plessen et al. (1994) G. von Plessen, T. Meier, J. Feldmann, E. O. Göbel, P. Thomas, K. W. Goossen, J. M. Kuo, and R. F. Kopf, Phys. Rev. B 49, 14058 (1994).
- Diez et al. (1998) E. Diez, F. Domínguez-Adame, and A. Sánchez, Microelectron. Eng. 43-44, 117 (1998).
- Horak et al. (1998) P. Horak, J.-Y. Courtois, and G. Grynberg, Phys. Rev. A 58, 3953 (1998).
- Diener et al. (2001) R. B. Diener, G. A. Georgakis, J. Zhong, M. Raizen, and Q. Niu, Phys. Rev. A 64, 033416 (2001).
- Gavish and Castin (2005) U. Gavish and Y. Castin, Phys. Rev. Lett. 95, 020401 (2005).
- Ospelkaus et al. (2006) S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, K. Sengstock, and K. Bongs, Phys. Rev. Lett. 96, 180403 (2006).
- Wang et al. (2004) D.-W. Wang, M. D. Lukin, and E. Demler, Phys. Rev. Lett. 92, 076802 (2004).
- Wu and Niu (2001) B. Wu and Q. Niu, Phys. Rev. A 64, 061603(R) (2001).
- Holthaus (2000) M. Holthaus, J. Opt. B: Quantum Semiclassical Opt. 2, 589 (2000).
- Witthaut et al. (2005) D. Witthaut, M. Werder, S. Mossmann, and H. J. Korsch, Phys. Rev. E 71, 036625 (2005).
- Menotti et al. (2003) C. Menotti, A. Smerzi, and A. Trombettoni, New J. Phys. 5, 112 (2003).
- Modugno et al. (2004) M. Modugno, C. Tozzo, and F. Dalfovo, Phys. Rev. A 70, 043625 (2004).
- Schulte et al. (2008) T. Schulte, S. Drenkelforth, G. K. Büning, W. Ertmer, J. Arlt, M. Lewenstein, and L. Santos, Phys. Rev. A 77, 023610 (2008).
- Drenkelforth et al. (2008) S. Drenkelforth, G. K. Buning, J. Will, T. Schulte, N. Murray, W. Ertmer, L. Santos, and J. J. Arlt, New J. Phys. 10, 045027 (2008).
- Zener (1934) C. Zener, Proc. R. Soc. Lond. Ser. A 145, 523 (1934).
- Feit et al. (1982) M. D. Feit, J. A. Fleck Jr., and A. Steiger, J. Comput. Phys. 47, 412 (1982).
- Breid et al. (2006) B. M. Breid, D. Witthaut, and H. J. Korsch, New Journal of Physics 8, 110 (2006).
- Breid et al. (2007) B. M. Breid, D. Witthaut, and H. J. Korsch, New Journal of Physics 9, 62 (2007).
- et al (2009) B. Gadway et al., in preparation (2009).
- Catani et al. (2009) J. Catani, G. Barontini, G. Lamporesi, F. Rabatti, G. Thalhammer, F. Minardi, S. Stringari, and M. Inguscio, Phys. Rev. Lett. 103, 140401 (2009).
- Pertot et al. (2009) D. Pertot, D. Greif, S. Albert, B. Gadway, and D. Schneble, J. Phys. B: At. Mol. Opt. Phys. 42, 215305 (2009).