the effect of a magnetic field on an s=1/2 antiferromagnetic chain has been extensively investigated theoretically over many years .
the hamiltonian is written : @xmath0.\label{ham}\ ] ] an important conclusion was that the groundstate remains gapless right up to the saturation field .
the low - energy excitations can be described by bosonization which predicts gapless excitations at wave - vectors 0 and @xmath1 and also at the incommensurate wave - vectors @xmath2 where @xmath3 is the magnetization per site , @xmath4 .
the first detailed experimental study of such systems at large fields with @xmath5 of @xmath6 were only performed very recently,@xcite on cu benzoate .
this material has a relatively small exchange energy , @xmath7 mev so that @xmath8 for a field of 7 t. while these experiments verified , in detail , the expected field - dependent shift of the wave - vector at which a gap minimum occurs , they also discovered an unexpected result .
a non - zero gap appeared which seemed to scale as approximately @xmath9 , with strong dependence on the field orientation .
dender et al.@xcite suggested that this gap might arise from the staggered @xmath10 ( gyromagnetic)-tensor , associated with the low symmetry of the crystal structure and the presence of 2 crystallographically inequivalent cu sites on each chain .
thus the last term in eq .
( [ ham ] ) must be replaced by : @xmath11s^b_j.\ ] ] this results in the presence of an effective staggered field , @xmath12 , upon the application of a uniform field .
such a staggered field , which couples directly to the nel order parameter , is expected to produce an ordered antiferromagnetic moment and a gap which scale with field .
this idea was developed in detail in ref .
[ ] where it was found that a staggered dzyaloshinskii - moriya@xcite ( dm ) interaction also contributes a roughly equal amount to this effective staggered field
. this corresponds to an additional term in the hamiltonian : @xmath13    it was found that several aspects of the experiments could be explained in detail by this model .
these include the field orientation dependence of the gap and its scaling with field magnitude .
much of this work used the bosonization technique which maps the problem onto the sine - gordon model , for which various exact results are available .
the excitations observed in neutron scattering were identified with the soliton , antisoliton and `` breather '' ( soliton - antisoliton boundstate ) spectrum of the sine - gordon model .
additional results , further supporting this approach , were obtained by essler and tsvelik.@xcite the purpose of this paper is to provide more details and some extensions of the results in [ ] . while cu benzoate is the only example of such a system that we discuss in the present paper , it should be possible to apply our theory to other quasi - one dimensional system with similar crystallographic structure .    in the absence of a staggered field ,
the critical behaviour of the antiferromagnet is determined by 3 field dependent quantities : the magnetization , @xmath3 [ which determines the soft wave - vectors via eq .
( [ softk ] ) ] , the spin - wave velocity , @xmath14 and the boson compactification radius ( which determines the critical exponents ) , @xmath15 .
all three of these quantities can be determined very accurately by numerical solution of bethe ansatz equations .
furthermore , we derive an exact relationship between these three functions using field theory arguments . we also derive the logarithmic dependence of @xmath16 on @xmath17 , as @xmath18 , from the existence of a marginally irrelevant operator , using the renormalization group .
we determine the scaling of gap with field , order parameter with field and susceptibility with temperature .
the calculations are done including logarithmic corrections which arise from a marginally irrelevant operator and take into account both uniform and staggered fields .
furthermore , the exact amplitudes of the scaling functions are determined using a recent result of lukyanov and zamalodchikov.@xcite ( after this calculation was finished we received the preprint [ ] which gives the same result for the gap , without a discussion of logarithmic corrections . )
we also give some further discussion of the structure factors @xmath19 , measured in neutron scattering .
we discuss a hidden su(2 ) symmetry of the model .
we prove that the longitudinal structure factor ( for  @xmath20 corresponding to the uniform field direction ) gets contributions from only the soliton and anti - soliton intermediate states , in agreement with experiment . on the other hand
, the transverse structure function gets contributions only from the breathers . using the approximate @xmath21 symmetry
, we discuss the relative intensity of the various single - particle peaks in the neutron - scattering cross - section , taking into account the polarization dependent factors which arise from fourier transforming the dipole interaction between neutrons and spins which were omitted in ref .
[ ] . a comparison is made with experimental results .
in particular , the problem of determining a consistent value for the dm vector , @xmath22 is discussed .
the susceptibility of the sine - gordon model is calculated , using the integrability of the model , giving essentially the field and temperature dependence of the staggered susceptibility .    in section
ii we discuss the dm interaction and the mapping of the system into a heisenberg model with orthogonal uniform and staggered fields . in section iii
we treat this problem using conventional spin - wave theory . in section
iv we discuss bosonization in the presence of a uniform magnetic field . in section v we extend the bosonization approach to the case with staggered field and analyze the induced gap . in section vi we discuss structure factors and compare with the observed neutron scattering cross - section . in section vii we present estimates of the dm interaction based on several experimental results . in section viii
we discuss the magnetization and susceptibility .
the crystal structure of cu benzoate is shown in fig .
( [ fig : crystal ] ) and fig .
( [ fig : crystallarge ] ) .
the chain direction is the c - axis .
note that each cu atom is surrounded by 6 ligands with a local symmetry which is almost tetragonal .
however the principal axes for this tetragonal symmetry alternate along the chain , with the two inequivalent c - axes being rotated by @xmath23 relative to each other .
these corresponding c - axes are labelled i and ii in fig .
( [ fig : crystal2 ] ) .
the b - axes are the same for both cu sites .
neither of these sets of principal axes correspond to the crystal axes .
it is expected that the principal axes for the gyromagnetic tensor will also alternate , corresponding to the local tetragonal axes around each cu ion .
the principal axis for the anisotropic exchange interaction is expected to be the @xmath24 axes , the perpendicular bisector of i and ii axes . on the other hand , the principal axis for the dipole interaction , which is of roughly the same order of magnitude is essentially the c - axis .
combining these two types of contributions to the nearest neighbor spin - spin interaction , gives a principal axis which roughly bisects @xmath24 and @xmath25 , and is denoted @xmath26 in fig .
( [ fig : crystal3 ] ) .
it is convenient to refer the g - tensor to this @xmath27-@xmath28-@xmath26 coordinate system . from electron spin resonance ( esr )
measurements,@xcite it takes the form : @xmath29 with the @xmath30 referring to the 2 inequivalent cu sites . @xmath31 and @xmath32 are the uniform and staggered parts of the g - tensor .
this staggered g - tensor produces an effective staggered field , @xmath33 , while the uniform g - tensor produces an effective uniform field @xmath34 . in the special cases where the applied field is along the b axis or in the @xmath35 plane the effective staggered field is perpendicular to the applied field and also to the effective uniform field . for general directions of the applied field they are almost perpendicular ( to within a few % ) .
as discussed by dzyaloshinskii@xcite and moriya@xcite , in magnetic crystals of low symmetry an additional antisymmetric exchange interaction occurs , the dm interaction : @xmath36 the possible values of the dm vector @xmath37 can be limited by considering crystal symmetries of cu benzoate .
firstly , the compound is invariant under a translation along the c - axis by two sites .
this means the dm vectors are the same among the even ( or odd ) links , but even and odd dm vectors can be different .
secondly , there is a symmetry under rotation by angle @xmath1 about an axis parallel to the b - axis that passes through the mid - point of two neighboring sites ( @xmath38 and @xmath39 ) , along the chains ( c - axis ) .
as noticed by moriya,@xcite this implies the dm vector for the interaction between @xmath38 and @xmath39 must be orthogonal to the b - axis .
this can be shown as follows : assume we have the dm interaction @xmath40 .
now apply the rotation described above .
it acts on the spin operators as @xmath41 and @xmath42 .
thus , for the @xmath28-component of @xmath43 , the dm interaction would be inverted while it is unchanged for the @xmath44-component of @xmath43 implying that @xmath45 . finally , the crystal structure is invariant under the combined operation of one site translation along the chain ( c ) direction _ and _ reflection in the ac - plane .
considering the fact that the spin vector @xmath46 is an axial vector , the operation acts as @xmath47 and @xmath48 .
since the dm vector is orthogonal to the @xmath28-axis ( and thus one factor of @xmath49 always appears in the outer product ) , the dm interaction term is inverted by the combined operation : @xmath50 .
thus , the dm vector is alternating as in eq .
( [ hdm ] ) .
there are apparently no other restrictions that can be placed on the dm vector using symmetry alone . by considering a tight - binding model for the exchange interactions it was estimated@xcite that @xmath51 is of o(@xmath52 where @xmath53 is the deviation of g from twice the identity matrix .
apart from the antisymmetric dm interaction , the remaining exchange anisotropy is believed to be quite negligible ( about 1% of j ) and we will henceforth ignore it .
taking @xmath54 , we may write the hamiltonian : @xmath55 \nonumber \\ & &      + j \sum_j [ s^z_{2j-1 } s^z_{2j}+s^z_{2j}s^z_{2j+1}],\end{aligned}\ ] ] where @xmath56 . performing a rotation@xcite of the spins by an angle @xmath57 : @xmath58 where @xmath59 the hamiltonian
is transformed to the standard xxz model : @xmath60.\ ] ] with some assumptions this anisotropic exchange may cancel the pre - existing one . in any event , it is small and we will ignore it .    now consider an external magnetic field , approximating the @xmath10-tensor as 2 times the identity matrix .
the spin redefinition of eq .
( [ redef ] ) introduces an effective staggered field .
for example , for a uniform field in the x - direction : @xmath61.\ ] ] combining the actual form of the g - tensor with the dm interaction , we can obtain effective uniform and staggered fields corresponding to an arbitrary applied one . writing the rotation matrices by @xmath57 around @xmath22 as : @xmath62
the effective uniform and staggered fields are defined by : @xmath63\vec h \nonumber \\
h^s & \equiv & [ { \cal r}^s g^u + { \cal r}^ug^s]\vec h. \label{huaeff}\end{aligned}\ ] ] in general @xmath64 is nearly orthogonal to @xmath65 . for small @xmath32 and @xmath51
, the staggered field can be approximated as @xmath66 namely the sum of two contributions .
henceforth , since we are ignoring the small residual exchange anisotropy and assuming that @xmath67 , we will take @xmath65 to be in the z - direction and @xmath64 to be in the x - direction and refer to them as simply @xmath17 and @xmath68 respectively . also setting @xmath69 we arrive at the simple effective hamiltonian : @xmath70 , \label{heff}\ ] ] with @xmath71 .
[ note that we have switched the directions of the uniform and applied fields relative to our earlier paper@xcite which , unfortunately , contained some inconsistencies of notation . ]
in this section we summarize the results of spin - wave theory ( leading order 1/s expansion ) for the effective hamiltonian of eq .
( [ heff ] ) .
[ as far as we know , spin - wave theory results for this problem were first published in ref .
although this misses certain features caused by quantum fluctuations in one dimension , it is still quite instructive .
the classical groundstate is a canted antiferromagnetic structure , shown in figure ( [ fig : class ] ) .
the spins on both sub - lattices lie in the @xmath72 plane canted towards the @xmath73-axis by an angle @xmath74 from the @xmath30 x - axis .
the classical energy of this state is : @xmath75 this is minimized for @xmath76 the solution of : @xmath77 for @xmath78 , the solution is : @xmath79 in order to do a systematic 1/s expansion , it is convenient to regard @xmath17 and @xmath68 as being of o(s ) .
we assume that @xmath17 is less than the saturation field , @xmath80 .
the leading order spin - wave expansion for a spin pointing in the x - direction is : @xmath81\ ] ] here @xmath82 is a boson annihilation operator . to consider small fluctuations about the canted structure we simply write :
@xmath83 where @xmath84 is a rotation about the y - axis by @xmath85 : @xmath86 and @xmath87 is a rotation by @xmath1 about the z - axis : @xmath88 using the facts that : @xmath89 the hamiltonian may be written in a manifestly translationally invariant way : @xmath90 .
\label{hsw}\ ] ] it is the translational invariance of eq .
( [ hsw ] ) which motivated the somewhat peculiar looking choice of transformation matrices in eq .
( [ swt ] ) . substituting eq .
( [ swt ] ) into the hamiltonian of eq .
( [ hsw ] ) , we find that the term of @xmath91 is a c - number and the term of @xmath92 vanishes .
the term of @xmath93 is : @xmath94.\end{aligned}\ ] ] fourier transforming and performing a bogoliubov transformation that mixes @xmath95 with @xmath96 , we obtain a single band of spin waves in the paramagnetic brillouin zone , @xmath97 with dispersion relation : @xmath98 ^ 2\nonumber \\ & & -[js(1+\cos 2\theta ) \cos k]^2\bigr\}^{1/2}.\end{aligned}\ ] ] the above transformation has allowed us to obtain a single band in the paramagnetic brillouin zone
. we may equivalently fold the dispersion relation into the antiferromagnetic brillouin zone , @xmath99 .
this gives us 2 branches of spin - waves with dispersion relations : @xmath100 ^ 2\nonumber \\ & & -[js(1+\cos 2\theta ) \cos k]^2\bigr\}^{1/2}.\end{aligned}\ ] ] while this is a fairly simple and explicit formula for the energies in terms of @xmath76 , it must be borne in mind that @xmath76 is determined in terms of @xmath101 , @xmath17 and @xmath68 by eq . ( [ theta ] ) . in the special case @xmath78 , using eq . ( [ thetah0 ] ) ,
we obtain : @xmath102^{1/2}.\ ] ] note that at @xmath103 , or equivalently @xmath104 , @xmath105 the @xmath106 mode is the goldstone mode corresponding to a uniform precession about the z - axis .
a non - zero @xmath68 gives this mode a gap , pinning the spins along the x - axis .
we may calculate this gap , to lowest order in @xmath68 using , from eq .
( [ theta ] ) @xmath107 where @xmath108 to linear order the gap is given by @xmath109 thus : @xmath110}[1-(h/4js)^2]^{1/4}+o(h^{3/2}).\ ] ] note that @xmath111 is a singular function of the staggered field , exhibiting a mean field exponent of 1/2 . on the other hand
, it depends only weakly on the uniform field , being almost independent of @xmath17 up to @xmath17 of @xmath112 .
while this mean field exponent changes when one - dimensional quantum fluctuations are taken into account , it is reasonable to expect this weak dependence on @xmath17 to remain true .
the upper mode , @xmath113 depends strongly on h but only weakly on @xmath68 . for @xmath71 @xmath114 .
in the one - dimensional case , an exact picture of the low - energy behavior can be obtained using bosonization and rg arguments . here
we summarize the results for the case of a uniform field , but no staggered field .
we begin with the case where the uniform field also vanishes .
the low energy degrees of freedom of the quantum spin variables can be represented in terms of a free boson with lagrangian density : @xmath115 .
\label{lag}\ ] ] [ here @xmath116 is the spin - wave velocity which we will generally set equal to 1 . ] the boson field , @xmath117 can be separated into left and right moving terms : @xmath118 their difference defines the dual field : @xmath119 for @xmath120 there are low energy degrees of freedom at wave - vectors 0 and @xmath1 .
the spin operators can be approximated as : @xmath121 + c(-1)^je^{i2\pi r\tilde \phi}. \label{abbos}\end{aligned}\ ] ] ( the first constant above is universal .
the next three ( real ) constants are not , but @xmath122 has been recently determined using the integrability of the model@xcite and will be discussed in the next section . )    for the h=0 heisenberg model , @xmath123 . for the xxz model @xmath16 ( and @xmath122 )
vary with the anisotropy parameter . writing the hamiltonian : @xmath124,\label{hxxz}\ ] ] exact bethe ansatz results
determine : @xmath125 @xmath16 varies between @xmath126 and 0 along the xxz critical line , @xmath127 .    in order to understand the vicinity of the isotropic antiferromagnetic point , @xmath128
, it is convenient to use non - abelian bosonization : @xmath129 here @xmath10 is the su(2 ) matrix field of the wess - zumino - witten , k=1 non - linear @xmath130-model ( wzw model ) . @xmath131 and @xmath132 are the left and right moving conserved currents associated with the su(2 ) symmetry of the spin chain . by comparing eqs .
( [ abbos ] ) and ( [ nonabbos ] ) we may read off the correspondences between the wzw fields and the free boson fields .
these are : @xmath133 here @xmath134 .
the bosonized spin - chain lagrangian contains , in addition to the free boson lagrangian , interaction terms : @xmath135 .
\label{marg}\ ] ] for the isotropic heisenberg model , @xmath136 .
including small anisotropy , @xmath137 .
these obey the rg equations : @xmath138 the rg trajectories are hyperbolas as shown in fig .
( [ fig : kt1 ] ) . for @xmath139
, they end at the @xmath140 axis , corresponding to the xxz critical line .
for @xmath141 they lead towards strong coupling , corresponding to the easy axis ordered phase . reverting to abelian bosonization
, we see that the fixed point lagrangian contains the extra term : @xmath142.\ ] ] since this is proportional to the free lagrangian , we can eliminate it by a rescaling of @xmath117 .
this corresponds to a rescaling of the parameter , @xmath16 , decreasing it by an amount of @xmath143 in agreement with the bethe ansatz result of eq .
( [ delta ] ) .
we now consider the heisenberg model with a uniform external field ( but no staggered field ) .
the extra term in the lagrangian becomes , upon bosonization , @xmath144 this term can be eliminated by a redefinition of the boson field : @xmath145 this leaves the free lagrangian unchanged .
however , it does effect the interaction term and the bosonization formulae .
the interaction term is changed due to the shift of the @xmath30 components of the currents : @xmath146 ( note that the phases add , rather than cancel in the interaction term @xmath147 . ) the effect of these phases on the rg equations can be determined from a consideration of the operator product expansion ( ope ) .
one of the ope s gets shifted while the other one does not : @xmath148 the one - loop rg equations can be conveniently derived using an ultraviolet cut off on the distance between any pairs of insertions of the interaction lagrangian in perturbation theory .
( see for example [ ] . ) as the position space ultraviolet cut off is increased from @xmath20 to @xmath149 , we integrate over that range of separation of the two points , using the ope .
the net effect is that , when the cut - off is small compared to @xmath150 the phase factor in the ope is nearly constant and can be ignored .
however , when the cut - off is large compared to @xmath150 the phase factor produces rapid oscillations which tend to cancel out the term from the effective renormalization . reverting to an energy cut - off , this means that the rg equations of eq .
( [ rge ] ) are essentially correct for @xmath151 , but for @xmath152 the right hand side should be replaced by 0 in the first equation .
that is , @xmath140 ceases renormalizing , at @xmath153 of o(h ) whereas @xmath154 continues to renormalize as before .
thus the rg trajectories are essentially 2 straight lines .
for @xmath151 the couplings renormalize along the isotropic separatrix , @xmath155 , but for @xmath152 @xmath140 is constant and @xmath154 renormalizes to 0 .
see fig .
( [ fig : kt2 ] ) .
thus , in order to determine @xmath156 we simply need to calculate @xmath157 using the isotropic rg equations of eq .
( [ rge ] ) . for @xmath158 ,
these give : @xmath159 this argument may seem rather naive and it certainly does not give the correct rg trajectory for e of o(h ) .
however , due to the weak , logarithmic , dependence of @xmath160 on e we expect that this argument gives the correct behaviour of @xmath156 with h for @xmath158 .
this argument then determines the dependence of r on h for @xmath158 : @xmath161 precisely this result was obtained from the bethe ansatz  @xcite in the limit @xmath158 .
( [ fig : r(h ) ] ) . for larger values of h higher order terms in the rg equations would be needed and additional interactions would have to be considered , some of which break lorentz invariance .
the net effect is that the h - dependence of r becomes more complicated at larger h and also the spin - wave velocity also change with h. we have so far set it equal to @xmath162 .
it is known to have the value @xmath163 from the bethe ansatz , for the heisenberg model at @xmath120 .
it can be determined numerically from the bethe ansatz integral equations  @xcite , and is given in fig .
( [ fig : v ] ) .    apart from the functions @xmath14 and
@xmath15 we will also be interested in the behavior of the magnetization , @xmath164 from bosonization we obtain : @xmath165 as @xmath18 .
one way of obtaining this result is from calculating the zero field susceptibility : @xmath166    bosonization leads to an exact relation between the three functions @xmath14 , @xmath15 and @xmath3 .
this follows from calculating the susceptibility at arbitrary field , @xmath17 , using eq .
( [ abbos ] ) with the corresponding value of @xmath15 : @xmath167 in particular , in the limit of small @xmath17 this predicts a logarithmic correction to the magnetization : @xmath168 . \label{mh}\ ] ] eqs .
( [ rh ] ) , ( [ mrv ] ) and ( [ mh ] ) are universal
. they should remain true for generic half - integer spin isotropic antiferromagnets in the gapless phase . for the particular case of the nearest neighbour s=1/2 heisenberg model we can set @xmath169 in eq .
( [ mh ] ) .
( [ mrv ] ) agrees very well with our numerical solution of the bethe ansatz equations .
[ see figs .
( [ fig : r(h ) ] ) , ( [ fig : v ] ) , ( [ fig : mag ] ) .
] as far as we know the first ( numerical ) calculation of these quantities for @xmath170 were ref .
[ ] for @xmath3 , refs .
[ , ] for @xmath15 and ref .
[ ] for @xmath14 .
( numerical calculations in refs .
[ , ] were based on the bethe ansatz integral equations of ref .
once @xmath171 , @xmath172 and @xmath16 are determined from the bethe ansatz , all low energy properties of the system are determined by bosonization . from eq .
( [ abbos ] ) the magnetization can be written : @xmath173 therefore the exact formula for the field - induced shift in @xmath117 is : @xmath174 all low energy green s functions are then determined from eq .
( [ abbos ] ) after shifting @xmath175 .
in particular we see that @xmath176 has the soft wave - vectors @xmath177 and @xmath178 whereas @xmath179 has the soft wave - vectors @xmath180 and @xmath1 .
now we consider the effective hamiltonian of eq .
( [ heff ] ) with both uniform and staggered fields non - zero and @xmath71 .
we begin by using the results of the previous section to obtain the @xmath78 theory with the shifted and rescaled boson field characterized by @xmath15 , @xmath3 and @xmath14 . for @xmath78 , upon making these tranformations , the lagrangian density is simply the free boson one of eq .
( [ lag ] ) . from the bosonization formulae of eq .
( [ abbos ] ) the staggered field adds the interaction term : @xmath181 noting the duality transformation between @xmath117 and @xmath182 : @xmath183 we may also write the free lagrangian in terms of @xmath182 : @xmath184.\ ] ] hence we have the standard sine - gordon field theory .
an impressive array of conjectured exact results are available on this model , which can be brought to bear on the spin - chain problem .
the interaction term is sometimes written as @xmath185 , so we see that we have @xmath186 . in the next section
we will discuss details of the excitation spectrum of this model . in this section
we discuss the dependence of the gap on the uniform and staggered field .
the renormalization group scaling dimension of this operator is @xmath187 from which it follows that the gap scales as : @xmath188 for some function @xmath189 .
since , for small uniform fields , @xmath190 , it follows that the exponent is approximately 2/3 , as found in the experiment . note that this formula is valid for @xmath191 at fixed @xmath17 .
this is a reasonable order of limits for describing the experiments since @xmath68 is only a few per cent of @xmath17 .
the exponent is determined by @xmath15 , given in figure ( [ fig : r(h ) ] ) .
we wish to improve on this result in two ways .
first of all , if we consider the case where @xmath17 is strictly 0 , then this formula is modified to : @xmath192 we will calculate exactly the amplitude @xmath193 . this is not of direct relevance to the experiments however since they are in the opposite limit @xmath194 .
more importantly , we can determine the amplitude function @xmath195 in eq .
( [ delta ] ) for @xmath158 : @xmath196 we will also determine exactly the numerical factor , @xmath197 , which is different than @xmath193 in eq .
( [ gaph0 ] ) . while the logarithmic factors in eqs .
( [ gaph0 ] ) and ( [ ah ] ) are universal and follow from an rg treatment of the marginal interaction , the exact numerical coefficients are specific to the ordinary nearest neighbour heisenberg s=1/2 model and are obtained by using a remarkable exact conjecture made recently by lukyanov and zamalodchikov@xcite extended to the heisenberg point following the method in ref .
[ see also ref .
our calculations follow the notation of refs .
[ ] and [ ] .
we consider the renormalization group equation obeyed by the effective coupling constant @xmath198 multiplying the @xmath199 interaction in eq .
( [ lsg ] ) , with bare value @xmath200 for a bare cut off @xmath101 , taking into acount , to linear order , the effect of the marginal interactions of eq .
( [ marg ] ) .
this is : @xmath201 g .
\label{rgg}\ ] ] the anomalous dimension is given , to low order , by : @xmath202 @xmath140 obeys the rg equation of eq .
( [ rge ] ) . by the usual scaling arguments
, we determined the gap , @xmath111 by reducing the ultraviolet cut - off down to a scale @xmath111 such that @xmath203 is o(1 ) .
taking into account the dependence of the effective coupling constant , @xmath198 on the bare coupling constant , @xmath204 , then determines the dependence of the gap on @xmath68 . in the case
@xmath120 , the rg flow of @xmath205 is given by eq .
( [ rge ] ) for all @xmath153 . on the other hand , for finite @xmath17
, @xmath205 essentially stops renormalizing at a scale of order @xmath17 . integrating eq .
( [ rgg ] ) , gives : @xmath206 this integral can be conveniently evaluated by changing integration variables to @xmath207 using the second of eqs .
( [ rge ] ) .
this gives : @xmath208^{1/4}.\ ] ] determining the gap by the condition @xmath209 and setting @xmath210 , gives : @xmath211^{1/4 } , \label{hlambda}\ ] ] for some non - universal constant @xmath212 , of o(1 ) .    in the case
@xmath120 , the solution of the rg equation eq . ( [ rge ] )
is : @xmath213 thus , @xmath214^{1/6}[1+o(1/\ln ( j / h ) ) ] .
\label{deltah0}\ ] ]    for non - zero @xmath17 , the first of eq .
( [ rge ] ) is only valid for @xmath151 . at lower energies
, @xmath140 stops renormalizing .
its fixed value at low e determines @xmath15 : @xmath215 we can extend somewhat the accuracy of our results to larger @xmath17 , by expressing the subsequent results in terms of @xmath15 , determined numerically from the bethe ansatz , rather than by using the above asymptotic small @xmath17 result for @xmath15 . for @xmath216 ,
we use the second of eq .
( [ rge ] ) with @xmath140 fixed at @xmath156 as given by eq .
( [ lambdaz0 ] ) and the initial condition @xmath217 to obtain : @xmath218(\frac{e}{h})^{2[1 - 2\pi r(h)^2]}. \label{lambdar}\ ] ] setting @xmath219 and substituting into eq .
( [ hlambda ] ) gives : @xmath220^{1/4}.\ ] ] thus we obtain eq .
( [ delta ] ) with @xmath221^{1/4}\right\}^{-1/(2-\pi r^2)}.\label{ah2}\ ] ] note that in eq .
( [ delta ] ) and ( [ ah2 ] ) @xmath16 is a function of @xmath17 , shown in fig .
( [ fig : r(h ) ] ) and given approximately by eq .
( [ lambdaz0 ] ) . as @xmath18
, we may evaluate @xmath222 explicitly , from eq .
( [ lambdaz0 ] ) . using : @xmath223
this gives : @xmath224 note that the same numerical constant , @xmath212 , occurs in the @xmath120 case , eq .
( [ deltah0 ] ) and the @xmath225 case , eq .
( [ delta ] ) and ( [ ah0 ] ) .
however , in the latter case it gets multiplied by an extra factor of @xmath226 .
finally we wish to determine the dimensionless amplitude , @xmath212 , appearing in eq .
( [ hlambda ] ) and below .
this can be done using two remarkable results .
one of them is the exact proportionality constant in the bosonization formula for the staggered part of @xmath227 , that is the constant c in eq .
( [ abbos ] ) .
this fixes the coupling constant in the sine - gordon model , @xmath228 , in eq .
( [ lsg ] ) .
the other recent result is the exact relationship between the sine - gordon coupling constant and the mass of the lightest particle in the spectrum of the sine - gordon model , which is the gap , @xmath111 .
this determines the exact relationship between @xmath111 and @xmath68 .
this calculation was done in zero uniform field for the xxz s=1/2 antiferromagnet of eq .
( [ hxxz ] ) .
the calculation was performed for all @xmath229 along the critical line , @xmath230 . due to the logarithmic corrections at the isotropic point , @xmath128 , an additional calculation
is needed at that point .
this can be done using the rg .
we essentially just need to apply eq .
( [ hlambda ] ) to the case of @xmath120 but @xmath229 slightly less than 1 . comparing to the exact result for all
@xmath231 then determines the coefficient , @xmath212 .
for @xmath231 and @xmath120 we use the rg equations of eq .
( [ rge ] ) at low energies .
the rg flows are hyperbolas terminating on the positive @xmath140 axis .
these flows are conveniently labelled by : @xmath232= 2 \frac{\cos^{-1}{\delta}}{\pi}. \label{epsilon}\ ] ] the solution of eq .
( [ rge ] ) with this initial condition is : @xmath233}. \label{lambdaep}\ ] ] substituting into eq .
( [ hlambda ] ) gives : @xmath234}\right\}^{1/4}.\ ] ] note that in the isotropic limit , @xmath235 , we recover our previous logarithmic result of eq .
( [ gaph0 ] ) . on the other hand ,
taking @xmath236 with @xmath237 held fixed it gives : @xmath238 by comparing this to the exact result for arbitrary @xmath229 ( and hence @xmath237 ) we may extract the value of the amplitude , @xmath212 .    in ref .
[ ] the spin correlation function in the xxz antiferromagnet is shown to have the asymptotic behavior : @xmath239 with an exact expression determined for the amplitude @xmath240 . in the isotropic limit , @xmath241 : @xmath242 in ref .
[ ] the operator @xmath243 is normalized so : @xmath244\cos[2\pi r\tilde \phi ( r ) ] \rangle \to { 1\over 2}|r|^{-2\pi r^2},\ ] ] ( after accounting for a difference in normalization of the free boson lagrangian by a factor of @xmath245 ) .
this determines the exact proportionality constant in the bosonization formula of eq .
( [ abbos ] ) : @xmath246 . \label{prop}\ ] ] hence the coupling constant in the sine - gordon lagrangian is precisely @xmath247 .
the exact relationship between this coupling constant and the mass , @xmath111 , of the soliton of the sine - gordon model is : @xmath248^{2-\pi r^2}.\ ] ] here we have inserted , by dimensional analysis , the spin - wave velocity v. taking the isotropic limit on both sides of this equation , using @xmath249 , gives : @xmath250^{3/2}.\ ] ] we see that this is equivalent to the rg result of eq .
( [ hxxz ] ) with the amplitude determined to be : @xmath251^{3/2 } \label{b}\ ] ] from eq .
( [ deltah0 ] ) , for @xmath120 , the gap behaves as : @xmath252 with : @xmath253^{2/3 } { \gamma ( 1/6)\over \gamma ( 2/3)}\approx 1.77695.\label{a0ex}\ ] ] since we were not aware of the result [ ] at the time , in [ ] the behavior of the gap with staggered field ( for @xmath120 ) was estimated numerically by extrapolating lanczos results for lengths up to 22 sites .
a very good fit to eq .
( [ deltah02 ] ) was obtained with : @xmath254 considering the numerical difficulties related to logarithmic corrections this is remarkably good agreement with the exact result of eq .
( [ a0ex ] ) . in the experimentally relevant case ,
@xmath216 , the gap behaves as in eq .
( [ delta ] ) and , for small @xmath255 , the amplitude is given by eq .
( [ ah ] ) with , from eq .
( [ ah0 ] ) , @xmath256 thus our expression for the gap becomes : @xmath257^{1/6}(h
/ j)^{1/[2-\pi r(h)^2)]}. \label{deltaf1}\ ] ] for larger @xmath255 greater accuracy might be obtained by using eq .
( [ ah2 ] ) with @xmath212 given in eq .
( [ b ] ) . that is : @xmath258/2}[2(1 - 2\pi r(h)^2)]^{1/4 } \right\}^{-1/[2-\pi r(h)^2 ] } ( h / j)^{1/[2-\pi r(h)^2]},\ ] ] where @xmath15 is given , from the bethe ansatz , in fig .
( [ fig : r(h ) ] ) . inserting its asymptotic value at low @xmath17 , @xmath259 gives back eq .
( [ deltaf1 ] ) .
we note that , to actually fit the experimental data , we take @xmath260 for some constant of proportionality which depends on field direction but is generally of order a few % . thus the actual scaling of gap with field is not a pure power law .
for @xmath120 it is convenient to use non - abelian bosonization , eq .  ( [ nonabbos ] ) so that the interaction term is written : @xmath261 in this case the model has an su(2 ) symmetry . note
that when both uniform and staggered fields vanish and ignoring the marginal operator , the symmetry is actually su(2)@xmath262su(2 ) : @xmath263 these two independent su(2 ) s act on left and right - movers separately .
the ordinary su(2 ) symmetry of the spin chain is the diagonal subgroup with @xmath264 .
this symmetry is broken by the staggered field . however
, a different su(2 ) subgroup of the original su(2)@xmath262su(2 ) survives for which : @xmath265 we may redefine the field @xmath10 by @xmath266 in which case the interaction becomes @xmath267 , which has the diagonal su(2 ) symmetry .
in fact , this continuum limit interaction arises from a staggered heisenberg exchange interaction , as occurs in the spin - peierls problem .
the equivalence of the continuum limit of these two apparently very different problems , is a non - trivial consequence of the chiral symmetry which maps @xmath268 into @xmath269 .
if we assume that the only effect of the uniform field is to add a term to the hamiltonian : @xmath270 then the uniform field can be removed by the gauge transformation : @xmath271 or equivalently : @xmath272 [ here we have set @xmath273 . ]
this transforms the matrix field @xmath10 as : @xmath274 this gauge transformation leaves invariant the staggered field term tr@xmath275 .
thus the exact su(2 ) symmetry remains in this approximation .
however , additional irrelevant terms in the hamiltonian in the presence of a uniform field break the exact su(2 ) symmetry , as evidenced by the change in the parameter @xmath16 with field .
the su(2 ) symmetry is only present for @xmath123 .
the non - abelian bosonization formula for the staggered part of the spin operators ( for @xmath120 ) is : @xmath276 the magnetic field leads to the gauge transformation of eq .
( [ gauge ] ) .
the parts of the spin operators with wave - vectors near @xmath1 thus become : @xmath277 now making the transformation @xmath278 , this becomes : @xmath279 in general the spectrum of the sine - gordon theory consists of the soliton , anti - soliton and breathers ( soliton - anti - soliton boundstates).@xcite in the su(2 ) symmetric case , the excitation spectrum of the sine - gordon model with @xmath280 consists of a triplet composed of soliton , anti - soliton and lowest breather and a second breather , heavier by a factor of @xmath281 .
the degeneracy of the triplet is a result of the su(2 ) symmetry . due to the su(2 ) symmetry ,
the 3 elements of the triplet are produced by the 3 operators tr@xmath282 with equal intensity , whereas the singlet is produced by the operator tr@xmath10 .
thus setting @xmath120 ( but not @xmath68 ) the structure functions @xmath283 and @xmath284 would be equal .
note that , for @xmath120 , @xmath285 , creates the y - polarized member of the triplet .
this can be regarded as a linear combination of the soliton ( created by tr@xmath286 ) and the anti - soliton ( created by tr@xmath287 ) .
@xmath288 consists of 2 identical contributions from the @xmath289 and @xmath290 terms in eq .
( [ strans ] ) .
each contributes exactly @xmath291 .
the effect of @xmath17 is to split the soliton and anti - soliton contributions to @xmath284 into two separate contributions at different wave - vectors @xmath292 .
thus , ignoring the small @xmath21 symmetry breaking : @xmath293    it is also interesting to note that the staggered part of the energy density is given by : @xmath294 this operator couples to lattice displacements ( phonons ) and is used to describe raman scattering experiments . upon making the gauge transformation of eq .
( [ gauge ] ) and the redefinition of eq .
( [ redefg ] ) this becomes : @xmath295 thus this operator also creates the soliton and anti - soliton .
hence this theory predicts a single particle excitation observable in raman scattering at the same field - dependent wave - vector and frequency as the incommensurate mode observed in neutron scattering .    upon allowing for su(2 ) symmetry breaking the radius changes . after making the gauge transformation , a u(1 ) symmetry still survives , corresponding to shifting @xmath117 by a constant .
the triplet is now split , with the lowest breather having a different mass than the degenerate soliton antisoliton pair . since the operators @xmath296 have charge @xmath297 with respect to this @xmath298 symmetry
, we see that the soliton and anti - soliton are created by the @xmath299 fourier modes of @xmath300 respectively .
the breathers can be classified as even or odd with respect to the discrete symmetry @xmath301 .
the odd breathers are created by the @xmath302 component of @xmath303 and the even breathers by the @xmath302 component of @xmath304 .
it can be shown that even and odd breather alternate in the spectrum of the sine - gordon model .
furthermore , the number of breathers increases with decreasing @xmath16 .
a third breather drops below the soliton antisoliton ( @xmath305 ) continuum immediately as soon as @xmath16 decreases below the isotropic value , @xmath126 with another one dropping below the @xmath305 continuum each time @xmath306 passes through an integer .
the mass of the @xmath307 breather , expressed in terms of the soliton mass , @xmath308 , is : @xmath309 where : @xmath310 thus the odd - numbered breathers contribute single - particle poles to @xmath283 and the even - numbered ones to @xmath311 while the soliton and antisoliton contribute single - particle poles to @xmath284 .
in addition various multi - particle continua contribute to the three spectral functions .    for a field of 7 t.
we estimate @xmath312 .
there are 3 breathers at this point , with masses .79 m , 1.45 m and 1.87 m .
a resolution limited peak was observed at @xmath313 , at energy .22 mev .
we identify this with the soliton ( or anti - soliton ) contribution to @xmath284 ; hence m=.22mev . a resolution
limited peak is clearly observable in the neutron scattering data at @xmath302 and an energy of @xmath314mev=@xmath315 .
this agrees very well with the prediction for the first breather mass .
we may also test the su(2 ) prediction of eq .
( [ intrat ] ) .
the su(2 ) symmetry is broken by various small effects as exemplified by the fact that @xmath316 .
in particular , this implies that @xmath317 has a second peak , of very low intensity , corresponding to the third breather . ignoring these effects
we expect the intensity of the lowest breather peak in @xmath283 to be approximately twice the intensity of the soliton peak in @xmath284 .
however , before a comparison can be made with experiment it must be taken into account that the unpolarized neutron scattering cross - section contains an important direction dependence arising from the fourier transform of the dipole - dipole interaction between the neutron and the spins .
the cross - section can be written : @xmath318 where @xmath319 is a unit vector in the direction of @xmath320 and the function @xmath321 is slowly varying .
thus the soliton , even breathers and odd breathers are weighted by different factors @xmath322 , @xmath323 and @xmath324 , respectively
. we also have to consider the variation of @xmath321 in examining the relative intensity of solitons to breathers since they occur at different values of @xmath325 .
it must be recalled that @xmath326 here refers to the direction of the effective staggered magnetic field ( and @xmath327 the direction of the uniform field ) . in the neutron scattering experiments the field was along the b - axis .
we note that the @xmath27 axis is rotated by about @xmath328 from the crystal axis , @xmath20 .
[ we define the rotation angle in the @xmath329- ( @xmath35- ) plane so that @xmath25 axis is @xmath330 rotated from @xmath20-axis . ]    strictly speaking , we must take into account the effect of the redefinition of the spin operators , discussed in section ii , that was used to eliminate the dm interaction . letting @xmath331 to be the rotated spin operators , defined in eq .
( [ redef ] ) , and inverting the transformation , we may write the structure function for the original spin operators in terms of the structure function for the rotated operators , which we write as @xmath332 . in this way we obtain : @xmath333 @xmath334 remains 0 due to translation invariance and the @xmath335 are unaffected by the transformation . here
@xmath336 and @xmath337 refer to two axes orthogonal to @xmath22 ( _ not _ orthogonal to @xmath338 as in most of this paper . )
we expect the second terms in eqs .
( [ gredef ] ) to be negligible since @xmath339 is small .
furthermore , the @xmath340 are small for @xmath341 also making the second terms in eq .
( [ gredef ] ) smaller than the first for @xmath342 near @xmath1 .
henceforth we ignore this small correction and simply use @xmath343 .
the direction ( and magnitude ) of the effective staggered field depends on both the staggered part of the @xmath10-tensor and dm interaction . since the dm interaction in cu benzoate is unknown , the direction @xmath344 is not known at present .
( the dm interaction in cu benzoate may be estimated from various experimental results based on the present theory .
we will discuss this issue in the next section .
) however , the direction @xmath344 can be deduced from the polarization analysis@xcite of the neutron scattering experiment .
dender et al .
analysed the polarization of the neutron scattering at constant energy @xmath345 and various momentum transfers . for magnetic field
@xmath346 b and momentum transfer along the chain @xmath1 , this should probe the lowest ( @xmath347 ) breather . as discussed above , this odd breather is polarized orthogonal to the total effective staggered field . in ref .
[ ] it was claimed that the observed scattering is polarized in the @xmath27 direction .
also note that a misstatement of the crystal orientation occurred in ref .
[ ] so that the wrong sign appeared there for @xmath348.@xcite correcting this error , the polarization is @xmath349 from a - axis , or equivalently @xmath350 from @xmath27-axis .
this implies that , for an applied field in b - direction , the total effective staggered field @xmath351 ( @xmath352 ) from a - axis , or equivalently @xmath353 ( @xmath354 ) from @xmath27-axis .
[ there is actually another feature of the polarization analysis in ref .  [ ] which appears inconsistent with the theory presented here , namely the polarization analysis in zero field .
the strong dependence of the intensity on the a - component of the momentum is taken to indicate that @xmath355 .
our work ignores any anisotropy in the zero field limit and therefore predicts @xmath356 in that limit .
we do not understand the source of this discrepancy at present . ]
the constant-@xmath357 ( momentum transfer ) scan experiments , sensitive to the breather modes , were carried out at the fixed momentum transfer : @xmath358 where the lattice constants are @xmath359 , @xmath360 and @xmath361 .
here we have again corrected the misstatement@xcite of the crystal orientation in ref .
[ ] mentioned above .
note that the antiferromagnetic wave - vector , in the chain direction , is actually @xmath362 rather than the normal @xmath363 , because there are 2 cu atoms per unit cell along the c - axis . the crystal axes @xmath20 , @xmath28 and @xmath25 are essentially orthogonal .
thus , in the a - b - c system : @xmath364 @xmath319 is rotated @xmath365 from the @xmath20 axis , or @xmath366 from the @xmath27-axis .
we note that , this @xmath367 is almost parallel to the direction of the total effective staggered field estimated from the polarization analysis above .
the constant-@xmath357 scan for the soliton modes was done for a slightly different momentum transfer @xmath368 however , the direction of this momentum transfer @xmath367 is almost the same as the above and we will ignore the difference .    the fact that , an intense first breather peak is observed experimentally , supports the deduction that @xmath319 is nearly parallel to the total effective staggered field .
in fact , assuming that @xmath319 is completely parallel to the staggered field , the approximate su(2 ) prediction is that the lowest peak at @xmath104 ( from the first breather in @xmath283 ) should have twice the intensity of the lowest peak at @xmath369 ( from the soliton / antisoliton in @xmath284 ) .
this prediction is only a very approximate one due to the breaking of su(2 ) .
there could also be corrections from the function @xmath321 in eq .
( [ neut ] ) .
experimentally this ratio appears to be about @xmath370 .
this is perhaps satisfactory agreement .
this agreement is only worsened if @xmath367 deviates from the direction of the staggered field .
thus , analyses of polarization and of the scattering intensity of the lowest breather mode are consistent , and apparently lead to the conclusion that the total staggered field for h@xmath371b is almost parallel to the @xmath367 direction of eq .
( [ kdir ] ) used in the constant-@xmath357 scan .    on the other hand , in ref .
[ ] we discussed a feature in the experimental data at @xmath104 , and @xmath372 .
this is very close to the predicted mass of the second breather , which contributes to @xmath311 .
a recent calculation , based on integrability of the sine - gordon model , indicates that the relative intensity of the second breather in peak in @xmath311 should be roughly 1/2 of the intensity of the first breather in @xmath283 .
however , we expect this to have essentially zero intensity in the neutron scattering cross - section due to the factor of @xmath323 in eq .
( [ neut ] ) . [ in ref .
[ ] an apparently good agreement between theory and experiment was obtained because the factors of @xmath373 were not included . in ref .
[ ] intensities were not considered . ] a possible resolution of this disagreement is that this `` feature '' at @xmath374 , discussed by two groups of theorists , is just noise . as stated in ref .
[ ] `` given the quality of the data , this double - gap conjecture is highly speculative '' .
clearly more data is needed to determine whether or not there is really another sharp peak at this frequency .
an even more statistically insignificant feature in the data , at @xmath104 and @xmath375 , was discussed in ref .
[ ] where it was interpreted as the third breather peak in @xmath283 .
both the energy ( @xmath376 ) and the intensity ( very approximately 1/6 the intensity of the first breather ) agree with the predictions of the sine - gordon model .
note that in this case the factor of @xmath324 is common to first and third breathers so it does nt affect the intensity ratio .
however , once again , considerably more data is needed to determine if there is really a peak at this frequency .
another striking feature of the neutron scattering data is a second resolution limited peak at @xmath104 and @xmath377 .
this also has a possible interpretation in our field theory approach .
it is natural to assume that this peak actually comes from @xmath284 . at this wave - vector , from eq .
( [ strans ] ) , this is proportional to @xmath378 at wave - vector @xmath379 .
we expect the continuum limit to hold for some range of wave - vectors close to the gap minimum at wave - vector @xmath380 .
thus , at least for weak enough fields , we would expect the soliton to persist as a single - particle excitation in @xmath284 up to wave - vector @xmath1 .
its energy should obey the lorentz invariant formula : @xmath381 ( here @xmath342 is measured from the incommensurate wave - vector @xmath382 and @xmath172 is set equal to 1 . ) since @xmath383 this gives approximately @xmath384 , as seen in the experiment .
the intensity of this feature in @xmath284 can be easily calculated .
the result follows from the fact that the @xmath342 near @xmath1 parts of the spin operators are all lorentz scalars .
the matrix element between groundstate and a single particle excited state of a lorentz scalar operator is independent of the momentum of the particle , by lorentz invariance ( assuming a lorentz - invariant normalization of the state ) .
it then follows that the soliton and anti - soliton peaks in @xmath284 have an intensity that is proportional to @xmath385 .
the energy is approximately four times higher at @xmath104 .
we must also take into account that both soliton and anti - soliton are contributing at @xmath104 which increases the intensity by a factor of 2 .
thus , we expect the single - particle peak at @xmath386 , @xmath104 to have an intensity approximately 1/2 that of the peak at @xmath387 .
experimentally this ratio looks somewhat larger than 1/2 but it must be remembered that the peak at @xmath386 , @xmath104 is sitting on top of a background from @xmath283 . taking this into account
, the agreement looks fair .
there is actually a possible objection to this argument .
the same reasoning would seem to imply sharp peaks near @xmath384 at @xmath387 coming from @xmath283 and @xmath311 .
these were not observed experimentally ; at most a small shoulder , was observed beginning at @xmath384 for @xmath387 .
this may simply mean that the breathers have merged into the continuum by this wave - vector , ( due to non - relativistic effects not contained in the continuum limit field theory ) whereas the soliton has not .
we note that spin - wave theory fails to capture the one - dimensional quantum fluctuation dominated physics in various ways .
it predicts a single low energy mode with @xmath388 , instead of @xmath389 with soft wave - vectors @xmath1 and @xmath177 , missing the incommensurate shift .
it also predicts another single particle mode at energy approximately @xmath17 , at the same wave - vectors .
in the present framework , the only unknown parameters of cu benzoate are the dm vector , which has not been determined directly in previous studies . based on the present theory , we can in principle determine the dm vector from several experimental results .
actually , there seem to be no solution that can perfectly fit all the available experimental data , as explained below .
presumably , a precise error estimate on an experiment gives a permissible region for the dm vector , and such constraints from several experiments would give a region of possible values of the true dm vector .
however , it should be noted that we have been ignoring the interchain effects , irrelevant operators etc . which might be necessary in such a precise discussion .
firstly , as argued in sec .
ii @xmath22 must lie in the @xmath390 plane , leaving two free parameters .
the total staggered field is determined by eq .
( [ huaeff ] ) .
here we list the constraints on the dm vector from various experiments .      as observed in ref .
[ ] , the induced gap is strongly dependent on the direction of the applied uniform field . in the present theory ,
this is accounted by the field direction dependent constants of proportionality between the uniform and staggered magnetic fields . as discussed in sec .
ii , the proportionality constant should be given by the @xmath10-tensors and the dm vector @xmath43 .
since the @xmath10-tensors were obtained previously@xcite , the measured gap can be used to determine @xmath43 .
the gap is proportional to essentially the 2/3 power of this effective staggered field .
the specific heat measurements of ref .
[ ] were fit by the authors to the specific heat of a collection of free massive relativistic bosons .
the masses were found to scale approximately as @xmath9 with a direction dependent amplitude in the @xmath391 . in ref .
[ ] , this ratio was used to estimate the dm vector .
very recently@xcite the specific heat of the sine - gordon model was calculated from the thermal bethe ansatz and fit to the cu benzoate data . again
a good fit of the gap to @xmath9 was obtained for fields in the @xmath28 or @xmath26 directions with a somewhat different amplitude ratio @xmath392 .
( the velocity , @xmath14 is another parameter in the fit .
this may also be determined from bethe ansatz for the s=1/2 chain in a uniform field
. a slightly better fit to the specific heat data was obtained in [ ] by letting @xmath14 be a free parameter . )    a reasonable fit was not obtained for the field in the @xmath27 direction where the specific heat data is nearly linear .
this suggests that , for h@xmath393 the apparent gap structure was either due to some sort of experimental error or due to other mechanisms@xcite than the effective staggered field . in any case , it seems that the effective staggered field for h@xmath393 is rather close to zero .
this implies the cancellation of the effective staggered field coming from the staggered @xmath10-tensor and the dm interaction .
it is not quite unnatural , because for the applied field in ac - plane , both the staggered field generated by the staggered @xmath10-tensor and the dm interaction point to b - direction .
thus there is a direction in the ac - plane where the cancellation occurs , for a wide range of parameters .
actually , the cancellation of the staggered field is also consistent with the electron spin resonance ( esr ) result , which will be explained later .    assuming the cancellation of the staggered field for h@xmath393 , and that observed gap for h@xmath371b and h@xmath394
are entirely due to the staggered field , the ratio of the proportionality constants between staggered and uniform fields for h@xmath393 , h@xmath371b and h@xmath394 are @xmath395 .
the ratio gives a constraint on the dm vector through  ( [ huaeff ] ) .
based on several exact results on the sine - gordon field theory and on the @xmath396 heisenberg antiferromagnetic chain , we have determined the magnitude of the gap for a given staggered field @xmath68 .
thus , comparing this with the gap estimated from the specific heat measurement , we can fix the proportionality constant @xmath25 between the staggered field @xmath68 and the applied field @xmath17 ( @xmath397 ) .    here
we use only the result for h@xmath394 , which is presumably most reliable . for h@xmath394 ,
the gap is very well fit by the power law @xmath398 , without introducing the logarithmic correction .
the proportionality constant is given by @xmath399 if @xmath111 and @xmath17 are measured in units of kelvin and tesla , respectively . by comparing this with eqs .
( [ deltah02 ] ) and ( [ a0ex ] ) , where we assume the logarithmic factor to be close to unity for the present case , we obtain @xmath400 .
this also gives a constraint on the dm vector through eq .
( [ huaeff ] ) .
an anomalous broadening of esr , which is strongly dependent on the field direction , at low temperatures was observed@xcite in cu benzoate .
the mechanism of this broadening was left unexplained .
however , we have recently developed a field - theory approach to esr on quantum spin chains at low temperature.@xcite according to the theory , the contribution of the staggered field to the esr linewidth is given by @xmath401 this diverges at lower temperature , in agreement with the experiment .
the direction dependence can also be explained by the direction - dependent proportionality constant between the effective staggered field and applied uniform field .
actually , this is consistent with the previous discussion on the field - induced gap , at least qualitatively .
in particular , for h@xmath393 , the low - temperature anomalous part of the esr linewidth vanishes .
this implies the cancellation of the staggered field for h@xmath393 . while this appears to contradict the apparent gap found in ref .
[ ] , it is rather consistent with more refined analysis discussed in the last subsection .
in addition , h@xmath394 gives the largest linewidth , which is consistent with the larger field - induced gap for h@xmath394 .
on the other hand , the ratio of the staggered field is not quantitatively consistent with the specific heat measurement .
the estimate of the staggered field is somewhat subtle because there are also contributions to the esr linewidth from other sources ( most importantly exchange anisotropy / dipolar interaction ) .
the low - temperature anomalous part , which is related to staggered field , appears to be approximately @xmath402 for h@xmath371b and h@xmath394 .
this gives the ratio of the staggered field @xmath403 for h@xmath371b and h@xmath394 .
this is smaller than expected from the specific heat analysis .
as we have discussed in section vi , the analyses on the polarization and intensity of the first breather suggests that the total effective staggered points to @xmath351 ( @xmath352 ) from a - axis , or equivalently @xmath353 ( @xmath354 ) from @xmath27-axis , if the external field is applied in b - direction .
this gives another constraint on the dm vector .
there are several experimental data which give some constraints on the dm vector , and they are not perfectly consistent .
the estimate of dm vector is also sensitive to the assumed form of the g - tensor , extracted from esr measurements.@xcite the experimental data are presumably subject to several errors which has not been identified precisely .
we hope that more experimental data will be available in the future to make more precise comparison with the theory .
the esr linewidth was measured for various directions of the applied field .
thus it is perhaps quite reliable that the staggered field precisely cancels at h@xmath393 .
as we have discussed , this is rather consistent with the refined specific heat result .
this gives a single constraint on the dm vector . in the linearized approximation
, it reads @xmath404    the estimated ratio of the staggered field for h@xmath405b and h@xmath394 was inconsistent between the specific heat and esr .
however , it should be recalled that each analysis has its own problem . in the specific heat analysis , an apparent gap structure , which is unrelated to the staggered field , was observed for h@xmath393 .
whatever the origin of this gap structure , it is natural to expect similar contributions also for other field directions . unfortunately
, we do not know how to estimate these effects at present . on the other hand ,
there are also contribution to the esr linewidth from other sources than the staggered field , and the subtraction causes some uncertainty , in addition to the estimate of the linewidth itself .    the constraint from various experimental result on the dm vector is summarized in fig .
( [ fig : dmest ] ) . for the case of neutron scattering polarization
we have included an estimated error bar from the polarization analysis of ref .
we have not attempted to estimate error bars in the other cases .
we see that a candidate dm vector @xmath406 , which satisfies the most reliable requirement of the cancellation for h@xmath393 , is roughly consistent with all the constraints except for the ratio of the gap between h@xmath371b and h@xmath394 .
this suggests that the observed ratio of the gap was wrong , or was affected by other factors than the staggered field .
more experimental data are needed to draw a reliable conclusion .
in the first subsection we consider the staggered susceptibility , resulting from the application of a staggered magnetic field , using the mapping onto the sine - gordon model . in the next subsection
we combine this with the standard uniform magnetization of the s=1/2 heisenberg model to obtain the total physical susceptibility .
we first discuss the susceptibility of the sine - gordon model . writing the sine - gordon lagrangian in the form : @xmath407 and adopting units where the velocity is set to 1 ,
we define the sine - gordon susceptibility as : @xmath408 where @xmath409 is the free energy .
the groundstate energy ( i.e. the @xmath410 free energy ) is expressed in terms of the gap as:@xcite @xmath411 this determines the @xmath410 susceptibility using the exact relationship between the coupling constant , @xmath412 and the mass @xmath111 : @xmath413 where @xcite @xmath414^{2/3}\left[{\gamma ( 1/6)\over \gamma ( 2/3)}\right ] \approx 4.82764 .
\label{atildedef}\ ] ] this gives the @xmath410 susceptibility : @xmath415 the high - temperature susceptibility is given by : @xmath416 this green s function is normalized to @xmath417 ( at @xmath418 ) . at finite @xmath419 and @xmath420
we have : @xmath421^{1/2}}.\label{gft}\ ] ] here @xmath422 .
thus , the susceptibility becomes : @xmath423^{1/2}}={1\over t}\left[{\gamma ( 1/4)\over \gamma ( 3/4)}\right]^2 \approx { 8.75376 \over t}.\label{chisght}\end{aligned}\ ] ] ( this result for the integral can be obtained , by analytic continuation , from the general results of schulz.@xcite ) an integral equation determining the sine - gordon free energy at finite @xmath420 was given in ref .
[ ] . we may determine the susceptibility by differentiating twice with respect to @xmath412 .
note however that this integral equation actually determines @xmath424 so we must add the zero - temperature part of the susceptibility , given above .
the resulting susceptibility , @xmath425 is plotted versus @xmath426 in fig .
( [ fig : susc ] ) .
as expected , it agrees quite well with high temperature result @xmath427 down to @xmath428 .
it has a maximum at about @xmath429 .
the @xmath410 value is given by eq .
( [ chi0sg ] ) , ( [ atildedef ] ) .
up to a multiplicative factor and logarithmic corrections , the sine - gordon susceptibility essentially gives the staggered susceptibility of the s=1/2 chain , i.e. its response to a staggered field : @xmath430 in order to determine this factor and estimate the logarithmic corrections we first consider the @xmath410 staggered magnetization of the s=1/2 chain .
we refer to the staggered magnetization as @xmath431 : @xmath432 in the continuum limit , @xmath433 since this operator has scaling dimension @xmath187 , a standard rg scaling argument gives the scaling of the staggered magnetization with staggered field : @xmath434 for some function @xmath435 . for weak fields
the exponent is approximately 1/3 . in a similar way to our analysis of the gap in the previous section , by combining the exact results for the xxz model with an rg analysis of the marginal operator we may determine the scaling of magnetization with staggered field in the case of zero uniform field : @xmath436^{1/3 } \label{magscale2}\ ] ] and determine the behavior of @xmath437 in eq .
( [ magscale1 ] ) for small field : @xmath438^{1/3}.\label{cas}\ ] ] @xmath431 obeys a standard rg equation relating a change in the cut - off energy scale , @xmath153 , to a change in the coupling constant , @xmath160 : @xmath439m=0.\ ] ] working to linear order in the marginal couplings , as before , we set @xmath440 and use eq .
( [ gamma ] ) . using eq .
( [ rge ] ) , and lowering the cut - off scale to the gap , @xmath111 , this gives : @xmath441^{-1/4 } , \label{mlambda}\ ] ] for some constant , @xmath409 .
( [ hlambda ] ) , and ( [ a0ex ] ) ) this can be written : @xmath442^{-1/3},\ ] ] where @xmath193 is defined in eq .
( [ gaph0 ] ) and determined in eq .
( [ a0ex ] ) . as in the previous sections we obtain the various formulae from the different asymptotic scaling of @xmath443 in the 3 cases : @xmath120 , @xmath225 and @xmath120 with exchange anisotropy @xmath237 , defined in eq .
( [ epsilon ] ) .
( [ lambdaiso ] ) we obtain eq .
( [ magscale2 ] ) with : @xmath444 using eq .
( [ lambdar ] ) we obtain eq .
( [ magscale1 ] ) with : @xmath445^{-1/3 }      \left ( \frac{h}{j}\right)^{-2(1 - 2\pi r^2)/3}.\ ] ] using eq .
( [ lambdaz0 ] ) we obtain eq .
( [ cas ] ) with : @xmath446    using eq .
( [ lambdaep ] ) , in the limit @xmath236 , in eq .
( [ mlambda ] ) we obtain an expression for @xmath431 in terms of @xmath111 with exchange anisotropy : @xmath447 we may determine the constant @xmath409 , by comparing to the exact result of ref . [ ) ] . from eq .
( [ prop ] ) : @xmath448>,\ ] ] with the exact formula for @xmath449 give in eq .
( [ a(r ) ] ) .
this expectation value is given in terms of the soliton mass , @xmath111 in [ ] , eq .
( 15 ) , with : @xmath450 inserting a power of the spin - wave velocity , @xmath172 by dimensional analysis , and taking the limit of small @xmath237 , we obtain : @xmath451 this gives a result consistent with eq .
( [ manis ] ) and determines the constant , @xmath409 , to be : @xmath452^{-3/2}=2a_0^{3/2}/(3\sqrt{3}\pi ) , \ ] ] where the constant @xmath193 is defined in eq .
( [ gaph0 ] ) and its value is given in eq .
( [ a0ex ] ) .
here we have used the exact identity @xmath453 .
hence the amplitude of eq .
( [ magscale2 ] ) is given by : @xmath454 we thus obtain the @xmath410 staggered susceptibility of the s=1/2 chain by differentiating @xmath431 with respect to @xmath68 : @xmath455 for @xmath120 . comparing to eqs .
( [ atildedef ] ) and ( [ chi0sg ] ) , we see that : @xmath456    the susceptibility for @xmath457 ( and @xmath120 ) follows from eqs .
( [ gft])-([chisght ] ) .
here we use the exact result for the @xmath410 correlation function of the s=1/2 chain:@xcite @xmath458 this differs from the correlation function of the sine - gordon model by a factor of @xmath459 .
[ note the factor of 4 difference in the susceptibilities due to the factor of 2 in the interaction term of the sine - gordon lagrangian of eq .
( [ sgnorm ] ) . ] upon going to finite @xmath420 and fourier transforming at zero frequency and wave - vector , we expect the logarithmic factor to become : @xmath460 . thus : @xmath461 comparing to eq .
( [ chia0 ] ) suggests the heuristic formula : @xmath462\over 2(2\pi ) ^{3/2}}\chi_{sg}(t,\delta ) .\ ] ] @xmath463\approx \chi_s$ ] is plotted in fig .
( [ fig : susc ] ) .    including a small uniform field , @xmath17 , only makes unimportant changes in these formulas . at @xmath410 , the power of @xmath111 in eq .
( [ chia01 ] ) changes by a small amount ; the argument of the logarithm changes to @xmath464 and the amplitude by a factor of @xmath465 . for @xmath457 we must distinguish two regimes , depending on the relative magnitude of @xmath17 and @xmath420 . for @xmath466 , the power of @xmath420 changes .
on the other hand , for @xmath467 we expect to obtain eq .
( [ chialt ] ) .      above we considered the staggered susceptibility resulting from ( independent ) staggered and uniform fields . to make any comparison with experiments we must take into account that the effective staggered field is proportional to the uniform field , @xmath468 where the constant of proportionality
, @xmath25 is strongly dependent on field direction .
the physical susceptibility is conveniently obtained from its thermodynamic definition : @xmath469 where @xmath409 is the free energy .
writing @xmath409 as a function of uniform and staggered fields , we must set @xmath260 before taking the @xmath17-derivative .
we may calculate the free energy in the rotated spin basis of eq .
( [ redef ] ) , used throughout this paper .
noting that the first derivative of @xmath409 with respect to @xmath17 or @xmath68 gives the uniform magnetization , @xmath470 and staggered magnetization @xmath431 respectively ( in the rotated basis ) , we obtain : @xmath471 where we have used : @xmath472 we may ignore the dependence on @xmath68 of the first term and use the standard result for the uniform susceptibility of the @xmath396 chain , @xmath473 . for low fields and temperatures
this gives : @xmath474 independent of field and temperature .
the second term in eq .
( [ chiphys ] ) is larger than the third so we approximate : @xmath475 where @xmath476 is the staggered susceptibility discussed in the previous subsection .
thus , when measuring the physical susceptibility of the present system , one actually probes also the staggered susceptibility  @xcite .
while the first term , the standard result for the susceptibility of the s=1/2 chain , goes to a finite constant , at @xmath420 an @xmath18 , the second term , resulting from the effective staggered field , is highly singular . at zero field
it blows up at @xmath477 as @xmath478 .
thus , although it is multiplied by the small constant , @xmath479 , it eventually dominates for low enough @xmath420 . for any finite field
, the divergence of the second term is cut off , at essentially the gap energy , @xmath480 , at a value of approximately @xmath481 , as shown in fig .
( [ fig : susc ] ) . at low fields
the behavior looks quite similar to a paramagnetic impurity contribution .
it can be distinguished from that , however , by its very strong field - direction dependence .
the effect is largest for the field in the c - direction when the parameter @xmath482 .
the experimental susceptibility@xcite of cu benzoate shows very peculiar behavior at low fields and temperatures .
as the temperature is lowered , in low fields , the susceptibility grows .
this effect is highly direction dependent with the biggest effect occurring for fields in the c - direction .
this effect is cut off by the application of a field .
qualitatively , the experimental results look similar to our predictions .
however , at a more quantitative level there are many differences .
the anomalous low field low t part of @xmath483 is much larger in the experiments than in the theory .
for instance , at zero field , a temperature of @xmath484 , and @xmath485 , the value that we find for h@xmath371c , the second term in eq .
( [ chiphys ] ) is still smaller than the first by a factor of about 1/2
. the anomalous contribution could be somewhat larger due to the third term in  ( [ chiphys ] ) and the logarithmic correction .
however , the anomalous term in @xmath486 observed in cu benzoate at this zero field and low t is about six times _ larger _ than the @xmath473 .
this is perhaps too large to be explained by our theory which is purely one - dimensional .
furthermore , the detailed experimental dependence on t and h is quite complex .
the maximum susceptibility occurs at a field of about 30 g . and t=1.5k .
at lower fields or temperatures the susceptibility decreases somewhat . at low fields , two extra low temperature peaks
are observed in the susceptibility as a function of t ( in addition to the normal one for an s=1/2 chain at @xmath487 ) . while it is tempting to try to identify the higher-@xmath420 peak with the peak in @xmath476 shown in fig .
( [ fig : susc ] ) , and the lower-@xmath420 peak with the onset of nel order , both susceptibility peaks are broad , unlike what would be expected from a phase transition .
( indeed , no evidence for magnetic order was found from neutron scattering,@xcite down to t=@xmath488k . ) a strong frequency dependence of the low t and h susceptibility was also observed .
we expect that a proper theoretical description of the low field and temperature susceptibility of cu benzoate will require the inclusion of inter - chain coupling effects .
the interchain super - exchange paths look complicated and it is not even clear what is the sign of the interchain coupling .
we remark that a perfectly one - dimensional s=1/2 antiferromagnet with su(2 ) symmetry broken only by the dm interaction has a disordered groundstate .
this follows from our spin redefinition in sec .
ii which maps the system into an _ easy plane _ xxz model which is well - known to have a disordered groundstate with power - law correlations .
we may think of the spins as fluctuating primarily in the plane perpendicular to @xmath22 with a tendency for the cross - product of neighboring spins to be parallel to @xmath22 .
this system has a u(1 ) symmetry .
interchain coupling would normally be expected to produce long - range order with both an anti - ferromagnetic component and a possible perpendicular ferromagnetic component .
for example , for @xmath489 , @xmath490 see fig .
( [ fig : dmclass ] ) .    the standard mean field treatment of the interchain interactions
@xcite would suggest a critical temperature of order the interchain coupling ( assuming that it is antiferromagnetic ) . on the other hand , when a magnetic field is applied ( not parallel to @xmath22 )
, the u(1 ) symmetry is broken and the phase transition should disappear . @xmath431 and
@xmath470 now becomes non - zero even in the purely one - dimensional system and are no longer order parameters for spontaneous symmetry breaking .
thus the application of a magnetic field smoothes out the phase transition in this system .
surprisingly , neutron scattering experiments in zero field have failed to detect a nel transition .
we would like to thank y. ajiro , c. broholm , w. hardy and d. reich for helpful discussions .
m.o . thanks t. kato for his help on drawing crystal figures .
this research was supported in part by nserc of canada and by the ubc killam fund .
dender , p.r .
hammar , d.h .
reich , c. broholm and g. aeppli , phys .
lett . * 79 * , 1750 ( 1997 ) .
dender , ph.d .
thesis , johns hopkins university , 1998 .
m. oshikawa and i. affleck , phys .
lett . * 79 * , 2883 ( 1997 ) .
i. dzyaloshinskii , j. phys .
solids * 4 * , 241 ( 1958 ) .
t. moriya , phys . rev .
* 120 * , 91 ( 1960 ) .
essler and a.m. tsvelik , phys . rev .
* b57 * , 10592 ( 1998 ) .
essler , cond - mat/9811309 .
k. oshima , k. okuda and m. date , j. phys .
. jpn . * 44 * , 757 ( 1978 ) .
l. shekhtman , o. entin - wohlman and a. aharony , phys .
* 69 * , 836 ( 1992 ) .
t. sakai and h. shiba , j. phys .
japan * 63 * , 867 ( 1994 ) .
cardy , _ scaling and renormalization in statistical physics _ , cambridge university press , 1996 .
s. lukyanov and a. zamalodchikov , nucl .
phys . * b493 * , 571 ( 1997 ) .
i. affleck , j. phys . *
a31 * , 4573 ( 1998 ) .
s. lukyanov , preprint cond - mat/9712314 .
i. affleck , d. gepner , h.j .
schulz and t. ziman , j. phys .
* a22 * , 511 ( 1989 ) . n.  m. bogoliubov , a.  g. izergin and v.  e. korepin , nucl .
* b275 * , 687 ( 1986 ) ; v.  e. korepin , n.  m. bogoliubov and a.  g. izergin , _ quantum inverse scattering method and correlation functions _ , cambridge university press ( 1993 ) .
griffiths , phys .
rev . * 133 * , a768 ( 1964 ) .
a. fledderjohan et al .
phys . rev . * b54 * , 7168 ( 1996 ) .
cabra et al .
phys . rev . * b58 * , 6241 ( 1998 ) .
hammar et al .
phys . rev .
* b59 * , 1008 ( 1999 ) .
r. dashen , b. hasslacher and a. neveu , phys . ref . * d11 * , 3424 ( 1975 ) .
d. reich , private communication .
k. okuda , h. hata and m. date , j. phys .
. jpn . * 33 * , 1574 ( 1972 ) m. oshikawa and i. affleck , preprint cond - mat/9904199 .
baxter , _ exactly solved models in statistical mechanics _ .
london : academic press ( 1982 ) . c. destri and h. de vega , nucl .
b * 358 * , 251 ( 1991 ) .
schulz , phys . rev . *
b34 * , 6372 ( 1986 ) .
m. fowler and x. zotos , phys .
b * 25 * , 5806 ( 1982 ) .
m. matsuura and y. ajiro , j. phys .
. jpn . * 41 * , 44 ( 1976 ) .
dender , d. davidovi , d.h .
reich and c. broholm , phys . rev . * b53 * , 2583 ( 1996 ) .
see for example , i. affleck , m. gelfand and r. singh , j. phys . * a27 * , 7313 ( 1994 ) ; i. affleck and b. halperin , j. phys .
* a29 * , 2627 ( 1996 ) .