One-dimensional sawtooth and zigzag lattices for ultracold atoms.

We describe tunable optical sawtooth and zigzag lattices for ultracold atoms. Making use of the superlattice generated by commensurate wavelengths of light beams, tunable geometries including zigzag and sawtooth configurations can be realised. We provide an experimentally feasible method to fully control inter- (t) and intra- (t') unit-cell tunnelling in zigzag and sawtooth lattices. We analyse the conversion of the lattice geometry from zigzag to sawtooth, and show that a nearly flat band is attainable in the sawtooth configuration by means of tuning the lattice parameters. The bandwidth of the first excited band can be reduced up to 2% of the ground bandwidth for a wide range of lattice setting. A nearly flat band available in a tunable sawtooth lattice would offer a versatile platform for the study of interaction-driven quantum many-body states with ultracold atoms.

Recent experimental and theoretical advancements in optical lattices hold out a promise of engineering the band structures, the amplitude and the phase of inter-site tunnelling energies, and on-site interactions in the Hubbard regime12. Furthermore, the highly tunable optical lattices345 would realise unconventional band structures that are not easy to study in solid-state materials, including flat bands with completely quenched kinetic energies or topological bands. Flat band occurs in a variety of condensed-matter systems, ranging from the Landau levels of an electron gas6, edge states of graphene7 to unconventional superconductors8. Such flattened band structures significantly enhance interaction effects and show a wide range of interaction-driven many-body phenomena from supersolidity9 to flat-band ferromagnetism10. Intriguingly, the divergence in the density of states at the energy of the flat band prevents natural ordering in the system1112, and the question still remains whether or not Bose-Einstein condensation is stable in a flat band1314. In addition, the interaction-induced localised state reveals crystalline phases in fractional-filled bosons and fermions in a flat band9 and more interestingly topological flat bands with nonzero Chern number open a new way for realising a fractional quantum hall state without Landau levels121516.

In the ultracold atomic system, a bulk flat band appears though the whole system when the hopping matrix elements are destructively interfered such as in kagome5, Lieb17, sawtooth918 lattices. However, it has not been experimentally explored in great depth so far how the flat dispersion affects the ground state of bosons and fermions except the very recent work in which the decay of the bosons from the flat excited band has been studied19. This is because the flat band is usually not a ground band in a typical optical lattice with negative tunnelling energy20. Indeed, the surface flat band has been also theoretically investigated in optical checkerboard lattices when the system is topologically nontrivial21. Intriguingly, a nearly flat band can be also obtained by including the spin degrees of freedom in the presence of the spin-orbit coupling222324 or the dipolar interactions25. We note similar flattened quartic dispersion is also available in the shaken lattice26.

In this work, we describe a one-dimensional tunable optical zigzag lattice2728 and then analyse the geometric conversion from a zigzag into a sawtooth optical lattice that provides a nearly flat band918. Especially, a sawtooth lattice has attracted significant interest due to the flat band structure but the creation of tunable sawtooth geometry has not been proposed to our knowledge. Our approach relies on a tunable superlattice generated by commensurate wavelengths that are red-detuned to the principle optical transition of atoms. We first demonstrate a method to fully control tunnelling matrix elements in the zigzag configuration. The full control of zigzag lattice is the critical step to access a nearly flat band in a sawtooth lattice which is converted from the original zigzag configuration. Only for certain lattice configuration satisfying the relation between inter-(t) and intra-(t′) unit-cell tunnelling energies, a flat band shows up in the sawtooth geometry. We analyse the band structure of the sawtooth and zigzag lattices for various lattice settings and provide a recipe for obtaining a nearly flat band in the sawtooth configuration. Besides the sawtooth lattice containing a flat band, a one-dimensional zigzag chain2728 also reveals particularly rich physics due to the interplay between frustration imposed by lattice geometry and two (three)-body interactions. Indeed, a chiral bosonic superfluid is expected in a zigzag lattice, which may become a Mott insulator even at small interactions2728.

Results

Creation of optical zigzag and sawtooth lattices

An array of one-dimensional zigzag lattices can be created in two-dimensional superlattice geometry. Here, the superlattice is generated by the incoherent superposition of bichromatic lattice beams. We assume that both short-wavelength (SW) and long-wavelength (LW) lattices are red-detuned from the principle optical transition of the atom and commensurate with each other. For example, the wavelength of 532 nm (SW) and 1064 nm (LW) lights are employed for ytterbium atoms with the principle optical transition of 399 nm2930. The SW lattice consists of the three intersecting beams of the same optical frequency along x- and y- directions (see Fig. 1). All lights are vertically polarised perpendicular to the lattice plane, which results in the SW lattice potential (say V). Additionally, the LW lattice is added to decouple a y-directional zigzag chain from the adjacent one. We note other geometries such as a triangular or a square lattice can be also generated as described in Fig.1(b), but we focus on the zigzag and sawtooth lattices for the rest of the discussion. Finally, additional square lattice generated by the 1064 nm light is incoherently added with the maximum trap depth of V in order to convert the zigzag to the sawtooth lattice. Here, the frequency difference (e.g. at least 10 MHz or larger) between and V laser fields is asuumed to minimise the interference effect.

Figure 1

Schematics of the optical superlattice generating zigzag and sawtooth lattices.

(a) Three bichromatic light beams with potential depth of V, V and form a tunable zigzag optical lattice with the two-dimensional potential shown in (b). Adding a two-dimensional square lattice with the depth of V, the lattice is converted to the sawtooth geometry. (b) Two-dimensional triangular and square lattices are also attainable without LW lattice when V = V and 16V = V respectively. A typical zigzag or sawtooth lattice is shown in (b) with the setting of (V, V, , V) = (1, 1, 1, 0) or (1, 1, 1, 1) × E. In this case, the time phase is set to ϕ = 0, ϕ = 0, and ϕ = −π.

With all beams in hand, the total optical potential is given as:

where time-phase of the short-wave lattice is , ϕ the phase of the light at the position of atoms, ϕ the accumulated optical phase for the round trip from the atom where i ∈ x, y. Here, the lattice depth V where is the depth of the single-beam normalised by the recoil energy for the short-wavelength where the wavenumber for λ = 532 nm. The time-phase of the lattice beam can be actively stabilised by means of controlling optical paths or the laser frequency3519. For the rest of the discussion, we shall assume all optical phases are fixed as ϕ = 0,  = π/2, ϕ = 0, and ϕ = −π, but only each lattice depths (, V) are varied.

Zigzag and sawtooth lattices

To achieve a sawtooth lattice with a nearly flat band, we begin by considering an alternative lattice structure which is easier to obtain, a zigzag lattice. In our optical potential (Eq. (1)), if we set V = 0, then a zigzag lattice is at hand. The zigzag lattice has the two lowest energy s-orbital bands corresponding to two basis lattice sites, labeled by A and B, within the unit cell (see Fig. 2). The inter-unit-cell tunnelling between neighbouring A (or B) wells is defined as t (or t″), and the intra-unit-cell tunnelling between A and B is t′. We employ a large potential depth of the LW lattice to isolate a one-dimensional array in vertical direction, and therefore the coupling between adjacent one-dimensional zigzag lattices is negligible. With symmetric zigzag lattice configuration, we have t″ = t. In this case, the band structure is determined by inter- (t, t″) and intra- (t′) unit-cell tunnelling energies as shown in Fig. 2(a) for variable tunnelling ratio t′/t.

Figure 2

Zigzag and sawtooth lattices.

(a) A one-dimensional zigzag lattice consists of two sub-lattice sites A and B in a unit cell (gray in (a)) with tunnelling eneergies (t, t′ and t″). The dispersion of zigzag lattice in the case of ε = ε and t″ = t is shown for various cases: , , and respectively (from left to right). (b) The ratio between inter-(t) and intra-(t′) unit-cell tunnelling is calculated as a function of Vx and Vy in the zigzag lattice. (c) The zigzag lattice is converted to a sawtooth lattice when the tunnelling t″is suppressed. The square lattice with the depth of V increases the barrier between the potential wells B but decreases the one between the potential wells A. The lattice potential of (V + V), described by the solid line, exhibits sawtooth configuration.

To search for a flat band in the sawtooth lattice, let us first consider the particular ratio of t′/t required to support a flat band in a tight-binding model. A sawtooth lattice with a flat band will be obtained if we make t′′ = 0, and the relation is satisfied9. These two conditions can be achieved by varying the parameters of the optical potential such as . To achieve the first condition of t″ = 0, we introduce the square lattice of the potential depth V as described below, and then we vary the ratio between V and V to satisfy the second requirement of in the sawtooth lattice.

Here we first demonstrate the way to fully control the tunnelling ratio t′/t in the zigzag case before we discuss the transform of the lattice into the sawtooth configuration. As it is difficult to obtain an obvious relation connecting the potential parameters V and V with the tunnelling energy t and t′, we first calculate the Bloch state of the two-dimensional superlattice using the plane-wave expansion method with variable lattice settings3132. Then, we perform least-square fitting of the band structure along the y-direction by the one-dimensional tight binding model of the zigzag lattice. The Hamiltonian for the zigzag lattice in the second quantised form is given:

where is the annihilation operator at the site s ∈ {A, B} in the i-th unit cell. Then in the momentum-space basis , the Hamiltonian is diagonalised giving the tight-binding matrix H of the zigzag lattice.

where ε denotes the lowest eigen-energy of the quantum well A, B and a the lattice constant.

Since ε = ε and t″ = t are naturally satisfied in the zigzag lattice we consider here, the tight binding dispersion for a zigzag lattice is then simplified to:

This dispersion is quite sensitive to the ratio t′/t as shown in Fig. 2(a). Therefore fitting the band structure obtained by the plane-wave expansion method with the Eq.(4), we obtain the ratio of t′/t from the potential parameters V and V as shown in Fig. 2(b). To fulfil the condition , the ratio of V/V should be extremely high around the order for about ~100. In this case, the hopping integrals are quite sensitive to the overlap between Wannier functions. Indeed, they are exponentially dependent on space coordinates, and therefore the hopping elements are very sensitive to the separation between the associated potential wells. For current form of the potential well, the separation between two neighbouring A wells that is about times larger than the separation between neighbouring A and B wells, and therefore it will generate an extremely large t′/t ratio. To reduce this ratio down to t′/t ~ O(1), we have to modify the shape of the potential wells reducing t′ and enhancing t. So it’s necessary to increase V/V ratio, making the potential wells longer and narrower. This method reduces the overlap between A and B wells, but enhances the overlap between neighbouring A wells explaining the reason for such a large value of V/V.

A nearly flat band in a sawtooth lattice

After we tune the zigzag lattice with various values of t′/t, the next step is to transform a zigzag lattice into a sawtooth lattice. To suppress the tunnelling energy t″ appropriately, we now introduce a two-dimensional square lattice potential with the depth of V. The potential maximum (minimum) of the V lattice is aligned to the middle point of neighbouring B (A) wells as described in Fig. 2(c). Then, the square lattice with the depth V will generate an additional barrier between neighbouring B wells reducing t′′ significantly. Because the barrier between neighbouring A wells are reduced by V, and therefore the Vx/Vy should get modified to take this effect into account. From our analysis, making t′/t ~ 10 (in the zigzag lattice) is necessary to balance the effect of V.

Indeed we investigate the dependance of the flatness of the first excited band on the lattice setting (V, V, V). To measure the flatness of the first excited band, we define the ratio between the bandwidths of the ground and the first excited band as R, and numerically calculate R as a function of Vx, Vy, and V as shown in Fig. 3. Since the ratio V/V must be large to obtain t′/t ~ O(1) in the sawtooth lattice, we scan the potential depth V around V~E setting the potential depth V from 28E to 46E. As V decreases, the ratio t′/t increases (see Fig. 2(b)). Therefore, the relatively deep square lattice potential of V is required to keep the small value of t′/t for achieving a nearly flat band. From our numerical analysis, we find the flatness of R > 50 can be easily achieved for a wide range of the lattice setting. In Fig. 3, we provide a set of parameters (V, V, V) where the flatness R is at least over 50 from V = 28E and V = 46E.

Figure 3

Where a flattened band emerges in a sawtooth lattice.

The ratio R between the bandwidths for the ground and the first excited band is calculated for variable lattice depths V, V and V. In the diagram, the solid curve represents the lattice setting (V, V, V) in which the flatness ratio R is larger than 50. The maximum flatness R decrease as a function of the lattice depth V as indicated by the colour scale. For each value of V, the bar around the date point indicates the characteristic range of (V, V) in which the flatness R is larger than 50. The solid and dashed lines are guides to the eye.

For a real flat band in a sawtooth lattice, however, it’s necessary to set the tunnelling t′′ to be exactly zero. This requires the barrier between neighbouring B sites to be as large as possible, and therefore we should increase V. However the increase of V also has consequence in reducing the barrier between neighbouring A sites. So with the presence of the square lattice, A and B potential wells are no longer identical to each other. This will give rise to slight difference in eigenvalues of the two potential wells ε and ε causing that the first excited band becomes dispersive (see Eq.(3)). Therefore the presence of the square lattice with V suppresses the tunnelling t″ but introduces the finite dispersion of the first excited band at the same time. Those limitations are prominent at the small value of V → 28E in Fig. 3, and thus the bandwidth ratio between the ground and the first excited bands has a maximum value of R < 400. For a larger value of V → 46E, the first excited band become flattened at the small value of V. In this case, the tunnelling t″ cannot be completely suppressed causing a small dispersion (or R˜50) of the band. This limitation can be readily solved by employing the lattice shaking technique which allows us to control the tunnelling t and t′ independently3334.

Finally we demonstrate how the dispersive band in the zigzag lattice is adiabatically changed into a nearly flat band in the sawtooth lattice. As an example in Fig. 4(a), the two lowest energy bands are shown as V increases with the lattice setting of . The first excited band becomes flatter as the lattice geometry is changed from the zigzag to the sawtooth limit. Indeed, the bandwidth of the first excited band becomes bosonic atoms is shown for zigzag and sawtooth lattices, which is obtained from the calculation of the single-particle ground state.

Figure 4

Adiabatic transformation from the zigzag to the sawtooth lattice as V is increased.

At the lattice setting of , the energies of the lowest two bands are calculated for V = {0, 0.25, 0.5, 0.85, 1.75} × E (left to right). (b) In this setting, the first excited band becomes flattened with bandwidths (triangle) from h × 50 Hz to h × 0.2 Hz where the maximum flatness R (circle) is ~300. Here the bandwidth is calculated for ytterbium atoms as an example. (c) The real space optical potential and its corresponding momentum distribution are shown. The momentum distribution is calculated from the Fourier transform of the single-particle ground state in the real space. In the momentum distribution, the area of the black dot reflects the fractional population at the reciprocal site and the G and a denote the reciprocal and the lattice vectors of the one-dimensional zigzag lattice respectively.

In optical lattices, a flat band often emerges in either the second or third lowest energy band for sawtooth (or Lieb) and kagome lattices respectively. In general, one may shake the optical lattice to flip the sign of the tunnelling energy by which atoms can be directly loaded in the flat band3334. For bosons, atoms can be coherently transferred at the particular quasi-momentum (say k = 0 point) by manipulating the energy level of the sub-lattice site19. The bosonic atoms in the flat band may decay into the lowest band as observed in the Lieb lattice19, but the lifetime of atoms in the flat band can be lengthened by increasing the band gap19. The effect of the flat band for fermions, however, should be observable with minimal decay into the lowest band when the Fermi energy is controlled near the energy of the flat band.

Conclusion

In this work, we provide a recipe to realise a one-dimensional sawtooth optical lattice with ultracold atoms. The highly tunable superlattice proposed here would offer a flexible experimental platform for the study of a flat band which is attainable in the sawtooth geometry. By simply adjusting lattice parameters and therefore controlling intra- and inter-unit-cell tunnelling, the first excited band becomes flat up to 2% of the ground bandwidth for a wide range of the lattice setting. Our superlattice yields the zigzag and sawtooth lattices for alkali-earth-like bosons and fermions (e.g. 174,173Yb and 84,87Sr isotopes), in which various interaction-driven phenomena can be addressed in a flat band.

Additional Information

How to cite this article: Zhang, T. and Jo, G.-B. One-dimensional sawtooth and zigzag lattices for ultracold atoms. Sci. Rep. 5, 16044; doi: 10.1038/srep16044 (2015).