in feynman s statistical mechanics textbook , it is stated that the many - body ground state wavefunctions of bosons are _ positive - definite _ in the coordinate representation provided no external rotation is applied and interactions are short - ranged @xcite .
the proof is very intuitive : due to the time - reversal ( tr ) symmetry , the ground state wavefunction @xmath0 can be chosen as real . if it is not positive - definite , _
i.e. _ , it has nodes , the following surgery can be done to lower its energy .
we first take its absolute value of @xmath1 , whose energy expectation value is exactly the same as that of @xmath0 .
however , such a wavefunction has kinks at node points . further smoothing the kinks results in a positive - definite wavefunction , and the kinetic energy
is lowered by softening the gradients of the wavefunction .
although the single body potential energy and the two - body interaction energy increase , they are small costs of a high order compared to the gain of the kinetic energy . we can further conclude that the ground state is non - degenerate because two degenerate positive - definite wavefunctions can not be orthogonal to each other .
this `` no - node '' theory is a very general statement , which applies to almost all the well - known ground states of bosons , including the superfluid , mott - insulating , density - wave , and even super - solid ground states .
it is also a very strong statement , which reduces the generally speaking complex - valued many - body wavefunctions to be positive - definite distributions .
this is why the ground state properties of bosonic systems , such as @xmath2he , can in principle be exactly simulated by the quantum monte - carlo method free of the sign problem .
furthermore , this statement implies that the ordinary ground states of bosons , including bose - einstein condensations ( bec ) and mott - insulating states , can not spontaneously break time - reversal ( tr ) symmetry , since tr transformation for the single component bosons is simply the operation of the complex conjugation
.    it would be exciting to search for exotic emergent states of bosons beyond this `` no - node '' paradigm , whose many - body wavefunctions can be complex - valued with spontaneous tr symmetry breaking . since properties of complex - valued functions are much richer than those of real - valued ones , we expect such states can exhibit more intriguing properties than the ordinary ground states of bosons .
for this purpose , we have recently made much progress with two different ways to bypass feynman s argument , including the meta - stable states of bosons in the high orbital bands of optical lattices @xcite and multi - component bosons with spin - orbit coupling linearly dependent on momentum @xcite .
clearly , the `` no - node '' theorem is a ground state property which does not apply to the excited states of bosons . the recent rapid development of optical lattices with ultra - cold bosons provides a wonderful opportunity to investigate the meta - stable states of bosons pumped into high orbital bands . due to the lack of dissipation channels
, the life time can be long enough to develop inter - site coherence @xcite .
we have shown that the interaction among orbital bosons are characterized by the `` _ orbital hund s rule _ '' @xcite , which gives rise to a class of complex superfluid states by developing the on - site orbital angular momentum ( oam ) moments . in the lattice systems ,
the inter - site hoppings of bosons lock the oam polarization into regular patterns depending the concrete lattice structures .
furthermore , the `` no - node '' theorem even does not apply to the ground states of bosons if their hamiltonians linearly depend on momentum , i.e. , the gradient operator .
a trivial example is the formation of vortices in bose condensates with the external rotation in which the coriolis force is represented as the vector potential linearly coupled to momentum .
however , in this case tr symmetry is explicitly broken .
a non - trivial example is that bosons with spin - orbit ( so ) coupling , whose hamiltonians also linearly depend on momentum and are tr invariant .
the invalidity of the `` no - node '' theorem can also give rise to complex - valued ground state wavefunctions .    although @xmath2he is spinless and most bosonic alkali atoms are too heavy to exhibit the relativistic so coupling in their center - of - mass motion , so coupling can be important in exciton systems in two dimensional quantum wells .
we have shown that the rashba so coupling in the conduction electron bands induces the same type of so coupling in the center - of - mass motion of excitons . in a harmonic trap
, the ground state condensate wavefunction can spontaneously develop half - quantum vortex structure and the skyrmion type of the spin texture configuration @xcite . on the other hand , effective so coupling in boson systems can be induced by laser beams , which has been investigated in several publications by other groups in refs @xcite .    in the following
, we will review our work in both directions outlined above including many new results unpublished before .
in section [ sect : boson ] , we explain the characteristic feature of interacting bosons in high orbital bands , the `` orbital hund s rule '' , and the consequential complex - superfluid states with the ordering of on - site oam moments .
the ordering of oam moments in the mott - insulating states is also investigated . in section [
sect : spinorbit ] , we review the tr symmetry breaking states of spinful bosons with spin - orbit coupling .
interesting properties including half - quantum vortex , and the skyrmion - like spin - textures are studied .
conclusions are made in section [ sect : conclusion ] .
due to the limit of space , we will not cover interesting related topics of orbital bosons , such as the nematic superfluid state @xcite and the algebraic superfluid state @xcite , and the orbital physics with cold fermions @xcite , which has also aroused much research attention recently .
in this section , we will review the `` complex - condensation '' with tr symmetry breaking , which is a new state of the @xmath3-orbital bosons .
this is an example of novel orbital physics in optical lattice with cold bosons .
below let us give a brief general introduction to orbital physics for the general audience .
orbital is a degree of freedom independent of charge and spin .
it plays important roles in magnetism , superconductivity , and transport in transition metal oxides @xcite .
the key features of orbital physics are orbital degeneracy and spatial anisotropy .
optical lattices bring new features to orbital physics which are not easily accessible in solid state orbital systems .
first , optical lattices are rigid lattices and free from the jahn - teller distortion , thus orbital degeneracy is robust .
second , the meta - stable bosons pumped into high orbital bands exhibit novel superfluidity beyond feynman s `` no - node '' theory @xcite .
third , @xmath3-orbitals have stronger spatial anisotropy than that of @xmath4 and @xmath5-orbitals , while correlation effects in @xmath3-orbital solid state systems ( e.g. semiconductors ) are not that strong .
in contrast , interaction strength in optical lattices is tunable .
we can integrate strong correlation with spatial anisotropy more closely than ever in @xmath3-orbital optical lattice systems @xcite .
recently , orbital physics with cold atoms has been attracting a great deal of attention@xcite . for orbital bosons , a series of theory works
have been done @xcite , including the illustration of the ferro - orbital nature of interactions @xcite , orbital superfluidity with spontaneous time - reversal symmetry breaking @xcite , the nematic superfluidity@xcite .
the theory work on the @xmath3-orbital fermions is also exciting including the flat band and associated strong correlation physics in the honeycomb lattice @xcite , orbital exchange and frustration @xcite , and topological insulators @xcite .
the experimental progress has been truly exciting @xcite .
_ @xcite have realized the meta - stable @xmath3-orbital boson systems by using the stimulated raman transition .
the spatially anisotropic phase coherence pattern has been observed in the time of flight experiments .
sebby - strabley _
_ @xcite have successfully pumped bosons into excited bands in the double - well lattice .
in addition , @xmath3-orbital bose - einstein condensation ( bec ) has also been observed in quasi one - dimensional exciton - polariton lattice systems @xcite .    below we illustrate the important feature of interactions between orbital bosons , the `` orbital hund s rule '' , which results in tr symmetry breaking .
the most remarkable feature of interacting bosons in high orbital bands is that they favor to maximize their onsite orbital angular momentum ( oam ) as explained below @xcite .
let us illustrate this hund s rule type physics through the simplest example : a single site problem with two degenerate @xmath6-orbitals filled with two spinless bosons . because bosons are indistinguishable , the hilbert space for the two - body states only contains three states , which can be classified according to their oam : an oam singlet as @xmath7 depicted in fig [ fig : pboson ] .
a , and a pair of oam doublets @xmath8 , the interaction energy of the former is calculated as @xmath9 while that of the latter is @xmath10 with the definition of u = gd r^2 |_p_x ( r)|^4= gd r^2 |_p_y ( r)|^4 . in the oam singlet state , bosons occupy polar ( real ) orbitals ( e.g. @xmath11 ) , whose angular distribution in real space is narrower than that of the axial ( complex ) orbitals ( e.g. @xmath12 ) for the oam doublets . by occupying the same axial orbital and therefore maximizing oam , two bosons enjoy more room to avoid each other .
this `` ferro - orbital '' interaction is captured by the following multi - band hubbard hamiltonian for the @xmath3-orbital bosons as h_int=_r \{n^2_r- l_z^2 } , [ eq : bshamint ] where @xmath13 and @xmath14 .
the first term is the ordinary hubbard interaction and the second term arises because of the orbital degree of freedom . in three dimensional systems in which all of the three @xmath3-orbitals are present
, we only need to replace @xmath15 with @xmath16 and @xmath17 defined as @xmath18 , @xmath19 .
when more than two bosons occupy a single site , bosons prefer to go to the same single particle state by their statistical properties . again going into the same axial state and thus maximizing oam
can reduce their repulsive interaction energy .
this is an analogy to the hund s rule of electron s filling in atomic shells .
the first hund s rule of electrons maximizes electron total spin to antisymmetrize their wavefunction , and the second hund s rule further maximize their oam to extend the spatial volume of wavefunction .
the key feature is that electrons want to avoid each other as possible as they can . for the orbital physics of spinless bosons
, the same spirit applies with the maximizing of oam .
now we have understood the single site physics in which bosons develop rotation . from the symmetry point of view
, it is similar to the @xmath20 superconductors .
the @xmath3-orbital hamiltonian within the tight - binding approximation can be written as h_0= t__i j [ p^_i , e_ij p_j , e_ij+h.c.]-t__i j [ p^_i , f_ij p_j , f_ij + h.c .
] , [ eq : bandsquare ] where the unit vector @xmath21 is along the bond orientation between two neighboring sites @xmath22 and @xmath23 and @xmath24 is perpendicular to @xmath21 .
@xmath25 and @xmath26 are the projections of @xmath3-orbitals along ( perpendicular to ) the bond direction respectively as defined below p_e_ij= ( p_x e_x + p_y e_y ) e_ij ,  p_f_ij= ( p_x e_x + p_y e_y ) f_ij .
the @xmath27-bonding @xmath28 and the @xmath29-bonding @xmath30 describe the hoppings along and perpendicular to the bond direction as depicted in fig .
[ fig : pboson ] c and d , respectively .
the opposite signs of the @xmath27 and @xmath29-bondings are due to the odd parity of @xmath3-orbitals .
@xmath30 is usually much smaller than @xmath28 in strong periodical potentials .
[ eq : bandsquare ] exhibits degenerate band minima in both the square and triangular lattice@xcite . in the square lattice , the brillouin zone ( bz )
is a square with the edge length of @xmath31 where @xmath32 is the length between the nearest neighbors .
the spectra read @xmath33 and @xmath34 .
as depicted in fig .
[ fig : band ] .
a , the band minima are located at @xmath35 for the @xmath36-band and @xmath37 for the @xmath38-band , respectively .    in the triangular lattice
, we set the @xmath29-bonding to zero which does not change the qualitative band structures .
we take unit vectors from one site to its six neighbors as @xmath39 with @xmath40 , @xmath41 .
the brillouin zone takes the shape of a regular hexagon with the edge length @xmath42 .
the energy spectrum of @xmath43 is e(k)= t_\ { f_k } , with f_k= _
i=1 ^ 3 ( k e_i ) ,  g_k = _
3i > j1 ( k e_i ) ( k e_j ) .
the spectrum contains three degenerate minima located at the three non - equivalent middle points of the edges as @xmath44 , the factor @xmath45 takes the value of @xmath46 uniformly in each horizontal row but alternating in adjacent rows .
if the above pattern is rotated at angles of @xmath47 , then we arrive at the patterns of @xmath48 .
each eigenvector is a 2-component superposition vector of @xmath36 and @xmath38 orbitals .
the eigenvectors at energy minima are @xmath49 .
@xmath50 can be obtained by rotating @xmath51 at angles of @xmath52 respectively as depicted in fig .
[ fig : band ] b.      if there were no interactions , any linear superposition of the above band minima would be a valid condensate wavefunctions .
however , interactions will select a particular type combination which exhibits tr symmetry breaking and agrees with the above picture of `` orbital hund s rule '' .      in the square lattice
, we take the condensate wavefunction as @xmath53 under the constraint of @xmath54 .
any choice of @xmath55 minimizes the kinetic energy .
however , the interaction @xmath56-terms break the degeneracy .
the degenerate perturbation theory shows that the ground state values of @xmath55 take @xmath57 or its equivalent tr partner of @xmath58 .
the mean field condensate can be described as \ { ( ^_k_p_x + i ^_k_p_y)}^n_0 |0 , [ eq : square ] where @xmath59 is the particle number in the condensate .    for better insight
, we transform the above momentum space condensate to the real space .
the orbital configuration on each site reads e^i_r(|p_x+i_r |p_y ) , where the @xmath60 phase @xmath61 is specified at the right lobe of the @xmath3-orbital .
the ising variable @xmath62 denotes the direction of the oam , and is represented by the anti - clockwise ( clockwise ) arrow on each site in fig .
[ fig : square ] b. each site exhibits a nonzero oam moment and breaks tr symmetry .
the condensate wavefunction of eq .
[ eq : square ] describes the staggered ordering of @xmath63 .
we check that the phase difference is zero along each bond , and thus no inter - site bond current exists .
the condensate described in eq .
[ eq : square ] breaks both tr reversal and translational symmetries , and thus corresponds to a bec at non - zero momenta .
this feature should exhibit itself in the time - of - flight experiments which have been widely used to probe the momentum distribution of cold atoms .
the coherence peaks of eq .
[ eq : square ] are not located at integer values of the reciprocal lattice vectors but at @xmath64 , and @xmath65 as depicted in fig .
[ fig : square ] a. furthermore unlike the @xmath66-orbital condensate the @xmath3-wave wannier function superposes a _ non - trivial profile _ on the height of density peaks . as a result
, the highest peaks are shifted from the origin  a standard for the @xmath66-wave peak  to the reciprocal lattice vectors whose magnitude is around @xmath67 where @xmath68 is the characteristic length scale of the harmonic potential of each optical site .
a detailed calculation of form factors is given in ref .
@xcite .
the elementary excitations in @xmath3-orbital becs consist of both the gapless phonon mode and the gapped orbital mode .
the former corresponds to the goldstone mode related to the @xmath60 symmetry breaking , and the latter corresponds to the flipping of the direction of oam moments . in the following
, we will take the 2d square lattice as an example .
we assume the condensate as @xmath69 $ ] with @xmath70 and @xmath71 .
the boson operators take the expectation values as = ( -1)^r_x ,  = i(-1)^r_y . from minimizing the onsite part of the free energy respect to the condensate order parameter @xmath72
f =- n+ ( n^2- l_z^2 ) , we have @xmath73 respect to the band minima . the fluctuation around the expectation value is defined as p_x(r)=+p_x ,  p_y(r)=+p_y
. then the interaction terms in eq .
[ eq : bshamint ] are expanded as h_int&= & ||^4 + [ p^_x p_x + p^_y p_y ] ||^2 + ^*,2 [ p_x p_x -p_y p_y ] + & + & ^2 [ p_x^p_x^-p_y^p_y^ ] + ( -)^r_x+r_y [ ^2
p_x^p_y^- ^*,2 p_x p_y ] . combined with the free part
, we arrive the mean field hamiltonian as h_mf=_k , a , b^^(k)_a m_ab(k ) _
b(k ) where @xmath74 $ ] , @xmath75 , and the summation is only over half of the brillouin zone .
the matrix kernel reads m(k)= (    cccc _ x(k)+||^2 & 0 & ^2 & i ^2 + 0 & _ y ( k + q)+ ||^2 & i ^2 & ^2 + u ^*,2 & - u ^*,2 & _ x(k)+||^2&0 + - u ^*,2 & u ^*,2 & 0 & _ y ( k + q)+ ||^2    ) .
[ eq : kernel ] the spectra are reduced to generalized eigenvalue problem of x^-1 \{(1,1,-1,-1 )  m ( k ) } x = ( e_1,e_2,-e_3,-e_4 ) , where @xmath76 are excitation eigenvalues , and @xmath77 contains the eigenvectors .    let us consider the excitation close to the condensation wavevector @xmath78 , and then @xmath79 . at small value of @xmath80 , we obtain the excitation spectra as e_1,3(k)&= & ,   e_2,4(k)= |(q ) + u ||^2 .
\approx \frac{t_\pp + t_\perp}{2 } ( q_x^2 + q_y^2)$ ] .
clearly , the gapless mode with the linear dispersion relation describes the superfluid phase fluctuations ; the gapped mode describes the orbital excitations corresponding to the flipping of orbital angular momenta .      in the triangular lattice ,
the oam moments instead form a stripe ordering , i.e. , the oam moments along one row polarizes along the @xmath82-axis and those along the neighouring rows polarizes with the opposite direction .
this can be intuitively understood as follows . in the superfluid state ,
the oam moments behave like vortices whose interactions are long range .
the above stripe configuration of positive and negative vortices is the optimal configuration to minimize the globe vorticity .
let us first examine the weak coupling limit .
again we write a general form for condensation wavefunction as a linear superposition of the three band minima _ c(r)= ( c_1 _ k_1 + c_2 _ k_2 + c_3 _
k_3 ) , where @xmath83 are the locations of band minima defined in sect .
[ subsect : band ] . without loss of generality , we set @xmath84 , and define @xmath85 , @xmath86 and @xmath87 and @xmath88 , then the normalization factor @xmath89 .
the interaction energy persite is calculated as e_int&= & \ { ( x_1 ^ 2+x_2 ^ 2-y_1 ^ 2-y_2 ^ 2)^2 + 8 x_1 ^ 2 + 8 y_1 ^ 2 } + 1 + & - & \ { 4 ( x_1 y_2 -x_2 y_1)^2 + 2x_2 ^ 2 + 2 y_2 ^ 2}. the terms in the first line are from the density - density interaction which can be minimized by setting @xmath90 and @xmath91 .
this means only one of @xmath92 is nonzero and purely imaginary . in this case
, the particle number on each site is uniform .
the terms in the second line can also be minimized at this condition with a further requirement of @xmath93 . if we take @xmath94 nonzero , then @xmath95 .
thus the mean field condensate can be expressed as @xmath96 with @xmath97 the vacuum state and @xmath59 the particle number in the condensate .
this state breaks the @xmath60 gauge symmetry , as well as tr and lattice rotation symmetries , thus the ground state manifold is @xmath98 .
this state also breaks lattice translation symmetry , which is , however , equivalent to suitable combinations of @xmath60 and lattice rotation operations .
again we transform the above momentum space condensate to the real space , whose orbital configuration takes e^i_r(|p_x+i_r |p_y ) [ eq : triangular ] with @xmath99 as @xmath100 .
the general configuration of @xmath101 is depicted in fig .
[ fig : triangular ] b for later convenience . at @xmath100 , @xmath6
are not equally populated , and the moment per particle is @xmath102 .
this does not fully optimize @xmath103 which requires @xmath104 .
however , it fully optimizes @xmath43 which dominates over @xmath103 in the weak coupling limit .
we check that the phase difference is zero along each bond , and thus no inter - site bond current exists .
interestingly , as depicted in fig .
[ fig : triangular ] b , oam moments form a stripe order along each horizontal row .
this stripe ordering in the weak coupling limit is robust at small values of the @xmath29-bonding @xmath30 because it does not change the location of the band minima and the corresponding eigenfunctions of @xmath105 at all .
the driving force for this stripe formation in the sf regime is the kinetic energy , i.e. , the phase coherence between bosons in each site .
by contrast , the stripe formation in high t@xmath106 cuprates is driven by the competition between long range repulsion and the short range attraction in the interaction terms @xcite .
this stripe phase should manifest itself in the time of flight ( tof ) signal as depicted in fig .
[ fig : triangular ] a. in the superfluid state , we assume the stripe ordering wavevector @xmath107 , and the corresponding condensation wavevectors at @xmath108 . as a result , the tof density peak position after a fight time of @xmath109 is shifted from the reciprocal lattice vectors @xmath110 as follows _
t & & _ g\ { |_2 ( , k)|^2 ^2(k - k_2- g ) + |_3 ( , k)|^2 ^2(k - k_3-g ) } , where @xmath111 ; @xmath112 is the fourier transform of the wannier @xmath3-orbital wavefunction @xmath113 , and @xmath114 $ ] with @xmath115 integers .
thus bragg peaks should occur at @xmath116 $ ] . due to the form factors of the @xmath3-wave wannier orbit wavefunction @xmath117 , the locations of the highest peaks
is not located at the origin but around @xmath118 . due to the breaking of lattice rotation symmetry
, the pattern of bragg peaks can be rotated at angles of @xmath52 .      in this subsection ,
the ordering of the oam moments in the strong coupling superfluid regime is examined .
we will employ the lattice gauge theory formalism developed by moore and lee in ref .
@xcite in the context of the @xmath119 josephson junction array systems .    in the superfluid regime
, each site @xmath22 is denoted by a @xmath60 variable @xmath120 which denotes the superfluid phase along the @xmath121-direction , and by an ising variable @xmath122 for the direction of the oam moment as depicted in fig .
[ fig : strong ] . along a general direction of the azimuthal angle @xmath123
, the superfluid phase is @xmath124 .
the intersite @xmath27 and @xmath29-bonding become the intersite josephson coupling as h_ij =- nt_(_i-_j - a_(i , j ) ) -nt_(_i-_j - a _ ( i , j ) ) , [ eq : josephson ] where @xmath125 is the average particle number per site .
the phase differences in eq .
[ eq : josephson ] takes into account the geometric orientation of the bond @xmath126 as captured by the gauge fields @xmath127 and @xmath128 with the definition that @xmath129 , @xmath130 , and @xmath131 are the azimuth angles relative to the @xmath121-axis defined in fig .
[ fig : strong ] .
let us first consider the case of @xmath132 and write down an effective hamiltonian for the ising variables @xmath133 .
in order to minimize both the josephson couplings of the @xmath27 and @xmath29-bonds , we need @xmath134 and thus @xmath135 . otherwise , if @xmath136 , then @xmath137 .
this costs an energy of @xmath138 , and leads to an effective antiferro - orbital interaction between the ising variables as h_eff = n t__ij _ z , i _ z , j .
thus the antiferro - orbital ordering of oam moments in the square lattice remains valid in the strong coupling limit at non - vanishing @xmath30 , which enforces that both the @xmath27 and @xmath29-bonds achieve phase coherence as depicted in fig .
[ fig : square ] .
b.    on the other hand , in the limit of a vanishing @xmath30 , the leading order effect involves multiple site interaction around a plaquette .
we follow the method described in ref .
@xcite , and perform the duality transformation for the @xmath60 phase variables @xmath72 under the background of the geometric gauge potential @xmath127 .
we separate the contributions from the phonon part and the vortex part as z&=&z_ph _
m^\{- _ x ( m(x ) - _ x)^2 + & + & nt__xx^ ( m(x ) -_x ) ( m ( x^ )
-_x^ ) } , [ eq : vortex ] where @xmath139 is the phonon contribution in the vortex - free configuration ; @xmath121 marks the dual lattice site ( or the plaquette index in the original lattice ) ; @xmath140 is the vortex charge ; @xmath141 is the external flux through the plaquette defined as _
x=_ij a _ , ij .    in the square lattice , the geometric gauge flux around each plaquette @xmath121
is calculated as _
x=- ( _ i + _ i+e_x+_i+e_x+e_y + _ i+e_y ) .
if @xmath141 is integer - valued , it can be absorbed by shifting the zero of the vortex charge as @xmath142 which remains integer - valued , thus there is no cost for energy .
the feature can be captured by the effective hamiltonian in ref .
@xcite as h_eff=-
k _ ijkl _ i _ i _ k _
l , [ eq : z2 ] where @xmath143 are four sites around a plaquette centered at @xmath121 ; @xmath144 is the energy scale of the @xmath27-bonding .
this model has a sub - extensive @xmath145 symmetry investigated in ref .
@xcite , i.e. , flipping the sign of @xmath133 along each row or column leave eq .
[ eq : z2 ] invariant .
thus it can not develop the ordering of oam at any finite temperature .
let us come back to the superfluid sector @xmath72 , its most relevant topological defect is the half - quantum vortices because @xmath141 can take the values of @xmath146 .
the kosterlitz - thouless transition is associated with the unbinding of half - quantum vortices . as a result
, the low temperature phase is the quasi - long range ordering of pairing of bosons @xcite .
next we move to the strong coupling theory in the triangular lattice .
we start from the limit of @xmath147 .
the geometrical gauge flux becomes _ x=- ( _ 1 + _
2+_3 ) , [ eq : triagflux ] where @xmath148 are three sites around a triangular plaquette . there is no way to form an integer flux .
the smallest vorticity per plaquette is @xmath149 which corresponds to either two @xmath150 s and one @xmath151 , or two @xmath151 s and one @xmath150 to minimize the vortex core energy .
in such a dense vortex configuration , eq . [ eq : vortex ] rigorously speaking does not apply because its validity relies on the assumption of small vortex fugacity .
however , the structure of interactions among vortices still implies that vortices form a regular lattice with alternating positive and negative vorticity . the dual lattice ( the center of the triangular plaquette )
is the bipartite honeycomb lattice .
it is tempting to assign @xmath149 alternatively to each plaquette , but actually it is not possible due to the following reason .
consider a plaquette with vorticity @xmath152 , thus its three vertices are with two @xmath150 s and one @xmath151 .
the neighboring plaquette sharing the edge with two @xmath150 s must have the same vorticity , and merges with the former one to form a rhombic plaquette with vorticity @xmath153 as depicted in fig .
[ fig : strong ] b. thus the ground state should exhibit a staggered pattern of rhombic plaquette with vorticity of @xmath154 .
this arrangement precisely corresponds to the stripe order of the ising variables .
this configuration breaks both the lattice rotational and translational symmetries , which is six - fold degenerate .
if we turn on the small @xmath30 term , it results in an antiferromagnetic ising coupling with nearest bonds , the stripe configuration also satisfies its ground state requirement . thus we believe that @xmath30 term does not change the stripe configuration in the strong coupling limit either .
having established the stripe order for the oam moments , it is straightforward to further optimize the @xmath60 phase variable @xmath72 . the result is depicted in fig .
[ fig : triangular ] b with @xmath72 marked on the right lobe of the @xmath155 orbit .
each horizontal bond has perfect phase match , while each of tilted bond has a phase difference of @xmath156 .
thus around each rhombic plaquette , the phase winding is @xmath157 and gives rise to the josephson supercurrent along the bonds as j= , [ eq : pqcurrent ] with @xmath158 and the directions specified by arrows in fig .
[ fig : triangular ] b. in other words , in addition to the stripe order of the ising variables which corresponds to the onsite oam moments , there exists a staggered plaquette bond current order
. this feature does not appear in the square lattice @xcite .
the orbital plaquette current ordering has been studied in the strongly correlated electron systems , such as the circulating orbital current phase @xcite and the d - density wave states in the high t@xmath106 compounds @xcite .
it is amazing that in spite of very different microscopic mechanism and energy scales , the two completely different systems exhibit a similar phenomenology .
we have shown that in the triangular lattice the stripe ordering of the oam moments exist in both weak and strong coupling superfluid states .
the orbital configurations are slightly different : the orbital mixing angle @xmath101 defined in eq .
[ eq : triangular ] equals to @xmath159 in the weak coupling limit , while it equals to @xmath160 in the strong coupling limit .
this arises from the competition between the kinetic energy and the onsite hund s rule . below
we will see that as @xmath161 goes from small to large , the stripe ordering remains the same with a smooth evolution of @xmath101 from @xmath159 to @xmath160 as depicted in fig .
[ fig : triangular ] b.    we use a gutzwiller type mean field theory for a @xmath162 lattice under the periodic boundary condition .
we explicitly checked for three sets of parameters @xmath163 marked as points of @xmath164 and @xmath165 in fig .
[ fig : triang - phase ] b. the stripe - ordered ground state with a @xmath166 unit cell depicted in fig .
[ fig : triangular ] a is found stable against small random perturbations in all the three cases .
then we further apply the gutzwiller type theory assuming a @xmath166 unit cell and obtain the phase diagram depicted in fig .
[ fig : triang - phase ] a which includes both the stripe ordered superfluid phase and mott - insulating phases .
so far , we have only discussed the orbital ordering in the superfluid states .
since the oam moment is a different degree of freedom from the superfluid phase , we expect that its ordering can survive even inside the mott - insulating phases . in this subsection
, we continue to study the exchange physics of orbital bosons and the related ordering of oam moments in the absence of the superfluidity order . for simplicity
, we will use the triangular lattice as an example .
we consider the mott - insulating phases with @xmath125 spinless bosons per - site and two degenerate orbitals of @xmath178 and @xmath38 .
we define the tr doublets of all the particles in the states of @xmath179 as the eigenstates of the ising operator @xmath133 with the eigenvalues of @xmath46 .
the ising part of the effective exchange hamilton occurs at the level of the second order perturbation theory , while the ising variable flipping process occurs at the @xmath180-th order perturbation theory .
we consider the large-@xmath125 case in which the system is deeply inside the ising anisotropy class , and no orbital - flip process occurs at the leading order . in the following , we will study the physics of the @xmath181-site exchange processes .
let us consider the virtual hopping processes in the mott - insulating states along the bond depicted in fig .
[ fig : strong ] a. both the @xmath27-bonding @xmath182 and @xmath29-bonding @xmath30 are kept . with the definition of @xmath183
we can express the hopping as h_t&=&-_^ \
{ t_p^_(j ) p_^ ( i ) e^-i(_j-^_i ) + t_p^_(j ) p_^ ( i ) e^-i(-^ ) ( _ i + ) + h.c . }
+ & = & _ \{p^_j , p_i+h.c . }
+ _ \{p^_j , p_i , - e^-2 i _ i+h.c . }
, [ eq : motthop ] where the definition of angles @xmath184 follows the convention depicted in fig .
[ fig : strong ] a.    we calculate the energy shifts in both the ferro - orbital configurations ( @xmath136 ) and the antiferro - orbital configurations ( @xmath185 ) within the second order perturbation theory .
the energy difference between these two configurations is e_fo&= & - - + e_afo&= & - - , where @xmath186 and @xmath187 .
the antiferro - orbital configuration has lower energy , which arises from the fact that the orbital - flip process has a larger amplitude than that of the orbital non - flip process described in eq .
[ eq : motthop ] .
this is because that the @xmath27-bonding and @xmath29-bonding amplitudes have a @xmath29-phase shift .
the ising coupling @xmath188 reads h=
j_afo _ ij _ z ( i ) _ z ( j ) [ eq:2ptcl ] where @xmath189 with only the leading order contribution kept .
if @xmath30 is set to zero , @xmath188 varnishes .
this is clear from the fact that we can flip the sign of the @xmath3-orbit component perpendicular to the bond direction .
this operation changes the value of @xmath133 , but has no effect on the energy . in this case , we need to further study the multi - site virtual hopping processes .
the @xmath27-bonding term by itself gives non - zero ring exchange terms for the multi - site processes , thus we neglect the contribution from the @xmath29-bonding part . in the following ,
we only keep the leading order virtual process proportional to @xmath190 .
we consider a triangular plaquette with three sites @xmath191 each of which is denoted by the particle number and the ising variables as @xmath192 , respectively .
there are 12 different ring - hopping processes whose contribution is at the order of @xmath190 .
we enumerate 4 of them explicitly as & & ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) ( n-1,_z , i ) ( n+1,_z , j ) ( n,_z , k ) + & & ( n-1,_z , i ) ( n,_z , j ) ( n+1,_z , k ) ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) + & & ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) ( n+1,_z , i ) ( n-1,_z , j ) ( n,_z , k ) + & & ( n,_z , i ) ( n-1,_z , j ) ( n+1,_z , k ) ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) + & & ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) ( n-1,_z , i ) ( n+1,_z , j ) ( n,_z , k ) + & & ( n,_z , i ) ( n+1,_z , j ) ( n-1,_z , k ) ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) + & & ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) ( n+1,_z , i ) ( n-1,_z , j ) ( n,_z , k ) + & & ( n-1,_z , i ) ( n,_z , j ) ( n+1,_z , k ) ( n,_z , i ) ( n,_z , j ) ( n,_z , k ) .
the other @xmath193 processes can be obtained by a cyclic permutation @xmath194 .
after sum over all these processes , the corresponding energy shift reads e ( _ i , _ j , _ k)6 ( ) ^3 ( n ( n+1))^3/2 & & = -j_3 ( 2_3 ) , [ eq:3ptcl ] where @xmath195 and @xmath196 ; @xmath197 is defined as before in the superfluid case . again for each plaquette , @xmath197 takes the value of @xmath198 which correspond to @xmath133s taking the values of two @xmath199s and one @xmath151 or two @xmath151s and one @xmath199 .
. [ eq:3ptcl ] plays a similar role to vortex core energy in eq .
[ eq : vortex ] .
[ eq:3ptcl ] by itself has a huge ground state degeneracy because the requirement is the same as that of the antiferromagnetic ising model .
we need to add the four - site process to lift the degeneracy .
the four - site process mimics the interaction between two adjacent vortices . by the spirit of perturbation theory ,
the 4-particle process is further suppressed by a factor @xmath200 , thus the interaction between vortices should be of short range .
this makes sense in the mott insulating state due to the condensation of vortices and the resulting screening effect .
again we only consider the effect from the @xmath27-bonding term , and only keep the leading order contribution at the level of @xmath201 .
similarly , we have 48 different virtual hopping processes . for a plaquette with vertices @xmath202 ,
the energy shift is & & e(_i,_j,_k,_l ) -j_4 ( 2_4 ) , [ eq:4ptcl ] where @xmath203 . in the square lattice ,
@xmath204 . if @xmath147 , eq .
[ eq:4ptcl ] is the leading order contribution , thus the effective hamiltonian looks the same as that in the superfluid case but with a reduced coupling constant . if @xmath132 , then the two - site exchange of eq .
[ eq:2ptcl ] results in the long range antiferro - orbital order .    in the triangular lattice
, the @xmath205 of the four - site plaquette @xmath202 just equals the sum of the flux of the two triangular plaquettes .
since each triangular plaquette takes the flux value of @xmath149 , the rhombic four - site plaquette can only take the flux value of @xmath206 , @xmath172 .
the stripe ordering pattern gives the maximum number of the zero - flux four - site plaquettes , and thus is selected as the ground state .
the stripe ordering appears at temperatures @xmath207 and disappears at @xmath208 in which the three - site exchange process dominates the physics .
in this section , we discuss another route for `` complex - condensation '' of bosons , i.e. , spinful bosons with spin - orbit ( so ) coupling .
the proof of feynman s `` no - node '' theorem crucially relies on the fact that the kinetic energy of bosons only contains even power of momentum , thus so coupling which linearly depends on momentum invalidates feynman s proof . in the following , we consider the simplest example of the two - component bosons with the rashba - like so coupling and investigate their condensate in the inhomogeneous harmonic traps .
a related work on the bec of two - component bosons has also been studied by stanescu _
@xcite .
we consider bosons with two - components denoted as pseudospin up and down .
the many - body hamiltonian for interacting bosons is written as h&=&d^2 r   ^_\ { -- + v_ext(r ) } _ + & + & _ r ^ _ ( -i _ y _ x + i _ x _ y)_+^ _ ^ _ _ _ ,   [ eq : mnbdy ] where @xmath209 is the boson operator ; @xmath101 refers to boson pseudospin @xmath210 and @xmath211 ; @xmath212 is the external potential ; @xmath213 describes the @xmath66-wave scattering interaction which is assumed to be spin - independent for simplicity ; @xmath214 is the rashba so coupling strength .
although eq .
[ eq : mnbdy ] is of bosons , it satisfies a suitably defined tr symmetry like fermion systems as @xmath215 with @xmath216 where @xmath217 is the complex conjugate and @xmath218 operates on the pseudospin doublet . in the homogeneous system ,
the single particle states are the helicity eigenstates of @xmath219 with the dispersion relations of @xmath220 where @xmath221 .
the energy minima are located at the lower branch along a ring with radius @xmath222 .
the corresponding two - component wavefunction @xmath223 with @xmath224 can be solved as @xmath225 , where @xmath226 is the azimuth angle of @xmath227 .    instead of discussing the condensation in the homogeneous system in which frustrations occur due to the degeneracy of the single particle ground states , we consider a 2d system with an external harmonic trap with @xmath228 .
we define the characteristic so energy scale for the harmonic trap as @xmath229 where @xmath230 is the length scale of the trap , and correspondingly a dimensionless parameter as @xmath231 .
the characteristic interaction energy scale is defined as @xmath232 where @xmath59 is the total particle number and the dimensionless parameter @xmath233 .
the gross - pitaevskii equation can be obtained as the saddle - point equation of eq .
[ eq : mnbdy ] as \ { - & + & _ r ( -i _ y _ x , + i _ x _ y , ) + g ( ^ * _ _ ) + m ^2_t r^2 } _ ( r , ) + & = & e _ ( r , ) .
[ eq : gp ] due to the 2d rotational symmetry , the ground state condensate wavefunctions can be denoted by the total angular momentum @xmath234 , which can be represented in the polar coordinate as     c f(r ) + g(r ) e^i    ) ,  |_-= (    c -g(r ) e^i + f(r )    ) , where @xmath235 and @xmath72 are the radius and the azimuthal angle , respectively .
both @xmath236 and @xmath237 are real functions .
[ eq : gp ] is numerically solved in the harmonic trap with parameters of @xmath238 and @xmath239 .
we choose the condensate as one of the tr doublet @xmath240 , and present its radial density profiles of both spin components @xmath241 and @xmath242 at @xmath238 in fig .
[ fig : trap ] a. further in this strong so coupling case with @xmath243 , the condensate wavefunction has nearly equal weight in the spin up and down components , i.e. @xmath244 , thus the average spin moment along the @xmath82-axis equals to zero .
the total angular momentum per particle @xmath245 is mainly from the orbital angular momentum polarization , i.e. , one spin component stays in the @xmath66-state and the other one in the @xmath246-state .
this is an example of half - quantum vortex configuration @xcite , thus spontanously breaking tr symmetry .
clearly this is a `` complex - valued '' ground state wavefunction beyond the `` no - node '' theorem .
this condensate wavefunction exhibits interesting spin density distributions in real space as skyrmion - like spin textures .
the radial wavefunction in pseudo - spin up and down components @xmath236 and @xmath237 exhibit oscillations with an approximate wavevector of @xmath247 , which originates from the ring structure of the low energy states in momentum space , thus are analogous to friedel oscillations in fermion systems .
the pseudo - spin up component is @xmath66-wave like , thus @xmath236 reaches the maximum at @xmath248 ; while the down component is of the @xmath3-wave , thus @xmath249 at @xmath248 . in other words , approximately there is a relative phase shift of @xmath250 between the oscillations of @xmath236 and @xmath237 . in fig .
[ fig : trap ] a , @xmath241 and @xmath242 are plotted .
the spin density distribution @xmath251 can be expanded as s_z ( r , ) & = & ( |f(r)|^2 -|g(r)|^2 ) ,  s_x ( r , ) = f(r ) g(r ) , + s_y ( r , ) & = & f(r ) g(r ) . along the @xmath121-axis
, the spin density lies in the @xmath82-@xmath121 plane as depicted in fig .
[ fig : trap ] b. because the circulating supercurrent is along the tangential direction , the spin density distribution is along the radial direction and exhibits an interesting topological texture configuration which spirals in the @xmath82-@xmath121 plane at the pitch value of the density oscillations . the distribution in the whole space can be obtained through a rotation around the @xmath82-axis .
this spin texture configuration is of the skyrmion - like .
next we discuss the possible realization of the above exotic bec in exciton systems .
excitons are composite objects between conduction electrons and valence holes @xcite .
in particular , the recently progress on the indirect exciton systems greatly enhances the lift - time @xcite , which provides a wonderful opportunity to investigate the exotic state of matter of the exciton condensation @xcite .
the ordinary bosons are too heavy to exhibit the relativistic so coupling in their center - of - motion . due to the small effective mass of excitons ,
so coupling can result in important consequences , including anisotropic electron - hole pairing @xcite , spin hall effect of the center - of - mass motion of excitons @xcite , and the berry phase effect on exciton condensation @xcite .
we will show that the rashba so coupling in the electron band also survives in the center - of - mass motion of excitions .
we begin with the hamiltonian of indirect excitons h_e&=&- ( _ e , x^2 + _ e , y^2 ) + i _
r , el ( _ e , x _ y - _ e , y _ x ) , [ eq : hamrashba ] + h_hh&=&- ( _ h , x^2 + _ h , y^2 ) , [ eq : hamhh ] + h_e - hh&= & - , where @xmath252 is the effective mass of conduction electrons ; @xmath253 is the rashba so coupling strength of the conduction electron , @xmath254 is the dielectric constant ; @xmath4 is the thickness of the barrier .    for small exciton concentrations ,
we only need to consider the heavy hole ( @xmath255 ) band with the effective mass @xmath256 and @xmath257 , which is separated from the light hole band with a gap of the order of 10 mev .
we consider the center - of - mass motion in the bec limit of excitons , which can be separated from the relative motion in @xmath258 and @xmath259 . similarly to ref .
@xcite , the effective hamiltonian of the 4-component @xmath255 excitons denoted as @xmath260 can be represented by the matrix form as [ eq : hamexciton ] where @xmath227 is the center - of - mass momentum ; @xmath261 is the total mass of the exciton ; @xmath262 ; @xmath263 ; @xmath264 is the exchange integral and @xmath265 has the @xmath4-wave structure as @xmath266 , both of which are exponentially suppressed by the tunneling barrier for indirect excitons , and will be neglected below .
consequentially , @xmath267 becomes block - diagonalized .
we consider using circularly polarized light to pump the exciton of @xmath268 , and then focus on the left - up @xmath269 block of the @xmath267 matrix with heavy hole spin @xmath270 .
if we use the electron spin number @xmath210 and @xmath211 as the exciton component , we will arrive the eq .
[ eq : mnbdy ] with the renormalized so coupling strength @xmath271 .
we next justify the above choice of the values of parameters based on experimental situations .
the effective masses of electrons and holes in gaas / algaas quantum wells are @xmath272 and @xmath273 in ref .
so coupling strength @xmath274 can reach @xmath275 ev.m in ref .
thus we can estimate a reasonable value of @xmath276 @xmath277 . for a harmonic trap with @xmath278 , @xmath279 and @xmath280 mk . in two dimensional harmonic traps , the critical condensation temperature @xmath281 @xcite .
if we take the exciton density @xmath282 @xmath283 and the effective area @xmath284 , we arrive at @xmath285 mk , which is an experimentally available temperature scale @xcite .
the average interaction energy per exciton in the typical density regime of @xmath286@xmath277 is estimated around 2mev in ref .
@xcite , thus we take the interaction parameter @xmath239 in the calculation above .
the spatial periodicity of the spin texture is about @xmath287 and , thus , is detectable by using optical methods .
we have reviewed the exotic condensations of bosons whose many - body wavefunctions are complex - valued in the coordinate representation , thus they go beyond the well - known paradigm of the `` _ no - node _ '' theorem .
we studied two possible ways to bypass this theorem to achieve unconventional condensations , i.e. , meta - stable states of bosons in the high orbital bands in optical lattices , and spinful bosons with so coupling .
the first mechanism of orbital bosons is essentially an interaction effect which is characterized by the hund s rule . with the orbital degeneracy
, bosons favor to enlarge their spatial extension to reduce the inter - particle repulsion , which results in the maximization of their onsite oam .
we reviewed the ordering of the oam moments in both the square and triangular lattices due to the inter - site coupling in both the weak and strong coupling limits .
the low energy excitations include both the gapless phonon mode and the gapped orbital - flip mode .
the survival of the oam ordering in the soft mott - insulating regime is also discussed .
the second mechanism employing so coupling , which linearly depends on momentum , is a kinetic energy effect . in this case , taking the absolute value of a non positive - definite wavefunction will change its energy and thus invalidate feynman s proof .
we have shown that the condensate wavefunction in a harmonic trap becomes a half - quantum vortex and develops skyrmion - like spin textures . in both cases ,
tr symmetry is spontaneously broken which are not possible in the conventional bec .
there are a number of interesting open issues for future study .
for example , the ordering of oam moments in various different lattices and in three dimensions generates a variety of challenging problems of the lattice gauge theory @xcite .
more importantly , in order to facilitate the communication between the theory work with cold atom experiments , detailed calculation and optimization of the life time of orbital bosons in different lattices are highly desired .
furthermore , the dynamics of bosons in high orbital bands in another challenging topic . taking into account the tremendous progress of cold atom physics
, it would be great if this research can stimulate general interests on the unconventional condensates beyond the `` no - node '' theorem .
i thank l. balents , d. bergman , h. h. hung , w. c. lee , j. moore , i. mordragon - shem , w. v. liu , s. das sarma , v. stajonovic , and s. z. zhang for fruitful collaborations on this topic . in particular , i am grateful to w. v. liu for his introducing me to this research direction , and s. das sarma for his guidance in the research .
i also thank i. bloch , l. butov , h. deng , l. m. duan , m. fogler , j. hirsch , t. l. ho , n. kim , z. nussinov , l. sham , y. yamamoto , s. c. zhang and f. zhou for helpful discussions .
this work is supported by the nsf grant no .
dmr-0804775 and sloan research foundation .