coupled spin 1/2 chains have attracted much attention lately due to the large number of their experimental realizations , as well as to the variety of theoretical techniques , both analytical and numerical , available to study the relevant models .
spin 1/2 chains , owing to the jordan - wigner transformation@xcite , show properties remarkably similar to those of interacting one dimensional fermions with their low energy properties described by an effective luttinger liquid theory@xcite .
recently it has been realized that these properties are drastically modified when the spin 1/2 chains are coupled together forming ladder systems@xcite . in this case , in a way very similar to haldane s spin - s problem @xcite , a gap is found to open for an even number of chains while the system remains gapless if the number of chains is odd .
this phenomenon has been thoroughly investigated both analytically @xciteand numerically @xcite , and corresponding experimental systems were identified .
typical examples of two - chain spin 1/2 ladders exhibiting a gap are @xmath1@xcite and @xmath2@xcite . on the other hand ,
an example of a gapless three - chain ladder is @xmath3@xcite . from the theoretical point of view
, the difference between odd and even number of legs@xcite or odd and even spin@xcite has been understood qualitatively through a generalization of the lieb - schultz - mattis theorem .
the theorem states that if the ground state is unique the system will necessarily be gapless .
this is indeed the case for an odd number of coupled spin 1/2 chains .
more recent theoretical work has emphasized the role of boundary conditions in the transverse direction in the formation of spin gaps .
the results quoted above - the opening of a gap for an even number of coupled chains , gaplessness if the number is odd - are valid for open boundary conditions ( obc ) .
it turns out that when periodic boundary conditions ( pbc ) are imposed on the transverse direction a gap opens for both cases of even and odd number of chains although the underlying reasons and the natures of the gaps are different . in the case of an even number of coupled chains ,
the reason for the gap is the formation of spin singlets along the transverse direction , similarly to the case of chains with open boundary conditions . when the number of coupled chains is odd a two fold degenerate dimerized ground state
is obtained in the case of pbc@xcite - in contrast to its uniqueness in the obc - allowing for a gap in the spin excitations .
the degeneracy of the ground state in the case of pbc can be understood as a consequence of the fact that pbc are frustrating for an odd number of legs . to date , no experimental system described by coupled antiferromagnetic spin 1/2 chains with periodic boundary conditions and an odd number of legs has been reported .
when , further , leg - leg biquadratic interactions are included new states emerge .
a system of two spin-1/2 spin chains enters a spontaneously dimerized phase with a gapped spectrum exhibiting non - haldane spin - liquid properties@xcite .
the elementary excitations are neither spinons nor magnons , but pairs of propagating triplet or singlet solitons connecting two spontaneously dimerized ground states . subsequently more non - haldane spin - liquid models have been proposed @xcite .
recently , these models of spin ladders with biquadratic exchange have been advocated@xcite as possible models for the formation of a spin gap in @xmath0 and @xmath4 .    in this paper
we study the low energy physics of the three - leg ladder with periodic boundary conditions ( see fig .
[ fig : spin_tube ] ) .
this model is called the spin - tube model in the following . by analyzing in detail its effective low energy hamiltonian ( leh ) , which consists of two _
non - equivalent _ coupled xxz chains of spins and chiral degrees of freedom , in the presence of a biquadratic exchange@xcite , we will show that the dimerized ground state of this model falls in the universality class of the non - haldane spin liquids . in the case of xxz chains with a biquadratic coupling and for @xmath5 , we find that the ising antiferromagnetic order can compete with the dimer order , and we will describe the resulting phase diagram .
the organization of the paper is as follows . in sec.[model ] we introduce the spin tube model .
we sketch the derivation of the strong interchain coupling effective hamiltonian and discuss its relation with two spin 1/2 chains coupled by a biquadratic interaction .
we then discuss the symmetries and recall the results already known on the two chains with biquadratic interactions , in particular the prediction of a spin gap .    in sec .
[ bosonization ] we discuss the bosonization treatment of two non - equivalent xxz chains coupled by a biquadratic interaction .
we show that this hamiltonian contains a term @xmath6 that is responsible for the formation of a spin - gap and singlet order , as well as terms @xmath7 that cause antiferromagnetic order . to analyze the competitions of these terms ,
we derive renormalization group equations . using these equations
, we estimate in which part of the phase diagram we should expect a competition of antiferromagnetism and dimer order .
we also estimate the equation of the phase boundary between the singlet and the antiferromagnetic order .    in sec.[sec :
dimerized ] we analyze in details the dimerized phase .
the dimerized ground state is two - fold degenerate and is formed of alternating singlets of spins and pseudospins .
the elementary excitations above the ground state carry a spin @xmath8 as well as a pseudospin @xmath8 .
they can be seen as formed by the introduction of a spin and a pseudospin in the pattern of alternating singlets .
similar excitations were obtained in ref .   and it was shown that they lead to response functions very different from those of a haldane spin liquid .
numerical evidence for such excitations was also obtained in ref . .
in section [ sec : gap_closure ] , we analyze the effect of a magnetic field strong enough to close the gap of the dimer phase . when both chains are magnetic , the situation is similar to the one obtained with two coupled spin ladders @xcite . in the more physical situation , when only one of the two chains is magnetic , the ground state properties are those of a two component luttinger liquid . in the case of spin - orbital models , spin and orbital degrees of freedom
decouple completely but the presence of orbital modes should affect the specific heat . in the case of the spin - tube , we show that in contrast to spin - orbital models , the spin - correlation functions are affected by the presence of auxiliary gapless modes .
we also show that for non - equivalent chains ( a case that is realized in spin - orbital models ) the exponent of the spin correlations is _ non - universal _ at the transition .
finally , in sec.[conclusions ] we will give some concluding remarks .
we wish to study the three - chain ladder hamiltonian ,    @xmath9    where @xmath10 ( resp .
@xmath11 ) is a chain ( resp .
site ) index , @xmath12 is the coupling along the chain and @xmath13 the transverse coupling .
we impose periodic boundary conditions along the rungs by identifying @xmath14 .
we call the resulting model the spin - tube ( see fig .
[ fig : spin_tube ] ) since it can be realized by placing 3 spin 1/2 chains forming an equilateral triangle . in order to investigate the low - energy physics of the spin - tube
, we consider the limit of strong interchain coupling ( @xmath15 ) .
this is the appropriate starting point that yields a good effective description of the properties of the spin tube in the whole range of @xmath16 .
( starting from the opposite limit @xmath17 , and treating @xmath18 as a perturbation , one finds@xcite that it gives rise to relevant terms in the hamiltonian . as a result ,
the initial @xmath18 grows until it is of order of @xmath12 at which point the weak coupling bosonization scheme is no more valid . )    to derive the effective low - energy hamiltonian , let us first consider the case @xmath19 .
the system is a collection of independent rungs , each described by the following hamiltonian @xmath20    the ground state of each rung is four - fold degenerate , composed of two doublets of spin 1/2 excitations , corresponding to the left and right chirality ( -/+ ) , with energy @xmath21@xcite , and the excited states form a spin 3/2 quadruplet with energy @xmath22 .
turning on @xmath23 allows the rungs to exchange spins . in the limit @xmath15
, the quadruplet of s=3/2 excitations can be neglected , and only the degenerate low energy subspace of spin 1/2 states needs to be taken into account . in this subspace
the hamiltonian transforms into an effective hamiltonian with a biquadratic coupling between the spin and chirality degrees of freedom@xcite ,    @xmath24    where @xmath25 is the total spin , and @xmath26 are operators exchanging left and right chiralities .
the original spin operators can be expressed in terms of the effective spin @xmath25 and the chirality @xmath27 in the following way : @xmath28 the spin tube model has already been investigated numerically@xcite using dmrg in the case of coupled xxx chains .
the system was shown to exhibit a spin gap @xmath29 , with exponentially decaying correlation functions @xmath30 and @xmath31 .
this behavior was shown to be related to the formation of a dimer order .
a qualitative discussion of the origin of this dimer order can be found in ref .
@xcite . in ref .
the dispersion relation of excitations having @xmath32 and @xmath33 as a function of momentum was obtained numerically for a system of 12 sites showing only gapped excitations .    in ref .
, a generalization of the effective hamiltonian eq .
( [ heff ] ) to coupled anisotropic spin chains has been derived ( see eq .
( [ eq : heff_anisotropic ] ) .
numerical diagonalizations were performed .
a gap to @xmath34 excitations was obtained for@xmath35 with results in good agreement with those of ref . for @xmath36 .
it was also shown that the ground state was degenerate in agreement with the dimer order picture of ref . .
the dispersion relation of excitations having @xmath37 and total spin @xmath38 or @xmath39 was obtained@xcite .
it was shown to be the bottom of a two particle continuum .
the fundamental particle , the spinon , was conjectured to have @xmath40 .
the spinon dispersion relation was obtained numerically by considering a system with an odd number of sites@xcite , showing that the spinons were massive for all momentum .
the hamiltonian ( [ heff ] ) is part of the class of the hamiltonians consisting of two non - equivalent coupled xxz chains in the presence of a biquadratic exchange , @xmath41 where h_0 & = & _ i _
= 1,2 j _ ( s_,i^xs_,i+1^x+s_,i^ys_,i+1^y ) + j_^z s_,i^zs_,i+1^z [ h0 ] + h_b&= & _ i ( s_1,i^xs_1,i+1^x+s_1,i^ys_1,i+1^y+ _ 1 s_1,i^zs_1,i+1^z ) ( s_2,i^xs_2,i+1^x+s_2,1^ys_2,i+1^y+ _ 2 s_2,i^zs_2,i+1^z).[hbiq ] in the case of the spin - tube , @xmath42 corresponds to spin @xmath43 , and @xmath44 is associated with the chiral degrees of freedom @xmath26 .
another way of writing a class of hamiltonian generalizing the spin tube model ( [ heff ] ) is :    @xmath45    the effective hamiltonian for the spin tube is obtained for ,    @xmath46    for a tube made of xxz chains , one has , instead@xcite : @xmath47 the parameters @xmath48 and @xmath49 in eq .
( [ ham ] ) measure the xxz anisotropy for spin and chirality , respectively . when both of them are equal to 1 , the hamiltonian eq .
( [ ham ] ) is @xmath50 symmetric . for @xmath51 , @xmath52 ,
@xmath53 the two chains are equivalent , and can be parameterized as :    @xmath54    hamiltonians of the type eq .
( [ ham ] ) can be mapped onto hamiltonians of the type eq .
( [ eq : basic_hamiltonian ] ) .
the correspondence is given by : @xmath55 with the identification @xmath56 , @xmath57 .
since writing the hamiltonian in the form ( [ eq : basic_hamiltonian ] ) is less restrictive than in the form ( [ ham ] ) ( i.e. the former includes the case @xmath58 ) , we focus on the hamiltonian ( [ eq : basic_hamiltonian ] )
. hamiltonians of the type eq .
( [ eq : basic_hamiltonian ] ) are also encountered in a different context than the spin - tube model .
in particular , a hamiltonian of the type ( [ eq : basic_hamiltonian ] ) has been proposed by mostovoy and khomskii as a model for the spin gap formation in the @xmath0 @xcite and @xmath4 @xcite compounds . in that case , the @xmath59 spins correspond to the real spin of the system whereas the @xmath60 spins are pseudospins associated with _
orbital _ degrees of freedom .
these spin - orbital models can be derived from a multiband hubbard model .
the derivation is reviewed for instance in ref . .
let us discuss first the case @xmath61 , @xmath62 .
the hamiltonian describes two coupled heisenberg chains with a biquadratic coupling preserving the @xmath63 symmetry .
actually , the full symmetry group is larger than su(2 ) .
one has : @xmath64=0 \hspace{0.5 cm } [ h,\vec{s}_{2,\text{tot.}}]=0\ ] ] and the full symmetry group is therefore @xmath50 rather than the su(2 ) symmetry that follows from @xmath65=0 $ ] . as a result
the spectrum consists of @xmath66 multiplets@xcite . for @xmath67 ,
the hamiltonian ( [ eq : basic_hamiltonian ] ) has been shown to have an even larger @xmath68 symmetry and to reduce to an integrable @xmath68 spin chain @xcite .
the spectrum of the su(4 ) spin chain has been obtained by the bethe ansatz @xcite , and the correlations functions have been obtained by non - abelian bosonization techniques@xcite , identifying the low energy effective theory describing the spin chain as the @xmath69 wznw model .
the integrable @xmath68 spin - chain has also been intensively studied numerically using quantum monte carlo ( qmc ) @xcite or density matrix renormalization group ( dmrg ) and lanczos exact diagonalization ( ed)@xcite in the context of the spin - orbital models .
the ground state energy , and excited state energy were obtained in good agreement@xcite with analytical calculation using the bethe ansatz@xcite .
the numerical calculation of the correlation functions @xcite reproduces results of the the continuum field theoretical treatment @xcite . at @xmath70 perturbations
are generated which lower the @xmath68 symmetry to @xmath71 and render the hamiltonian non - integrable .
these perturbations have been recently studied @xcite by describing the chain away from the integrable point as a perturbed @xmath69 wzw model .
it was found that for @xmath72 , a gapless phase is obtained while for @xmath73 , a gap is formed . a different field theoretical treatment in the limit
@xmath74 also leads to the appearance of a gapped phase .
for @xmath75 , the ground state wavefunction could be obtained exactly in matrix product form@xcite , of singlet states along the legs of the ladder .
this picture is in good agreement with the predictions of the field theoretical treatment .
the dependence of the gap on the coupling for the range @xmath76 was obtained by dmrg calculations@xcite .
it was found that for @xmath77 , the gap increases proportionally to @xmath78 in agreement with the weak coupling bosonization treatment , and vanishes for @xmath79 , which is the su(4 ) symmetric point as predicted .
the dmrg calculations of ref . showed however a power law gap opening .
such power law gap opening can only be explained by the presence of a relevant operator in the continuum description .
however , no such operator was obtained in the bosonization treatment@xcite . moreover ,
if a relevant operator was present in the continuum theory , the absence of a gap at @xmath80 would be the result of the coefficient of this operator vanishing precisely at @xmath80 .
but then , in contrast to what is observed in numerical calculations , a gap would also obtain for @xmath81 .
a solution to this puzzle has been suggested recently@xcite . in ref . , the gap has been calculated by assuming that the first excited state was in the subspace @xmath82 .
this assumption was shown to be incorrect in the gapped phase in which the first excited state lies in the subspace @xmath83 .
when corrected@xcite , a slow increase with @xmath84 is found , compatible with an exponential gap opening .
the correlation functions @xcite do not show incommensuration , in agreement with the field theoretic approach @xcite .
taking @xmath85 and @xmath62 in ( [ eq : basic_hamiltonian ] ) preserves the @xmath71 symmetry .
this case has been investigated numerically@xcite with @xmath86 .
it was shown that a gapless
gapped phase transition obtained as @xmath87 was increased .
this work was followed by analytical investigation based on the perturbed @xmath69 wzw continuum theory@xcite .
the analytical investigations established the existence for a given @xmath78 of an extended gapless region in the plane @xmath88 that contains the line @xmath89 previously discussed .
in the present work , we consider the general case of @xmath90 and @xmath91 .
this case includes in particular the spin tube model .
we will focus on the regime @xmath92 and apply methods similar to those of ref . .
in this section , we derive the phase diagram of two xxz spin chains weakly coupled by a biquadratic exchange .
we first recall the bosonization of a spin chain in sec . [
sec : single_chain ] .
this section can be skipped by readers familiar with bosonization .
then , we discuss the bosonization of the two coupled chains system in sec . [ sec : biquadratic_coupling ] .
this allows us to derive renormalization group ( rg ) equations for the coupling constants of the problem . finally , in sec .
[ phasediagram ] , we discuss the phase diagram deduced from the analysis of rg equations .
in this section , we recall briefly the derivation of the bosonized hamiltonian and spin operators that describe an isolated xxz chain .
we follow the well known abelian bosonization procedure for spins @xcite .
the xxz spin chain is described by the hamiltonian : @xmath93 the xxz spin chain hamiltonian is first transformed into an interacting fermionic system on the lattice by expressing the spin operators @xmath94 , @xmath95 , @xmath96 in terms of fermion operators @xmath97 , @xmath98 , using the jordan - wigner transformation@xcite :    @xmath99    this transformation turns the xxz hamiltonian into a model of spinless fermions with nearest neighbor interaction described by the hamiltonian , @xmath100 for @xmath101 , the hamiltonian ( [ spinlessf ] ) describes non - interacting fermions and is easily diagonalized .
to proceed , we restrict our attention to the low - energy sector of the theory , captured by the continuum theory .
introduce the left ( right ) chiral fermion fields @xmath102 containing momenta close to the fermi points @xmath103 , ( @xmath104 with @xmath98 the lattice spacing ) , @xmath105 the continuum fermi fields are then reexpressed in terms of bosonic field@xcite , as follows :    @xmath106    where the pair of conjugate fields , @xmath107 satisfy the following commutation relation : @xmath108=\imath
\delta(x - x^\prime),\ ] ] and the field @xmath109 , dual to @xmath110 , is defined as : @xmath111    the spin operators eqs .
( [ jwp])([jwz ] ) , can be expressed in the continuum limit as : @xmath112 \mbox { , } \mbox { } s^z(x)=\frac{s^z_n}{a}=-\frac { \partial_x \phi } \pi + \frac{e^{\imath \pi \frac x a } } { \pi a } \cos 2 \phi,\end{aligned}\ ] ]    introducing normal ordering with respect to the fermion vacuum , one has : @xmath113    where @xmath114 indicates normal ordering .
the average @xmath115 is the ground state energy of the chain and can be found by the bethe ansatz@xcite .
however , as long as one is only interested in the correlation functions of the single chain , one can simply drop this contribution and focus on the normal ordered terms that describe the excitations above the ground state .
our task is thus to derive a bosonized expression of the normal ordered product in eq .
( [ eq : normal_order ] ) .
some care is needed in order to obtain correct results@xcite , and one finds : @xmath116 + \frac { e^{\frac { \imath \pi x } a}}{\pi } \sin 2 \phi \nonumber \\ : s_i^zs_{i+1}^z+ s_i^z s_{i+1}^z : & = & \frac 2 { \pi^2 } ( \partial_x \phi)^2 + \frac 2 { ( \pi a)^2 } e^{\imath \frac{\pi x } a } \sin 2 \phi + \frac { 2 \cos 4\phi}{(2\pi a)^2 } + \text{irrelevant terms \ldots}\end{aligned}\ ] ] the oscillating terms in eq .
( [ eq : bosonized_quad ] ) are dropped from the hamiltonian after integration over @xmath117 .
the hamiltonian @xmath118 , eq . ( [ spinlessf ] ) then becomes :    @xmath119-\frac { 2\delta}{(2\pi a)^2 } \int dx \cos(4\phi),\ ] ]    where for @xmath120 , u&=&aj(1+)^1/2 + k&=&(1+)^-1/2 + & = & j_za.[deltaequiv ]    thus , the bosonized form of @xmath118 reduces to a sine - gordon hamiltonian , where the cosine terms come from the intrachain umklapp process@xcite .
the renormalization group treatment shows@xcite that in the vicinity of the xy point , the cosine terms are irrelevant , so that asymptotic properties are described by a free scalar field with renormalized @xmath121 . since the xxz chain is integrable , it can be shown that the gapless spectrum extends to the whole region @xmath122 .
moreover , it is possible to obtain an analytic expression for the renormalized @xmath121 from the exact solution@xcite .
one finds : @xmath123 @xmath124 the isotropic point @xmath36 ( antiferromagnetic heisenberg model ) corresponds , in the bosonization description , to @xmath125 and @xmath126 . at this point @xmath127
is marginally irrelevant . for @xmath128
, the exact solution shows that a gap opens in the excitations of the system@xcite , and that an ising order of the spins along the z axis is obtained .
this result can also be found from a bosonization procedure valid in the vicinity of the isotropic point .
we now proceed to derive , using the results reviewed in the previous subsection , the bosonized hamiltonian of two non - equivalent coupled xxz chains . to bosonize the hamiltonian @xmath129 , eq .
( [ h0 ] ) describing two decoupled xxz chain we introduce two pairs of dual fields ( one pair for each chain ) @xmath130 ( @xmath131 ) as defined in sec .
[ sec : single_chain ] .
the bosonized form of the hamiltonian @xmath129 is then :    @xmath132-\frac{2 \delta_\alpha}{(2\pi a)^2}\int dx \cos 4\phi_\alpha \right\},\end{aligned}\ ] ]    with the fields @xmath133 and @xmath134 having _ a prori _ different velocities and luttinger couplings , @xmath135 the bosonization formulas for the spins ( [ spin - operators ] ) are unchanged except for the obvious introduction of a chain index .    in order to have the full bosonized hamiltonian , we now have to derive the bosonized form of the biquadratic exchange ( [ hbiq ] ) .
the first step is to normal order using eq .
( [ eq : normal_order ] ) .
this step is important since it leads to non trivial contributions to the quadratic part of the hamiltonian .
we introduce the terms : @xmath136 where @xmath131 that measure the strength of these corrections .
having expressed the hamiltonian in a normal ordered form , we can apply eq .
( [ eq : bosonized_quad ] ) to obtain the bosonized expression .
the oscillating terms of eq .
( [ eq : bosonized_quad ] ) give rise to an interchain coupling term , the strength of which is given by :    @xmath137    collecting everything , one obtains the following bosonized biquadratic exchange : @xmath138 \nonumber \\ & + & a_2 \int \frac{dx}{\pi } \left[(\pi \pi_1)^2 + ( 1 + 2\delta_1/\pi ) ( \partial_x \phi_1)^2 \right ] + \text{irrelevant terms \ldots},\end{aligned}\ ] ] the terms in @xmath139 form a mean - field like interchain interaction .
for @xmath140 , these terms merely produce a renormalization of the velocities and luttinger liquid exponents with respect to the decoupled chains . however , in the case where @xmath141 or @xmath142 , these terms become crucial since they give a finite velocity to the @xmath133 ( resp .
@xmath134 excitations ) excitations .
such a case is realized in the spin - tube problem where the exchange constant of the pseudospins @xmath27 is zero .
we see that the mean - field like contribution of the spin fluctuations contributes in that system to provide an exchange constant and thus a finite velocity to the pseudospin excitations . of course , in such case ,
the bosonization procedure is not really justified since interactions are of the order of magnitude of the bandwidth of the spin or pseudospin excitations .
however , it is usual in quasi - one dimensional systems to have a continuity between the weak and the strong coupling regime .
moreover , a mean field theoretical treatment in the xy limit ( @xmath143 ) leads to similar results to the bosonization treatment .
details can be found in app .
[ app : mean - field ] .
we will therefore assume that although not fully justified in the spin tube case , bosonization nevertheless leads to qualitatively correct results concerning the phases of the system and the overall behavior of correlations in these various phases .
a quantitative treatment ( in particular of the phase boundaries ) requires numerical simulations that are beyond the scope of this paper .
the present treatment shows us that in weak coupling the two chains are described by a bosonized hamiltonian : @xmath144}\nonumber \\ & + &   \frac { 2g}{(2 \pi a)^2}\int dx \sin 2\phi_1 \sin 2\phi_2   \end{aligned}\ ] ]    we shall analyze the model perturbatively , with results valid up to a given value of @xmath78 : it is known that in the strong coupling regime of the two chains with biquadratic exchange there is a special value of the interchain coupling at which the model has @xmath68 symmetry . at this special point ,
the model develops a gapless phase described by a @xmath69 wznw model .
such an effect is non - perturbative : the resulting critical point has a conformal anomaly @xmath145 whereas the original unperturbed model has @xmath146 .
this implies by zamolodchikov s @xmath147theorem@xcite that the transition to the @xmath69 can not be predicted by a rg calculation that always leads to a decrease of @xmath147 . beyond this special value of @xmath78
, the weak coupling theory would lead to incorrect predictions and an alternative approach such as the one of lecheminant and azaria is needed . on the other hand ,
it is expected to give qualitatively correct predictions when the coupling is smaller than the critical value . in the remainder of the paper we shall thus work in the weak coupling regime where the weak coupling theory is valid . in the following section
, we will discuss the rg treatment of the weak coupling model .      in the following analysis
, we will neglect the velocity difference between chains 1 and 2 since usually velocity differences do not play an important role in the derivation of the phase diagram by rg techniques .
however , for the sake of completeness , we have given in app .
[ app : momentum_shell_rg ] the rg equations for non - equal velocities derived using a momentum shell integration technique@xcite .
when velocity differences are neglected , the renormalization group equations for @xmath148 can be easily derived from the hamiltonian ( [ non - equivalent - chains - ni])([hbiqbos ] ) using operator product expansions ( ope ) @xcite .
the renormalization group equations for @xmath149 neglecting velocity differences are : @xmath150 while the renormalization group equations for @xmath151 are : @xmath152    we now proceed to deduce the weak coupling phase diagram of the model
( [ eq : scaling_rge ] ) indicates that @xmath153 is a relevant variable when @xmath154 , while the variables @xmath155 become relevant when respectively , @xmath156 , @xmath157 .
there are therefore four cases to distinguish . in the first case : @xmath158 , there are no relevant operators and the system is in a gapless state . in the second case : @xmath159 , @xmath160 , there is a single relevant operator ; for @xmath161 ( or equivalently @xmath162 ) , there are two , while for @xmath163 there are three relevant operators .
these four different regions are shown in fig .
[ fig : phase_diagram ] .    in the presence of relevant operators , rg equations cease to be valid as soon as the dimensionless coupling become @xmath164 , and we have to determine the nature of the strong coupling fixed points in order to predict the phase diagram .
we see that @xmath151 are driven to zero by the flow of the rg when there is a relevant operator so that the fields become classical . as a result ,
@xmath165 are locked at average values @xmath166 that minimize the ground state energy and a gap opens in the excitations of these fields@xcite .
when these fields are locked , it is also known that the exponentials of the dual field have exponential decay@xcite .
let us consider first the case with @xmath153 the only relevant operator .
then , the minimization of the ground state energy requires @xmath167 , i.e. @xmath168 .
it is then clear that @xmath169 as well as @xmath170 , which are the staggered component of @xmath96 and @xmath94 , will display an exponential decay .
it can also be shown that the uniform component of the spins also present an exponential decay . as a result
we have a spin - liquid phase that presents only short range order in all its spin correlations . as we shall see , this spin - liquid phase is a dimerized phase whose properties are discussed extensively in sec . [ sec : dimerized ] . in the case
@xmath171 ( decoupled chains ) the analysis is even simpler .
then the two chains remain decoupled , and we are left with the analysis of the usual sine - gordon model .
it is then known that when @xmath172 is relevant , the field @xmath173 is locked .
the analysis of the resulting strong coupling fixed point can be found for instance in ref . .
it is found that the strong coupling fixed point is associated with the ising antiferromagnetic phase of the xxz chain at @xmath128 in which the staggered component of @xmath96 has a non zero expectation value in the ground state .    when @xmath153 and at least one of the @xmath172 are relevant , there are two possible candidates for the ground state
. considering eq .
( [ eq : scaling_rge ] ) , we see that the effect of the biquadratic interchain interaction is to reduce the effective @xmath174 .
physically , this means that the tendency to form singlets competes with the tendency to form an ising antiferromagnet .
two scenarios are possible .
one is that there is a well defined phase boundary between a pure spin - liquid state and a pure ising antiferromagnet state .
the second scenario is that there is a crossover between the two pure states as @xmath175 is varied . in such case
, increasing @xmath78 would lead to a gradual disappearance of antiferromagnetic order leaving a purely singlet state as @xmath176 .
since in both phases it is the same field that orders , there is _ a priori _ no reason to exclude a mixed spin - liquid antiferromagnet order .
thus , the first scenario appears extremely unlikely .
a numerical investigation of the crossover could be very interesting as a toy model of a crossover from spin - liquid to antiferromagnetism .
it is interesting to remark that if @xmath177 , and @xmath178 , the equations ( [ eq : scaling_rge ] ) can be integrated analytically .
two phases are obtained , separated by a line @xmath179 . in the first one ,
@xmath180 and @xmath181 which corresponds to a ground state with singlet order . in the second one , @xmath182 and @xmath183 , which corresponds to antiferromagnetic order .
the rg equations cease to be valid for @xmath184 or @xmath185 .
if when this scale is reached @xmath153 and @xmath186 are of the same order of magnitude , there is a possibility of obtaining a mixing of antiferromagnetism and dimer order . note that even on the line @xmath187 , there is a finite correlation length .
this is a further evidence for a progressive crossover from dominant antiferromagnetic order to dominant dimer order .
it is also possible to give a purely classical treatment for @xmath177 by simply minimizing the ground state energy with respect to @xmath188 and @xmath189 . in the case @xmath190
, one finds that there are three different regimes . for @xmath191
one obtains @xmath192 corresponding to a purely dimerized phase . for @xmath193 ,
one obtains @xmath194 and @xmath195 .
this corresponds to persistence of dimerization in the chain with the smallest tendency to antiferromagnetic order , whereas the chain with the strongest tendency to antiferromagnetism is found in state with mixed dimer and antiferromagnetic order .
the antiferromagnetic order parameter in that chain , @xmath196 , then assumes the value @xmath197 .
finally in the region @xmath198 , both chains display antiferromagnetism with @xmath199 .
these results are summarized on figure [ fig : classical_phase_diag ] . for @xmath200 , the region with mixed antiferromagnetic and dimer
order shrinks to a single point
. it would be interesting to see how quantum fluctuations affect the present picture and in particular determine whether sharp transitions are preserved or if they evolve into crossovers .      in this section
, we try to estimate the position of the crossover between the spin - liquid and the antiferromagnet .
this can be done roughly by comparing the correlation lengths in the antiferromagnet and in the spin - liquid phase . using the rg equations , for small @xmath201
, we can neglect the renormalization of @xmath151 .
this leads to :    @xmath202    when any of these quantities become of the order of the energy cutoff @xmath203 , the rg equations cease to be valid and the phase that is obtained is determined by minimizing the ground state energy .    for @xmath153
, the strong coupling is obtained at a length scale : @xmath204 this is the correlation length of spin fluctuations in the spin - liquid phase . for @xmath205,@xmath206 the strong coupling is obtained respectively at length scale : @xmath207 it is clear that the shortest length corresponds to the first operator to attain strong coupling .
therefore , the phase that is obtained is the one with the shortest correlation length . for @xmath177 ,
this is in agreement up to a constant with the criterion derived from the rg equations ( [ eq : scaling_rge])([eq : second_order_rge ] ) .
the comparison of correlation length allows to draw a rough phase boundary between the antiferromagnet and the dimerized phase that could also be obtained by numerically integrating the rg equations ( [ eq : scaling_rge])([eq : second_order_rge ] ) starting from weak coupling and any @xmath151 .
the equation of the phase boundary is in the case of equivalent chains @xmath178 , @xmath208 : @xmath209 let us note that in the isotropic spin - tube case , the operators causing antiferromagnetic order are ( marginally ) irrelevant so that there is only singlet order . however , in the case of an anisotropic spin tube ( 3 coupled xxz spin chain with @xmath210 ) , such competition becomes possible .
completely decoupled chains exhibit for @xmath210 an ising antiferromagnetic phase . introducing a strong enough biquadratic interchain coupling favors on the other hand a spin liquid phase .
the competition of the two should produce a crossover from the ising antiferromagnet to the spin liquid of the type discussed in the preceding section .
in this section , we discuss the properties of the spin liquid phase . in the case of equivalent chains , further progress can be made by using symmetric and antisymmetric modes allowing in particular a refermionization of the problem and the calculation of some correlation functions@xcite .
this will be discussed in sec .
[ sec : equivalent_chains ] . in the general case with inequivalent chains ,
such decoupling is no longer possible .
however , it is still possible to present a simple semiclassical picture of the nature of excitations above the ground state
. this will be the subject of section [ sec : non_equivalent_chains ] .
we will focus on region @xmath211 , @xmath159 , in which the only relevant operator is the biquadratic exchange , @xmath212 .
the ground state shows long - range order of the fields @xmath133 and @xmath134 .
the expectation values of the ordered fields are : @xmath213 and as a result , @xmath214 . using eq .
( [ eq : bosonized_quad ] ) , this implies that a dimerized order develops both in spin variables and pseudospin variables , i.e. @xmath215 ( @xmath131 ) . in parallel with that
, we have @xmath216 , so that by eq .
( [ spin - operators ] ) @xmath217 , and the correlation functions @xmath218 as @xmath219 .
it is also well known that when the fields @xmath220 are ordered , the correlation functions of the disorder operators @xmath221 decay exponentially at large distances .
using again the bosonization formulas for the spins , eq .
( [ spin - operators ] ) , this implies an exponential decay of all correlation functions : @xmath222 and @xmath223 where @xmath224 .
therefore , the dimerized phase appears as a spin liquid state formed of singlets of spins on both chain 1 and chain 2 .
such a conclusion was reached previously in a numerical investigation of the spin tube @xcite and by considering the equivalent isotropic chains at the solvable point@xcite @xmath225 where the ground state wavefunction can be obtained exactly in matrix product form .
our results show that such mechanism of spin liquid ground state formation does not require su(2 ) symmetry .
this mechanism is somewhat reminiscent of the spin - peierls transition@xcite , the pseudospins playing here the role of the phonons .
this is the essence of the mostovoy - khomskii model@xcite for the `` spin - peierls '' transition at @xmath226k in @xmath0 .    in the case of the spin - tube the same picture of the ground state obtains , with the ground state of the spin - tube formed by _
singlet of spins _ on even bonds and _ singlet of chiralities _ on odd bonds or by _ singlet of spins _ on odd bonds and _ singlet of chiralities _ on even bonds ( see fig .  [
fig : ground_state ] ) .    for the moment
, we have only been able to discuss the nature of the ground state .
however , it is also important to discuss the nature of the excitations as well as the various correlation functions . in order to do that
, it is worth to restrict first to the simple case of two equivalent chains , in which the physical picture is the clearest .
when the two chains are equivalent , we introduce the fields@xcite : @xmath227 and their conjugate fields : @xmath228    so that the total hamiltonian can be completely decoupled into the symmetric and antisymmetric parts ,    @xmath229    @xmath230 - \frac{g}{(2\pi a)^2 } \int dx \cos \sqrt{8}\phi_s\ ] ]    @xmath231 + \frac{g}{(2\pi a)^2 } \int dx \cos \sqrt{8}\phi_a,\ ] ]    where the magnetic field couples only to @xmath232 , and only the most relevant operators have been taken into account .
the elementary excitations can be discussed in terms of solitons of two decoupled sine - gordon models .
the solitons of the hamiltonian @xmath233 carry a total spin @xmath234 whereas those of hamiltonian @xmath235 carry a spin 0 as they do not couple to the magnetic field .
these solitons are represented in fig .
[ fig : solitons ] .
it is also convenient to use the canonical transformation @xmath236 , @xmath237 followed by a refermionization . introducing the fermion operators : @xmath238    where @xmath239 , one can finally rewrite @xmath240 in the form : @xmath241 where @xmath242 and @xmath239 .
the couplings are given by : @xmath243    at the isotropic point ( @xmath244 ) , one has @xmath245 , so that @xmath246 is a free fermion hamiltonian @xcite .
similarly to the spin ladder , where @xmath247=0 $ ] , the excitation spectrum can be split into a singlet and a triplet with spin @xmath248 .
however , in contrast to the spin ladder case , the triplet and the singlet here have the same mass .
this is the signature of a larger symmetry group , @xmath249 .
the correlation functions can be obtained from mapping the free fermion hamiltonian onto two non - critical ising model @xcite exhibiting the @xmath250 symmetry .
remarkably , although the system has a spin gap , these correlation functions are very different from those of a spin-1 chain @xcite or a spin ladder @xcite . in the chain with biquadratic exchange
, the response functions do not show any particle - like delta function peak in their imaginary part , but only a two particle continuum@xcite even in the vicinity of @xmath251 .
this is to be contrasted with the spin ladder @xcite which shows a delta function peak associated with a single particle excitation at @xmath251 .
the two chains with biquadratic interactions thus form a `` non - haldane '' spin liquid .    away from the isotropic point
, the hamiltonian can still be refermionized but the fermions ( solitons ) have interactions which preclude a mapping on a non - critical ising model .
however , in the anisotropic case , when the antiferromagnetic intrachain interaction is irrelevant one has @xmath252 .
this implies that no bound states of solitons can be formed since there is no coherent propagation of two solitons@xcite , hence the absence of a single particle peak .
a calculation of correlation functions is in principle possible using as form factors approach , but this is far beyond the scope of the present paper . since @xmath235 and @xmath233 remain decoupled , the excitations having @xmath234 do not interact with the excitations having @xmath253 .
this can be understood as a consequence of the @xmath254 symmetry of the problem and the resulting separate conservation of @xmath255 and @xmath256 .      in the case of two non - equivalent chains
, one can not decouple the hamiltonian into two sine - gordon hamiltonians .
therefore , determining the nature and the quantum numbers of elementary excitations is not as straightforward .
the preceding canonical transformation leads to the following hamiltonian : @xmath257 + \frac{(u_1/k_1+u_2/k_2 ) } 2 [ ( \partial_x \phi_s)^2 + ( \partial_x \phi_a)^2]\right .
\nonumber \\   & + & \left .
\frac { u_1 k_1 -u_2 k_2}{2 } \pi^2 \pi_s \pi_a + \frac { u_1 /k_1 - u_2 /k_2 } 2 \partial_x \phi_s \partial_x \phi_a \right]+ \frac{2g}{(2\pi a)^2 }
\int dx [ \cos \sqrt{8 } \phi_s -\cos \sqrt{8 } \phi_a ] \end{aligned}\ ] ] shifting @xmath258 , we see that the resulting hamiltonian has a @xmath259 invariance under @xmath260 .
a rescaling @xmath261 , @xmath262 , and a refermionization brings the hamiltonian to the form : @xmath263    \nonumber \\ & + &   g(\rho_{r , s}(x)\rho_{r , a}(x)+\rho_{l , s}(x)\rho_{l , a}(x))+ \tilde g_1 \sum_{r\ne r'}\rho_{r , r}(x)\rho_{l , r'}(x ) + \tilde g_2 \sum_r \rho_{r , r}(x)\rho_{l , r}(x)\end{aligned}\ ] ] as a result , there is now an interaction between the excitations of spin @xmath264 ( the @xmath98 fermions ) and the excitations of spin @xmath265 ( the @xmath266 fermions ) .
we have : @xmath267\\ g&=&\pi \left [ u_1\left(k_1+\frac 1 { 4k_1}\right ) -        u_2\left(k_2+\frac 1 { 4k_2}\right)\right ] \\ \tilde g_1&= & \pi\left [ u_1\left(\frac 1 { 4k_1}-k_1\right ) +        u_2\left(k_2-\frac 1 { 4k_2}\right)\right ] \\   \tilde g_2&= & \pi\left [ u_1\left(\frac 1 { 4k_1}-k_1\right ) -        u_2\left(k_2-\frac 1 { 4k_2}\right)\right ] \\
m & = & \frac{g}{4 \pi a } \end{aligned}\ ] ]    the fermionic version of the model is a generalization of the massive thirring model .
the fields carry spin and the interactions break spin rotation symmetry .
it can be checked that the interaction of @xmath98 with @xmath266 fermions disappears only for @xmath268 , @xmath269 .
not much can be said of the elementary excitations of hamiltonian ( [ eq : fermion_non_equivalent ] ) due to the absence of an exact solution . in order to gain some insight into the elementary excitations of the dimerized state in the case of non - equivalent chains , we resort to semi - classical approximations . clearly , we can search for a semi - classical minimum of ( [ eq : rotated_non_equivalent ] ) with either @xmath270 or @xmath271 .
this corresponds to a single soliton in the system . in the general case ,
such excitations are associated with a single fermion , either @xmath98 or @xmath266 . in order to obtain a physical picture for this type of elementary excitations ,
let us calculate the average `` magnetization '' for the fields @xmath133 and @xmath134 when there is a soliton connecting the two degenerate dimerized ground - states . in the case
@xmath270 , @xmath133 decreases from @xmath272 to @xmath273 while @xmath134 increases from @xmath273 to @xmath272 .
we immediately get :    @xmath274    another solution is obtained by reversing the sign of @xmath220 leading to @xmath275 .
the case @xmath271 corresponds to @xmath276 . in each case , the elementary excitations are formed by breaking a singlet on neighboring sites on each chain .
such objects then propagate coherently . in the case of the spin tube
, this corresponds to having one unpaired spin associated with one unpaired chirality pseudospin .
such an excitation has @xmath277 , @xmath278 .
it is the spinon of ref . .
thus the elementary excitations of the model can be easily visualized as an unpaired spin and an unpaired chirality forming a triplet or a singlet diagonal bond ( see fig .
[ fig : singlet_triplet_ne ] ) .
this is the soliton in the dimer order .
this picture generalizes from isotropic equivalent spin chains @xcite to the case of inequivalent and anisotropic chains .
an open question is whether in the case of inequivalent chains bound states of these elementary excitations can be obtained in contrast to the case of equivalent chains .
a necessary condition is that there exists attractive interactions between @xmath98 and @xmath266 fermions .
the study of such bound states could be of interest in relation with light scattering experiments on @xmath0 .
besides semi - classical approximation , another approximate treatment is possible . in the case where @xmath268 , @xmath279 and @xmath280 , the hamiltonian ( [ eq : full_bosonized_hamiltonian ] ) is the double sine gordon model@xcite , or the bukhvostov - lipatov model@xcite .
the model is known to be bethe ansatz solvable@xcite but its elementary excitations and @xmath25 matrix have only been obtained recently @xcite and shown to be identical to those of the double sine - gordon model on its integrable line .
the spectrum of the model @xcite is consists of four massive particles carrying two quantum numbers @xmath281 , and the @xmath25 matrix is@xcite : @xmath282 where @xmath283 is the @xmath25 matrix of a sine gordon model having the parameter @xmath284 and @xmath285 . as a result ,
one of the sine gordon models is in the attractive regime and has a spectrum made of fermions and bound states of fermions whereas the second sine - gordon model is is the repulsive regime and has a spectrum containing only fermions .
the mass of the fundamental fermion is @xmath286 . for @xmath287 ,
one recovers the equivalent chain and the mass of the bound states is then : @xmath288    since for @xmath289 , the operator @xmath290 is the most relevant , it is reasonable to expect that neither the marginal perturbations due to @xmath291 nor the less relevant perturbations @xmath292 change the gapped nature of the spectrum . in this regime
, one should not observe any bound state , since the @xmath293 bound state only exists for @xmath294 . in terms of the original spin chain
, the double sine - gordon regime should be accessible if one chooses to have one chain with @xmath210 and the other chain with @xmath128 in such way that @xmath280 and @xmath295 .
however , the double sine - gordon regime would not be observed in the spin tube case , in which we should have @xmath296 .
we expect that in the spin tube case , the system has more quantum fluctuations than in the double sine - gordon regime .
therefore , no bound states of spinons will form .
this heuristic argument agrees with the numerical calculations on the spin tube that show no bound state of spinons@xcite .
in this section , we discuss the effect of the application of a magnetic field on the dimerized phase . in general , the application of a magnetic field to a gapped one dimensional spin liquid system results in the closure of the gap for a magnetic field of the order of magnitude of the gap .
below the critical magnetic field the magnetization is zero . above the critical magnetic field
the magnetization increases as @xmath297 , the system is in a single component luttinger liquid state and has incommensurate spin correlations@xcite . here
, we have to consider two _ a priori _ different cases . in the first one and most academic ,
the two coupled chains carry real spins that couple to the magnetic field . in the second case ,
only one chain carry a spin that couple to the magnetic field , the other one carrying only a pseudospin degree of freedom that does not couple to a magnetic field .
this last case corresponds to the spin - orbital models@xcite and to the spin tube model@xcite .
the behavior of the spin - orbital model in a magnetic field has been discussed in refs . and . both references study models with @xmath50 symmetry in the absence of the magnetic field .
it is shown that in the vicinity of the su(4 ) point , the spin - orbital model becomes equivalent to the @xmath298 gross neveu model describing the gapped modes plus a @xmath299 cft that describes the gapless magnetic modes . also , in ref . , the weak coupling case has been discussed .
the case where both chains carry spin will be dealt with in sec .
[ sec : both_w_spin ] and the case of a single chain carrying spin in sec .
[ sec : one_w_spin ] .
if both chains carry spin , the coupling to the magnetic field ( taken parallel to the z axis ) is : @xmath300 in that case , we can use the decoupling of sec . [
sec : equivalent_chains ] , eqs .
( [ hs ] ) and ( [ ha ] ) .
the coupling to the magnetic field is : @xmath301      we consider the case where the two chains are equivalent .
this problem has been discussed in the context of spin ladders@xcite .
the gap in the antisymmetric modes @xmath220 is not affected by the presence of the magnetic field . on the other hand , for a sufficiently large magnetic field @xmath302
, the gap in the symmetric modes closes and the magnetization @xmath303 behaves as @xmath304 for @xmath302 .
moreover , it can be shown that the low energy modes of @xmath232 are described by the following hamiltonian : @xmath305\ ] ] and the exponent @xmath306 takes the universal value 1/2 at @xmath307 .
correlation functions in the incommensurate phase can be obtained in the isotropic case by a calculation similar to the spin - ladder case @xcite .
this time , @xmath308 is finite so that @xmath309 has exponentially decaying correlations . as a result
, one has : @xmath310 a simplified expression for the operator @xmath311 is : @xmath312 and the correlations of @xmath96 are : @xmath313 there are subleading power law corrections at @xmath314 . as shown by furusaki and zhang ,
these corrections are missed if one naively neglects the band curvature after refermionization @xcite .
these corrections can also be obtained@xcite by using the haldane expansion of the spin operators and retaining the terms up to @xmath315 as we did here .      in the case of non - equivalent chains
, the problem gets more complicated .
the magnetic field still couples only to @xmath232 .
however , when the magnetic field @xmath316 exceeds the field @xmath317 needed to close the gap in the @xmath266 modes , the appearance of a non - zero @xmath318 creates an effective magnetic field that couples to @xmath319 . if this magnetic field is not strong enough to close the gap in the @xmath220 modes , the system remains in a one - component luttinger liquid for fields @xmath316 not much stronger than the critical field @xmath317 . for @xmath316 sufficiently large , the generated effective field can close the gap in the @xmath98 modes leading to a two component luttinger liquid . since no experimental ladder system with a biquadratic exchange much larger than the quadratic exchange much larger that the quadratic exchange and made of two non - equivalent chains is presently available , this two - step transition to a two component luttinger liquid is unlikely to be observed experimentally .
this case includes the spin tube and the mostovoy - khomskii model and is therefore the most relevant physically .
this case has been discussed for two coupled @xmath71 symmetric chains in ref . .
in the case where only one chain , say chain @xmath320 to fix notations , carries spin , the interaction with the magnetic field is given by a term : @xmath321 in the case of equivalent chains , we can use the same decoupling of sec .
[ sec : equivalent_chains ] , eqs .
( [ hs ] ) and ( [ ha ] ) as in the preceding section
. however , there is an important difference .
now , we have the coupling to the magnetic field in the form : @xmath322 as a result , now the magnetic field couples to both @xcite @xmath232 and @xmath220 .
moreover , the strength of the couplings is exactly the same .
we start with the discussion of equivalent chains .
then , we discuss non equivalent chain .
we will show that in both case , the closure of the gap leads to a two - component luttinger liquid behavior in contrast with the usual spin - liquid systems that lead to a single component luttinger liquid behavior@xcite
. we will also give the expression of correlation functions in the luttinger liquid phase .
as a result of the symmetry between @xmath220 and @xmath232 , the gap in the symmetric and the antisymmetric mode close simultaneously , leading to a two component luttinger liquid ground state under strong enough magnetic field .
contrarily to the case where the magnetic field couples to both 1 and 2 spins , there is no intermediate single component luttinger liquid phase .
we expect that the magnetization @xmath323 behaves as @xmath324 close to the threshold .
it is also important to note that , the two sine - gordon model being equivalent , one has for the fixed point hamiltonian @xmath325 and @xmath326 . as a result , the fixed point hamiltonian is invariant under any rotation in the @xmath327 plane .
we can thus write the fixed point hamiltonian as : @xmath328\ ] ] where @xmath329 and @xmath330 .
moreover , in the case of @xmath50 symmetry , we have @xmath331 .
it can the be shown easily , using the refermionization procedure that @xmath125 for any magnetic field .
equation ( [ eq : fp_mf_equiv ] ) has important consequences for the correlation functions , which are of the form ( for equivalent chains ) :    @xmath332   \nonumber \\ \langle ( s^z_{\alpha}(x , t)-m_\alpha ) ( s^z_{\alpha}(0,0)-m_\alpha ) \rangle   =    \cos(\pi x ( 1 - 2m_\alpha))(x^2-t^2)^{-k_\alpha}+ \text{constant } \mbox { } \frac{1}{4\pi^2 } \left ( \frac{x^2+t^2}{(x^2-t^2)^2 } \right ) , \nonumber \end{aligned}\ ] ]    where the index @xmath333 indicates the spin in chain 1 or 2 , @xmath334 is the magnetization , with @xmath335 ( total magnetization ) and @xmath336 . from the above expressions we deduce the following : the correlation function parallel to the field , @xmath337 , has a staggered part shifted from the wave vector @xmath251 to @xmath338 , while the correlation function perpendicular to the field , @xmath339 , has an unshifted staggered mode and the uniform magnetization mode shifted to @xmath340 .
the correlation functions for the spin of type 2 , instead , are completely unaffected by the presence of an external magnetic field .
in this section , we consider the case where the two chains are not necessarily equivalent .
in particular , this is the case that is realized in the spin - tube under a magnetic field .    in this case ,
is convenient to use the rotation , @xmath341 to bring the hamiltonian to the form : @xmath342 + \frac{(u_1/k_1+u_2/k_2 ) } 2 [ ( \partial_x \phi_s)^2 + ( \partial_x \phi_a)^2]\right .
\nonumber \\   & + & \left . \frac { u_1 k_1 -u_2 k_2}{2 } \pi^2 \pi_s \pi_a + \frac { u_1 /k_1 - u_2 /k_2 } 2 \partial_x \phi_s \partial_x \phi_a \right]+\nonumber \\ & + & \frac{2g}{(2\pi a)^2 } \int dx [ \cos \sqrt{8 } \phi_s -\cos \sqrt{8 } \phi_a ] -\frac h { \sqrt{2 } } \int dx \partial_x(\phi_s+\phi_a)\end{aligned}\ ] ]    shifting @xmath343 , renders hamiltonian invariant under the interchange of @xmath232 and @xmath220 . as a consequence ,
the gaps in @xmath232 and @xmath220 are _ identical _ and have to close simultaneously .
therefore , as in the case of identical chain , there is a transition from a gapped phase into a gapless two - component luttinger liquid phase .
the most general hamiltonian for a two - component luttinger liquid is : @xmath344.\end{aligned}\ ] ] in the case of coupled xxz chain , the effective hamiltonian ( [ eq:2component_ll ] ) can be simplified by making use of symmetry .
we know that the original hamiltonian , eq .
( [ eq : rotated_mf ] ) is invariant under the transformation : @xmath345 this symmetry must be preserved by the renormalized hamiltonian .
therefore , one must have @xmath346 and @xmath347 . as a result
, the effective hamiltonian ( [ eq:2component_ll ] ) is diagonalized by returning to the original variables @xmath348 .
the effective hamiltonian is therefore the following : @xmath349\ ] ] the absence of coupling between @xmath133 and @xmath134 in this hamiltonian is somewhat surprising .
this can however be understood by the fact that the hamiltonian even in the presence of a magnetic field is invariant by a rotation by @xmath350 @xmath351 .
thus , @xmath352 for any @xmath316 , implying that @xmath353 terms can not appear in the two component luttinger liquid hamiltonian .
in contrast to the spin ladder case@xcite , the exponents @xmath354 in the hamiltonian eq .
( [ eq : effective_gapless ] ) are _ non - universal_. this can be shown in the following way .
let us consider the case where @xmath355 and @xmath356 are small compared to @xmath357 , namely where the couplings in the two chains are nearly identical . in that case , we can neglect in a first approximation the terms @xmath358 and @xmath359 in the hamiltonian [ eq : rotated_mf ] .
we are thus left with two decoupled sine - gordon models under a magnetic field .
these sine - gordon models undergo a commensurate - incommensurate transition at a critical magnetic field .
it is well known@xcite that the exponent at the transition assumes a universal value that renders scaling dimension of the cosine term equal to one .
therefore , in our case , the universal value of the exponent at the transition is : @xmath125 . as a result , close to the transition
the effective hamiltonian is : @xmath360 +   \int \frac{dx}{2\pi } \left [ \frac { u_1 k_1
-u_2 k_2}{2 } \pi^2 \pi_s \pi_a +   \frac { u_1 /k_1 - u_2 /k_2 } 2 \partial_x \phi_s \partial_x \phi_a \right].\end{aligned}\ ] ] returning to @xmath133 and @xmath134 , one obtains an hamiltonian of the form eq .
( [ eq : effective_gapless ] ) with : @xmath361 indicating that except in the case of equivalent chains , one should not expect universal exponents at the transition .
this should be the case in particular for the spin - tube @xcite .
it may provide an experimental test for the spin - orbital model@xcite of @xmath0 since the exponent @xmath362 controls the temperature dependence of the nmr relaxation rate@xcite .
however , since the transition temperature between the gapped and the gapless phase in @xmath0 is @xmath363 , the magnetic field needed to close the gap should be of order @xmath364 which could make the experiment impossible . in @xmath365 , with @xmath366 , the situation is even worse .
the spin and pseudospin operators are , @xmath367 \nonumber \\
\tau^z(x)&=&-\frac{\partial_x \phi_2 } { \pi } + \frac{e^{\imath \frac{\pi x } a}}{\pi a } \sin 2\phi_2 , \mbox { } \mbox { } \;\;\;\;\;\;\;\;\;\;\;\;\;\ ; \tau^+(x)=\frac{e^{\imath \theta_2}}{\sqrt{\pi a } } \left [ e^{\imath { \pi x } a } + \sin 2\phi_2   \right],\end{aligned}\ ] ] where @xmath368 is the total magnetization . in the case of the spin - orbital model ,
the spin - spin correlation functions are therefore given by the usual formulas@xcite .
the situation is however more interesting in the case of the spin tube . using the formula ( [ eq : spin_stot_chirality ] ) , one has : @xmath369,\ ] ] where @xmath370 , and @xmath10 is the chain index .
explicitly , @xmath371    this correlation function , which enters in particular in the calculation of the nmr relaxation rate , contains power - law divergences at wave vector @xmath372 but also @xmath373 due to the fluctuations of chiralities .
we have presented a field - theoretical analysis of the low - energy physics of the anisotropic spin - orbital model and the three - leg ladder with periodic boundary conditions ( the spin tube ) in the strong interchain coupling limit .
we gave a derivation of the field theoretical model from the lattice hamiltonian and then analyzed the phase - diagram using renormalization group equations .
the system is found to exhibit a gapless phase , a spin - liquid phase or an ising antiferromagnetic phase depending on the microscopic couplings .
the spin liquid ground state is two - fold degenerate , formed either by singlets of spins on even bonds and singlets of orbital pseudospins ( chirality in the spin - tube case ) on odd bonds or the other way round .
the antiferromagnetic phase competes with the spin - liquid and we have discussed this competition briefly .
the spin liquid phase obtains in the case of the spin tube with @xmath63 symmetry , and we discussed the nature of excitations above this spin liquid ground state .
these excitations have spin @xmath374 and pseudospin @xmath375 .
they are formed by introducing a free spin as a defect in the spin singlet pattern as well as a pseudospin in the pseudospin singlet pattern .
these excitations lead to a kind of non haldane spin liquid analogous to the one discussed by nersesyan and tsvelik@xcite .
an interesting consequence is the absence of a magnon peak at @xmath376 in the spin - spin correlation functions of the spin tube
. this behavior could be tested in numerical simulations .
we have investigated the effect of an applied magnetic field @xmath316 on the dimerized phase .
a strong enough magnetic field @xmath302 causes the closure of the gap and the disappearance of dimer order .
the resulting gapless phase is a _ two _ component luttinger liquid in contrast to the one component luttinger liquid that is observed in spin ladders .
the exponents appeared to be _ non - universal _ at the transition point , in contrast with the spin ladder case@xcite
. it would be interesting to obtain numerically the luttinger liquid exponents for the spin tube or the anisotropic spin orbital model .
for the spin tube , we have also shown that new soft modes appeared in the spin - spin correlation functions above @xmath317 by comparison with the soft modes of the single chain .
this is the result of the presence of soft chirality modes .
this may be tested numerically .
e. o. acknowledges support from nsf under grant nsf - dmr 96 - 14999 .
we thank p. lecheminant and h. saleur for discussions and useful comments on the manuscript .
@xmath377 permanent address : cnrs - laboratoire de physique thorique de lecole normale suprieure 24 , rue lhomond 75231 paris cedex 05 france
we start from the hamiltonian : @xmath378 after the usual jordan wigner fermionization , this gives : @xmath379    the mean field approximation is obtained by taking : ( recall , @xmath380 is a chain index ) . @xmath381
the mean field equation are then : @xmath382 \nonumber \\ t=-\frac{2}{\pi}(1+\lambda t)\sqrt{1 + \left(\frac{\delta}{1+\lambda
t}\right)^2}\frac{\mathbf{e}\left(1-\left(\frac{\delta}{1+\lambda t}\right)^2\right)}{1-\left(\frac{\delta}{1+\lambda t}\right)^2}\end{aligned}\ ] ] where @xmath383 are the complete elliptic integrals of the first and second kind respectively . for small @xmath186 , the mean field equations ( [ eq : mean_field_symmetric ] ) reduce to : @xmath384    and one obtains : @xmath385    this result from the mean field theory agrees with the prediction of bosonization since bosonization would predict the same essential singularity in @xmath186 at small @xmath78 , the interchain coupling being marginal .      in this section , we consider a problem with a hamiltonian of the form : @xmath386 after a jordan wigner transformation , the hamiltonian becomes : @xmath387    the mean field approximation is now : @xmath388    the mean field equations are now : @xmath389 \nonumber \\
t_2=-\frac{2}{\pi}(1+\alpha t_2)\sqrt{1 + \left(\frac{\delta_1}{t_1}\right)^2}\frac{\mathbf{e}\left(1-\left(\frac{\delta_1}{t_1}\right)^2\right)}{1-\left(\frac{\delta_1}{t_1}\right)^2 } \nonumber \\
\delta_2=-\frac{2\alpha \delta_1}{\pi}\sqrt{1 + \left(\frac{\delta_1 } { t_1}\right)^2}\left [ 2\mathbf{k}\left(1-\left(\frac{\delta_1}{t_1}\right)^2\right)-\frac{\mathbf{e}\left(1-\left(\frac{\delta_1}{t_1}\right)^2\right)}{1-\left(\frac{\delta_1}{t_1}\right)^2 } \right ] \\\end{aligned}\ ] ]    for small @xmath333 , the mean field equations for @xmath390 and @xmath205 reduce to : @xmath391 one sees that a finite bandwidth @xmath392 is produced for the @xmath27 spin waves at least for small @xmath333 .
the gap equations admit solutions at small @xmath186 .
one expects at weak coupling an essential singularity in @xmath393 .
therefore , @xmath394 should go to zero for @xmath395 .
this implies that one can write : @xmath396 one can solve the mean field equations for @xmath205 and @xmath206 .
one obtains : @xmath397    we obtain therefore self - consistently a gap much smaller than the smallest bandwidth . as a consequence ,
the correlation length is much larger than the lattice spacing , which justifies a continuum approximation .
the bosonization treatment is therefore valid for @xmath395 and close to the xy limit .
in this section , we derive rge for the quantum sine gordon model using the method of knops and den ouden @xcite .
the quantum sine gordon model is defined by the following lattice hamiltonian : @xmath398 - \frac{2g}{(2\pi a)^2 }   \cos 4\phi_i \right\}\end{aligned}\ ] ] the cutoff in this model is only on space and not on time .
the euclidean action of the quantum sine gordon model is obtained as : @xmath399=\int_{-\infty}^{\infty } \frac{d\omega}{2\pi } \int_{-\frac \pi a}^{\frac \pi a } \frac{dk}{\pi } \frac{u}{2\pi k } \left(k^2+\frac{\omega^2}{u^2 } \right ) |\phi(k,\omega)|^2 -\frac{2g}{(2\pi a)^2 } \int dx d\tau \cos 4\phi\end{aligned}\ ] ]    where @xmath400
. one can integrate from @xmath401 to @xmath402 using a cutoff function @xmath403 .
this action breaks rotation invariance in the @xmath404 space in contrast with the action of the classical sine - gordon model @xcite . instead of the sharp cutoff , one can use any cutoff function @xmath405 such that @xmath406 and @xmath407 .
an example is @xmath408 .        where : @xmath410}\frac { d \theta}{\pi } \nonumber
\\ f_2(\rho)=\int_0^{2\pi}\cos^2 \theta e^{16 k \left[v(\rho,\theta)+ \frac 1 2 \ln \rho\right]}\frac { d \theta}{\pi } \nonumber \\\end{aligned}\ ] ] and : @xmath411      the eqs .
( [ eq : rge_momentum_cutoff ] ) have to be contrasted with the usual rg equations for the sine gordon model with a cutoff isotropic in the @xmath414 space . in the latter case , the equations are such that @xmath415 , and only @xmath284 is flowing under renormalization . in the case
we have considered , both @xmath416 and @xmath284 are flowing under rg transformation .
however , for @xmath244 , the rg equations ( [ eq : rge_momentum_cutoff ] ) are considerably simplified . since @xmath417 and @xmath418 , one has the following simplified rg equations : @xmath419    these rg equations are identical to the ones obtained with a cutoff isotropic in @xmath414 space .
we conclude that in the vicinity of the bkt transition point , the anisotropy between space and time does not matter .                        where : @xmath429 - 2(k_1-k_2)\ln \rho } \nonumber \\ f_8(\rho)=\int_0^{2\pi } \frac{d\theta}{2\pi } e^{4\left[k_1 v(\rho,\theta ) - k_2 v(\rho,\theta,\frac{u_2}{u_1})\right]-2(k_2-k_1)\ln \rho } \nonumber \\
f_9(\rho)=\int_0^\infty e^{8k_1 \left[v(\rho,\theta)+\frac 1 2 \ln \rho\right]}\end{aligned}\ ] ]                                                    a. furusaki and s .- c .
zhang , dynamical spin correlations in the heisenberg ladder under a magnetic field and correlation functions in the so(5 ) ladder , to be published in phys .
60 , number 1 , 1999 , cond - mat/9807375 .