Arxiv:1209.2671V2 [Cond-Mat.Quant-Gas] 13 Sep 2012

Home , Spin density wave

arXiv:1209.2671v2 [cond-mat.quant-gas] 13 Sep 2012 ordinate er,i eae to related is metry, e a ecs opeeyaaoost htof that to analogous or- completely antiferromagnetic cast such be can of 1(a). der the theory Fig. with the in SDW, Interestingly, depicted exhibits pattern coupling, checkerboard weak well-known in- at half-filling, even has at example which, model quintessential Fermi-Hubbard that square-lattice A the phenomenon is decades. a for physics scientists matter magnetism, trigued physics. condensed quantum of many-body subjects is study central the to of avenues One new opened has etm rdcigteeitneo hl ra of array whole a and of SDW robust states, existence quite SDW mo- the angular is predicting higher analogy with mentum, This pairing particle-hole in to [3]. pairs extends theory Cooper of BCS particle-hole condensation the to of corresponding condensation pairs with – superconductivity eandeuie hs xtcsae ipa modula- a vector display spin states of tion exotic These elusive. remained o pnsnltcag est ae (CDW waves density charge spin-singlet for is ()ad1b.Seicly h re aaee of parameter order the Specifically, SDW 1(b). and 1(a) Figs. nFg.1c n () sopsdt ouaino on- of modulation to opposed SDW conventional as in 1(d), variables depicted site and arrows red 1(c) of Figs. pattern in checkerboard the as such pt o h elzto fsi-rpe SDW spin-triplet of realization the now to up cates hresco,asmlraaoypeit h xsec of existence the predicts analogy generalized similar a sector, charge oantcodro i.1a osiue SDW constitutes 1(a) Fig. of order romagnetic 4 p SDW on unu nttt,Ntoa nttt fSadrsa Standards of Institute National Institute, Quantum Joint wv uecnutr n superfluid and superconductors -wave 2 h deto lr-odaoi n oeua gases molecular and atomic ultra-cold of advent The hl eea addt ytm aebe discussed been have systems candidate several While etu ¨rOtsh unetcnlge n ntttfü Institut and Quantentechnologien für Optische Zentrum p p ℓ etrn atcehl ar with pairs particle-hole featuring , stepril-oeaaou ftespin-triplet the of analogue particle-hole the is 0 = 1 r colo hsc,AtooyadCmuainlSine,G Sciences, Computational and Astronomy Physics, of School i , ν − 1 ftedplroinain n ics o hs Dsarise numbers: PACS SDWs these how phases. wave discuss density and dipola charge o of orientation, and defined diagram dipolar is phase The the moreover, quantum rich and, of lattice. the space, optical determine spin on We in 6] sites. quantity [5, vector molecules a or is [4] atoms of gases eeso htsi-rpe atcehl odnae wit condensates functi particle-hole unbiased spin-triplet an that Using show here [3]. momenta systems angular electron non-zero dimensional with waves density counterparts, ayucnetoa uecnutr ihCoe ar ffi of pairs Cooper with superconductors unconventional many xml 2,cnb ecie sacnest fparticle-ho of condensate a as described ℓ be can [2], example wv hredniywv (CDW wave density charge -wave , ,aaoost odnaeo atcepril ar in pairs particle-particle of condensate a to analogous 0, = r 2 h ovninlsi est ae(D)pae[] sfound as [1], phase (SDW) wave density spin conventional The .G Bhongale G. S. j .. , 3 ˆ = novninlSi est ae nDplrFriGases Fermi Dipolar in Waves Density Spin Unconventional eateto hsc n srnm,Uiest fCalifor of University Astronomy, and Physics of Department epciey[] h aiirantifer- familiar The [3]. respectively ν S S S hr h label the where , y S endo h od ftelattice, the of bonds the on defined η aldsacsaei nt fthe of units in are distances (all = h a ˆ i,α † 1 σ ˆ , 4 αβ η .Mathey L. , s a ˆ n CDW and j,β ν i = 3 p o eaieco- relative for e[] nthe In [7]. He y ,p ,.. d, p, s, obtlsym- -orbital ν s Dtd aur 7 2018) 17, January (Dated: s 6= 2 ν whereas , hnWnTsai Shan-Wen , s hw in shown ). ν 3 8], [3, ) 6= s -wave s indi- has dTcnlg nvriyo ayad atesug D20 MD Gaithersburg, Maryland, of University & Technology nd aepyi,Uiesta abr,271Hmug Germa Hamburg, 22761 Universität Hamburg, Laserphysik, r 2] n ododrsld[8]. supersolid solid bond-order [18–20], and crystalline [23], liquid [13–17], superfluid futacl emoi ehahmlclsof molecules Feshbach formation fermionic the Recently, ultra-cold lattices. of optical into loaded and unu eeeay[] h emoi oa molecule polar fermionic The magne- Bohr 40 [4]. below well 10 degeneracy cooled of quantum and trapped moment successfully dysprosium magnetic been have of large ton, atoms a Fermionic quan- with the 161, in limit. accessible degenerate become tum have fermions dipolar of long- between anisotropic competition interaction. the the dipolar the and ranged in of inter-atomic regime result short-ranged BCS a the the as and diagram, CDW phase the between emerge hssaepeitd hredniywv [9–12], wave density charge quantum predicted: gases Numerous body are Fermi groups. many phases many dipolar side, by explored theory (spinless) been have the single-species On of [6]. physics achieved been has ih nrdet ostabilize the have to ingredients interaction right dipole-dipole dominant with lattice ie r umdoe) eeˆ Here in- over). repeated summed and are paper, dices this throughout constant lattice cie yo-ieodrparameters order on-site by scribed h al arcs n ˆ and matrices, Pauli the h in(rain prtr o pseudo-spin for operators (creation) tion n odaoal poiests e.g., sites, opposite diagonally ( to with similarly ing defined is parameter der tended site ncnrs,tecnetoa SDW conventional the contrast, In n ˆ K i eew osdratocmoet(suosi 1 (pseudo-spin two-component a consider we Here nanwgnrto feprmns lr-odgases ultra-cold experiments, of generation new a In h e nih fti ae sta emosi 2D a in fermions that is paper this of insight key The i 87 h 3 i = hre .Clark W. Charles , h SDW The . bhsbe oldna unu eeeay[5] degeneracy quantum near cooled been has Rb ℓ nlrnraiaingopaayi,we analysis, group renormalization onal P og ao nvriy ara,V 22030 VA Fairfax, University, Mason eorge mregnrclyi ioa Fermi dipolar in generically emerge 1 = s episwt eoaglrmomentum, angular zero with pairs le re aaee fteeeoi SDWs exotic these of parameter order aeol enhpteie ntwo- in hypothesized been only have , wv n a and -wave atc od ahrta nlattice on than rather bonds lattice n ovninlsprodcos While superconductors. conventional α emosa affiln safunction a as half-filling at fermions r h mdtcmeiinwt superfluid with competition amidst nite n ˆ i,α i,Rvrie A92521 CA Riverside, nia, natfroantcmtlfor metal antiferromagnetic in i ℓ = s aebe icvrd their discovered, been have + P d hw nFg ()cnan nex- an contains 1(d) Fig. in shown α d h xy a ˆ 4 i,α † a ra Zhao Erhai , i,α ( aecmoet n t or- its and component, wave † a ˆ ) p i,α and - stefrincannihila- fermionic the is σ i η respectively. with s d wv Ds They SDWs. -wave n CDW and h 1 r a ˆ , η 4 i ,j i, i,α † − { ∈ α correspond- ) σ ˆ r αβ η j {↑ ∈ ,y z y, x, a ˆ ˆ = s 23 i,β r de- are p Na , x i -wave ↓} } ˆ + 899 40 and ny / are at 2) y K . 2

FIG. 1: Schematic spin and charge order for pseudo-spin 1/2 dipolar fermions in two dimensions. The central image illustrates the onsite interaction U between opposite spins (red arrows) and the dipolar interaction Vr, here assumed 2 3 to be spin-independent. The characteristic scale for Vr is Vd = d /a , with a the lattice spacing. Vr depends also sensitively on the orientation of the dipoles labelled by angles (θ,φ). (a) The conventional antiferromagnetic spin density wave (SDWs). (b) Checkerboard charge density wave (CDWs), where “charge” is defined as the total density. (c) An example of p-wave spin density waves (SDWp) with modulation of y bond variables. (d) An example of mixed (extended) s- and d-wave spin density waves (SDWs+d). Red arrows in (c) and (d) indicate the direction of the spin vector S defined on the bonds (yellow ellipsoids). dipolar Fermi gas [24–26] in an optical square lattice at To give a heuristic argument about possible orders of half filling. The two pseudo-spin states can be two hy- the system, we consider a simplified version of model perfine states of Dy atoms, or two rovibrational states (1) retaining only the nearest and next-nearest neigh- of KRb molecules. This provides a tunable platform for bor dipolar interactions, denoted Vxˆ(ˆy) and Vxˆ+ˆy respec- quantum simulation of interacting fermions with long- tively, see Fig. 1. First, for dˆ=z ˆ, dipolar interactions are range interactions [21, 22], beyond the Fermi-Hubbard purely repulsive. For U ≫ Vd, the Hamiltonian reduces model. The system is described by the Hamiltonian to the Fermi-Hubbard model, implying a ground state with SDWs order at half-filling, Fig. 1(a). For U ≪ Vd, U Hˆ = − taˆ† aˆ + nˆ nˆ + V nˆ nˆ . (1) the dipolar energy is reduced by placing same-spins on X j,σ i,σ 2 X i,σ i,−σ X ij i j hi,ji,σ i,σ i6=j diagonally opposite sites, while opposite spins share the same site with only a small energy cost U. This implies a The lattice is aligned along the x- and y-directions, with checkerboard modulation of the total density ni = hnîi, nearest neighbor hopping t and on-site interaction U. U i.e. CDWs, shown in Fig. 1(b). contains contributions from the bare short range inter- As dˆ is tilted away fromz ˆ towards thex ˆ-direction, ⊥ action, and the on-site dipolar interaction Vii , defined there exists a region of tilting direction for which the below. We assume that all dipoles are aligned in the nearest neighbor interaction Vxˆ becomes attractive while ˆ same direction d = dd = (d, θ, φ) by an external mag- Vyˆ and Vxˆ+ˆy remain repulsive. For instance, for φ = 0, netic (or electric) field. In general, the off-site dipole- this region is bounded by two critical values of θ: ϑc1 = −1 ◦ −1 ◦ dipole interaction can be decomposed into equal- and cos (p2/3) ≈ 35 and ϑc2 = sin (p2/3) ≈ 54 . In unequal-spin components, labeled by k and ⊥, respec- the simpler case of spinless dipolar fermions, a checker- k ⊥ ˆ tively, Vij nîσnˆj,σ + Vij nîσnˆj,−σ, and depends on d and board bond order solid is formed [8] in this region. Then k(⊥) k(⊥) k(⊥) it is plausible that for the spin 1/2 case, unconventional r = r − r via Vr (dˆ) ≡ V (dˆ)= hij|V (dˆ)|iji = i j ij dd SDWs of non-s wave symmetry may be stabilized by k(⊥) ˆ 2 3 ⊥ ˆ Vd [1 − 3(ˆr · d) ]/r . We will mostly assume Vd (d)= interaction-induced correlated hopping either along the k ˆ ˆ Vd (d) ≡ Vd(d), as in Eq. 1, which arises naturally when xˆ,y ˆ, or the diagonalx ˆ +y ˆ direction. The spatial sym- the two states are associated with the same hyperfine metry of these SDWs depends on the value of φ. This ⊥ ˆ k ˆ manifold. The Vd (d) 6= Vd (d) case will be discussed scenario is illustrated in Figs. 1(c) and 1(d). briefly further down. Finally, for large dipole tilting angles, e.g., θ > ϑc2 3

FIG. 2: Phase diagram of dipolar fermions on a square lattice at half filling obtained from FRG. It is shown on the surface of a sphere as a function of the dipole orientation angle θ and φ for fixed interactions (a) Vd = 0.5, U = 0.1; and (b) Vd = 0.5, U = 0.5 in units of t. Unconventional SDWs (SDWp and SDWs+d) are sandwiched between the CDW and BCS superfluid phase. The FRG eigenwavefunction corresponding to a representative point (marked by “∗”) in each phase is shown in the kx − ky polar plots with matching colors. Note that FRG predicts a mixed s + d-wave (rather than pure d-wave) SDW. As the on-site interaction is increased from U = 0.1 in (a) to U = 0.5 in (b), the SDWs+d phase expands and squeezes out the neighboring phases. The s-wave component of SDWs+d increases with U, while the d-wave component diminishes. for φ = 0, the dominant dipolar interaction is attractive. ing momenta. (3) Finally we identify the most dominant The leading instability is towards formation of Cooper instability of the Fermi surface from the most divergent pairs. Again, the precise orbital symmetry of the BCS eigenvalue of the interaction matrix. The correspond- phase is determined by the value of φ. ing eigenvector provides information about the orbital We now determine the phase diagram of these intricate symmetry of the incipient order parameter. The phase competing orders in the weak coupling limit, {U, Vd}

SDWs+d takes the form of Fig 1(d): atoms delocalize across a plaquette, in the diagonal direction perpendic- ular to the dipole tilting direction. In contrast to the triplet BCSp phase at small φ, the BCSs+d phase is a superfluid of singlet Cooper pairs with mixed orbital symmetry, ℓ =0, 2. Next we illustrate how the phase diagram changes as the model approaches the repulsive Fermi-Hubbard model (U > 0, Vd = 0). We calculate the FRG flows for increased on-site interaction, U = 0.5, while keep- ing Vd fixed at 0.5. The phase diagram is shown in Fig. 2(b). Since the on-site interaction U favors antiferromagnetism, the SDWp phase shrinks, while the SDWs+d phase extends to cover a broader region, in- cluding that previously occupied by BCSs+d. Note that the d-wave component of SDW, even though diminished, is always present since the dipole interaction is kept finite. When U is further increased such that U ≫ Vd, only the isotropic component (c0) will survive, indicat- ing the SDWs+d crosses over to SDWs, the conventional antiferromagnetic ordering of spins in Fig. 1(a). To corroborate the FRG prediction of the unconventional spin density waves, we use self-consistent mean field theory. For a square lattice of finite size N × N, we impose periodic boundary conditions and retain the dipole interactions up to a distance of 12 lattice con- stants. We define the various normal and anomalous av- † FIG. 3: Mean field results for the SDWp phase. (a) erages, ρi,σ,j,σ′ = haˆ ′ aî,σi and mi,σ,j,σ′ = haî,σaˆj,σ′ i. S j,σ Preferred direction of the spin polarization vector S as a The corresponding mean field Hamiltonian is solved self- k ⊥ function of the ratio Vd /Vd . It is along the z-direction for consistently by starting from an initial guess of the gen- ⊥ k ⊥ k Vd < Vd and lies in the x-y plane for Vd > Vd . (b) The eralized density matrix, and iterating until desired con- Sz energy difference between the mean field states with and vergence is reached. At each step the chemical poten- Sx order. The inset shows the magnitude of the correspond- k tial is tuned to maintain half filling. The results are ing order parameter. The parameters are U = 0.1, Vd = 1.0, ◦ checked to be size-independent by varying N > 24. In θ = 47 , φ = 0. All energies are in units of t. (c) The mean Fig. 3, N is set to 28. Although mean field results are field energy gap of the SDWp phase, in units of the Fermi energy, as a function of the inter-spin long-range dipolar in- only suggestive, they provide an independent confirma- ⊥ ◦ teraction Vd for θ = 47 and φ = 0. For the SU(2) symmetric tion of the FRG results and unveil the real space pat- case plotted in (c), the energy gaps for different vector polar- terns of S directly in the SDWν phases. They can also izations are degenerate. be used to investigate the direction of S for the general- ⊥ ˆ k ˆ ized model with Vd (d) 6= Vd (d). We search for unconventional SDW phases with homogeneous spin density, and BCSs d, all of which contain a dxy-wave as well as + ni,σ = 1/2. In and around the SDWp region predicted s-wave components. The contributions of the isotropic by FRG, we indeed find solutions with order parameter s-wave, extended s-wave, and d-wave components are η † η S = haî,ασˆαβaˆj,βi, rj − ri =y ˆ. Further, the mean field inferred by fitting the FRG wavefunctions using func- energy for Sx order is identical to that for Sy and Sz, due tion c + c cos k cos k + c sin k sin k , with {c ,c ,c } k 0 1 x y 2 x y 0 1 2 to the SU(2) symmetry of Hˆ imposed by V ⊥(dˆ)= V (dˆ). as fitting parameters. As a general trend, for increas- d d This degeneracy is lifted for V ⊥(dˆ) 6= V k(dˆ). In Fig. 3, ing θ, the magnitude of isotropic s-wave c0 reduces, d d we compare the mean-field energies of the SDWp solu- while the magnitudes of c1 and c2 are comparable and tion with order parameter Sz and Sx. The z(x)- polar- increase. The CDWs+d phase can be viewed as the z x k ⊥ natural continuation of the CDWs as c1 and c2 be- ized order S (S ) is energetically favored for Vd > Vd k ⊥ x come appreciable. The two representative points shown (Vd < Vd ). However, the degeneracy between S and y in Fig. 2(a) for the SDWs+d and BCSs+d are fit by S remains. The mean field results support our inter- 0.05 − 0.16 cos kx cos ky − 0.18 sin kx sin ky and 0.01 + pretation of the SDWp order as schematically shown in 0.23 cos kx cos ky − 0.19 sin kx sin ky, respectively. Since Fig. 1(c). A similar analysis can be performed for the c0 is small, the real space modulation pattern for such SDWs+d phase. 5

In conclusion, we have established the emergence of [11] Parish, M. M. & Marchetti, F. M., Density Instabilities in unconventional spin density wave orders, SDWp and a Two-Dimensional Dipolar Fermi Gas. Phys. Rev. Lett. SDWs+d, along with other exotic phases with non-zero 108, 145304 (2012) angular momentum, within ultra-cold spin-1/2 dipolar [12] Gadsbolle, A. -L. & Bruun, G. M., Dipolar fermions fermions on the square lattice. These phases occupy in a two-dimensional lattice at non-zero temperature. arXiv:1206.4974 a sizable region of the phase diagram mapped out via [13] Baranov, M. A., Dobrek, L., & Lewenstein, M., Superflu- the functional renormalization group approach. Further- idity of Trapped Dipolar Fermi Gases. Phys. Rev. Lett. more, the self-consistent mean field estimation for the 92, 250403 (2004). energy gaps of SDWp, shown in Fig. 3(c), indicates a [14] Baranov, M. A., Mar’enko, M. S., Rychkov, V. S., & critical temperature Tc ∼ 0.08TF for Vd/t =1.5. Consid- Shlyapnikov, G. V. Superfluid pairing in a polarized dipo- ering the currently reported temperature of degenerate lar Fermi gas. Phys. Rev. A 66, 013606 (2002). dysprosium T ≈ 0.2T [4], this suggests experimental [15] Bruun, G. M. & Taylor, E., Quantum Phases of a Two- F Dimensional Dipolar Fermi Gas. Phys. Rev. Lett. 101, accessibility of these emergent phases in the near future. 245301 (2008). Our study thus articulates the dipolar Fermi system as [16] Cooper, N. R. & Shlyapnikov, G. V., Stable Topological an intriguing and novel test bed for exotic many-body ef- Superfluid Phase of Ultracold Polar Fermionic Molecules. fects. They provide a fresh perspective on systems with Phys. Rev. Lett. 103, 155302 (2009). competing orders. [17] Zhao, C.,Jiang, L., Xunxu, L., Liw, W. M., Zou, X., & We thank Benjamin Lev for invaluable discussions. SB Pu, H., Hartree-Fock-Bogoliubov theory of dipolar Fermi and EZ are supported by AFOSR (FA9550-12-1-0079) gases. Phys. Rev. A 81, 063642 (2010). [18] Fregoso, B. M., Sun, K., Fradkin, E., Lev, B. L., Biax- and ONR (N00014-09-1-1025A). LM acknowledges sup- ial nematic phases in ultracold dipolar Fermi gases, New port from the Landesexzellenzinitiative Hamburg, which Journal of Physics 11, 103003 (2009). is financed by the Science and Research Foundation [19] Quintanilla, J., Carr, S. T., & Betouras J. J., Metane- Hamburg and supported by the Joachim Herz Stiftung. matic, smectic, and crystalline phases of dipolar fermions SWT acknowledges support from NSF under grant DMR- in an optical lattice. Phys. Rev. A 79, 031601(R) (2009). 0847801 and from the UC-Lab FRP under award number [20] Sun, K., Wu, C., & Sarma, S. D., Spontaneous inhomo- 09-LR-05-118602. geneous phases in ultracold dipolar Fermi gases. Phys. Rev. B 82 075105 (2010). [21] Micheli, A., Brennen, G. K., & Zoller, P., A toolbox for lattice-spin models with polar molecules. Nature Physics 2, 341 (2006). [22] Gorshkov, A. V., Manmana, S. R., Chen, G. Ye, J., [1] Overhauser, A. W., Spin Density Waves in an Electron Demler, Lukin, M. D. & Rey, A. M., Tunable Super- Gas, Phys. Rev. 128, 1437 (1962). fluidity and Quantum Magnetism with Ultracold Polar [2] Fawcett, E., Spin-density-wave antiferromagnetism in Molecules. Phys. Rev. Lett. 107, 115301 (2011). chromium, Rev. Mod. Phys. 60, 209 (1988) [23] He, L. & Hofstetter, W., Supersolid phase of cold [3] Nayak, C., Density-wave states of nonzero angular mo- fermionic polar molecules in two-dimensional optical lat- mentum. Phys. Rev. B 62, 4880 (2000). tices. Phys. Rev. A 83, 053629 (2011). [4] Lu, M., Burdick, N. Q., & Lev, B. L., Quantum degen- [24] Fregoso, B. M. & Fradkin, E., Ferronematic Ground State erate dipolar Fermi gas. Phys. Rev. Lett. 108, 215301 of the Dilute Dipolar Fermi Gas. Phys. Rev. Lett. 103, (2012). 205301 (2009). [5] Ni, K. -K., Ospelkaus, S., de Miranda, M. G. H., [25] Sogo, T., Urban, M., Schuck, P., & Miyakawa, T., Spon- Pe’er, B., Neyenhuis, B., Zirbel, J. J., Kotochigova, S., taneous generation of spin-orbit coupling in magnetic Julienne, P. S., Jin, D. S. & Ye, J., A High Phase-Space- dipolar Fermi gases. Phys. Rev. A 85, 031601(R) (2012). Density Gas of Polar Molecules. Science 322, 231 (2008). [26] Li, Y. & Wu, C., Spin-orbit coupled Fermi liquid theory [6] Wu, C. -H., Park, J. W., Ahmadi, P., Will, S., & Zwier- of ultracold magnetic dipolar fermions. Phys. Rev. B 85, lein, M. W., Ultracold Fermionic Feshbach Molecules of 23 40 205126 (2012). Na K. arXiv:1206.5023. [27] Shankar, R., Renormalization-group approach to inter- [7] Leggett, A. J., A theoretical description of the new 3 acting fermions. Rev. Mod. Phys. 66, 129 (1994). phases of liquid He, Rev. Mod. Phys. 47, 331414 (1975). [28] Zanchi, D. & Schulz, H. J., Weakly correlated electrons [8] Bhongale, S. G., Mathey, L., Tsai, S. -W., Clark, C. W., on a square lattice: Renormalization-group theory. Phys. & Zhao, E., Bond Order Solid of Two-Dimensional Dipo- Rev. B 61,13609 (2000). lar Fermions. Phys. Rev. Lett. 108, 145301(2012). [29] Mathey, L., Tsai, S. -W., & Castro Neto, A. H., Com- [9] Yamaguchi, Y., Sogo, T., Ito, T. & Miyakawa, T., peting Types of Order in Two-Dimensional Bose-Fermi Density-wave instability in a two-dimensional dipolar Mixtures. Phys. Rev. Lett. 97, 030601 (2006); Fermi gas. Phys. Rev. A 82, 013643 (2010). [30] Mathey, L., Tsai, S. -W., & Castro Neto, A. H., Exotic [10] Mikelsons, K. & Freericks, J. K., Density-wave patterns superconducting phases of ultracold atom mixtures on for fermionic dipolar molecules on a square optical lat- triangular lattices. Phys. Rev. B 75, 174516 (2007). tice: Mean-field-theory analysis. Phys. Rev. A 83, 043609 (2011).