# Fulde-Ferrell-Larkin-Ovchinnikov superfluidity in one-dimensional optical lattices

## Abstract

Spin-polarized attractive Fermi gases in one-dimensional (1D) optical lattices are expected to be remarkably good candidates for the observation of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase. We model these systems with an attractive Hubbard model with population imbalance. By means of the density-matrix renormalization-group method we compute the pairing correlations as well as the static spin and charge structure factors in the whole range from weak to strong coupling. We demonstrate that pairing correlations exhibit quasi-long range order and oscillations at the wave number expected from FFLO theory. However, we also show by numerically computing the mixed spin-charge static structure factor that charge and spin degrees of freedom appear to be coupled already for small imbalance. We discuss the consequences of this coupling for the observation of the FFLO phase, as well as for the stabilization of the quasi-long range order into long-range order by coupling many identical 1D systems, as in quasi-1D optical lattices.

###### pacs:

03.75.Ss, 03.75.Lm, 71.10.Pm## I Introduction

Multi-component attractive fermionic systems with unequal masses, densities or chemical
potentials have attracted continued interest for many decades in several
fields of physics ranging from high-energy (1); (2) to
condensed matter (2); (3) and, more recently, atomic
physics (4); (5); (6).
The interplay between pairing and density imbalance of the different fermion
species leads to a rich scenario, which includes the possibility of
various exotic superconducting states (7).
In this context, the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase (8) has recently
attracted a great deal of interest from both the experimental and the
theoretical community (1); (2); (3); (4); (5); (6).
In the FFLO phase Cooper pairing occurs between a
fermion with momentum and spin and a fermion with momentum
() and spin .
As a result, the superconducting order parameter becomes spatially dependent.
Originally, the most favorable systems for the observation of the FFLO phase were predicted to be
clean superconducting films in the presence of an in-plane (*i.e.* Zeeman) magnetic field,
above the so-called Clogston-Chandrasekhar limit (9). Nevertheless, despite the fact that
the original prediction dates back to more than thirty years ago,
the FFLO phase has been very elusive to detect.

The experimental realization of interacting trapped Fermi gases with population imbalance (4); (5) has renewed the hope of observing the FFLO, thus stimulating an intense theoretical activity (6); (10). So far, most of the theoretical analysis has focused on 3D cold atomic systems. However, as in the case of solid-state superconductors, the region of phase diagram where the FFLO phase has been found to be stable is quite small (2); (6). On the other hand, quasi-one dimensional or strongly anisotropic systems (such as coupled chains, heavy-fermion, organic, high-, and superconductors) are believed to be good candidates for the realization of the FFLO phase (2); (3); (11); (12). Since the dimensionality of cold atomic systems can be easily tuned, and indeed cold atoms have already been successfully trapped in 1D geometries (13), it seems natural to consider these low dimensional systems as the ideal candidates to observe non-homogeneous superconductivity of the FFLO type.

Many important results are available on the properties of spin-polarized 1D Fermi systems with attractive interactions, which have been obtained by different methods and techniques. These include the Bethe-Ansatz solutions of certain exactly solvable models like the Hubbard or Gaudin-Yang models (17); (16); (15); (14), as well as different types of numerical approaches such as density-matrix renormalization-group (DMRG) (18); (19); (20) and quantum Monte Carlo (QMC) (21), or field theoretical techniques like bosonization (12); (18). Very recently Orso (14) and Hu et al. (15) have studied the phase diagram of harmonically-trapped 1D polarized Fermi gases by combining the exact solution of the Gaudin-Yang model with a local-density approximation. Mean-field theory has also been applied by Liu et al. (22), although it is known that it has a number of limitations (23) in 1D, particularly as far as paring correlations are concerned. DMRG has been employed by Feiguin et al. (19) and Tezuka and Ueda (20), and QMC by Batrouni et al. (21) to investigate the pairing correlations in the spin-polarized ground state of the attractive Hubbard model in the presence of harmonic trapping. Previously, Yang (12) used bosonization to study the pairing correlations and the phase diagram of a single 1D Fermi system as well as an array of weakly-coupled 1D Fermi systems in the presence of a Zeeman field.

Yang’s analysis is valid only close to a continuous magnetic-field-driven transition from a uniform BCS phase and an FFLO phase, which he assumed to belong to the commensurate-incommensurate universality class (24). Another important assumption in Ref. (12) is that charge and spin degrees of freedom are decoupled at low energies for small polarization. However, this scenario does not apply to the Hubbard model away from half-filling (17); (26); (25) nor to the Gaudin-Yang model (16), which are the relevant models for current 1D cold atomic systems. As we show in this work, charge and spin degrees of freedom are indeed coupled already for small polarization, which leads to important differences as compared to the scenario described by Yang. We also numerically demonstrate that the pairing correlation function exhibits prominent oscillations with a wave number equal (up to finite size corrections, see below) to the difference of Fermi wave numbers, , as predicted by FFLO theory in 1D (12) and in agreement with a number of Luttinger-type theorems (27); (28). Thus, the finite-wave-number oscillations in the pairing correlation function can be regarded as due to the excess of unpaired majority-spin fermions (12). This is because in 1D the Fermi wave number and the density are proportional to each other: . On the other hand, this relationship is no longer linear in dimensionality higher than one, and in this case the oscillations in the order parameter are related to the center-of-mass momentum of the Coopers pairs. These observations seem to indicate that in the 1D case it is not entirely clear whether there is a strict close parallelism with higher dimensional FFLO, and in some respects the system can be also understood as a a coupled Bose-Fermi mixture of spin-singlet pairs (the bosons) and unpaired fermions (26); (15).

## Ii The model

We consider a two-component mixture with a total of
fermionic atoms loaded in a 1D optical lattice of
sites (the lattice constant is taken to be unity). The fermions are assumed
to interact via attractive on-site interactions, whose strength can be tuned
*e.g.* by means of a Feshbach resonance. Sufficiently away from resonance(s),
this system is modeled by the attractive Hubbard model:

(1) |

where is the hopping parameter, () is the creation (destruction) fermion operator in the -site (), the pseudospin- index (in experiments this labels the two different atomic hyperfine states of the mixture), is the strength of the on-site Hubbard attraction, , and . In order to simulate the effect of an external trapping potential, open-boundary conditions (OBC) breaking translational symmetry will be used (these are indeed the most suitable conditions for the DMRG treatment of the above model (30)). Our calculations are performed in the canonical ensemble, and the results apply only to lattices away from half-filling, that is, when . In the calculations, the spin polarization was varied by decreasing while keeping constant the number of “background” up-spin atoms , from (the unpolarized case, i.e. ) all the way down to (the fully polarized case, i.e. ).

In the unpolarized case (), all fermions pair into spin singlets due to the attractive on-site interaction. This yields a gap to all spin excitations and therefore spin-spin correlations decay exponentially with distance. Singlet superconducting and charge-density wave correlations exhibit a slower decay (of power-law type in the ground state of a thermodynamically large system), being singlet superconducting correlations the ones that dominate at long distances in systems away from half-filing (29). The aim of this work is to study the nature of superfluidity for as a function of the dimensionless ratio (in the fully polarized case, where , describes a system of noninteracting fermions). The expectation values of all operators below are understood to be taken over the ground state of .

## Iii Numerical results and discussion

Due to the OBC (or, in general, to any external potential that breaks the Bloch translational invariance of the lattice) the spin-resolved site occupation profiles, , exhibit Friedel oscillations. This is illustrated in Fig. 1. In the unpolarized case the Friedel oscillations in are in-phase with those in giving rise to large-amplitude atomic-density waves (31); (32) in the total site occupation . As it is clear from the top panel of Fig. 1, in the general case the total site occupation displays maxima associated with the formation of spin-singlet pairs that are delocalized over the lattice. In the bottom panel of Fig. 1, we show the local spin polarization (which could be measured through phase-sensitive optical imaging (5)). For small (see e.g. the plot for and ) the local spin polarization displays maxima corresponding to the number of fermions that are left unpaired. With increasing though the spatial dependence of the local spin polarization becomes more complicated: the amplitude of the oscillations in the bulk becomes indeed smaller, thus making it hard to clearly identify maxima. These, however, are not distinctive and unambiguous signals of FFLO pairing.

We thus proceed below to present a study of pairing correlations: the model described in Eq. (1), in fact, cannot sustain any true long-range order (29) in 1D, i.e. the ground-state expectation value of the pairing operator is zero.

In the unpolarized case and for an extended system, the correlation function of the pairing operator decays with a power law at large distances, where is an interaction-dependent Luttinger-liquid dimensionless parameter (29). In the top panel of Fig. 2 we illustrate our DMRG results for the spin-polarization dependence of at , which measures real-space superfluid correlations between the site (the center of the trap) and all the other sites. For the power-law decay of the for is clearly visible. For finite , instead, the superfluid correlator is characterized by a distinctive oscillatory character (12) and a very simple nodal structure with exactly zeroes. We have carefully checked that the long-distance decay of is still power-law, signaling quasi-long range superfluid behavior also at finite .

A careful analysis of the oscillatory character of can be done by means of the Fourier transform of the pairing correlator

(2) |

where [with , ] are the eigenstates of the hopping term in Eq. (1). The mode with zero wave number is excluded from the allowed values due to the OBC. The lowest energy mode corresponds to . The diagonal part of the matrix will be simply denoted by . In the bottom panel of Fig. 2 we plot the difference between and its value in the noninteracting gas (i.e. at ), (33).

At possesses a very narrow peak at (see inset in the bottom panel of Fig. 2), signaling quasi-long-range superfluid order of the conventional BCS type. For a finite , instead, has a local minimum at and a single well-defined peak appears at a wave number , where are the spin-resolved Fermi wave numbers. The peak at in the Fourier transform of the pairing correlator, which is a direct consequence of the simple real-space nodal structure illustrated in the top panel of Fig. 2, is a clear-cut signal of FFLO pairing.

The DMRG data shown in Fig. 2 refer only to a single value of . We now turn to illustrate the dependence of the superfluid correlation functions on . In the top panel of Fig. 3 we illustrate the dependence of on the interaction strength for a fixed spin polarization . On decreasing the quasi-long-range FFLO order (i.e. the height of the peak at ), which is emphatically strong for large , survives all the way down to the weak coupling regime. This can be quantified better by analyzing the size of the anomaly at , which is measured by the difference between left and right (discrete) derivatives of evaluated at ,

(3) |

In the bottom panel of Fig. 3 we plot as a function of . In this range decreases in a smooth fashion to its noninteracting value (i.e. ) as is decreased to zero. In other words, for every finite , tends uniformly and smoothly to its noninteracting value as is decreased towards zero. For sufficiently large values of the FFLO phase can also be characterized by the peak visibility defined by

(4) |

This quantity is plotted in an inset to the bottom panel of Fig. 3.

Before concluding, we would like to illustrate the behavior of the density-density, spin-spin, and mixed density-spin static structure factors, , , and . These are defined by the sum over all frequencies of the corresponding dynamic structure factors (34) that can be in principle measured through Bragg spectroscopy or Fourier sampling of time-of-flight images (35). In practice, , , and are calculated from the following equations:

(5) |

where and . In Figs. 4-5 we show the dependence of , , and on for a slightly asymmetric system with and (). We remind the reader that in the unpolarized case has a peak at that signals real-space atomic-density waves (31); (32). In the spin-polarized case, this peak splits into two peaks at and . This is clearly visible in Fig. 4 in the static density-density structure factor (top panel), which presents a double-peak structure slightly below ( and for the system parameters in this figure). This double-peak structure is not so visible in the magnetic structure factor , most likely because magnetic correlations near and are still quite suppressed by the superfluid correlations, at least in the weakly polarized case (in the unpolarized case they are completely suppressed by the pairing gap). From Fig. 5 we note that is non-zero even at small , thus indicating that spin and charge degrees of freedom are coupled at long wavelength even for a small imbalance.

### Experimental signatures of the FFLO phase

The most direct way to detect FFLO pairing would be to measure the pairing correlation function . We would like to remark here that this correlation function is, at least in principle, measurable via interferometric schemes (36) in which two atomic wave packets are coherently extracted from the gas at different positions and then are mixed by a matter-wave beam splitter. The atom counting statistics in the beam splitter output channels has been shown (36) to reflect the spatial dependence of .

The oscillations of the pairing correlations will also leave a detectable signature in the noise correlations (37); (38), , where measures the number of fermions with momentum and spin in a time-of-flight experiment. With increasing spin-polarization, in fact, the peak at (39); (40) [here is the direction along the axis of the 1D system] will shift to a finite relative momentum (see e.g. the work by Yang in Ref. (10) and the very recent DMRG calculation by Lüscher et al. (41)).

However, it is worth pointing out that the strength of noise signal in a strictly-1D system will be strongly affected by finite-size and temperature effects. This is because in 1D the order is not long-range but quasi-long range, and therefore the slowest decay exhibited by correlations (like the pairing correlations) is a power law. Thus, in order to enhance the strength of the experimental signal for FFLO, it would be desirable to couple many 1D systems, as in a tight 2D optical lattice (arrays of “atomic quantum wires”) (42), so that the quasi-long range FFLO order can become true long-range order. The phase diagram of many coupled 1D systems has been worked out in Ref. (12), where the author found that, at small polarization , true long-range 3D FFLO order will occur when the Luttinger liquid parameter for the charge excitations, , is larger than . In such a case the low temperature properties of the system are dominated by hopping of pairs (i.e. Josephson coupling) rather than by single-particle hopping (the latter would turn the system into an anisotropic Fermi liquid, which could in turn become unstable to the FFLO state under appropriate conditions (43)).

However, the analysis of Ref. (12) assumed the transition from the unpolarized to the polarized case to belong the the commensurate-incommensurate universality class. This assumption implicitly neglects the coupling between charge and spin degrees of freedom at low energies, which is known to modify the behavior of physical observables at the transition (26); (25). In this work this coupling has been demonstrated to exist also at long wavelengths by an explicit numerical evaluation of the mixed static structure factor in a weakly polarized system (see Fig. 5). Thus, the phase diagram depicted in Fig. 1 of Ref. (12) seems not appropriate for coupled 1D Hubbard (or Gaudin-Yang (16)) models, and the transition to long-range order will not take place in general for and may in general depend on the system parameters (i.e. and the lattice filling for the Hubbard model).

To the best of our knowledge, a quantitative phase diagram of coupled 1D systems lacking spin-charge separation has not yet been calculated. Furthermore, it is worth noticing that in cold atomic systems with short range interactions, another important factor must be taken into account, namely the relative strength of the pair hopping when compared to the single-particle hopping. The strength of the latter is given by (, where is the Fermi energy, for the analysis based on coupled Tomonaga-Luttinger liquids to hold), where is the hopping amplitude between two neighboring 1D systems. However, in absence of long-range interactions pair hopping can be only generated by (virtual) single-particle tunneling events, which at lowest order, yield a (Josephson) coupling strength of order , where is the spin gap. The phase diagram predicted in Ref. (12) is the result of a calculation which only compares the scaling dimensions of the pair hopping and single-particle hopping operators, and thus does not take into account the microscopic details of the coupling between 1D systems. Thus, the stabilization of FFLO long-range order by a weak coupling between 1D systems in the FFLO phase (the extreme anisotropic limit that could not be accessed by the authors of Ref. (43)) does not seem easily achievable. In turn, the most likely scenario for arbitrary polarization is that the single-particle hopping will control the physics at low temperatures, and the system will behave as a spin-polarized normal Fermi liquid, which, in turn, could become unstable towards 3D FFLO ordering (43) under appropriate conditions.

## Iv Conclusions

In summary, we have shown how ultracold spin-polarized two-component Fermi gases confined in 1D optical lattices are FFLO superfluids whose pairing correlation functions are characterized by a power-law decay and a simple nodal structure. However, we have also shown that charge and spin degrees of freedom appear to be coupled already for a small value of the spin polarization. Finally, we have commented on the impact of this coupling on the detectability of true long-range FFLO order arising from Josephson coupling between 1D systems.

###### Acknowledgements.

We acknowledge useful discussions with P. Calabrese, I. Carusotto, M. Köhl, and A.H. MacDonald. M.P. gratefully acknowledges the hospitality of the Donostia International Physics Center and of the Department of Physics of the Zhejiang Normal University during the final stages of this work. M.A.C. thanks A. Nersesyan for drawing attention to Refs. (26). This work was partly supported by PRIN-MIUR and by the Academy of Finland (project number 115020). The DMRG calculations have been made using the DMRG code released within the “Powder with Power” Project (www.qti.sns.it).### References

- K. Rajagopal and F. Wilczek, in Frontier of Particle Physics, edited by M. Shifman (World Scientific, Singapore, 2002); also available as arXiv:hep-ph/0011333v2.
- R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 76, 263 (2004).
- K. Yang, in Pairing in Fermionic Systems: Basic Concepts and Modern Applications, edited by M. Alford, J. Clark, and A. Sedrakian (World Scientific, Singapore, 2006); also available as cond-mat/0603190.
- M.W. Zwierlein, A. Schirotzek, C.H. Schunck and W. Ketterle, Science 311, 492 (2006). G.B. Partridge, W. Li, R.I. Kamar, Y. Liao and R.G. Hulet, ibid. 311, 503 (2006); M.W. Zwierlein, C.H. Schunck, A. Schirotzek and W. Ketterle, Nature 442, 54 (2006).
- Y. Shin, M.W. Zwierlein, C.H. Schunck, A. Schirotzek and W. Ketterle, Phys. Rev. Lett. 97, 030401 (2006); G.B. Partridge, W. Li, Y.A. Liao, R.G. Hulet, M. Haque and H.T.C. Stoof, ibid. 97, 190407 (2006).
- M.M. Parish, F.M. Marchetti, A. Lamacraft and B.D. Simons, Nature Physics 3, 124 (2007).
- G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963); H. Müther and A. Sedrakian, Phys. Rev. Lett. 88, 252503 (2002); W.V. Liu and F. Wilczek, ibid. 90, 047002 (2003); P. F. Bedaque, H. Caldas, and G. Rupak, ibid. 91, 247002 (2003).
- P. Fulde and R.A. Ferrell, Phys. Rev. 135, A550 (1964); A.J. Larkin and Y.N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 47, 1136 (1964) [Sov. Phys. JEPT 20, 762 (1965)].
- In a superconductor in the presence of Zeeman splitting, the Clogston-Chandrasekhar limit determines the critical field below which a uniform BCS order parameter can exist. Above this field the electron system is fully polarized [B.S. Chandrasekhar, Appl. Phys. Lett. 1, 7 (1962); A.M. Clogston, Phys. Rev. Lett. 9, 266 (1962)].
- See for example R. Combescot, Europhys. Lett. 55, 150 (2001); T. Mizushima, K. Machida, and M. Ichioka, Phys. Rev. Lett. 94, 060404 (2005); K. Yang, ibid. 95, 218903 (2005); J. Carlson and S. Reddy, ibid. 95, 060401 (2005); D.E. Sheehy and L. Radzihovsky, ibid. 96, 060401 (2006); P. Pieri and G.C. Strinati, ibid. 96, 150404 (2006); H. Hu and X.-J. Liu, Phys. Rev. A73, 051603 (2006); M. Haque and H.T.C. Stoof, ibid. 74, 011602 (2006); M. Mannarelli, G. Nardulli, and M. Ruggieri, ibid. 74, 033606 (2006); Y.-P. Shim, R.A. Duine, and A.H. MacDonald, ibid. 74, 053602 (2006); J. Kinnunen, L. M. Jensen, and P. Törmä, Phys. Rev. Lett. 96, 110403 (2006); T.N. De Silva and E.J. Mueller, ibid. 97, 070402 (2006); C.-C. Chien, Q. Chen, Y. He, and K. Levin, ibid. 97, 090402 (2006); K. Machida, T. Mizushima, and M. Ichioka, ibid. 97, 120407 (2006); C.-C. Chien, Q. Chen, Y. He, and K. Levin, ibid. 98, 110404 (2007) and other related works cited in these papers.
- T. Pruschke and H. Shiba, Phys. Rev. B46, 356 (1992).
- K. Yang, Phys. Rev. B63, 140511(R) (2001).
- See e.g. B. Paredes, A. Widera, V. Murg, O. Mandel, S. Fölling, I. Cirac, G.V. Shlyapnikov, T.W. Hänsch and I. Bloch, Nature 429, 277 (2004); T. Kinoshita, T. Wenger, and D.S. Weiss, Science 305, 1125 (2004).
- G. Orso, Phys. Rev. Lett. 98, 070402 (2007).
- H. Hu, X.-J. Liu, and P.D. Drummond, Phys. Rev. Lett. 98, 070403 (2007).
- X.W. Guan, M.T. Batchelor, C. Lee, and M. Bortz, Phys. Rev. A76, 043617 (2007).
- T.B. Bahder and F. Woynarovich, Phys. Rev. B33, 2114 (1986).
- G. Roux, E. Orignac, P. Pujol, and D. Poilblanc, Phys. Rev. B75, 245119 (2007).
- A.E. Feiguin and F. Heidrich-Meisner, arXiv:0707.4172v3.
- M. Tezuka and M. Ueda, arXiv:0708.0894v1.
- G.G. Batrouni, M.H. Huntley, V.G. Rousseau, and R.T. Scalettar, arXiv:0710.1353v1.
- X.-J. Liu, H. Hu, and P.D. Drummond, Phys. Rev. A76, 043605 (2007).
- F. Marsiglio, Phys. Rev. B55, 575 (1997).
- G.I. Japaridze and A.A. Nersesyan, JETP Lett. 27, 334 (1978); G.I. Japaridze and A.A. Nersesyan, Phys. Lett. 85 A, 23 (1981); G.I. Japaridze and A.A. Nersesyan, J. Low Temp. Phys. 47, 91 (1983); V.L. Pokrovsky and A. L. Talapov, Phys. Rev. Lett. 42, 65 (1979).
- H. Frahm and T. Vekua, arXiv:0710.44981v1.
- F. Woynarovich, Phys. Rev. B43, 11448 (1991); K. Penc and F. Woynarovich, Z. Phys. 85, 269 (1991).
- F.D.M. Haldane, Proceedings of the International School of Physics “Enrico Fermi”, Course CXXI “Perspectives in Many-Particle Physics”, edited by R.A. Broglia and J.R. Schrieffer (North-Holland, Amsterdam, 1994), pag. 5; also available as arXiv:cond-mat/0505529v1.
- M. Yamanaka, M. Oshikawa, and I. Affleck, Phys. Rev. Lett. 79, 1110 (1997); M. Oshikawa, M. Yamanaka, and I. Affleck, ibid. 78, 1984 (1997); P. Gagliardini, S. Hass, and T.M. Rice, Phys. Rev. B58, 9603 (1998).
- T. Giamarchi, Quantum Physics in One Dimension (Clarendon Press, Oxford, 2004).
- U. Schollwöck, Rev. Mod. Phys. 77, 259 (2005).
- Gao Xianlong, M. Rizzi, M. Polini, R. Fazio, M.P. Tosi, V.L. Campo, Jr., and K. Capelle, Phys. Rev. Lett. 98, 030404 (2007).
- F. Karim Pour, M. Rigol, S. Wessel, and A. Muramatsu, Phys. Rev. B75, 161104(R) (2007).
- is calculated from the noninteracting real-space correlator , with .
- G.F. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid (Cambridge University Press, Cambridge, 2005).
- L.-M. Duan, Phys. Rev. Lett. 96, 103201 (2006).
- I. Carusotto and Y. Castin, Phys. Rev. Lett. 94, 223202 (2005).
- E. Altman, E. Demler, and M.D. Lukin, Phys. Rev. A70, 013603 (2004).
- M. Greiner, C.A. Regal, J.T. Stewart, and D.S. Jin, Phys. Rev. Lett. 94, 110401 (2005).
- L. Mathey, E. Altman, and A. Vishwanath, arXiv:cond-mat/0507108v3.
- A. Lüscher, A.M. Läuchli, and R.M. Noack, Phys. Rev. A76, 043614 (2007).
- A. Lüscher, R.M. Noack, and A.M. Läuchli, arXiv:0712.1808v2.
- F. Schreck, L. Khaykovich, K.L. Corwin, G. Ferrari, T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. Lett. 87, 080403 (2001); A. Görlitz, J.M. Vogels, A.E. Leanhardt, C. Raman, T.L. Gustavson, J.R. Abo-Shaeer, A.P. Chikkatur, S. Gupta, S. Inouye, T. Rosenband, and W. Ketterle, ibid. 87, 130402 (2001); M. Greiner, I. Bloch, O. Mandel, T.W. Hänsch, and T. Esslinger, ibid. 87, 160405 (2001); H. Moritz, T. Stöferle, M. Köhl, and T. Esslinger, ibid. 91, 250402 (2003); T. Stöferle, H. Moritz, C. Schori, M. Köhl, and T. Esslinger, ibid. 92, 130403 (2004); H. Moritz,T. Stöferle, K. Günter, M. Köhl, and T. Esslinger, ibid. 94, 210401 (2005).
- M.M. Parish, S.K. Baur, E.J. Mueller, and D.A. Huse, Phys. Rev. Lett. 99, 250403 (2007).