the kondo necklace ( kn ) type models are useful to discuss the quantum phase transitions between kondo singlet and antiferromagnetically ordered states such as found in heavy fermion compounds @xcite .
they have been originally proposed by doniach @xcite for the one - dimensional case as a simplified version of the itinerant kondo lattice ( kl ) models @xcite . thereby the kinetic energy of conduction electrons is replaced by an intersite exchange term . for a pure xy - type intersite exchange
this may be obtained by a jordan - wigner transformation .
however in higher d the replacement can not be justified easily .
the intuitive argument is that at low temperatures the charge fluctuations in the kondo lattice model are frozen out and the remaining spin fluctuation spectrum can be simulated by an antiferromagnetic inter - site interaction term of immobile @xmath1 spins coupled by a kondo interaction to the local noninteracting spins * s*. recent exact diagonalisation studies on finite clusters for both kondo lattice and kondo necklace models have indeed found that the competition between on - site kondo singlet formation and af inter - site correlations are very similar in both models @xcite .
a more formal way to get rid of charge fluctuations in the kl model is an inclusion of a hubbard term for conduction electrons ( klu model ) which leads to the isotropic kn model in the large u limit at half filling
.    nevertheless one should consider the kondo necklace model for d@xmath22 as a model in its own right which is suitable for studying quantum phase transitions between a kondo singlet ( ks ) and an antiferromagnetic ( af ) phase . in its original form
the local kondo exchange is isotropic while the intersite exchange is of xy - type .
this model has u(1 ) symmetry .
later more general models with arbitrary anisotropy of kondo as well as intersite exchange terms have been considered @xcite . indeed compounds which exhibit the ks to af quantum phase transition have mostly uniaxial symmetries .
a full account of the influence of uniaxial anisotropies of both terms in the kondo necklace hamiltonian on the quantum critical point has been given in ref . .
the d=2 kn model without any anisotropy may also be understood as a special case of a bilayer heisenberg model @xcite where the intersite bonds are cut in one of the layers .
a reintroduction of holes in this case leads to the kntj model which is related to the klu model away from half filling @xcite .    in the general anisotropic kn model the quantum phase transition
is achieved by varying the ratio of hopping t , i.e. , the intersite interactions of @xmath1 spins to the on - site kondo interaction j. in practice this is achieved by varying pressure ( hydrostatic or chemical ) .
an alternative way to arrive at the qcp is to apply an external magnetic field which breaks the local kondo singlets and leads to a field dependence of the critical t@xmath0 . starting from a noncritical or above critical t at zero field the system
may then be tuned to the qcp by varying the field strength .
this is indeed a practical method frequently applied @xcite . to investigate field - induced quantum phase transitions from kondo singlet to
af ordered state or vice versa we have extended our previous work to include the effect of the magnetic field within the kondo necklace model .
however , this introduces an additional parameter , namely the ratio of effective g- factors for the local kondo spins and the interacting spins
. they can be different due to the different strength of spin - orbit coupling and crystalline electric field effects involved in the formation of the @xmath1 and * s * ( pseudo- ) spins . a mainly numerical study of the itinerant isotropic kl model in a magnetic field has been given previously @xcite .    in sec .
[ sec2 ] we define the anisotropic kn model in a magnetic field and in sec .  [ sec2a ] briefly discuss the local state space of the kondo term in an external field . in sec .
[ sec3 ] we perform the transformation from spin to bosonic variables and in sec .  [ sec4 ] derive the selfconsistent equations for the mean field boson condensate amplitudes in the singlet and antiferromagnetic phases .. the influence of higher order terms in the hamiltonian is discussed in sec .
[ sec4a ] . in sec .
[ sec5 ] we investigate the numerical solutions and discuss the resulting quantum critical t - h phase diagram .
a discussion is provided in sec .
[ sec5a ] and sec .
[ sec6 ] finally gives the conclusions .
to investigate field induced quantum critical behaviour we start from the anisotropic kondo necklace ( kn ) model with u(1 ) symmetry @xcite where the field is applied along the anisotropy ( z ) direction : & = & _ t+_j+_z + & = & t_<n , m > ( _ n^x _ m^x + _ n^y _ m^y +
m^z ) + j_n(_n^x s_n^x+_n^y s_n^y + _
n^z s_n^z ) + _ is_n^z+_i_n^z , [ ham ] here , the summation over nearest neighbor ( n.n . )
is indicated by brackets and @xmath3 is the @xmath4-component ( @xmath4 = x , y , z ) of the itinerant electron spin at site @xmath5 whereas @xmath6 is the @xmath4-component of localized spins at position @xmath5 . for the exchange coupling between the itinerant and localized spins we generally use @xmath7 as reference energy scale in all figures except when stated otherwise .
the local anisotropy parameter @xmath8 is defined by the relation j@xmath9=@xmath8j@xmath10 between the z - axis and in - plane ( xy ) local exchange .
the hopping parameter of the itinerant electrons is proportional to @xmath11 with the anisotropy in the z - direction given by @xmath12 .
the present model has three control parameters : @xmath13 and the anisotropy parameters ( @xmath12,@xmath8 ) . in the zeeman term we defined @xmath14 and @xmath15 with h = @xmath16h where h is the strength of the applied field .
furthermore @xmath17 and g@xmath18 are the gyromagnetic ratios for localised ( @xmath6 ) and itinerant ( @xmath3 ) spins respectively .
they are determined by the combined effect of spin - orbit coupling and crystalline electric fields which depends on the degree of localisation or itineracy of electrons .
therefore g@xmath19 and g@xmath18 in general need not be equal .
we shall consider two extreme cases , namely ( @xmath20)=(2,0 ) and ( @xmath20)=(2,2 ) .
the former is more realistic , since in real heavy fermion compounds most of the magnetic response is due to the localised electrons with pseudo - spin * s*. as in previous work @xcite the present study of field induced quantum phase transitions in the anisotropic kn model is based on the bond operator formulation .
its hilbert space is spanned by local singlet - triplet states of ( @xmath21 ) spin dimers represented by bosonic degrees of freedom . in applying this technique to the finite field case
we will largely use the same or similar notations as in the previous zero field case @xcite for consistency .
before we perform the transformation to bosonic variables it is useful to have a clear understanding of the local singlet-triplet level scheme as function of anisotropy and magnetic field because certain tendencies of the quantum critical behaviour are correlated with the splitting of the ground state and the first excited state of a local bond .
therefore we first diagonalise the local hamiltonian @xmath22 .
this leads to eigenvalues @xmath23=e@xmath24/@xmath25 ( i=1 - 4 ) given by _ 1,2&=&_+ , + _
3,4&=&(4+_-^2)^- , + _ & = & -(g_sg _ ) , [ llev ] with @xmath26 and @xmath27 . defining s@xmath28=@xmath29 + s@xmath9 then @xmath30 are the energies of the triplet states @xmath31 with s@xmath32 1 respectively which exhibit the linear zeeman effect .
furthermore @xmath33 correspond to the triplet @xmath34 state @xmath35 and singlet state @xmath36 respectively .
they show level repulsion in a magnetic field , except for g@xmath19=g@xmath18 when the total s@xmath28-component is conserved and commutes with @xmath37 .
it is obvious that the spectrum of eigenstates in eq.([llev ] ) is invariant under the transformation ( g@xmath19,g@xmath18 ) @xmath38 ( g@xmath18,g@xmath19 ) .     ) on the local anisotropy @xmath8 ( xy - side ) or @xmath39=1 - 1/@xmath8 ( ising side ) at zero field .
numbers in parentheses denote the degeneracy of each level .
bottom : field dependence ( h=@xmath16h ) of local energy level for three extreme cases .
note that in the xy - case the excitation energy from ground state to first excited state vanishes asymptotically ( h@xmath40j@xmath10 ) whereas in the other cases it reaches a constant ( haf ) or is equal to a constant ( ising ) .
here we used ( g@xmath19,g@xmath18)=(2,0).,title="fig:",width=302 ] +   ) on the local anisotropy @xmath8 ( xy - side ) or @xmath39=1 - 1/@xmath8 ( ising side ) at zero field .
numbers in parentheses denote the degeneracy of each level .
bottom : field dependence ( h=@xmath16h ) of local energy level for three extreme cases .
note that in the xy - case the excitation energy from ground state to first excited state vanishes asymptotically ( h@xmath40j@xmath10 ) whereas in the other cases it reaches a constant ( haf ) or is equal to a constant ( ising ) .
here we used ( g@xmath19,g@xmath18)=(2,0).,title="fig:",width=302 ]    the variation of zero - field energy levels with @xmath8 is shown in fig .
[ fig1 ] ( top ) from the xy - u(1 ) limit ( @xmath8=0 ) via the heisenberg su(2 ) point(@xmath8=1 ) to the ising z@xmath41 limit ( @xmath8=@xmath42 ) .
the field dependence for a few selected @xmath8-values and g - factors ( g@xmath19,g@xmath18)=(2,0 ) are shown in fig .
[ fig1 ] ( bottom ) . for ( g@xmath19,g@xmath18)=(2,2 ) ( not shown in fig .  [ fig1 ] ) the ground state and first excited state levels cross at h@xmath43/j@xmath10=0.25 , 0.50 for ( @xmath44)=(0,0 ) and ( 1,1 ) respectively .
in such a case one may expect that the quantum critical t@xmath0(h@xmath43 ) for the transition from kondo singlet to af phase vanishes .
we apply the transformation from spin variables ( @xmath45,@xmath46 ) to bond variables ( s , t@xmath47 ) ( @xmath4 = x , y , z ) to the hamiltonian in eq .
( [ ham ] ) .
it is given by @xcite : s_n , = ( s^_n t_n , + t_n , ^ s_n -i _ t_n , ^t_n , ) , + _
n , = ( -s^_n t_n , - t_n , ^ s_n -i _ t_n , ^t_n , ) , [ bond ] where @xmath48 is the fully antisymmetric tensor . by construction
the singlet ( s ) and triplet ( @xmath49 ) operators generate the local eigenstates of @xmath37 for zero field and no anisotropy @xcite : s^|0&=&(|-|);t_x^|0=-(|-| ) , + t_z^|0&=&(|+|);t_y^|0=(|+| ) . here
the first and second arrows indicate the z - component of @xmath1 and @xmath50 spins respectively .
the singlet and triplet operators satisfy the usual bosonic commutation relations according to [ s_n , s_n^]=1 ; [ t_n , , t_n , ^]= _ , ; = 0 .
all other commutators vanish .
the physical states have to satisfy the local constraint @xmath51 .
this transformation leads to an effective hamiltonian in terms of singlet ( s ) and triplet ( t@xmath47 ) bosons .
it has been argued in @xcite that one may restrict to the terms which are bilinear in the triplet bosons and we will later in sect .
[ sec4a ] discuss to which extent this is justified .
for the moment we restrict to the bilinear contribution . it can be written as the sum of a local term @xmath37 and an inter - site interaction term @xmath52 . here _
l=_j+_z&= & _ n ( -(2+)s_n^s_n + ( 2-)t_n , z^t_n , z + ( t_n , x^t_n , x+t_n , y^t_n , y ) ) + & & + _ -_n(s^_nt_nz+t^_nzs_n ) -i_+_n(t^_nxt_ny - t^_nyt_nx ) , [ haml ] contains the on - site kondo interaction ( @xmath53 ) and the zeeman term ( @xmath54 ) .
the bilinear interaction term is given by ^(2)_1&=&_<n , m > _ = x , y ( s_n^ t_n , ( s_m^ t_m , + t_m , ^s_m)+ h.c . ) + & & + _ < n , m>(s_n^ t_n , z ( s_m^ t_m , z+t_m , z^s_m)+ h.c . ) .
[ ham1 ] the physical constraint on the local hilbert space is enforced by adding a lagrange term at each site with an associated chemical potential @xmath55 , leading to = _
l+_n _ n ( s_n^s_n+ _ = x , y , z t_n , ^t_n , -1)+^(2)_1 _ 0+^(2)_1 .
[ hammu ] this hamiltonian is diagonalised within a mean field approximation for the bond - operator singlet and triplet bosons .
we assume that in general there are three independent bosonic amplitudes which characterise the phases in the ( t , h)-plane : singlet @xmath56 denoting the strength of local singlet formation , staggered triplet @xmath57 which determines the af order parameter ( polarised along x ) and homogeneous triplet @xmath58 which determines the magnetisation caused by the external field along z - direction . for technical reasons
it is of advantage to transform to circular polarised transverse triplet coordinates ( u@xmath59,d@xmath59 ) with respect to the field direction ( z axis ) .
the transformation and its inverse are given by u_n&=&-(t_nx - it_ny ) ; t_nx=-(u_n - d_n ) , + d_n&=&(t_nx+it_ny ) ; t_ny=-(u_n+d_n ) .
[ circ ] in circular triplet coordinates the hamiltonian may then be written as _ 0&=&_n\ { [ -(2+)+_n]s_n^s_n+ [ ( 2-)+_n]t_nz^t_nz+ + & & [ + _ + + _ n]u_n^u_n+ [ -_++_n]d_n^d_n+ _ -(s_n^t_nz+t_nz^s_n ) -_n } , [ ham0 ] and ^(2)_1&=&-_n , m[(s_n^s_m^(u_nd_m+d_nu_m)+h.c . ) -(s_n^s_m(u_nu_m^+d_nd_m^)+h.c . ) ] + & & + _ n , m[(s_n^s_m^t_nzt_mz+h.c . ) + ( s_n^s_m t_nzt^_mz+h.c . ) ] .
[ ham1ci ] separating the mean values and the corresponding fluctuations the fourier components of singlet and triplet operators are then given by s_k&=&|s , + u_k&=&|u_k , q + _ k , + d_k&=&|d_k ,
q + _ k , + t_k , z&=&|t_0_k,0+_k , z .[separate ] here we assumed that if af order appears it will be of the in - plane type ( @xmath60 ) .
therefore the case of af order will only be considered for @xmath61 . instead of
@xmath62 it will later be convenient to use @xmath63 which is defined by the relation @xmath64 . in the paramagnetic phase at zero field only @xmath65
will be different from zero with a value close to one @xcite .
as the field increases @xmath65 will decrease and simultaneously the triplet amplitude @xmath62 associated with the uniform induced moment increases . if an af transition occurs , @xmath65 will further decrease and u , d will increase accordingly with the field .
the constraint @xmath66 is always respected for any field strength . in the limit
h / j@xmath10 @xmath401 the amplitudes @xmath65 and @xmath62 become asymptotically equal . on the mean field level
this signifies that the ground state of the local level scheme becomes an equal amplitude mixture of the two states with s@xmath67=0 which show the level repulsion in fig .
therefore it is possible to use the zero - field singlet - triplet bosons in eq .
( [ bond ] ) as a basis even in the present finite field problem . on the mean field level
, the physical constraint automatically takes care of the change in the local ground state wave function . inserting the above expressions into the hamiltonian of eqs .
( [ ham0],[ham1ci ] ) one obtains a bilinear form in the triplet fluctuation operators ( @xmath68 , @xmath69 , @xmath70 ) which may be diagonalised with two separate bogoliubov transformations given for the z - polarisation by a_k=(_k , z)_k , z + ( _ k , z)_-k , z^ , + a_-k^=(_k , z)_k , z + ( _ k , z)_-k , z^ , [ bglz ] and for the two circular polarised triplets by [ bglper ] a_k&=&(_k,)_k + ( _ k,)_-k^ , + b_-k^&=&(_k,)_k + ( _ k,)_-k^ , + b_k&=&(_k,)_k + ( _ k,)_-k^ , + a_-k^&=&(_k,)_k + ( _ k,)_-k^.the transformation angles @xmath71 are obtained from the diagonalisation conditions 2_k , z=;2_k,=- [ angles ] where the longitudinal auxiliaryy functions f@xmath9 and d@xmath9 are defined by f_z(k)&= & ( k ) ; d_z(k)=++ ( k ) .
[ fdz ] and for the in - plane case the auxiliary functions f@xmath72 and d@xmath72 are given by f_(k)&= & ( k ) ; d_u , d(k)=++ ( k)_+ , + d_(k)&=&(d_u+d_d)=++ ( k ) . [ fdp ] here @xmath73 with @xmath74 denotes the n.n .
structure factor in dimension @xmath75 ( @xmath76 = coordination number of the simple cubic lattice ) .
it should not be confused with the zeeman energies @xmath77 defined in sect .
[ sec2a ] .
the bogoliubov transformations in eq .
( [ bglz ] ) and eq .
( [ bglper ] ) yield the diagonalised bilinear hamiltonian _ mf = e_0+_k [ _ a(k)a^(k)a(k ) + _ b(k)b^(k)b(k)+_z(k ) a_k^a_k ] .
[ hmf ] where the triplet mode frequencies @xmath78 ( @xmath4=a , b , z ) in the mean field hamiltonian are given by [ modefreq ] _
a(k)&=&_a(k)+_+ , + _ b(k)&=&_b(k)-_+ , + _
z(k)&=&_z(k ) , + _
a(k)=_b(k)&=&[d_(k)^2 - 4f_(k)^2]^ _ ( k ) , + _
z(k)&=&[d_z(k)^2 - 4f_z(k)^2]^ .the excitation energies @xmath78 depend on the field both explicitly through @xmath79 ( for @xmath4=a , b ) and implicitly via @xmath80 ( k ) which is determined by the field dependent singlet and triplet amplitudes . in the nonmagnetic phase ( @xmath81=0 ) the singlet - triplet excitation gap
is given by the minimum excitation energy at the incipient ordering wave vector .
for d=2 this is at q=(@xmath82 ) , explicitly @xmath83 . in the approach from the nonmagnetic side the quantum critical line t@xmath0(h )
is then defined by the vanishing of e@xmath84 .
the ground state energy e@xmath85 is a function of three control parameters ( t / j@xmath10,@xmath8,@xmath12 ) , four singlet - triplet expectation values ( @xmath86 ) and the chemical potential @xmath87 .
writing the transverse mean values explicitly in terms of amplitudes and phases according to |u = u(i_u ) ,    e_0 ( , , ; , , , ) & = & n[-(2j_x+j_z)^2+^2-+ + & & ( + -zt^2)(u^2+d^2)+_+(u^2-d^2 ) + t^2ud(_u+_d ) + & & + ( ( 2j_x - j_z)++zt^2)^2 + 2_- ] + & & + _ k[_k^-d _ ] .as a first step we determine the triplet condensate phases ( @xmath88 ) by minimization . since only one term depends on the phases we obtain the extremal condition @xmath89)=0 .
the minimum must also satisfy @xmath90 .
this is achieved for @xmath91 with n an odd integer , i.e. the sum is only determined modulo 2@xmath92 . since there is no condition for the difference of phases
one of them is arbitrary .
we then choose ( @xmath88)=(0,@xmath92 ) .
the remaining continuous degeneracy with respect to the phase difference is a signature of the goldstone mode which is present throughout the af phase .    ) and triplet ( @xmath93 ) amplitudes and their associated uniform ( @xmath94 ) and staggered ( @xmath95 ) moments as function of field strength normalised to the critical field in the ks and af region . here
at the critical field @xmath96 . this case with t / t@xmath0(@xmath44)=0.6 ( t / j@xmath10 = 0.517 ) corresponds to the the upper curve of fig .
[ fig4 ] ( top ) . here
( @xmath97)=(2.0).,width=302 ]
the minimisation of the ground state energy e@xmath85 in eq .
( [ gsen ] ) leads to selfconsistent coupled equations for the condensate amplitudes @xmath98 and the chemical potential @xmath87 .
their structure is slightly different in the nonmagnetic ( @xmath99=0 ) and magnetic ( @xmath100 ) case , therefore we write them explicitly for both . for a convenient expression of the selfconsistency equations we define the quantity @xmath63 by the relation @xmath101 .
furthermore we introduce the brillouin zone integrals f@xmath47 and g@xmath47 ( @xmath4=a , b , z ) given by f_&=&_k , + g_&=&_k_k^. [ fg ] which appear in the extremal conditions when differentiating the last term in e@xmath85 with respect to @xmath87 or @xmath65 respectively . here
we have defined @xmath102 and @xmath103 .
note that even in the case of finite field ( @xmath104 ) it is the frequency @xmath105 and not the mode frequency @xmath106 which appears in the expressions for f@xmath47 and g@xmath47 . +    ( @xmath12,@xmath8 ) for quantum phase transition between kondo singlet ( t@xmath107t@xmath0 ) and antiferromagnetic ( t@xmath108t@xmath0 ) phases at zero field . for ( @xmath12,@xmath8)=(0,@xmath8 ) or ( @xmath12,1 )
t@xmath0 is obtained from a closed analytical expression ( eq.(25 ) in ref .  )
( full lines ) , for the other cases it is calculated from the zero - field version of eq .
( [ para ] ) ( dashed lines ) .
top : t@xmath0-dependence on anisotropy @xmath8 of local spins . in the ising case
( @xmath109 1 ) t@xmath0 for @xmath12=0 does not vanish because the af order is due to the mixing with the doublet separated by a finite excitation gap .
bottom : t@xmath0-dependence on anisotropy @xmath12 of interacting spins .
the value t@xmath0(0,1)/j@xmath10=0.7 agrees with results from mc simulatons for the genuine kondo lattice model in 2d .
@xcite , title="fig:",width=302 ] + ( @xmath12,@xmath8 ) for quantum phase transition between kondo singlet ( t@xmath107t@xmath0 ) and antiferromagnetic ( t@xmath108t@xmath0 ) phases at zero field . for ( @xmath12,@xmath8)=(0,@xmath8 ) or
( @xmath12,1 ) t@xmath0 is obtained from a closed analytical expression ( eq.(25 ) in ref .  )
( full lines ) , for the other cases it is calculated from the zero - field version of eq .
( [ para ] ) ( dashed lines ) .
top : t@xmath0-dependence on anisotropy @xmath8 of local spins . in the ising case
( @xmath109 1 ) t@xmath0 for @xmath12=0 does not vanish because the af order is due to the mixing with the doublet separated by a finite excitation gap .
bottom : t@xmath0-dependence on anisotropy @xmath12 of interacting spins .
the value t@xmath0(0,1)/j@xmath10=0.7 agrees with results from mc simulatons for the genuine kondo lattice model in 2d .
@xcite , title="fig:",width=302 ]    _ nonmagnetic kondo spin singlet phase in external field _
+ in this case the transverse triplet amplitudes vanish , i.e. , u = d=0 and the minimization of e@xmath85 with respect to ( @xmath110 ) leads to the selfconsistent set of equations [ para ] ^2&= & , + & = & - , + & = & [ ( j_x+j_z)-_g_-2_--zt^2 ^ 2 ] .here @xmath111 are found by iteration which determines @xmath63 and hence @xmath62 completely via the second equation above . for zero field ,
i.e. @xmath112 , the induced @xmath63 vanishes and then eq .  (
[ para ] ) reduces to the paramagnetic case of ref .
+     = singlet gap ) and antiferromagnetic af(m@xmath19 = staggered moment ) phases as function of increasing field strength h are shown : ( i ) field induced af with sequence ks -af .
( ii ) destruction of af and field induced gap opening corresponding to the sequence ks - af .
( iii ) reentrant behaviour corresponding to ks - af - ks .
top : singlet - triplet excitation gap e@xmath84 and af order parameter m@xmath19 as function of external field for the two cases @xmath113 = ( 0,0 ) ( circles ) and @xmath113 = ( 1,1 ) ( diamonds ) as function of the external field . a subcritical scaled value t/ t@xmath0(@xmath44 )
= 0.6 is used in both cases ( t@xmath0(0,0)/j@xmath10=0.350 , t@xmath0(1,1)/j@xmath10=0.862 ) . here
( g@xmath19,g@xmath18)=(2,0 ) is used .
this plot corresponds to the field induced af case ( i ) .
bottom : similar plot for @xmath113 = ( 0,1 ) and an above critical t/ t@xmath0(@xmath44 ) = 1.1 ( t / j@xmath10=0.77 ) and subcritical t/ t@xmath0(@xmath44 ) = 0.94 ( t / j@xmath10=0.658 ) . in the former case ( ii )
one has the sequence af - ks of phases ( @xmath114 ) and in the latter ( iii ) a reentrant situation with ks - af - ks sequence is observed ( @xmath115).,title="fig:",width=302 ] +   = singlet gap ) and antiferromagnetic af(m@xmath19 = staggered moment ) phases as function of increasing field strength h are shown : ( i ) field induced af with sequence ks -af .
( ii ) destruction of af and field induced gap opening corresponding to the sequence ks - af .
( iii ) reentrant behaviour corresponding to ks - af - ks .
top : singlet - triplet excitation gap e@xmath84 and af order parameter m@xmath19 as function of external field for the two cases @xmath113 = ( 0,0 ) ( circles ) and @xmath113 = ( 1,1 ) ( diamonds ) as function of the external field .
a subcritical scaled value t/ t@xmath0(@xmath44 ) = 0.6 is used in both cases ( t@xmath0(0,0)/j@xmath10=0.350 , t@xmath0(1,1)/j@xmath10=0.862 ) . here
( g@xmath19,g@xmath18)=(2,0 ) is used .
this plot corresponds to the field induced af case ( i ) .
bottom : similar plot for @xmath113 = ( 0,1 ) and an above critical t/ t@xmath0(@xmath44 ) = 1.1 ( t / j@xmath10=0.77 ) and subcritical t/ t@xmath0(@xmath44 ) = 0.94 ( t / j@xmath10=0.658 ) . in the former case ( ii )
one has the sequence af - ks of phases ( @xmath114 ) and in the latter ( iii ) a reentrant situation with ks - af - ks sequence is observed ( @xmath115).,title="fig:",width=302 ]    _ antiferromagnetic phase in external field _
+ the additional minimization with respect to u and d amplitudes leads to a direct relation for the chemical potential @xmath87 in terms of the singlet amplitude @xmath65 which is a generalised version of the zero - field af case @xcite : ( h)=zt^2-+[(zt^2)^2+_+^2]^ . at this stage
it is convenient to introduce an auxiliary parameter @xmath116 which controls the effect of the external field on the boson amplitudes ; it is defined by ( h)=[1+()^2]^ + .
[ kappa ] obviously @xmath116(0)=1 and @xmath117 since @xmath118 . using the expression for the zero - field chemical potential @xmath119=@xmath120 we may also write ( h)&=&_0+t^2 , + _ 0&=&t^2-j_z .
we define a total transverse triplet amplitude @xmath99 by @xmath121 which should not be confused with the hopping matrix element t of the hamiltonian .
the minimization with respect to u , d leads to a relation between the two amplitudes from which we obtain u^2+d^2=2 ^ 2 ; u^2-d^2=2 ^ 2 ( ) , [ ud1 ] or equivalently u^2= ; d^2==^2u^2 .
[ ud2 ] the remaining set of singlet and triplet amplitudes ( @xmath122 ) is again determined by the solution of three coupled equations : [ af ] ^2&= & , + & = & - -2_-_0 , + ^2&=&[2-(2j_x+j_z)+2_-+ _
g_]+^2 ^ 2 .these amplitudes are related to the spin expectation values by s_z&=&+(u^2-d^2 ) + _ z&=&-+(u^2-d^2 ) [ spinav ] for zero field ( @xmath63=0 ) the system of equations in eq .
( [ af ] ) may be shown to reduce to the one already discussed in ref .
[ ] . in this case @xmath116=1 and u
= d , i.e. , circular polarised triplet modes with equal amplitudes .
finally we discuss the uniform magnetisation m@xmath85 and the staggered magnetisation m@xmath19=m@xmath123=-m@xmath124 associated with af order . by definition m_0&=&g_ss_z+g__z + m_s&=&g_ss^b_x+g_^b_x .
[ unsta ] here the staggered moment is assumed to be polarised along x. using the expressions in eq .
( [ spinav ] ) for @xmath125 , @xmath126 and similar ones for @xmath127 , @xmath128 we obtain with the help of eqs .
( [ kappa],[ud1 ] ) m_0&=&2g_-2g_-m_0 , + m_s&=&2g_-2 g _ -m_s .
[ magnet ] here we used @xmath27 and defined = .
when both staggered and uniform component are nonzero the total magnetic moment is canted .
the canting angle @xmath4 , counted from the xy plane is then given by [ canting ] =( ) .
( h ) curves ) for various cases of anisotropies ( @xmath12,@xmath8 ) and g - factors ( g@xmath19,g@xmath18 ) . here
circles ( @xmath129 ) correspond to ( @xmath12,@xmath8)=(0,0 ) stars ( @xmath130 ) to ( @xmath12,@xmath8)=(1,1 ) and triangles ( @xmath114 ) to ( @xmath12,@xmath8)=(0,1 ) .
full symbols are obtained by m@xmath19=0 from the af side , open symbols by e@xmath84=0 from the kondo singlet side . for unequal g - factors the heisenberg case shows saturation at @xmath131 0.4 at larger fields , in the ising case t@xmath0 continues decreasing . for the mixed case ( @xmath12,@xmath8)=(0,1 ) ( the genuine kn model in 2d ) and unequal
g - factors the t@xmath0(h ) is nonmonotonic . in each case
the ks phase is below and the af phase above the corresponding t@xmath0(h ) curve .
a decreasing t@xmath0(h ) leads to the phase sequence ks - af ( field- induced af ) whereas an increasing t@xmath0(h ) entails the opposite sequence af - ks . in the nonmonotonic region
the reentrance behaviour ks - af - ks is observed .
for comparison see fig .
for equal g - factors in all cases t@xmath0 decreases monotonically and vanishes at rather small fields.,width=302 ]    the terms @xmath132 are higher order corrections , e.g. @xmath133 for m@xmath19 . for moderate fields
they may be neglected . in sec .
[ sec5 ] ( fig .
[ fig4 ] ) we only show the main contribution @xmath134 in accordance with ref . .
( [ canting ] ) a similar small field approximation may be employed for the canting angle leading to @xmath135 . for equal g - factors g@xmath19=g@xmath18 ( g@xmath136=0 and @xmath62=0 )
one has a special case : m@xmath19 = 0 due to moment compensation of local and itinerant spins , although this is still a phase with an af order parameter m@xmath19=@xmath137 .
the uniform moment m@xmath85 is now given by @xmath138 , oriented parallel to the external field ( @xmath139 ) .
+ _ the ground state energy : condensation vs. fluctuation _ + it is instructive to rewrite the ground state energy in a different way which renders a clearer interpretation of its individual contributions . the total mean field hamiltonian in eq .  ( [ hmf ] ) may be rewritten as _
mf=_0+_=z , a , b_k_(k)[n_+ ] , [ hmftil ] where n@xmath47 is the occupation number of bosons for each mode . the total ground state energy is then given by e_0=_0+_k_(k ) , [ gsentot ] where @xmath140 consists of three parts according to [ gsentil ] _ 0/n&=&-[(2j_x+j_z)(^2+)+(-^2 ) ] + & & + t(-1)^2 ^ 2 ^ 2+((2j_x - j_z)++t^2)^2 + & & + 2_-.here the first part ( @xmath141 ) in eq.([gsentot ] ) is the condensate contribution .
the second always positive term is due to triplet quantum fluctuations . according to eq.([gsentil ] ) @xmath141 is composed of three parts : the first ( negative ) contribution is due to the singlet formation .
the second ( positive ) part is due to the field induced triplet polarisation .
the last term is the zeeman energy contribution . for zero field
only the first term in eq.([gsentil ] ) is present @xcite .
the ground state with @xmath142 will be determined by the competition of these terms .
for example when h = 0 the quantum critical point where af order appears will be determined by the balance of two terms : the negative singlet formation energy and the positive energy of triplet quantum fluctuations .
the transformation of the hamiltonian in eq .
( [ ham ] ) to bond operator variables also creates third and fourth order terms in the triplet operators .
they have been neglected in the previous analysis based on eqs .
( [ haml],[ham1 ] ) .
in fact it was argued in @xcite that third order terms do not contribute to the ground state energy and fourth order terms are quantitatively negligible .
the latter was also found for the related bilayer heisenberg model @xcite .
hence higher order terms have no influence on the zero - field quantum critical properties . in this section
we investigate to what extent this is still justified in the presence of an external field . indeed in mean field approximation used here
the higher order terms lead to field induced effective bilinear triplet terms which have to be added to the genuine bilinear term in eq .
( [ ham1ci ] ) which is present already at zero field .
the quantitative influence of higher order terms on t@xmath0(h ) or h@xmath0(t ) is controlled by the ratio @xmath143 at the critical field h@xmath0 where @xmath144 is the field induced triplet amplitude .
if it is still moderate compared to the singlet amplitude @xmath145 then higher order terms have negligible influence . the field dependence of these amplitudes as obtained from eqs .
( [ para],[af ] ) is shown in fig .
[ fig2 ] for a typical case .
the contribution of third order in t@xmath146 ( n= site , @xmath4=x , y , z ) to the interaction term is given by @xcite : ^(3)_1&=&_n , m [ ( s_n^t_nx+t_nx^s_n)(t^_my t_mz - t^_mz t_my ) + ( s_n^t_ny+t_ny^s_n)(t^_mz t_mx - t^_mx t_mz ) + & + & ( s_n^t_n , z+t_n , z^s_n)(t^_m , x t_m , y - t^_m , y t_m , x ) ] .
[ ham3 ] which has to be added to eq .
( [ hammu ] ) .
we are only interested in the influence of this term on the quantum critical lines h@xmath0(t ) or t@xmath0(h ) to be discussed in sect .
as explained there this may be achieved both from the paramagnetic ( ks ) and af side of the critical line .
here we consider only the former for reasons explained at the end of the section . in this case
the circular triplet amplitudes @xmath147 .
the mean field approximation to h@xmath148 then contains only terms proportional to @xmath149 and no constant contribution to e@xmath85 appears .
transforming these terms to circular triplet coordinates one finally obtains another bilinear contribution with ^(3)_1mf = -z_n(u^_nu_n - d^_nd_n ) + _ nm(u^_nu_m - d^_nd_m ) .
[ ham3mf ] this contribution is field induced since for small fields @xmath150 . the first single site term in eq .
( [ ham3mf ] ) has the same structure as the @xmath79 part of the zeeman term and may be accommodated by a simple ( nonlinear ) rescaling of the applied field such that @xmath79 is replaced by @xmath151 according to [ hcscale ] _
+ ( h)=-g_+hf_s(h ) ; f_s(h)=1+t here @xmath152 is the rescaling function for the applied field .
note that for @xmath153 one has @xmath154 , i.e. no rescaling will occur in this case .
the second contribution in eq .
( [ ham3mf ] ) is an interaction to which the same ( transverse ) bogoliubov transformation as before may be applied .
these terms then have a simple effect : in eq .
( [ fdp ] ) one has to replace @xmath155 .
however only the average @xmath156 of the auxiliary functions , which is unchanged , enters the expression for @xmath157 in eq .
( [ modefreq ] )
. therefore @xmath157 is also unchanged by the third order contribution , as a result the total energy e@xmath85 in eq .
( [ gsen ] ) will be exactly the same as before .
also the selfconsistent equations eq .
( [ para ] ) will be unchanged .
note that in these equations @xmath158 is not rescaled .
the only effect of the third order terms is therefore the above rescaling of the external field in the @xmath79 zeeman term leading to the modified transverse mode frequencies _ a , b(k)=_a , b(k)_+ .
therefore the critical field h@xmath0 of the quantum phase transition which is defined as the field where one of the above modes ( @xmath159 for @xmath160 ) vanishes at the af wave vector * q*=@xmath161 will be changed by the scaling factor @xmath162 .
if h@xmath163 is the critical field without third order terms then approximately @xmath164 is the critical field with the effect of third order terms included .
this is only a quantitative modification which we will discuss in sec .
[ sec5 ] in connection with fig .
the qualitative topology of the phase diagram will not be changed by the third order term .
we note again that in the special case @xmath165=@xmath166 the third order term does not have any effect at all because @xmath167 .
furthermore for all cases g@xmath19=g@xmath18 there is no field induced triplet amplitude , i.e. @xmath168 according to eq .
( [ para ] ) leading again to @xmath167 .
we conclude that the influence of third order terms is not important for the field induced quantum critical behaviour .
+ now we discuss the effect of terms which are of fourth order in the triplet operators .
they are given by @xcite ^(4)_1&=&_n , m((t^_ny t_nz - t^_nzt_ny ) ( t^_my t_mz- t^_mzt_my ) + ( t^_nx t_nz- t^_nzt_nx)(t^_mx t_mz - t^_mzt_mx ) + & + & ( t^_nx t_ny- t^_nyt_nx ) ( t^_mx t_my- t^_myt_mx ) ) . the mean field approximation to @xmath169 for the nonmagnetic case ( @xmath170 ) does not produce a constant term but leads to an effective bilinear contribution in circular triplet coordinates ^(4)_1mf=^2_nm [ u_n^u_m + d_n^d_m + ( u_nd_n+d_nu_m+u_n^d_m^+d_n^u_m^ ) ] .
again this is a field - induced bilinear contribution with @xmath171 for small fields . due to
@xmath172 it has no longitudinal part .
it may be diagonalised by the same transverse bogoliubov transformation in eq .
( [ bglper ] ) as before .
the new mode frequencies including fourth order contributions are simply obtained by the replacement @xmath173 and @xmath174 in eq .
( [ fdp ] ) .
this will lead to the modified transverse frequency _ a , b(k)^2=_a , b(k)^2 [ 1+^2 ] , and new mode frequencies @xmath175 at the quantum critical point where @xmath176 this may be approximately written as _
a , b(k)_a , b(k ) [ 1+()^2_c(k ) ] , [ fourthc ] where the index c denotes the values of singlet and triplet amplitudes at the qcp . from fig .
[ fig2 ] we can estimate that the prefactor in eq.([fourthc ] ) is of the order @xmath177 since at the qcp the triplet amplitude is still quite small and the singlet amplitude basically unchanged from the zero - field value .
therefore the fourth order terms lead to corrections of the order one per cent in the frequencies @xmath157 which determine the last term in the ground state energy of eq .
( [ gsen ] ) .
indeed the correction to @xmath178 is even smaller since the momentum summation over @xmath157 in the last term of eq .
( [ gsen ] ) leads to a large amount of cancellation because @xmath179 is positive in one half of the brillouin zone and negative in the other half .
therefore the selfconsistent equations for @xmath65 and @xmath62 are nearly unchanged and we can conclude that the field induced quantum critical behaviour is not influenced by the inclusion of fourth order terms .
our analysis shows that the quantum critical lines h@xmath0(t ) are only marginally influenced by the presence of higher order terms in one example and have no effect in other cases .
however for fields with @xmath180 one would expect that the higher order terms should lead to non - negligible corrections , at least in the case of third order contributions . from fig .
[ fig2 ] one can see that for @xmath181 the ratio @xmath143 approaches 0.5 and the bilinear approximation might become inadequate .
in fact , by its very construction the bond - operator method is a strong coupling theory which assumes the dominance of the singlet state .
it can not be expected to be quantitatively correct for high fields when saturation of moments is approached , i.e. when singlet and triplet amplitudes become equal . in this limit
it is more appropriate to start from the polarised canted state and perform a conventional spin wave expansion for the two types of spins . + as a consequence of the above analysis we will only consider the genuine bilinear hamiltonian eq .  ( [ ham1ci ] ) as used in the previous sections for the following numerical calculations .
this restriction has an additional reason : although the critical @xmath182 or t@xmath0(h ) may be calculated by searching for the vanishing of the excitation gap e@xmath84 from the nonmagnetic singlet side , it is more satisfactory to search also from the af magnetic side for the vanishing of the staggered magnetisation m@xmath19 and confirm the agreement .
this can not be done easily when third and fourth order terms are included .
on the af side the latter lead to a mixing of the longitudinal and transverse modes and a closed analytic form of the bogoliubov transformation can not be found .
therefore it is better to include only the bilinear terms and be aware of the trivial effect of field rescaling by the third order terms and the tiny effect of fourth order terms .
this standpoint will be adopted for the following numerical analysis .
we shall now discuss the numerical solutions of the selfconsistent eqs .
( [ para],[af ] ) which describes the af - ks quantum phase transitions , i.e. , the field dependence of the kondo singlet gap e@xmath84 and the field dependence of the staggered magnetisation m@xmath19 .
the former was defined below eq .
( [ modefreq ] ) and the latter in eq .
( [ magnet ] ) .
we shall only discuss the results for d=2 and @xmath61 , i.e. , the easy xy - plane situation with @xmath183 c. in this case @xmath184 . as shown before the qualitative behaviour of d=3 is similar to d=2 @xcite .
the main purpose is to study the dependence of the quantum critical point t@xmath0(h ) on the external field , or equivalently the dependence of the quantum critical field h@xmath0(t ) on the hopping strength t. to check consistency the qcp has been obtained both by variation of t and of h and both from the paramagnetic ( e@xmath84=0 ) as well as the antiferromagnetic ( m@xmath19=0 ) side of the qcp . as mentioned in the previous section this requires restriction to the genuine bilinear hamiltonian terms .
+     of total moment moment ( counted from the xy - plane ) in the af phase of two different cases ( cf . fig .
[ fig4 ] ) with ( g@xmath19,g@xmath18)=(2,0 ) .
full curve is obtained from the exact expression eq .
( [ canting ] ) and dashed curve from the approximate expression @xmath185 . at the quantum critical field h@xmath0 the staggered moment vanishes and the total moment is aligned with the external field ( @xmath139 ) . here
h@xmath0(0,1)/j@xmath10=1.23 and h@xmath0(1,1)/j@xmath10=0.55.,width=302 ]    before discussing the field dependent results we show the behaviour of the quantum critical @xmath186 as function of the anisotropies in the zero field case as a starting point of our analysis .
the calculation of @xmath186 has been described in ref .
[ ] . here we present the results for the full range of anisotropies ( fig .
[ fig3 ] ) from the u(1 ) -xy case via heisenberg- su(2 ) to ising - z@xmath41 symmetry .
note the two different scales on the left ( j@xmath10 ) and right ( j@xmath9 ) half of each part of fig .
the quantum critical t@xmath0 for the ks - af transition reaches a singular maximum for the heisenberg case and generally vanishes for the ising case , except for @xmath12=0 .
this may be understood by comparing with fig .
[ fig1 ] . for @xmath187 ( i.e. @xmath188 ) the two singlets with s@xmath28=0 become degenerate and an arbitrary small interaction t leads to the af state .
however for @xmath12=0 the s@xmath28=0 states are not connected by a matrix element of @xmath189 and hence af can occur only via matrix elements between the singlet ground state and the excited doublet which requires a finite t@xmath0/j@xmath10 = 0.35 for af order to occur . + the effect of the external field for ( g@xmath19,g@xmath18 ) = ( 2,0 ) is shown in fig .
[ fig4 ] ( top ) for the two cases ( @xmath190 ) = ( 0,0 ) ( @xmath129 ) and ( 1,1 ) ( @xmath191 ) .
their corresponding t@xmath0 values may be taken from fig .
[ fig3 ] . in fig .
[ fig4 ] ( top ) we choose a subcritical value t / t@xmath192 = 0.6 meaning a nonmagnetic singlet phase with finite e@xmath84 exists for zero field . when the field is increased the gap is gradually reduced until it is closed at t@xmath0(h@xmath0 ) .
for the ( normalised ) t value of 0.6 obviously the difference in critical fields h@xmath0 is not significant .
however , the qualitative behaviour of the spin gap e@xmath84 in the two cases for h@xmath107h@xmath0 is different with a much smaller slope of e@xmath84 in the xy - case as compared to the heisenberg case . in both cases for h@xmath108h@xmath0 a field induced easy plane
af order parameter appears , i.e. , a phase sequence ks - af for increasing field .
the opposite behaviour is observed in the genuine kn case ( @xmath190 ) = ( 1,0 ) with ( g@xmath19,g@xmath18 ) = ( 2,0 ) in fig .
[ fig4 ] ( bottom ) . for an above- critical value t / t@xmath193 = 1.1 ( t / j@xmath10=0.77 ) the af order parameter
is gradually suppressed until it vanishes at h@xmath0/j@xmath10 = 1.23 and the kondo singlet phase appears , i.e. , the opposite phase sequence af - ks is realised . finally , for slightly subcritical t / t@xmath193 = 0.94 ( t / j@xmath10=0.658 ) an interesting reentrance sequence ks - af - ks of phases as function of field is observed .
it exists only in a narrow range of subcritical values @xmath194 . as discussed in sect .
[ sec4a ] this behaviour is robust because third order triplet contributions are exactly zero for @xmath12=0 .
+ we may collect the data of quantum critical fields h@xmath0 from numerous calculations such as presented in fig .
[ fig4 ] and similar ones for quantum critical t@xmath0 at fixed h in a h - t phase diagram .
it is shown in the form of t@xmath0(h ) curves in fig .
[ fig5 ] for two choices of the g - factors . for ( g@xmath19,g@xmath18 ) = ( 2,0 ) and the cases ( @xmath8,@xmath12)= ( 0,0 ) , ( 1,1 ) the scaled t@xmath195 is identical close to the zero - field qcp , i.e. , the slope does not depend on the anisotropy ( @xmath8,@xmath12 ) . for larger fields
one observes the plateau formation in the case ( @xmath8,@xmath12)=(1,1 ) and the further monotonic decrease for ( @xmath8,@xmath12)= ( 0,0 ) .
this behaviour may qualitatively be understood from the field dependence of the local energy levels of @xmath53 shown in fig .
[ fig1 ] . for increasing h the gap between the ground state and first excited state decreases and hence a smaller t@xmath0 is necessary to achieve the softening of the triplet excitation at q. for @xmath8=1 this effect eventually levels off because the splitting of the two lowest local levels becomes constant at large field and approaches 2@xmath8j@xmath10 so that t@xmath195 reaches a plateau at large fields for nonzero @xmath8 . on the other hand for @xmath8=0 the two lowest levels become asymptotically degenerate and
hence t@xmath0(h ) should approach zero for large fields . as already noted in the discussion of fig .
[ fig4 ] the intermediate case ( @xmath8,@xmath12 ) = ( 1,0 ) and ( g@xmath19,g@xmath18 ) = ( 2,0 ) behaves quite differently .
after an initial but much weaker decrease of t@xmath0(h ) it reaches a minimum at ( h / j@xmath10,t / j@xmath10)=(0.275,0.638 ) and then starts to increase again .
the region of the nonmonotonic behaviour with slightly subcritical t / t@xmath0(0,1 ) corresponds to the region where reentrance behaviour ks - af - ks is observed in fig .
[ fig4 ] . for above critical t / t@xmath0(0,1 ) , once t@xmath0(h)@xmath108 t the initial af phase is suppressed and we obtain the sequence af - ks , i.e. a field induced ks phase corresponding to fig ..  [ fig4 ] ( bottom ) . in all three cases with ( g@xmath19,g@xmath18 )
= ( 2,0 ) one has a finite t@xmath0(h ) for finite fields .
+ this is different if we chose equal g - factors ( g@xmath19,g@xmath18 ) = ( 2,2 ) . for this choice
the lowest two local levels cross at h=0.25 ( @xmath8=0 ) and at h=0.5 ( @xmath8=1 ) ( not shown in fig .
[ fig1 ] ) .
this is due to the fact that s@xmath28 commutes with the hamiltonian and hence there is no level repulsion of the two s@xmath67=0 states .
consequently at this point the quantum critical t@xmath0(h ) vanishes in all three cases of ( @xmath8,@xmath12 ) .
it decreases essentially linearly from its maximum value at h=0 . the phase diagram in fig .
[ fig5 ] shows a collection of the critical t@xmath0(h ) curves for the six cases considered .
the corresponding gapped ks and af ordered phases are below and above each curve respectively .
we comment again on the influence of third order terms discussed in sect .
[ sec4a ] on the phase diagram of fig .
remembering that all cases with @xmath196 or @xmath197 are unaffected by the third order contributions , only one of the critical field curves in fig .
[ fig5 ] ( @xmath130 with ( g@xmath19,g@xmath18 ) = ( 2,0 ) ) may be influenced . in this case
third order terms lead to an almost field independent and moderate rescaling of h@xmath0 by a factor @xmath1981.2 according to eq .
( [ hcscale ] ) and fig .
[ fig2 ] . + finally in fig .  [ fig6 ]
we show the dependence of the tilting angle of the total moment out of the xy - plane for two possible cases with af order .
it is calculated from eq .
( [ canting ] ) ( full line ) and a simplified ( low field ) expression given in the caption .
close to the quantum critical field the staggered moment vanishes and the total moment is then alligned with the external field ( @xmath199 z ) , meaning @xmath139 .
as we have mentioned in the introduction the interest for the kondo necklace ( kn ) model stems from two facts : ( i ) it is thought to be a simplified version of the original kondo lattice ( kl ) model with itinerant electrons and ( ii ) it is expected to describe some basic features of the observed quantum phase transitions in real kondo materials .
although there is no mapping or exact equivalence of the two models the agreement of their qualitative features by simply identifying the intersite interaction ( kn ) and hopping ( kl ) energies ( both given by t ) is generally quite good . as discussed in @xcite the 1d kondo necklace model does not exhibit quantum critical behaviour , this agrees with earlier numerical work on the kl model @xcite where the disordered spin gapped state always prevails @xcite .
+ we discuss the more interesting 2d case where quantum phase transitions in the kn model are possible @xcite .
one can compare the values of the quantum critical t@xmath0/j ( we set j@xmath200 j@xmath10 in this section to comply with literature conventions ) obtained for the kn model from fig .
[ fig3 ] with those for the kl model obtained in mc simulations obtained previously @xcite . in the spirit of doniach s replacement in 1d we should compare with the kn model for xy - type interactions , i.e. with the @xmath201 case in fig .
[ fig3 ] which has a t@xmath0/j=0.7 .
the mc simulations for the 2d su(2 ) fermionic kl model lead to quantum critical parameters ( average values ) t@xmath0/j= 0.71 @xcite , 0.69 @xcite and 0.68 @xcite , very close to the appropriate value in fig .  [ fig3 ] ( upper panel : center of full line ; lower panel : left corner of full line ) . note that the critical t@xmath0 s corresponding to full lines in fig .
[ fig3 ] are obtained from a simple analytical formula within the bond - operator mean field approach @xcite .
further support for this scenario comes from exact diagonalisation ( ed ) results for the kl model on small clusters @xcite using open boundary conditions .
the total average local moment @xmath202 calculated as function of j / t shows a reduction to half of the free moment size due to kondo singlet formation in the lattice for the value t@xmath0/j = 0.67 . if it is interpreted as the critical value in the thermodynamic limit this is again close to the above results .
this near equality of analytical kn and numerical kl results seems to suggest that even in 2d the su(2 ) kl model is better described by the u(1 ) xy - type kn model with ( @xmath44)=(0,1 ) since the su(2 ) kn model with ( @xmath44)=(1,1 ) has a considerably larger critical value @xmath203 .
+ the fermionic kl model is obtained from the periodic anderson model ( pam ) by a schrieffer - wolff transformation which eliminates charge fluctuations of f - electrons .
the resulting kondo coupling j is then given in terms of the anderson model parameters by @xmath204 where @xmath205 is the f - level position with respect to the fermi level at zero energy , u is the on - site f - electron repulsion and v their hybridisation strength with conduction electrons .
the periodic anderson model may be studied with the dynamical mean field theory ( dmft ) @xcite .
this method has been used more recently to investigate its quantum critical properties @xcite which should be related to those of the kondo lattice and hence also kn models .
the method is formally a d=@xmath42 approximation but has nevertheless been used for studying electron correlations in finite dimensional lattices such as the 3d pam in @xcite . using the value of j from the schrieffer - wolff transformation the extrapolated t=
0 phase diagram suggests that the critical hopping strength at the qcp is given by @xmath206 .
this is again in reasonable agreement with the mean - field bond operator result for the 3d kn model which predicts @xmath207 @xcite .
+ while the quantum critical behaviour as function of control parameter t / j is well presented , there are few microscopic investigations concerning the field induced quantum critical behaviour of kl / kn type models such as provided here .
the fermionic 2d kl model in an external field with equal g - factors corresponding to ( g@xmath19,g@xmath18 ) = ( 2,2 ) was studied in ref . by using variational and mc methods .
the numerical results for j / t =3 suggest that the transition between kondo and af ordered regime takes place in a field region ( b@xmath208 in notation of ref .  ) between b@xmath209/t = 1.0 and n b@xmath210/t = 2.25 .
in this regime both kondo - like features and transverse magnetic order coexist . from the present calculations of the 2d kn model again with ( @xmath44)=(0,1 ) , g - factors as above and using the data shown in figs .
[ fig3],[fig5 ] one obtains b@xmath0/t = 1.86 ( b@xmath0/j=0.62 ) which is close to the average 1.63 of the above b@xmath211 values . this value for the critical field
was also confirmed by exact diagonalisation ( ed ) results for the kl model on small clusters @xcite where the breaking up of on - site singlet formation was observed above the value b@xmath0/j @xmath212 .
we conclude that kl and kn models also seem to exhibit similar field induced quantum critical behaviour . another matter is the scaling behaviour of singlet gap and staggered magnetisation around the critical field .
the present mean field type theory does not give sufficient insight into that issue .
it has sofar been treated within continuum field theories involving order parameter fluctuations around the critical field @xcite . + the present bond operator mean field results for the kn model are apparently consistent with known numerical results for the fermionic kondo lattice models .
it is less clear whether they are useful for the interpretation of experimental results , especially for the field induced quantum phase transitions . since charge degrees of freedoms are eliminated in the kn model it is strictly speaking more appropriate for the kondo - insulator compounds which have a hybridisation gap due to half filled conduction bands .
however we ignore this subtlety in the following and also apply it to the magnetism of metallic kondo compounds . the zero field quantum critical behaviour as function of the control parameter
t / j can experimentally be mimicked by applying hydrostatic or chemical pressure ( by substitution of elements ) .
this changes mostly the hybridisation and hence the kondo coupling j. in the accessible pressure regime one may assume that t / j varies linearly with pressure . in this manner
it is feasible to drive af ordered heavy fermion systems to the quantum critical point where they become nonmagnetic heavy fermi liquids .
since hydrostatic ( positive ) pressure generally increases j it tends to suppress the af phase while with suitable substitution of elements ( negative ) chemical pressure may decrease j which favors af order . thus the af qcp may be approached from both sides .
there are many examples to be found especially among ce- compounds , for a review see @xcite . in most cases of ce - compounds , however , the af qcp may not be reached directly because it is enveloped by a dome of the superconducting phase and the critical pressure has to be obtained from extrapolation to t=0 . while the kl / kn type models discussed above naturally suggest the af qcp there is no direct way of experimental determination of the critical ( t / j)@xmath0 since its relation to the experimentally accessible critical pressure p@xmath0 ( or critical concentration of substituent ) is unknown .
in fact most of the interest on pressure induced qcp s focuses more on the scaling exponents of various quantities with respect to the distance ( p - p@xmath0 ) to the qcp .
this is extensively discussed in the reviews cited in the introduction .
+ for the field induced quantum phase transition much less experimental results are available .
a review of materials investigated is given in ref . .
a recent example of a field induced destruction of af order is the tetragonal ybrh@xmath41si@xmath41 compound @xcite . at ambient pressure and zero field
it has an af order which is indeed of the easy plane ( xy)-interaction type ( @xmath12=0 ) .
the af order is destroyed at 0.66 t for field along the hard c - axis and a nonmagnetic fermi liquid state appears above the critical field . if one assumes that the local kondo interaction which results from the schrieffer - wolff mechanism is more of the isotropic nature ( @xmath8=1 ) it qualitatively corresponds to the af - ks scenario shown in the lower panel of fig .
[ fig4 ] . as mentioned before the reentrant scenario shown also in this figure exists only in a narrow parameter range and it is perhaps understandable that no such example is known .
what is really surprising is the following observation : most of the parameter cases in fig .
[ fig5 ] would predict the ks - af sequence , i.e. the field - induced af order out of the kondo phase .
this is a natural consequence that in most cases a field reduces the singlet - triplet gap and supports the onset of magnetic order as shown in the upper panel of fig .
we would like to stress that this is also the only phase sequence found in the fermionic kl model @xcite .
however field induced af magnetic order is not easily found in heay fermion compounds .
for example ceni@xmath41ge@xmath41 at ambient pressure and zero field is in the kondo singlet phase and rather close to the af quantum phase transition which may be achieved by appropriate substitution of ni .
thus it is a complementary case to ybrh@xmath41si@xmath41 .
however in an external field there is no field induced transition to transverse af order as one might naively expect from the model discussion , rather it is driven further away from the af qcp as specific heat and resistivity measurements suggest @xcite .
there is however one known example of a heavy fermion system which exhibits field induced af order , though complicated by the appearance of superconductivity .
the compound cerhin@xmath213 is already an af at ambient pressure @xcite .
application of hydrostatic pressure moves it to the nonmagnetic ( ks ) side where it also becomes superconducting .
additional application of a magnetic field destroys the superconductivity and leads to the reestablishment of a field induced af order in a wide range of pressure form 1.4 to 2.4 gpa @xcite .
this is indeed the ks - af sequence of phases corresponding to fig .
[ fig4 ] ( top ) which appears in the majority of cases studied here .
experimentally , however , it seems to be the one most rarely encountered .
+ these examples suggest that while the kondo lattice or necklace type models have some qualitative and instructive properties which are relevant for magnetic quantum phase transitions they may not be able to explain some experimental observations in heavy fermion compounds .
it is well possible that the complete elimination of charge degrees of freedom is too radical by suppressing hybridisation effects of moments with conduction electrons and thus favoring the magnetic phases .
therefore the more fundamental pam model , including appropriate local and conduction electron zeeman terms may be a more promising starting point for future investigations of field induced quantum phase transitions in heavy fermion compounds .
in this work we have investigated the general anisotropic kondo necklace model and its quantum phase transition from paramagnetic singlet phase to af phase and vice versa in an external field .
we have derived a bosonic hamiltonian using the bond operator representation of local and interacting spins and diagonalised the bilinear part by a generalised bogoliubov transformation to circular polarised modes to include the effect of the external field .
higher order parts of the hamiltonian were found to be insignificant for the quantum critical behaviour . in mean field approximation one
obtains a total ground state energy that depends on the singlet amplitude @xmath65 , staggered and uniform triplet amplitudes @xmath214 and the chemical potential .
minimization leads to two sets of selfconsistent equations for kondo singlet and af phase , respectively .    from their numerical solution
we find that in most investigated cases the t@xmath0 value of the quantum critical point decreases monotonously with field strength leading to the kn - af phase sequence for increasing field . in the case of unequal g - factors and ( @xmath8,@xmath12 )
= ( 1,1 ) it levels off at a plateau value t@xmath195/t@xmath215 0.4 while for ( @xmath8,@xmath12 ) = ( 0,0 ) it continues to decrease at larger field .
however for the genuine kn case with ( @xmath8,@xmath12 ) = ( 1,0 ) the t@xmath0(h ) curve is nonmonotonic with a minimum at @xmath216 0.275 at small fields and continues to increase for larger fields .
depending on the size of t this implies two possibilities : for above critical t a suppression of af and the opening of a spin excitation gap for @xmath217 ( af - kn sequence ) or the reentrance behaviour ( kn - af - kn sequence ) of the gapped singlet phase for slightly subcritical t ( as compared to t@xmath0(0 ) ) . for equal g - factors
a crossing of local levels occurs due to a conserved s@xmath67 and the critical t@xmath0 reaches zero already at relatively small fields for all investigated cases of ( @xmath44 ) .
since our treatment is of mean - field type we will get mean - field exponents for the spin gap e@xmath84 and magnetisation m@xmath19 at the field - induced qcp .
this is apparent from fig.[fig4 ] and has not been discussed further .
improvement for the critical exponents will require a selfconsistent renormalisation theory for the triplet excitations close to the quantum critical point .
the authors would like to thank v. yushankhai for helpful suggestions .
a. l. would like to acknowledge the support of sharif university of technology .
l. also would like to acknowledge the max planck institute for the physics of complex systems - dresden where the initial part of this work was started during his visit .
99 a. langari and p. thalmeier , phys . rev .
b * 74 * , 024431 ( 2006 ) .