arxiv: v2 [cond-mat.quant-gas] 1 Feb 2018

Transcription

1 Critical behavior of a chiral superfluid in a bipartite square lattice arxiv:78.755v [cond-mat.quant-gas] Feb 8. Introduction Junichi Okamoto,, Wen-Min Huang 3,, Robert Höppner, and Ludwig Mathey, Zentrum für Optische Quantentechnologien and Institut für Laserphysik, Universität Hamburg, 76 Hamburg, Germany The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 49, 76 Hamburg, Germany 3 Department of Physics, National Chung-Hsing University, Taichung 47, Taiwan wenmin@phys.nchu.edu.tw Abstract. We study the critical behavior of Bose-Einstein condensation in the second band of a bipartite optical square lattice in a renormalization group framework at one-loop order. Within our field theoretical representation of the system, we approximate the system as a two-component Bose gas in three dimensions. We demonstrate that the system is in a different universality class than the previously studied condensation in a frustrated triangular lattice due to an additional Umklapp scattering term, which stabilizes the chiral superfluid order at low temperatures. We derive the renormalization group flow of the system and show that this order persists in the low energy limit. Furthermore, the renormalization flow suggests that the phase transition from the thermal phase to the chiral superfluid state is first order. Unconventional Bose-Einstein condensates (BECs) whose order parameter space is not simply the usual U() symmetry have been extensively studied. Examples from the field of ultracold atoms include Floquet engineered Bose gases [,, 3, 4] or spinor Bose gases [5], where the order parameter space has an additional Ising component or even more complex symmetry groups. The experimental realization of such systems in ultracold atomic systems are ideal to investigate phase transitions of those complex orders due to well-defined and tunable nature of these systems. Recently, BECs that break time-reversal (TR) symmetry have attracted increased attention from theorists [6, 7, 8, 9,,,, 3, 4, 5, 6, 7, 8, 9] and from experimentalists [,,, 3, 4]. According to Feynman s no-node theorem [5], such states cannot be a ground state of a conventional bosonic Hamiltonian with shortrange interactions, since breaking TR symmetry inevitably leads to a wave function with a node in real space. Therefore, to create BECs without TR symmetry, the assumptions of the no-node theorem have to be circumvented. One approach uses a long-lived

2 Critical behavior of a chiral superfluid in a bipartite square lattice metastable state of ultracold bosons in bands of higher orbitals []. Experimentally, a BEC in a p-band has been realized by first populating particles in a staggered pattern in a checkerboard lattice and then suddenly changing the potential shape. Due to the large anharmonicity in the energy spectrum, the life time of the metastable BEC is longer than ms [, ]. Since this BEC is not a ground state, breaking TR symmetry is not in contradiction with the no-node theorem. Indeed, a complex coherent superposition of the two p-band condensates (i.e., p x ± ip y order) that breaks TR symmetry has been realized by carefully tuning the lattice parameters [, 3, 4]. Since such a state hosts spatially staggered orbital currents, it is dubbed a chiral superfluid. Similar p x ± ip y paring has been proposed for the A phase of superfluid 3 H [6, 7] and for Sr RuO 4 [8]. In this paper, we investigate the stability of the chiral condensate and its critical behaviors by a renormalization group (RG) analysis. The analysis addresses the competition of chiral and non-chiral condensation, and the critical behavior. While the chiral superfluid has been confirmed experimentally, it is still important to know how stable and general the state is. In particular, since the energies of the competing non-chiral BEC and of the chiral BEC are close at the mean-field level, it is not trivial which of the two BECs becomes dominant at low temperatures. We find that the stable condition of the chiral BEC is always preserved at low energy scales, and the transition is expected to be first-order. The paper is organized as follows. In section we develop the field theoretical description of the mixed orbital model in a bipartite optical square lattice in the low energy limit. Section 3 is devoted to a mean-field analysis of the effective model. In section 4, we study the critical behavior of the model in the framework of a one-loop RG calculation. In section 5, we conclude. The details of calculations not covered in the main texts are summarized in appendix.. Effective field theory The system that we consider here is sketched in figure (a) and described by the Hamiltonian, { ( ) } H = dz b i (r, z) z µ 3D b i (r, z) + H xy, () m r,i with the tight-binding model of a bipartite optical square lattice H xy = J [ b (r, z)b (r + d, z) + b (r, z)b 3 (r + d, z) r ] b (r, z)b (r d, z) b (r, z)b 3 (r d, z) + h.c.. () Here r = a(n x, n y ) with n x,y Z, d / = (a/, ±a/ ). b i (r, z) with i =,, 3 represent the annihilation operators of bosons at the s, p x and p y orbitals respectively. We assume that bosons move freely along the z-direction. The hopping amplitudes between neighboring p-orbitals, J and J, are set to be zero for simplicity (see appendix

3 Critical behavior of a chiral superfluid in a bipartite square lattice 3 (a) z y x (b) (c) y p x p y J + J + + -J + J a x s Figure. (a) A bipartite optical square lattice is sketched in the xy plane with s-and p-orbitals arranged in a chequerboard pattern. Along the z-direction, bosons move freely. (b) The hopping amplitudes between orbitals on different sites, and the lattice constant a are illustrated. (c) The lowest band of the bipartite optical square lattice in units of E rec is plotted for the hopping amplitudes J =.3E rec, J = and J = ; these choices are motivated by the experimental setup in [3]. for a more general discussion). Converting the orbital representation into a band representation, we obtain three bands. The metastable BEC in experiments is loaded in the lowest band, whose dispersion is given by ( ) ( ) ɛ(k, k z ) = J cos kx cos ky + kz, (3) m where k = (k x, k y ) and we set a = in the following calculations. We illustrate the lowest band in momentum space in figure (c). We note that there are two energetic minima at k = (π/, ) and k = (, π/ ). These energetic minima are degenerate, thus giving rise to the Z symmetry of the noninteracting Hamiltonian. At low temperatures, bosons predominantly occupy momentum states near the two minima, and then condense below a critical temperature. To describe the critical behavior, we expand the bosonic operators near the two minima as b α (r, z) = e ik r u α (k)φ(k, z) N k a N j=, e ik j r u αj q j <Λ q d q j 4π eiq j r φ j (q j, z) a ψ j (r, z)u αj e ikj r, (4) j=, where ψ j (r, z) a N q j <Λ q d q j 4π e iq j r φ j (q j, z) with φ(k j + q j, z) φ j (q j, z). Λ q is the momentum cut-off, and N is the number of the unit cells. The kernel u α (k j ) = u αj represents the projection of the wave function of orbital α on the wave function of

4 Critical behavior of a chiral superfluid in a bipartite square lattice 4 the lowest band in the vicinity of the minimum k j. These are given by (u, u ) = ( i/, i/ ), (u, u ) = (/, /) and (u 3, u 3 ) = (/, /). We use the field decomposition to approximate the full Hamiltonian of equation () by an effective Hamiltonian with two components, H eff = j=, d 3 R ψ j (R) [ m R µ j ] ψ j (R), (5) where R = (r, z), µ j being the chemical potential of the j-th component and the effective mass being m = (m xym ) /3. As an example, we describe the experimental parameters of []; for 87 Rb atoms and for typical laser intensities that were used, we have m.m xy and m xy = / (λ L J) with the laser wavelength λ L = 64nm and J/E rec =.3. We note that to simplify the RG analysis we use an isotropic effective model, and the momentum cut-off in the field decomposition sets the energy cut-off of the effective Hamiltonian as ɛ Λ = Λ q/m. We further consider the on-site interaction; see [3], which gives the following terms, H I = dz { Us n s(r) [n s (R) ] r + U p n p(r) [n p (R) ] U p [ L z (R) n p (R) ] }, (6) where n s = b b, n p = b b +b 3b 3, and L z = i (b b 3 b 3b ) being an angular momentum operator. U s is the on-site interaction among s-orbitals. U p and U p are intra- and interorbital interactions among p-orbitals. In the tight-binding approximation, the strength of the on-site interactions can be calculated from the contact interaction and the Wannier functions [9]; for the details, see appendix. In the standard harmonic approximation, the on-site interactions follow U p = U p /3. The precise ratio between U s and U p depends on the depth of the optical potential since the harmonic frequency of the s-orbital sites are different from the one of the p-orbital sites. For a moderately deep potential, we find U s U p. Within the field-theory approximation (4), we represent the effective interaction as { HI eff = d 3 g R ψ (R)ψ (R)ψ (R)ψ (R) + g ψ (R)ψ (R)ψ (R)ψ (R) + g ψ (R)ψ (R)ψ (R)ψ (R) + g u [ ψ (R)ψ (R)ψ (R)ψ (R) + H.c.] }, (7) where g j (j =, ) is the intra-component interaction and g is the inter-component one. In this expression, there is an additional term with the interaction strength g u, which is an Umklapp term. This additional scattering process is not present in the previously studied triangular lattice system [9, 3, 3], which demonstrates that these two systems are in different universality classes. Our model is more general in the sense that three coupling constants flow independently under the renormalization equations. The Umklapp interaction allows interchange of bosons between the two components by lattice assisted collisions. In other words, the effective interaction only enforces

5 Critical behavior of a chiral superfluid in a bipartite square lattice 5 conservation of the total boson number, in stead of the boson number of each component as, for instance, in the frustrated triangular optical lattice [3]. Equivalently, we only have one global U() symmetry, instead of two U() symmetries for each component. The bare values of the coupling constants in terms of U s, U p and U p are: g / = U s + U p + 3U p, g = U s + U p U p, g u = 8 4 g. (8) The full symmetry of the effective action is U() Z Θ, where the U() symmetry corresponds to the invariance of the model under the global phase shift for both components, Z is the exchange of the two components, and Θ is the time-reversal symmetry. 3. Mean-field theory Before proceeding to the RG analysis, we study the ground state for the bare interactions within a zero-temperature approach. We assume that the bosons perfectly condense at the two energetic minima so that a many-body trial wave function is represented as, Ψ θ,φ = N! [cos θψ + e iφ sin θψ ] N,, (9) where N is the total number of bosons, and m, n stands for bosons occupation numbers m(n) at momentum k (k ) respectively. We will use the angles θ and φ as variational parameters. θ determines the relative population of the two minima and φ denotes the relative phase of the two-component condensates. Using the trial wave function, we compute the energy of the interacting effective Hamiltonian (7) as { H eff = N(N ) g I θ,φ cos4 θ + g sin4 θ + g cos θ sin θ } + g u cos θ sin θ cos φ. () First we note that for the frustrated triangular optical lattice, the Umklapp interaction does not occur ( g u = ), and we have g = g g, which follows from the common origin of these terms, i.e., the contact interaction between the atoms. In this case, the minimum of equation () occurs at θ = or π/. From equation (9), this means that bosons will condense in one of the energetic minima to break the Z symmetry [3]. However, in the square bipartite lattice that we study in this paper, the Umklapp interaction g u > exists due to the bare on-site repulsive interactions U s U p U p >. In this case, another energetic minimum may appear at (φ, θ) = (±π/, π/4) in equation [ N (); this corresponds to a chiral superfluid state Ψ = ψ ± iψ ], / N! given by a complex coherent superposition of two single particle states. This state breaks the time-reversal symmetry Θ, i.e., the chiral Z symmetry, in addition to the U() continuous symmetry of the phase (The situation is similar to the fully frustrated XY models [3]). Comparing this to the single condensate at θ =, we find that the chiral

6 Critical behavior of a chiral superfluid in a bipartite square lattice 6 superfluid state occurs when G g g + g u >, () where we set g = g = g. The above condition is also discussed in [9]. In terms of the interaction parameters in equation (6), we find g g + g u = U p/ >, and thus a chiral superfluid order occurs for any repulsive interaction. Even if we include non-zero values of J and J, we find that the condition is still satisfied as long as the energetic minima are located at k = (π/, ) and k = (, π/ ) (see appendix). However, the energy difference between the normal and chiral condensates is of the order of U p at the mean-field level, and at low temperatures the coupling constants get renormalized under the RG flow. Then it is nontrivial which superfluid order emerges at low temperatures. To study this problem more systematically, and to study the critical behavior in the low-energy limit, we apply the renormalization group method in the next section. 4. One-loop renormalization group method Following the RG method employed in [3] and ignoring quantum fluctuations, we calculate the RG equations at one-loop order as [33], dµ Λ = µ Λ T Λ F (µ Λ ) (g + g ), (a) dl dg = ɛg T Λ F (µ Λ ) ( 5g + g + g dl u), (b) dg = ɛg T Λ F (µ Λ ) ( 4g g + g + 4g dl u), (c) dg u = ɛg u T Λ F (µ Λ ) (g g u + 4g g u ), (d) dl where l = ln(λ q /Λ b ) is the logarithm of the ratio between the bare momentum cutoff Λ q and the running cutoff Λ b. Here ɛ = 4 d with d being the spatial dimension of the system, and ɛ = for our three-dimensional model. We have also defined dimensionless parameters, µ Λ = µ /ɛ Λ = µ /ɛ Λ, T Λ = k B T/ɛ Λ, g i = g i Λ 3 q/(π ɛ Λ ), and F (µ Λ ) = /( µ Λ ). If g u =, the flow equations become the ones of the twocomponent φ 4 -theory [3]. In the critical regime, µ Λ, where we can approximate F (µ Λ ), the RG equations exhibit four fixed points, except the trivial one, (µ Λ, g, g, gu) = (,,, ), ( ) (µ Λ, g, g, gu) = 5,,,, (3a) 5T ( Λ ) (µ Λ, g, g, gu) = 4,,,, (3b) 6T Λ 6T ( Λ ) (µ Λ, g, g, gu) = 5,,,, (3c) T Λ 5T Λ T Λ This corresponds to the lowest order of the ɛ-expansion [3, 34, 35, 36]. The effect of higher order contributions to the fixed points is discussed in appendix.

7 Critical behavior of a chiral superfluid in a bipartite square lattice (a) g.. (b) g.. 7. gu.. gu.. g. (c) g. g. (d)... g. gu.... gu. g... g Figure. (a) Three planes separating the RG flow in the g -g -gu space for gu >. The four fixed points are plotted as red dots. (b) The RG flow on the gu = plane, which are also introduced in [3]. (c) RG flows on the G = plane. (d) RG flows on the G = plane. (µ Λ, g, g, gu ) =,,, 5 5. (3d) All these fixed points are unstable, indicating that the system undergoes a first-order transition. In figure (a), we show the fixed points as red points for gu. 4.. Basic structures of RG equations While the full RG equations are complicated and can only be solved numerically, the basic structure of the equations give useful insight in the RG flow. In particular, we find three separatrix surfaces of the RG equations. RG flow cannot pass through these surfaces, and therefore the asymptotic behavior of the RG flow is severely constrained by the initial condition. The first separatrix is gu = ; the formal solution of equation (d) is Z l gu (l) = gu () exp dl F (µλ ) (g + 4g ), (4) which explicitly shows that gu does not change signs under the RG. In figure (b), we show the flow diagram on the gu = surface, which is also introduced in [3]. The two fixed points (g, g, gu ) = (,, ) and (/(5 ),, ) are unstable, and the one at (/(6 ), /(6 ), ) is marginally unstable.

8 Critical behavior of a chiral superfluid in a bipartite square lattice 8 The second separatrix surface is G g g + g u =. The RG equation for G is dg = G T Λ F (µ Λ ) G (5g + g 3g u ). (5) dl We emphasize that G > coincides with the mean-field stable condition for a chiral superfluid, (), and therefore, for general repulsive interactions that give G > as an initial condition, a chiral superfluid order always persists in the low-energy limit. In figure (c), we illustrate the flow diagrams on the G = plane for g u >. There are two unstable fixed points, (3b) and (3c), on the plane, in addition to the trivial one (g, g, gu) = (,, ). In particular the ray from (g, g, gu) = (,, ) to (/T Λ, /5T Λ, /T Λ ) separates the RG flows into two parts. For initial conditions g > g u on the G = plane, we find that system eventually flows into the fixed point (g, g, gu) = (/6T Λ, /6T Λ, ), where the system is reduced to the standard twocomponent φ 4 theory. However, if the initial conditions deviate from the G = plane by an arbitrarily small amount, g u will grow to a large positive value, and thus the fixed point (3b) is actually unstable. For initial conditions g < g u on the G = plane, the flow first approaches (/T Λ, /5T Λ, /T Λ ) and then runs away to larger positive values of g u. Finally the third separatrix surface is G g g g u =, and its RG equation is dg = G T Λ F (µ Λ ) G (5g + g + 3g u ). (6) dl The flow equation is similar to equation (5), and it guarantees that no flow passes through the G = plane. As illustrated in figure (d), we find that the Umklapp interaction g u always grows up on the G = plane. The increase of g u is also found when the initial couplings deviate from the G plane as we discuss below. 4.. RG flows for the effective model Now let us analyze the RG equations for the relevant parameter regime of our model. For our effective model, we have U s U p U p/3 >, and the RG flow is constrained to the space of g u >, G > and G <. Due to this constraint, as we will show below, the possible phases of our effective model are either the thermal gas phase (µ Λ ) or the chiral superfluid phase (µ Λ +). The phases are determined by the non-universal nature of the RG flows and initial conditions. First, when µ Λ remains positive, a typical flow of coupling constants behaves as in figure 3. For small positive initial interactions, a typical flow has monotonically increasing g u, which can also be seen from equation (4). For the evolution of g and g, there are three regimes that we can characterize: (i) < l < l (the initial regime): the linear terms in RG equations are dominant, and g and g gradually increase. (ii) l < l < l (the intermediate regime): the quadratic terms become more important with increasing g u, which eventually make g and g negative.

9 Critical behavior of a chiral superfluid in a bipartite square lattice 9 g l l l l g g u - Figure 3. A schematic typical RG flow for the three coupling constants when µ Λ remains positive. The characteristic length scales that separate the behavior of the flow are also plotted. (iii) l < l < l 3 (the asymptotic regime): the quadratic terms give asymptotically diverging behaviors, while the coupling constants are still smaller than the unity. Of course, we should stop the RG flow before any of the coupling constant becomes order of unity near l 3, above which the cubic terms become dominant. In the asymptotic regime, µ Λ approaches one, and the quadratic terms in the RG equations become dominant. Ignoring the linear terms, the flow can be analyzed by an ansatz g i (l) = kḡ i /( kl), where k is the inverse length scale at which these coupling constants diverge [37]. We find that the only asymptotic flow constrained in the space of g u >, G > and G < is given by ( (ḡ, ḡ, ḡ u ) = T Λ F (µ Λ ), 5, ). (7) This implies that G steadily increases and thus stabilizes the chiral superfluid order, (). At the same time, the quartic interactions g and g are renormalized to negative values. This indicates the breakdown of the quartic effective field theory, and we need to include higher order interactions such as g 6 d 3 R { ψ (R) 6 + ψ (R) 6}, (8) which is generated from three quartic vertices after one-loop renormalization. We note that at tree-level the g 6 contribution is marginal, and stabilizes the system [3]. At oneloop order, one contribution to the RG equation for g 6 is dg 6 /dl 4T Λ F (µ Λ ) g g 6. Since g flows to negative values, the above contribution further stabilizes the system. Second, when the chemical potential becomes negative, F (µ Λ ) = /( µ Λ ) gets suppressed, making the linear terms in the RG equations more dominant. Therefore, typical RG flows give a simple scaling behavior, (µ Λ, g i ) ( e l, e l). This corresponds to the thermal gas phase without condensation.

10 Critical behavior of a chiral superfluid in a bipartite square lattice Figure 4. Phase diagram for µ Λ =. and T Λ =.5. U s and U p are measured in units of Λ 3 q/(π ɛ Λ T Λ ). The phases are determined by the sign of the chemical potential, when the RG flow is terminated at g u =. Colors represent the final value of µ. The blue region is the condensed phase separated from the thermal phase with a first-order (discontinuous) phase transition. The pink area represents the thermal phase. As µ Λ increases or T Λ decreases, the condensed region becomes larger Phase diagram The RG analysis can be also used to study the critical behaviors of the phase transition between the two phases. In figure 4, we illustrate the phase diagram as a function of U s and U p. The pink region represents the thermal gas phase, in which µ Λ flow to negative values under RG transformations. Th blue region is the chiral p x ± ip y superfluid order with time-reversal symmetry breaking, where µ Λ approaches one. As µ Λ () increases or T Λ decreases, the phase boundary is shifted to enlarge the superfluid region. We note that along the asymptotic flow (7), the mean-field free energy can be written by two order parameters P ± ψ ± iψ as F Λ MF µ Λ ( P + P + ) ḡ( P 4 + P + 4 ) + O(P 6, P 6 +), (9) where ḡ > is the single coefficient characterizing the asymptotic flow (7), and the last term is the higher order correction stabilizing the system. This form suggests that the transition is first order. Another indirect support for this scenario is obtained by considering the strong coupling limit for g and g u, while g is set to. In this case, the model has effectively two XY spins on each space point that are orthogonal to each other due to the strong g u interaction. Such a model is known as Stiefel s V, model, and Monte Carlo studies show that it undergoes a first order transition [38, 39, 4]. For moderate interaction strength, we speculate that the first order nature becomes weaker. Clarifying if the transition remains first order for weak interactions seems to require an extensive numerical simulations, which are beyond the scope of the paper. The difficulty of determining the order of transitions is a common problem in the frustrated spin systems, whose effective Ginzburg-Landau models are similar to ours [3, 4].

11 Critical behavior of a chiral superfluid in a bipartite square lattice To detect this first order transition in experiments, we suggest two possibilities. The first one consists of measuring the condensate fraction as a function of temperature, preferably in a box potential [4, 4]. The measured temperature dependence will approach a non-zero jump at the condensation temperature for large systems. As a second approach we suggest to measure the spatial evolution of a phonon pulse in a condensate in a smoothly varying trap [43, 44, 45, 46]. Here, the phonon velocity will vary as the phonon pulse approaches the interface of the condensate and the thermal gas. For a second order transition the pulse velocity will smoothly approach zero, whereas for a first order transition the velocity will approach a non-zero value before the pulse is reflected, which gives a clear indication of a first order transition. 5. Conclusions In this paper we have investigated the critical behavior of unconventional Bose-Einstein condensates in the second band of an optical lattice. We have demonstrated that an Umklapp process between the two minima of the dispersion stabilizes a chiral superfluid state that breaks time reversal symmetry, first at the mean-field level and then within a renormalization group calculation. The latter shows that this stability is always persistent at low energy scales after integrating out thermal fluctuations. We obtain this result by identifying three separatrix planes in the RG flow, which constrain the low energy behavior to a stable regime. Furthermore, the RG flow suggests that the phase transition of the chiral superfluid state to the thermal state is of first order, in contrast to the usual second order transition of a conventional condensate. Acknowledgments We thank A. Hemmerich for helpful discussions on the experimental aspects of the system. J.O., R.H., and L.M. acknowledge the support from the Deutsche Forschungsgemeinschaft (through SFB 95 and EXC 74) and from the Landesexzellenzinitiative Hamburg, which is supported by the Joachim Herz Stiftung. W.M.H. especially acknowledges the support from Ministry of Science and Technology, Taiwan through Grant No. MOST 4--M-5-6-MY3. Appendix A. Effective interactions for general hopping amplitudes and interactions In this appendix, we start from a single particle picture in a bipartite optical lattice to derive the hopping amplitudes and interaction parameters of a Bose-Hubbard model. We then show that the condition for the chiral superfluid state G = g g + g u >

12 Critical behavior of a chiral superfluid in a bipartite square lattice in equation () is preserved as long as the band has two minima at k = (π/, ) and k = (, π/ ). Appendix A.. Derivation of Hubbard parameters We start from a following potential cos V (r) = V [k (x + y)/ ] [ + e iβ cos k (x y)/ ]. (A.) This has two local minima in a unit cell (see figure A(a)): one is a shallow local minimum hosting a s-orbital like Wannier state (denoted as A sites), and the other is a deep minimum hosting two p-orbital like Wannier sates (denoted as B sites). k = π/a with a being the distance between neighboring A and B sites (figure (b)). The energy difference between local minima at A sites and B sites is V = EA E B = 4V cos(β). We tune β so that the doubly degenerate first excited states in site B is close to the ground state in site A; in particular, we choose β so that the three bands are exactly degenerate at the Γ point in the following. After solving the single particle Schrödinger equation [ ] m + V (r) ψk(r) n = Ekψ n k(r), n (A.) we obtain a band dispersion as figure A(b). following tight-binding model H xy = J r A We fit the obtained dispersion by the [ b (r, z)b (r + d, z) + b (r, z)b 3 (r + d, z) ] b (r, z)b (r d, z) b (r, z)b 3 (r d, z) + h.c. J r B,ν=x,y [ b (r, z)b 3 (r + e ν, z) + b 3(r, z)b (r + e ν, z) + h.c. ] y (a) B B A x B B E(k) (b) Γ X M Figure A. (a) A bipartite optical potential V (r). (b) Numerically obtained band dispersions for the second composite bands (solid lines). Tight-binding fitting gives the dots.

13 Critical behavior of a chiral superfluid in a bipartite square lattice 3.. -ϵ A J -J -J V /E rec Figure A. The fitted tight-binding parameters as functions of the potential depth V in the logarithmic scale. J r B,ν=x,y i=,3 [ ] b i (r, z)b i(r + e ν, z) + h.c. + ɛ A b (r, z)b (r, z) + ɛ B b i (r, z)b i(r, z), r A r B i=,3 (A.3) where e = ( a, ) and e = (, a). The first term is already given in equation (). The degeneracy at the Γ point is achieved by setting ɛ B = ɛ A + 4J. The fitted band dispersion is plotted in figure A(b) as dots. The fitted parameters as functions of the potential depth V is given in figure A. We note that when the potential is relatively deep, the hopping between p-orbitals (J and J ) is much smaller than the one between s- and p-orbitals, J. Thus, in practice, we can ignore J and J. Now we use the Bloch wave functions for the obtained band dispersions to construct localized Wannier functions. Here we employ a simple projection approach [47]. These Wannier functions give the bare interactions of a Bose-Hubbard model as U s = g dz w z (z) 4 d r w (r) 4, (A.4) U p = g dz w z (z) 4 d r w /3 (r) 4, (A.5) U p = g dz w z (z) 4 d r w (r) w 3 (r), (A.6) where w z (z) is the Wannier function of a harmonic trap along the z-axis, and w i (r), i =,, 3 are the Wannier functions of the s-, p x and p y -orbitals on the xy-plane respectively. g is the contact interaction strength. The obtained values are plotted in figure A3. We find that U s U p U p/3 for moderately strong potential depth. Appendix A.. Effective interactions in field-theory approximations In this subsection, we show that the condition for the chiral superfluidity, (), is satisfied in general based on the Hubbard parameters determined above. With a general

14 Critical behavior of a chiral superfluid in a bipartite square lattice U s U p U p ' V /E rec Figure A3. Hubbard interactions obtained from localized Wannier functions as functions of the potential depth V. They are measured in units of E rec g dz w z (z) 4. The dashed lines are obtained under the harmonic approximation. dispersion in equation (A.3) and ɛ B = ɛ A + 4J, the projection of the Wannier orbitals to the two minima in the lowest band becomes (u, u ) = (λ, + λ λ ), (A.7) (u, u ) = (u 3, u 3 ) = (, ), (A.8) + λ (, ). (A.9) + λ with λ = i J J + J and J = J J. J With equation (8), we can show that the g g + g u = (A.) ( J 4U pj 4 [J + J [J ] J + J J ( ( J ) ) 4. (A.) 3 J + J J J For J J, J, the above quantity is always positive. Therefore, even when the system deviates from the simple limit of J = J =, the condition of the chiral superfluidity is still satisfied. )] Appendix B. The ɛ-expansion analysis of the fixed points In section 4, we have shown the fixed points that correspond to the lowest order of the ɛ-expansion [3, 34, 35, 36]. To estimate the high-order effect in the ɛ-expansion, we investigate the fixed points of the RG equations without expanding F (µ Λ ). Finding the

15 Critical behavior of a chiral superfluid in a bipartite square lattice 5 zeros of the right-hand sides of (a)-(d) leads to four fixed points, except the trivial one, (µ Λ, g, g, gu) = (,,, ), ( ) ɛ (µ Λ, g, g, gu) = 5 + ɛ, 5ɛ,,, (B.) (5 + ɛ) T ( Λ ) ɛ (µ Λ, g, g, gu) = 4 + ɛ, 8ɛ 8ɛ,,, (B.) 3(4 + ɛ) T Λ 3(4 + ɛ) T ( Λ ) ɛ (µ Λ, g, g, gu) = 5 + ɛ, 5ɛ 5ɛ 5ɛ,,, (B.3) (5 + ɛ) T Λ (5 + ɛ) T Λ (5 + ɛ) T ( Λ ) ɛ (µ Λ, g, g, gu) = 5 + ɛ, 5ɛ 5ɛ 5ɛ,,, (B.4) (5 + ɛ) T Λ (5 + ɛ) T Λ (5 + ɛ) T Λ Taking the lowest order in ɛ and setting ɛ = recovers the fixed points in (3a)-(3d). We emphasize that the above fixed points still lie on the three separatrix planes g u =, G = and G = in the g -g -g u space. This fact guarantees the persistence of the chiral superfluid in the low energy limit even when we include the higher order terms in ɛ in our one-loop RG study. [] Struck J, Ölschläger C, Targat R L, Soltan-Panahi P, Eckardt A, Lewenstein M, Windpassinger P and Sengstock K Science [] Struck J, Weinberg M, Ölschläger C, Windpassinger P, Simonet J, Sengstock K, Höppner R, Hauke P, Eckardt A, Lewenstein M and Mathey L 3 Nat. Phys [3] Parker C, Ha L and Chin C 3 Nat. Phys [4] Clark L W, Feng L and Chin C 6 Science [5] Kawaguchi Y and Uedaa M Phys. Rep [6] Isacsson A and Girvin S M 5 Phys. Rev. A [7] Liu W and Wu C 6 Phys. Rev. A [8] Kuklov A B 6 Phys. Rev. Lett [9] Wu C, Liu W, Moore J and Sarma S 6 Phys. Rev. Lett [] Lim L K, Smith C M and Hemmerich A 8 Phys. Rev. Lett. 34 [] Stojanovic V M, Wu C, Liu W V and Das Sarma S 8 Phys. Rev. Lett. 53 [] Wu C 9 Mod. Phys. Lett. B 3 [3] Lewenstein M and Liu W V Nat. Phys. 7 [4] Cai Z and Wu C Phys. Rev. A [5] Li X, Zhang Z and Liu W V Phys. Rev. Lett [6] Martikainen J P and Larson J Phys. Rev. A [7] Cai Z, Duan L M and Wu C Phys. Rev. A [8] Liu B, Yu X L and Liu W M 3 Phys. Rev. A [9] Li X and Liu W V 6 Reports Prog. Phys [] Ölschläger M, Wirth G and Hemmerich A Phys. Rev. Lett [] Wirth G, Ölschläger M and Hemmerich A Nat. Phys [] Soltan-Panahi P, Lühmann D S, Struck J, Windpassinger P and Sengstock K Nat. Phys. 8 7 [3] Ölschläger M, Kock T, Wirth G, Ewerbeck A, Morais Smith C and Hemmerich A 3 New J. Phys [4] Kock T, Ölschläger M, Ewerbeck A, Huang W M, Mathey L and Hemmerich A 5 Phys. Rev. Lett. 4 53

16 Critical behavior of a chiral superfluid in a bipartite square lattice 6 [5] Feynman R P 998 Statistical Mechanics: A Set of Lectures (New York: Avalon Publishing) [6] Leggett A J 975 Rev. Mod. Phys [7] Dobbs R Helium Three (Oxford: Oxford University Press) [8] Mackenzie A P and Maeno Y 3 Rev. Mod. Phys [9] Antonenko S A and Sokolov A I 994 Phys. Rev. B [3] Kawamura H 998 J. Phys. Condens. Matter 477 [3] Janzen P, Huang W M and Mathey L 6 Phys. Rev. A [3] Villain J 977 J. Phys. C Solid State Phys. 77 [33] Stoof H T C, Dickerscheid D B M and Gubbels K 9 Ultracold Quantum Fields (Dordrecht: Springer Netherlands) [34] Cardy J 996 Scaling and Renormalization in Statistical Physics (Cambridge: Cambridge University Press) [35] Wilson K G and Fisher M E 97 Phys. Rev. Lett. 8 4 [36] Domany E, Mukamel D and Fisher M E 977 Phys. Rev. B [37] Balents L and Fisher M P A 996 Phys. Rev. B [38] Kunz H and Zumbach G 993 J. Phys. A Math. Gen 6 3 [39] Loison D and Schotte K D 998 Eur. Phys. J. B [4] Itakura M 3 J. Phys. Soc. Japan 7 74 [4] Gaunt A L, Schmidutz T F, Gotlibovych I, Smith R P and Hadzibabic Z 3 Phys. Rev. Lett. 46 [4] Corman L, Chomaz L, Bienaimé T, Desbuquois R, Weitenberg C, Nascimbène S, Dalibard J and Beugnon J 4 Phys. Rev. Lett [43] Tey M K, Sidorenkov L A, Guajardo E R S, Grimm R, Ku M J H, Zwierlein M W, Hou Y H, Pitaevskii L and Stringari S 3 Phys. Rev. Lett [44] Sidorenkov L A, Tey M K, Grimm R, Hou Y H, Pitaevskii L and Stringari S 3 Nature [45] Weimer W, Morgener K, Singh V P, Siegl J, Hueck K, Luick N, Mathey L and Moritz H 5 Phys. Rev. Lett [46] Singh V P, Weimer W, Morgener K, Siegl J, Hueck K, Luick N, Moritz H and Mathey L 6 Phys. Rev. A [47] Marzari N, Mostofi A A, Yates J R, Souza I and Vanderbilt D Rev. Mod. Phys