A.I Bandstructure calculation

Experimental realization of strong effective magnetic fields in an optical lattice


We use Raman-assisted tunneling in an optical superlattice to generate large tunable effective magnetic fields for ultracold atoms. When hopping in the lattice, the accumulated phase shift by an atom is equivalent to the Aharonov-Bohm phase of a charged particle exposed to a staggered magnetic field of large magnitude, on the order of one flux quantum per plaquette. We study the ground state of this system and observe that the frustration induced by the magnetic field can lead to a degenerate ground state for non-interacting particles. We provide a measurement of the local phase acquired from Raman-induced tunneling, demonstrating time-reversal symmetry breaking of the underlying Hamiltonian. Furthermore, the quantum cyclotron orbit of single atoms in the lattice exposed to the magnetic field is directly revealed.

03.65.Vf, 03.75.Lm, 11.15.Ex, 73.20.-r

The application of strong magnetic fields to two-dimensional electron gases has led to the discovery of seminal quantum many-body phenomena, such as the integer and fractional quantum Hall effect (1). Ultracold atoms constitute a unique experimental system for studying such systems in a clean and well controlled environment and for exploring new physical regimes, not attainable in typical condensed matter systems (2); (3). However, charge neutrality of atoms prevents direct application of the Lorentz force with a magnetic field. An equivalent effect can be provided by the Coriolis force in a rotating atomic gas, which led to the observation of quantized vortices in a Bose-Einstein condensate (4). The regime of fast rotation, in which the atomic gas occupies the lowest Landau level, was achieved in Refs. (5) but the amplitude of the effective gauge field remained too small to enter the strongly correlated regime that requires a number of vortices on the order of the particle number (6); (2). An alternative route consists in applying Raman lasers to the gas in order to realize a Berry’s phase for a moving particle (7); (8). The effective gauge fields generated in such a setup resulted in the observation of a few vortices, but were still far from the strong-field regime.

In this Letter, we demonstrate the creation of strong effective magnetic fields for ultracold atoms in a two-dimensional optical lattice. Inspired by the proposal of Jaksch and Zoller (9) and subsequent work (10); (11); (12), our technique is based on atom tunneling assisted by Raman transitions [see Fig. 1(a)]. Due to the spatial variation of the Raman coupling, the wavefunction of an atom tunneling from one lattice site to another acquires a non-trivial phase, which can be interpreted as an effective Aharonov-Bohm phase. In our setup, the magnetic flux per four-site plaquette is staggered with a zero mean, alternating between and [see Fig. 1(b)] (13). We study the nature of the ground state in this optical lattice from its momentum distribution and show in particular that the frustration associated with the effective magnetic field can lead to a degenerate ground state for single particles, similar to the prediction of Ref. (14). We also study the quantum cyclotron dynamics of single atoms restricted to a four-site plaquette and obtain direct evidence for time-reversal symmetry breaking of the Hamiltonian.

Figure 1: Experimental setup. (a) The experiment consists of a 2D array of 1D potential tubes with spacing . While bare tunneling occurs along the direction with amplitude , it is inhibited along owing to a staggered potential offset . A pair of Raman lasers with wave vectors and frequency difference , induces a resonant tunnel coupling of magnitude whose phase depends on position. This realizes an effective flux per plaquette with alternating sign along . (b) Spatial distribution of the phase of the Raman-induced tunnel coupling realized in the experiment. The gray shaded area highlights the magnetic unit cell.
Figure 2: Momentum distribution measured after a time-of-flight of ms for (a) simple cubic lattice, (b) and (c) . The latter two are compared with theoretical profiles obtained by an exact diagonalization of Hamiltonian (1) on a lattice with harmonic confinement (15). Red squares in the theoretical profiles indicate the magnetic Brillouin zone and the crosses mark the center. (d) Projection along of the momentum peaks located at , as a function of . For the peaks are located at , while for the peaks are split due to the emergent ground-state degeneracy (see insets). The solid lines correspond to the minima of the lowest Bloch band for a translationally invariant system (15).

Our experimental setup consists of an ultracold gas of Rb atoms held in a two-dimensional square lattice, forming an array of 1D Bose gases. The lattice was created by two standing waves of laser light at nm (‘short’ lattices) and a third one with twice the wavelength (‘long’ lattice, ) to generate a staggered potential with amplitude as shown in Fig. 1(a). A pair of Raman lasers then induced tunneling along the staggered direction. Let us consider the Raman-assisted tunneling of an atom from a site of low energy at to a site of high energy at . Assuming , one obtains , where denotes the wave vector (15). The system is then effectively described by the non-interacting Hamiltonian


where the sign of the phase factor is positive for even sites of and negative otherwise.

The phase factors in can be interpreted as Aharonov-Bohm phases. For the propagation of the Raman beams shown in Fig. 1(a) along and and , we obtain a phase factor of . Therefore the phase accumulated on a closed path around a plaquette is equal to , alternating in sign along the direction. A different value of the flux could be achieved by choosing a different wavelength for the Raman lasers or by using a different angle between them, allowing for a fully tunable flux per plaquette in our setup.

Our experiment started by loading a Bose-Einstein condensate of about atoms into a staggered 2D optical lattice as shown in Fig. 1(a) with kHz, resulting in an array of tubes with no coherence along (15). We then switched on the Raman lasers on resonance with strength to restore the coherence. In the limit , the amplitude of Raman induced tunnel coupling is , with being the bare tunnel coupling along . For our experimental parameters, this yields a value of Hz, in agreement with an independent measurement of Hz (15). After holding the atoms in this configuration for  ms, we observed a momentum distribution with restored phase coherence as shown in Fig. 2(b). This can be attributed to a redistribution of entropy present in the random phases between the 1D condensates into their longitudinal modes.

To understand the momentum distribution, we calculated the bandstructure of this lattice in the tight-binding approximation according to Ref. (14). In the presence of the gauge field, the Hamiltonian remains periodic. However, the magnetic unit cell contains two non-equivalent sites, leading to a splitting of the tight-binding band structure into two subbands (15); (16); (17). The frustration introduced by the position-dependent phase factors in causes the phase of the atomic wavefunction to be non-uniform, leading to a condensation at non-zero momenta. In the case , we obtain a non-degenerate ground state. However, the wavefunction itself consists of two momentum components, therefore we observe two diffraction peaks within the first magnetic Brillouin zone, one being shifted by with respect to the minimum of the dispersion relation (15). Due to the unit cell containing more than one lattice site, the size of the magnetic Brillouin zone is reduced compared to the one of the square lattice [see Fig. 3(a)]. Therefore the ground-state momentum distribution exhibits several peaks in the short-lattice Brillouin zone, whose measured positions are in good agreement with the quasi-momenta of the Bloch states of lowest energy [see Fig. 2(b)].

Figure 3: Dispersion relation of the lowest Bloch band, calculated in the tight-binding approximation (15) for (a) and (b) and plotted for the first magnetic Brillouin zone. The red crosses mark the positions of the corresponding momentum peaks together with the ones shifted due to the intrinsic structure of the wavefunction. (c) Histogram of the measured fraction of atoms in peaks corresponding to the lower momentum state for . The measurement was performed over 172 identical experimental runs. (d)-(e) Spatial distribution of the phase and atomic density (color brightness) for the ground-state wave function (15). The vortices with different chirality in the phase distribution for (d) are illustrated by the rotation of the white arrows. While in this case the atomic density is uniform, it exhibits a charge density wave for (e). For the second degenerate ground state we observe similar behavior but the density pattern is shifted by one lattice site.

When varying the ratio by adjusting the -lattice depth, we observed the positions and number of peaks in the momentum distribution to remain unchanged for [see Fig. 2(d)]. Above this value, additional peaks appear [see Fig. 2(c)], which correspond to the population of two degenerate ground states. This behavior agrees with the bandstructure calculation, which shows a bifurcation at , above which the energy minimum is split into two degenerate ground states [see Fig. 3(b)]. We find the measured and predicted peak positions of these two ground states to be in good agreement [see Fig. 2(d)]. The nature of the bifurcation is identical to the one predicted in Ref. (14) where it is induced by a variation of the magnetic flux amplitude at a fixed value of . As shown in Fig. 3(d), the atomic density is uniform for , while in the case it is strongly modulated for both single-particle ground states [see Fig. 3(e)]. For the bare coupling dominates and the phases of the atomic wavefunction tend to align along the direction, thereby frustrating the phase relation imposed by the Hamiltonian. As a consequence, the density in every second stripe along is suppressed. Contrary to the case of a triangular lattice with frustrated hopping studied in (18), the atom fraction in each single-particle ground state does not fluctuate, as shown in Fig. 3(c), and we observe an equal population in both states as predicted for interacting (14) or finite size systems (19).

In order to directly reveal time-reversal symmetry breaking of the Hamiltonian, we probed the local structure of the lattice with artificial gauge field at the level of a four-site square plaquette, which allows us to isolate plaquettes with equal sign of the flux. This was achieved by applying superlattice potentials along both the - and -directions, and in order to avoid coupling to axial modes along the potential tubes, an additional lattice along the -direction (nm) was used in these measurements. The four sites of a single plaquette are denoted as [see inset Fig. 4(a)]. The relative phases of the long and short lattices were adjusted so that the plaquettes were symmetric along and tilted along , with an energy offset kHz. We first loaded single atoms in the ground state of the tilted plaquettes , and subsequently switched on the Raman lasers with in order to induce resonant coupling to the and sites. In the limit , the dynamics along would be suppressed and the initial state couples to the state , where the relative phase is induced by the Raman lasers. In our case () the evolution of the imprinted phase factors is more complex. We measured this evolution through the shape of the momentum distribution obtained after time-of-flight (15). This dynamics is a direct consequence of the complex phase factor, revealing the time-reversal symmetry breaking of the Hamiltonian [see Fig. 4(a)]. For the role of the Raman beams is exchanged leading to a sign reversal of the phase evolution.

We also investigated this dynamics in real space in order to exhibit the influence of the gauge field on the particle flow. By generalizing the site-resolved measurement technique performed in (20) for an array of double-wells to plaquettes, we measured the atom population per site () (15), thus obtaining the average atom positions and , with being the total atom number. In the initial state the atoms occupy the left wells and with equal weights. After switching on the Raman lasers we observe a coherent particle flow inside the plaquettes. Besides the particle current towards the right wells and [see Fig. 4(c)], we observe a deviation of the mean atom position along [see Fig. 4(d)]. This behavior is reminiscent of the Lorentz force acting on a charged particle in a magnetic field. As shown in Fig. 4(b), the mean atom position follows an orbit that is a small-scale quantum analog of the classical cyclotron orbits for charged particles. This coherent evolution is damped due to spatial inhomogeneities in the atomic sample. Having independently calibrated the values of and , we fit from the measured atom dynamics the value of the magnetic flux . The difference from the value expected for a homogeneous lattice stems from the smaller distance between lattice sites inside the plaquettes when separated. For the parameters used in Fig. 4(b)-(d) we calculate a distance yielding , which qualitatively explains the measured flux value. Residual deviations might be due to an angle mismatch between the Raman beams and the lattices beams.

Figure 4: Time reversal symmetry breaking and cyclotron orbits. (a) Phase evolution of the double-slit pattern along (integrated along ), as a function of time for (blue) and (gray). The inset shows the Fourier transformation for depicting two frequency components at kHz and kHz, in good agreement with theory (vertical lines). (b) Cyclotron orbits of the average particle position obtained from the mean atom positions (c) and (d) for kHz, kHz and kHz. Each data point is an average over three measurements. The inset in (b) shows the theoretical curve calculated for and a damping time of ms obtained from damped sine fits to and .

In conclusion, we have demonstrated a new type of optical lattice that realizes strong effective magnetic fields and breaks time-reversal symmetry. We have shown that the atomic sample relaxes to the minima of the magnetic bandstructure, realizing an analogue of a frustrated classical spin system. However, the spatial average of the magnetic flux is zero, hence the Bloch band is topologically trivial (21); (17); (10). By using a superlattice potential with more than two non-equivalent sites (10) or a linear tilt potential (9), it is possible to create a lattice with a uniform and non-zero magnetic flux. This system would realize the Harper Hamiltonian (22) and lead to the fractal band structure of the Hofstadter butterfly (23). In particular the lowest band would exhibit a Chern number of one and be analogous to the lowest Landau level (9); (12); (8); (24). Our work constitutes an important step towards the study of quantum Hall effect with ultracold atomic gases and the creation of strongly-interacting liquids such as the Laughlin state (25).

We acknowledge insightful discussions with N. Cooper. This work was supported by the DFG (FOR635, FOR801), the EU (STREP, NAMEQUAM, Marie Curie Fellowship to S.N.), and DARPA (OLE program).


Appendix A.I Bandstructure calculation

We calculated the bandstructure of our lattice according to Ref. (A1). For staggered fluxes a unit cell can be chosen as depicted in Fig. 1(b) of the main text which consists of two plaquettes with opposite sign of the flux. This cell contains two non-equivalent sites, which are denoted as even and odd. For the gauge used in our setup the following ansatz for the wavefunction can be made:


with the quasi-momentum in the ranges and . Substituting this ansatz into the Schrödinger equation using the Hamiltonian given in Eq. (1) of the main text, this reduces to a two-dimensional eigenvalue equation of the following form:


For the case we calculated the dispersion relation [see Fig. 3(a)] showing a single minimum located at . The wavefunction given in Eq. (S1), however, has two components at (even sites) and (odd sites). Therefore, in the time-of-flight measurements, we see two momentum components in the first magnetic Brillouin zone. The eigenvector in this case has equal weight on the two components, thus we see two equally populated peaks. The remaining peaks are shifted by replicas of the reciprocal lattice vectors given by and .

For the ground-state contains two minima as shown in Fig. 3(b). For the same reason described above we observe two additional peaks shifted by .

Appendix A.II Experimental sequence

The experimental sequence started by loading a Bose-Einstein condensate of about atoms with no discernible thermal fraction in 160 ms into a 2D optical lattice of depth and adjusted to a value, which corresponds to the desired bare tunnel coupling , where is the lattice recoil energy. We then ramped up the long-lattice along the -direction to during ms and subsequently decreased the short-lattice depth along to 9  in ms. The lattice potential along is then given by an alternating energy offset  kHz between neighboring sites, which was measured independently. The bare tunnel coupling Hz along was much smaller than , resulting in an inhibition of spontaneous tunneling along , therefore the initial state had no coherence along .

For the case of isolated plaquettes, the experimental sequence started by loading the condensate into a 3D optical lattice of depth , and , where . After that, a filtering sequence was applied to have at most one atom per lattice site (A2). By ramping up the short-lattice along the -direction during ms with the desired energy offset , a double well potential was created with single atoms localized on one side. Finally the short-lattice along the -direction was ramped up during ms to create a plaquette potential where the atom is in the initial state .

Appendix A.III Raman-assisted tunneling

In our experiments tunneling is inhibited along one direction of the lattice due to an energy offset . By using a pair of Raman lasers tunneling can be induced through an additional time-dependent potential , where is controlled by the Raman beam intensities, is the lattice vector, denotes the wave vector and the frequency difference of the Raman lasers. If we consider an atom tunneling from a site of low energy at to a site of high energy at , and assuming , the complex tunneling matrix element is given by the overlap integral:


Here is the Wannier function for and can be chosen real by convention. In the case of Raman-assisted tunneling from a high site to a low site, one obtains and the system can be described by the effective Hamiltonian given in Eq. (2) in the main text.

We also used Raman-assisted tunneling to characterize our system. To measure the value of the induced tunnel coupling , we performed an independent measurement in isolated double wells. The parameters of this system and are chosen to match the ones of the staggered potential. After loading single atoms on the lower site of the double wells we induced tunneling by the Raman beams. The atoms then undergo tunnel oscillations between the two sites of the double well whose frequency is equal to .

The energy offset can be measured starting from the same initial state by scanning the frequency difference of the two Raman beams. When one can observe particle transfer to the second site.

Appendix A.IV Data analysis for Phase measurement

In a system of isolated symmetric double wells, where we denote the two sites as and , the phase of an arbitrary single particle state


can be revealed from its momentum distribution after time-of-flight . The observed density distribution for a sufficiently large is proportional to times an envelope determined by the Wannier function, where , is the distance between the two sites and is the atomic mass.

For the case of isolated plaquettes we first integrate the momentum distribution along and then determine the phase along .

Appendix A.V Site-resolved detection

To detect the occupation numbers on the different sites of the plaquette we apply two mapping sequences along the and direction during which the populations are transferred to different Bloch bands analog to the technique described in Ref. (A2); (A3) for isolated double wells. A subsequent band-mapping technique allows us to determine the population in the Bloch bands by counting the atom numbers in different Brillouin zones. The colors used in Fig. A1(a) and (b) show the connection between the Brillouin zones and the corresponding lattice sites. A typical image obtained after ms of time-of-flight is shown in Fig. A1(c).

Figure A1: (a) Schematic of the four-site plaquette. (b) Brillouin zones of the 2D lattice. (c) Typical momentum distribution obtained after ms of time-of-flight.



  1. D. Tsui, H. Stormer and A. Gossard, Phys. Rev. Lett. 48, 1559 (1982); R. Laughlin, ibid 50, 1395 (1983).
  2. I. Bloch, J. Dalibard and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
  3. A. Fetter, Rev. Mod. Phys. 81, 647 (2009); N. Cooper, Adv. Phys. 57 (2008).
  4. K. Madison, F. Chevy, W. Wohlleben and J. Dalibard, Phys. Rev. Lett. 84, 806 (2000); J. Abo-Shaeer, C. Raman, J. Vogels and W. Ketterle, Science 292, 476 (2001).
  5. V. Schweikhard et al., Phys. Rev. Lett. 92, 40404 (2004); V. Bretin, S. Stock, Y. Seurin and J. Dalibard, ibid 92, 50403 (2004).
  6. N. Cooper, N. Wilkin and J. Gunn, Phys. Rev. Lett. 87, 120405 (2001).
  7. Y. Lin et al., Nature 462, 628 (2009).
  8. J. Dalibard, F. Gerbier, G. Juzeliūnas and P. Öhberg, arXiv:1008.5378 (2010).
  9. D. Jaksch and P. Zoller, New J. Phys. 5, 56 (2003).
  10. F. Gerbier and J. Dalibard, New J. Phys. 12, 033007 (2010).
  11. A. Kolovsky, Europhys. Lett. 93, 20003 (2011).
  12. E. Mueller, Phys. Rev. A 70, 041603 (2004).
  13. L. Lim, C. Smith and A. Hemmerich, Phys. Rev. Lett. 100, 130402 (2008); L. Lim, A. Hemmerich and C. Smith, Phys. Rev. A 81, 023404 (2010).
  14. G. Möller and N. Cooper, Phys. Rev. A 82, 063625 (2010).
  15. See Appendix.
  16. E. Blount, Phys. Rev. 126, 1636 (1962).
  17. Y. Wang and C. Gong, Phys. Rev. B 74, 193301 (2006).
  18. J. Struck et al., Science 333 996 (2011).
  19. In trapped systems the degeneracy is slightly lifted and the two lowest energy states are symmetric and anti-symmetric superpositions of the two states.
  20. S. Fölling et al., Nature 448, 1029 (2007).
  21. F. Haldane, Phys. Rev. Lett. 61, 2015 (1988).
  22. P. Harper, Proc. Phys. Soc. London, Sect. A 68, 874 (1955).
  23. D. Hofstadter, Phys. Rev. B 14, 2239 (1976).
  24. N. Cooper, Phys. Rev. Lett. 106, 175301 (2011); N. Cooper and J. Dalibard, arXiv:1106.0820 (2011).
  25. A. Sørensen, E. Demler and M. Lukin, Phys. Rev. Lett. 94, 86803 (2005); R. Palmer and D. Jaksch, ibid 96, 180407 (2006); M. Hafezi, A. Sørensen, E. Demler and M. Lukin, Phys. Rev. A 76, 023613 (2007).
  26. G. Möller and N. Cooper, Phys. Rev. A 82, 063625 (2010).
  27. S. Fölling et al., Nature 448, 1029 (2007).
  28. J. Sebby-Strabley et al., Phys. Rev. Lett. 98, 200405 (2007).
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.