much attention has been paid recently to the existence and properties of spin liquids in strongly correlated materials , both experimentally and theoretically@xcite .
such phases , first proposed by anderson , have been put forward in the context of high @xmath0 superconductivity as a possible precursor state , but they can occur in a much broader class of situations , both in real materials @xcite and in microscopic models @xcite .
their characterization is challenging  on the experimental side , it requires generally involved low - temperature techniques to prove the absence of long - range magnetic order down to temperature @xmath1 , and concerning the theory , either in hubbard - like models or in localized spin models , spin liquids usually result from frustrating interactions but compete with a variety of possible ordered phases , from usual antiferromagnetic ( af ) order to valence bond crystals , spin nematics etc.@xcite    clear evidence of a spin liquid ground state was found from an exact solution of the kitaev model @xcite of spins @xmath2 on the honeycomb lattice .
this liquid ground state is characterized by a gapless spectrum ( with low - energy excitations described in terms of majorana fermions ) and gapped vortices .
the kitaev model itself triggered a large variety of studies , focusing on topological properties , ( non)-abelian excitations , extensions to higher dimensions@xcite , effects of magnetic@xcite or non - magnetic impurities and of a magnetic field @xcite .
besides its mathematical beauty and the original properties mentioned above , the kitaev model was recently found @xcite to be relevant to orbitally degenerate systems with strong spin - orbit coupling such as layered iridates na@xmath3iro@xmath4 and li@xmath3iro@xmath4 . there
, around each ir atom , the octahedral environment [ see fig .  [ struc](a ) ] of oxygens results in a configuration with one hole in the @xmath5 manifold ; and the large spin - orbit coupling selects locally a low - energy doublet of states defining a pseudospin @xmath2 ; and finally , superexchange processes between ir ions driven by moderately strong coulomb interactions lead to an effective description by the kitaev - heisenberg model @xcite for pseudospins @xmath2 residing on the honeycomb lattice .
the study of this model showed the robustness of the spin liquid phase for finite af heisenberg interactions [ small compared to ferromagnetic ( fm ) kitaev interactions ] , and a quantum phase transition between this liquid phase and a stripe - ordered phase . in na@xmath3iro@xmath4 , where interaction parameters could be close to this transition , af order was found by magnetization and specific heat measurements@xcite and x - ray magnetic scattering@xcite .
the question of how vacancies influence the order and the magnetic response in the kitaev - heisenberg model is of both experimental and theoretical interest . indeed
the vacancy can act as a probe allowing to measure spin correlations in its vicinity ; and it can also have drastic effects on these correlations . for concreteness , in @xmath6 spin-@xmath2 antiferromagnets with long - range order , a non - magnetic impurity can enhance the staggered magnetization around the vacancy @xcite . in dimerized phases ,
the impurity is accompanied by a spinon , which is typically confined to the impurity in valence bond crystals and deconfined in resonating valence bond liquids @xcite , although the question of the confinement length can be subtle @xcite . in a nel - ordered phase , the impurity results in a curie- type contribution to the finite-@xmath7 susceptibility @xcite . in the spin - liquid phase of the kitaev model , a magnetic moment
is induced around a spinless vacancy ; interestingly , its contribution to the spin susceptibility diverges logarithmically for @xmath8 @xcite .
this makes it relevant to investigate the vacancy - induced magnetic response in a model interpolating between @xmath6 and kitaev - like interactions .
moreover , in the iridate samples studied up to now , a substantial amount of disorder might be present , e.g. , site - mixing effects with non - magnetic alkaline ions on the hexagonal ir - sublattice @xcite , which should be taken into account for proper understanding of their physical properties .
our aim is to study the effect of non - magnetic vacancies ( and of pairs of vacancies ) in the different magnetic phases of the kitaev - heisenberg model , at zero temperature .
we adopt similar notations as in ref .  , with a parameter @xmath9 interpolating between heisenberg and kitaev models , and a uniform magnetic field @xmath10 with amplitude @xmath11 : h= -2_i , j _ _
i^_j^+(1- ) _ i , j _ i _ j -_i h _ i , [ ham ] where @xmath12 labels simultaneously an axis in spin space , and ( in the interaction terms ) a bond direction of the honeycomb lattice . for convenience
we express @xmath13 in terms of pauli matrices @xmath14 .
we take as unit cell a _
@xmath15-bond _ [ see fig .
[ struc](b , c ) ] , on which kitaev interactions are of the form @xmath16
. each site is labeled by an index @xmath17 , in bijection with @xmath18 where @xmath19 denotes the position of the unit cell and @xmath20 the sublattice index , respectively .
elementary translations between neighbor cells are @xmath21 in cartesian coordinates .
we consider periodic clusters of @xmath22 sites ( this cluster can be represented with an hexagonal shape , see fig .
[ cormag ] ) , @xmath23 or @xmath24 sites ( both with the shape of a parallelogram ) .
these clusters respect all spatial symmetries of the model , i.e. translation and @xmath25-rotation - other symmetries combining spin- and spatial rotations @xcite are not used here .    in section [ vfs ] ,
we discuss several features of the model in a magnetic field .
we will show how a spin canting mechanism allows one to understand the phase diagram for a field oriented along an easy axis , as well as to understand a directional anisotropy of the static susceptibility .
also the perturbations due to heisenberg interactions and due to small magnetic fields on the spin correlation functions are explored . in section [ ces ]
, we study the effect of a single non - magnetic vacancy in the system , in various phases of the model , and analyze the magnetization pattern in the vicinity of this vacancy in a small field . in section [ ibv ] , we discuss situations where two vacancies , close enough to each other , effectively interact .
we focus mainly on the stripe phase and explain the underlying mechanism of local selection of an ordered pattern by a vacancy pair .
section [ scr ] provides a short summary and some concluding discussions .
we first consider some basic properties of the kitaev - heisenberg model in a magnetic field , and study the magnetic response to a field oriented along different crystallographic directions .
this comparison is important because of the symmetry of the model , which is clearly not @xmath6 invariant due to the kitaev interactions . in the context of iridates ,
this anisotropy is due to spin - orbit coupling which results in the easy spin axes @xmath26 $ ] , @xmath27 $ ] , @xmath28 $ ] [ corresponding to the octahedral axes , see fig .
[ struc](a ) ] , compared to the cubic diagonal @xmath29 $ ] direction .
the clearest anisotropic features are seen in the stripe - ordered phase , and we compare this situation to that of the nel phase .
we first want to characterize the various phases found in a magnetic field .
it is already known @xcite that at zero field , the spin liquid phase extends from @xmath30 to @xmath31 ( the kitaev limit ) ; the stripe - ordered phase , with spins pointing either in @xmath32 , @xmath33 or @xmath15 direction [ for the last case , the spatial structure , i.e. @xmath15-stripe pattern , is shown in fig .  [ stripe2](a ) ] , is located between @xmath34 and @xmath35 ; and a nel - ordered phase is found for @xmath36 . to determine whether these phases persist in a magnetic field and to understand the magnetic response of the system , one needs adapted structure factors which act as order parameters . for the nel phase , it reads : s_neel=_r , r _ ,  _ r , _ r ,  with @xmath37 being 1 for @xmath38 , and @xmath39 otherwise ; for the stripe - ordered phase , the structure factors corresponding to @xmath32- and @xmath15- stripe patterns are defined as follows : s^x(_x)&=&_r , r,, e^i _ x ( r -r ) _ ,
^x_r , ^x_r ,  , + s^z(_z)&= & _ , r,, e^i _ z ( r - r ) ^z_r , ^z_r ,  , with @xmath40 shown in fig .
[ struc](d ) ; for instance , @xmath41 which is equivalent to the more intuitive @xmath42 ( but which is not in the first brillouin zone ) . both wave vectors account for a phase factor @xmath39 between neighboring unit cells separated by @xmath43 in real space .
these structure factors also evidence the symmetries spontaneously broken in each of the 6-fold degenerate ground states in the stripe phase : indeed at zero field , in addition to translational symmetries , the hamiltonian possesses a @xmath44 symmetry - @xmath45 standing for the @xmath46 transformation and @xmath47 ( labeled as @xmath48 in ref .  ) for the invariance by a cyclic permutation of the @xmath32 , @xmath33 and @xmath15 spin components coupled to a spatial rotation by @xmath49 .
( color online ) ( a ) schematic picture of an octahedron ( 6 oxygens surrounding a ir@xmath50 ion ) with the octahedral axes @xmath32 , @xmath33 , @xmath15 - which are easy axes on identically - labeled bonds in the kitaev - heisenberg model .
( b ) a hexagon in the ir@xmath50 plane , indicating bonds labeled @xmath32 , @xmath33 , @xmath15 according to their local easy axis regarding kitaev interactions . ( c )
honeycomb ( direct ) lattice with the labeling @xmath51 of bonds according to their in - plane orientation .
a unit cell ( centered on a @xmath15-bond ) is represented in the dashed rectangle , as well as elementary translations @xmath52 .
( d ) first brillouin zone , with the position of stripe - order wave vectors @xmath53 ( @xmath12).,title="fig:",width=132 ]   ( color online ) ( a ) schematic picture of an octahedron ( 6 oxygens surrounding a ir@xmath50 ion ) with the octahedral axes @xmath32 , @xmath33 , @xmath15 - which are easy axes on identically - labeled bonds in the kitaev - heisenberg model .
( b ) a hexagon in the ir@xmath50 plane , indicating bonds labeled @xmath32 , @xmath33 , @xmath15 according to their local easy axis regarding kitaev interactions .
( c ) honeycomb ( direct ) lattice with the labeling @xmath51 of bonds according to their in - plane orientation . a unit cell ( centered on a @xmath15-bond )
is represented in the dashed rectangle , as well as elementary translations @xmath52 .
( d ) first brillouin zone , with the position of stripe - order wave vectors @xmath53 ( @xmath12).,title="fig:",width=132 ] +   ( color online ) ( a ) schematic picture of an octahedron ( 6 oxygens surrounding a ir@xmath50 ion ) with the octahedral axes @xmath32 , @xmath33 , @xmath15 - which are easy axes on identically - labeled bonds in the kitaev - heisenberg model .
( b ) a hexagon in the ir@xmath50 plane , indicating bonds labeled @xmath32 , @xmath33 , @xmath15 according to their local easy axis regarding kitaev interactions .
( c ) honeycomb ( direct ) lattice with the labeling @xmath51 of bonds according to their in - plane orientation . a unit cell ( centered on a @xmath15-bond )
is represented in the dashed rectangle , as well as elementary translations @xmath52 .
( d ) first brillouin zone , with the position of stripe - order wave vectors @xmath53 ( @xmath12).,title="fig:",width=136 ]   ( color online ) ( a ) schematic picture of an octahedron ( 6 oxygens surrounding a ir@xmath50 ion ) with the octahedral axes @xmath32 , @xmath33 , @xmath15 - which are easy axes on identically - labeled bonds in the kitaev - heisenberg model .
( b ) a hexagon in the ir@xmath50 plane , indicating bonds labeled @xmath32 , @xmath33 , @xmath15 according to their local easy axis regarding kitaev interactions . ( c )
honeycomb ( direct ) lattice with the labeling @xmath51 of bonds according to their in - plane orientation . a unit cell ( centered on a @xmath15-bond )
is represented in the dashed rectangle , as well as elementary translations @xmath52 .
( d ) first brillouin zone , with the position of stripe - order wave vectors @xmath53 ( @xmath12).,title="fig:",width=136 ] +     ( color online ) ( a ) representation of a @xmath15-stripe pattern ( stripe pattern with spins oriented along the @xmath54 $ ] axis - neighboring spins are aligned on @xmath15 bonds and antialigned on other bonds .
( b ) schematic representation , at the level of a single hexagon , of a canted stripe pattern : here a @xmath15-stripe pattern , with spins tilted of an angle @xmath55 towards a field along [ 100 ] direction.,title="fig:",width=158 ]   ( color online ) ( a ) representation of a @xmath15-stripe pattern ( stripe pattern with spins oriented along the @xmath54 $ ] axis - neighboring spins are aligned on @xmath15 bonds and antialigned on other bonds .
( b ) schematic representation , at the level of a single hexagon , of a canted stripe pattern : here a @xmath15-stripe pattern , with spins tilted of an angle @xmath55 towards a field along [ 100 ] direction.,title="fig:",width=158 ] +    both af ( nel and stripe ) ordered phases of the model persist under moderate fields , which tilt spins by an angle proportional to the small field amplitude ( _ canted nel phase _
@xcite or _ canted stripe phase _
@xcite ) . to allow for this tilting ,
the untilted spins ( for @xmath56 infinitesimal ) should not be collinear with the field @xcite , and the staggered magnetization of the selected ordered pattern has an angle with @xmath57 as close to @xmath58 as possible .
we first discuss the effect of a small field in the nel phase , which allows to precise what we will understand in the following by ( an)isotropy in the magnetic response . in the pure heisenberg model , both classical and quantum cases ,
spontaneous breaking of @xmath6 symmetry in the thermodynamic limit leads to an anisotropic zero - field susceptibility . on the other hand ,
the ground state of a finite cluster is known to be a singlet , which does not break @xmath6 symmetry . in order to perform a classical treatment of the heisenberg model in absence of spontaneous symmetry breaking , one can consider for instance a @xmath6 symmetric combination of classical nel states - such a state can mimic , from an experimental viewpoint , the response of a polycristalline sample .
the magnetic response of this state is isotropic in the sense that a small field results in a magnetization of amplitude @xmath59 along the field , and @xmath59 does not depend on the field direction .
when considering the kitaev - heisenberg model for small finite @xmath60 , with the classical approach above , one can find for any field direction a nel pattern with spins perpendicular to it ; within this approximation , the magnetic response is again isotropic , despite the absence of @xmath6 symmetry due to small kitaev interactions .
however , the combination of quantum fluctuations and small kitaev interactions have been found to select one of the cubic axes as an easy spin axis and open thereby a small spin gap @xcite ; this will lead to small anisotropies in the spin susceptibility .
in the stripe phase , even at the classical level only six ordered patterns are allowed at zero field . for a generic field
@xmath57 with @xmath61 , @xmath62 and @xmath63 taking three distinct values , the three allowed stripe orientations are energetically differentiated , depending on the field direction .
for instance , in a small field @xmath64 $ ] , the two @xmath15-stripe patterns are favored ( along with the @xmath33-stripe patterns ) , since they allow for the field to tilt spins by a small angle @xmath55 [ see fig .
[ stripe2](b ) ] , such that the uniform magnetization develops an @xmath32-component .
this effect is evidenced by structure factors characteristic of the ordered phase - in the case of fig .
[ diag ] ( with @xmath65 ) , in a small @xmath66$]-oriented field the @xmath32-stripe pattern becomes more favorable than the @xmath15-stripe pattern : @xmath67 jumps from approximatively @xmath68 at @xmath69 to @xmath2 at small finite @xmath63 , while @xmath70 ( not shown ) drops to almost @xmath71 .
as for the magnetization per site @xmath72 , it grows linearly with @xmath63 in this regime .
the tilting angle @xmath55 can be determined simply in the classical approach , where the energy per unit cell as function of @xmath55 , for given @xmath60 and @xmath63 , is : e_s(h_z,)=-(1+)-2h_z+4(1-)^2 , and its minimization w.r.t .
@xmath55 gives a magnetization per spin @xmath73 , in relatively good agreement with the linear behavior observed at small field ( see fig .
[ diag ] for the case @xmath65 ) ; also , the predicted @xmath60-dependence of the ( uniform ) susceptibility per spin , ^z=()_h=0= , [ chz ] is in good agreement with the numerical results in fig .
[ khiz1 ] .
note that the direction of the field is relevant in the stripe phase : a field @xmath74 $ ] differentiates stripe patterns , and selects two of them where spins are exactly orthogonal to the field as @xmath75 ; in contrast , in a @xmath76 $ ] field the three stripe directions remain equivalent ( thus @xmath77 for @xmath78 ) ; also , none of the uncanted stripe patterns is orthogonal to a @xmath29$]-field , so the magnetic response for this field direction necessarily differs from the response to a @xmath66$]-field .
( color online ) magnetization @xmath59 ( per site ) and @xmath79stripe structure factor @xmath67 versus @xmath63 for various values of @xmath60 and @xmath80 .
dotted lines : perturbative results for @xmath59 in the polarized phase .
inset : phase diagram in the @xmath81 plane ; continuous transition lines are interpolations of points from @xmath22 data and dashed lines are perturbative estimates.,width=283 ]     ( color online ) evolution of correlations @xmath82 with @xmath60 on a @xmath22 cluster at @xmath69 , for various relative positions @xmath83 with @xmath17 and @xmath84 indicated in the inset .
the red dot indicates the value @xmath85 of the @xmath86 correlation at @xmath31 , which coincides with an analytical result from ref .
.,width=294 ]     ( color online ) susceptibility per spin @xmath87 , evaluated as @xmath88 for small @xmath89 , and its equivalent @xmath90}$ ] for a @xmath29$]-oriented field , on clusters @xmath22 or @xmath23 .
shown also the quantity @xmath91 for @xmath23.,title="fig:",width=294 ] +    we now turn to the spin liquid regime , corresponding to the weakly perturbed kitaev model , and describe the effect of perturbations to the system .
our main focus is on the @xmath92 correlations . in the pure kitaev model ( @xmath31 , @xmath69 )
, these correlations are non - zero only when @xmath93 are nearest neighbors on a @xmath15-bond @xcite , where they take a value @xmath94 which we can see in fig .
[ diag2 ] ( the finite - size correction to this quantity , i.e. to the ground state energy , is negligible even for the @xmath22 sites cluster ) .
the gapless nature of excitations persists in a small axial field @xmath63 ; in contrast , a @xmath29$]-oriented field opens a gap resulting in a topologically ordered phase @xcite . yet even in the first case , the field causes peculiar @xmath92 correlations with spatial oscillations perpendicular to @xmath15-bonds and a power - law decaying envelope , as shown for @xmath31 in ref .  .
some of these trends appear in the spin - spin correlations displayed in fig .
[ cormag](a ) . due to the small size
, one can not identify a power - law behavior ; but one notices that correlations decrease with distance , much more slowly in the direction parallel to @xmath15 bonds than in the direction perpendicular to them ( i.e. parallel to @xmath95 ) . given the periodic boundary conditions on the @xmath22 lattice ,
this is in agreement with the period of three unit cells for spatial oscillations predicted in ref .  .
the liquid phase is stable with respect to small heisenberg interactions ( for @xmath96 at zero field @xcite ) .
these interactions induce correlations between further neighbors , of different type than those induced by a small field [ compare fig .
[ cormag](a ) and ( b ) ] .
the spatial alternation of sign on fig .
[ cormag](b ) reminds of the @xmath15-stripe pattern , with the difference that here these correlations are short - ranged .
a proof of this last point has been given in ref .
for @xmath69 , and indeed absolute values @xmath97 of correlations shown in fig .
[ diag2 ] seem to decrease rapidly  possibly exponentially  with distance , at least in the spatial direction orthogonal to @xmath15 bonds . on nearest neighbors ,
these correlations are proportional to @xmath98 ( not too close to @xmath35 ) . with both perturbations at play , i.e. when both @xmath98 and @xmath63 are finite and small , their effects should add up ; therefore , and since the correlations induced by heisenberg terms decay exponentially , one should observe a correlation pattern with a power - law tail as in the case of @xmath99 , but modulated as in a @xmath15-stripe pattern for short distances , at least for @xmath100 .
eventually , for any value of @xmath60 , a sufficiently high field drives the system into a _ polarized phase _ with spins parallel to @xmath101 . in the case of a @xmath66$]-field , by considering transverse interactions @xmath102 and @xmath103 as perturbations [ the unperturbed ground state being fully polarized for @xmath104 , within second order perturbation theory one finds m_0 = _ z 1- , in good agreement with numerical data for @xmath105 ( see dashed curves in fig .
[ diag ] ) .      when a magnetic field @xmath63 is applied to a system which is initially stripe- or nel ordered , aside from the spin canting in the af ordered state , the excitation energy of a competing polarized state decreases with @xmath63 so that this state becomes the ground state above a critical field @xmath106 .
this field can be evaluated by estimating the energies ( both per unit cell ) @xmath107 of the canted stripe phase and @xmath108 of the polarized phase , as a function of field - then the equation @xmath109 determines @xmath106 . for both quantities ,
we add to the classical value a term ( last terms in both lines of eq .  [ ecantpol ] ) accounting for quantum fluctuations , considered in second order perturbation theory , and evaluated respectively from the uncanted stripe- or the fully polarized state : e_s(h_z)&=&-(1+)- - , + e_fz(h_z)&=&-2h_z+3 - 5- .
[ ecantpol ] we compare this estimate of @xmath106 with numerical estimates obtained from the @xmath63-dependence of the stripe structure factor @xmath67 and the magnetization @xmath59 for @xmath22 and @xmath23 clusters  both quantities show very weak size - dependence ( see fig .  [ diag ] ) and their derivatives with respect to @xmath63 are maximized for values very close to each other . for @xmath65 , this numerical estimation
gives @xmath110 compared to a value @xmath111 found from the criterion @xmath109 .
both types of estimates are shown as a function of @xmath60 , in the phase diagram in the inset of fig .
[ diag ] . for @xmath112 , the magnetization curve @xmath72 on a finite cluster at fixed @xmath60 ( @xmath113 in fig .
[ diag ] ) shows a sharp increase in a field range centered on @xmath114 ; this value is used to locate the transition to the polarized phase ( our precision on @xmath115 is limited here by larger finite size effects in the liquid state than in ordered phases ) .
this can be understood by modeling the energy per unit cell for the liquid phase by e_l=3(1 - 3)^^_^(=1,h=0)-^z h_z^2 , and looking for the field value such that @xmath116 .
this gives a value @xmath117 , which increases slightly with @xmath98 .
similarly , the position of the stripe - liquid transition , which at zero field occurs at @xmath35 , is almost insensitive to a small field @xmath63 .
the schematic phase diagram in fig .  [ diag ] has a similar structure as its counterpart for a @xmath29$]-oriented field @xcite , except that we do not see a multicritical point at @xmath118 .
this may be connected with the absence of @xmath48 symmetry in the hamiltonian eq .
[ ham ] in a @xmath66$]-field ; moreover , the nature of the spin liquid for @xmath112 in a small axial field differs from that of the topologically ordered phase in a small @xmath29$]-field , mainly in the structure of the low - energy spectrum .
we analyze now in more detail the zero - field susceptibility : it has unusual properties due to a special spin symmetry of the model , which is lower than @xmath6 . in the stripe phase ( @xmath119 ) , if we consider the general case of a @xmath120$]-oriented field , that is @xmath57 parallel to @xmath121 , the corresponding susceptibility @xmath122}$ ] behaves as @xmath123 , where @xmath124 is a rational number which can be found from the classical picture of a spin - flop mechanism .
one can see in fig .
[ khiz1 ] that susceptibilities along different axes follow well this type of law in the stripe phase .
the coefficient @xmath124 is maximized for a @xmath66$]-field ( @xmath125 , see eq .
[ chz ] ) , but also for a @xmath126$]-field , and minimized for a @xmath29$]-field ( @xmath127 ) .
note here that the @xmath29 $ ] axis corresponds , in the context of susceptibility measurements in ref .
, to the @xmath128 axis of na@xmath3iro@xmath4 , i.e. perpendicular to the ir - honeycomb layers [ see fig .
[ struc](b ) ] . concerning the susceptibility in the @xmath129 plane
, this model predicts that  in the stripe phase
@xmath130 takes a value between @xmath131}=\frac{1}{4(1-\alpha)}$ ] and @xmath132}=\frac{1}{5(1-\alpha)}$ ] .
an average over the field orientations in the @xmath129 plane gives @xmath133 .
no such simple expressions for the susceptibility are available for the liquid phase ; but there , despite large finite - size effects ( compared to the stripe phase ) one notices that for a given axis susceptibilities tend to increase with @xmath60 , with again @xmath90 } \le \chi^{[001]}$ ] . as for the nel phase
, a classical approach does not predict there any anisotropy in the static susceptibility ( see related discussion in section [ phas ] ) , and the anisotropy resulting from quantum fluctuations at small finite @xmath60 should be smaller than in the stripe phase  but the nel phase is obtained for values of @xmath60 likely unrealistic for na@xmath3iro@xmath4 @xcite and also not supported by experiment@xcite .
the experimental data of ref .
reveal an anisotropy in the static susceptibility as well , but the measured susceptibility is significantly larger along the @xmath128 axis than in the @xmath129 plane , in contrast to our model prediction @xmath134 .
this suggests that other factors could contribute to the susceptibility anisotropy@xcite .
for instance , anisotropic contributions are expected to arise from the van vleck susceptibility @xmath135 , and also via the g - value anisotropy in the ground - state kramers doublet , i.e. , @xmath136 , due to lattice distortions lowering the octahedral symmetry .
these effects involve the inclusion of higher , excited spin - orbital , states .
nevertheless , we believe that the kitaev - heisenberg model describes some characteristic features of layered iridates : for instance , recent susceptibility measurements@xcite on polycrystals of li@xmath3iro@xmath4 and of na@xmath3iro@xmath4 showed that the zero - temperature susceptibility was about twice higher for the first ( li ) compound than for the second . assuming that both compounds can be described by the kitaev - heisenberg model and are in the stripe phase ,
the value of @xmath60 for the li compound should be significantly larger than for the na compound - probably @xmath137 if one compares the paramagnetic curie temperatures of both compounds and follows theoretical results of ref .  .
this is actually in agreement with the law @xmath138 which we find for any field direction ( and thus for polycrystalline samples ) : since by assumption @xmath139 , the observed ratio@xcite @xmath140 would imply @xmath141 .
this suggests that the li - compound is quite close to the stripe - liquid transition .    in the following
, we shall confine our study to the kitaev - heisenberg model and explore the role of non - magnetic vacancies , unavoidably present in real compounds , and which , in general , can significantly alter the properties of frustrated magnetic systems @xcite .
we now study the effect of a single non - magnetic vacancy on the magnetic response of the system , especially at low - field . as previously
, we will address separately the different phases of the model .
( color online ) ( a , b ) correlation functions @xmath92 with respect to site @xmath17 indicated by a black circle .
( c , d ) induced magnetization patterns @xmath142 due to the presence of a vacancy , on the @xmath22 cluster .
circle areas are proportional to the quantities shown .
parameters are respectively @xmath143 in ( a ) and ( c ) , while @xmath144 in ( b ) and ( d ) [ a small symmetry - breaking field @xmath145 was added in ( d)].,title="fig : " ]   ( color online ) ( a , b ) correlation functions @xmath92 with respect to site @xmath17 indicated by a black circle .
( c , d ) induced magnetization patterns @xmath142 due to the presence of a vacancy , on the @xmath22 cluster .
circle areas are proportional to the quantities shown .
parameters are respectively @xmath143 in ( a ) and ( c ) , while @xmath144 in ( b ) and ( d ) [ a small symmetry - breaking field @xmath145 was added in ( d)].,title="fig:",width=145 ] +   ( color online ) ( a , b ) correlation functions @xmath92 with respect to site @xmath17 indicated by a black circle .
( c , d ) induced magnetization patterns @xmath142 due to the presence of a vacancy , on the @xmath22 cluster .
circle areas are proportional to the quantities shown .
parameters are respectively @xmath143 in ( a ) and ( c ) , while @xmath144 in ( b ) and ( d ) [ a small symmetry - breaking field @xmath145 was added in ( d)].,title="fig:",width=137 ]   ( color online ) ( a , b ) correlation functions @xmath92 with respect to site @xmath17 indicated by a black circle .
( c , d ) induced magnetization patterns @xmath142 due to the presence of a vacancy , on the @xmath22 cluster .
circle areas are proportional to the quantities shown .
parameters are respectively @xmath143 in ( a ) and ( c ) , while @xmath144 in ( b ) and ( d ) [ a small symmetry - breaking field @xmath145 was added in ( d)].,title="fig:",width=141 ] +      we start by considering the liquid phase , and choose the field direction along the @xmath66 $ ] axis , for which features characteristic of the kitaev model are the clearest .
first , we show that the combination of an infinitesimal field and a spin vacancy induces a peculiar magnetization pattern around the vacancy .
the relevant quantity is the induced magnetization : ^z_j = ^z_j -m_0 , which gives , for a given field @xmath63 and site @xmath84 , the change of the magnetization at this site relative to its value ( @xmath59 ) in the vacancy - free system .
we relate the vacancy - induced magnetization pattern to the correlation function @xmath146 of the vacancy - free system .
the most striking feature of the magnetization patterns shown in fig .
[ cormag](c , d ) is the large magnetization on the @xmath15-neighbor of the vacancy ( referred to as site @xmath147 in fig .
[ difmz ] from now on ) .
the magnetization is significantly larger there than on any other site  and only on this site does it remain finite at vanishing field in the kitaev model : @xmath148 for @xmath31 and @xmath149 ( @xmath150 for @xmath22 ) . actually , in the thermodynamic limit ( tl ) the magnetization ( both the total one and @xmath151 ) should behave with field as@xcite @xmath152 ; our results are consistent with this prediction in the sense that @xmath153 decreases with increasing @xmath80 , so that it should vanish in the tl as expected . as for the logarithmic field - dependence ,
we stress that this feature pertains to the tl , while for finite systems we can only notice weak non - linearities in the magnetization curves , discussed later in more detail .    the magnetic response localized on the @xmath15-neighbor of the vacancy
is clearly connected with the fact that the correlation function @xmath154 is non - zero only between nearest - neighbors sharing a @xmath15-bond @xcite : if one spin is removed at site @xmath155 , the term @xmath156 obviously disappears from the hamiltonian and thus the spin @xmath147 is easily polarizable by the @xmath15-oriented field ( one can also interpret this behavior in terms of dangling majorana fermions @xcite ) .    yet , we already saw in section [ phas ] that a small field induces @xmath154 correlations beyond nearest neighbors , it makes sense to focus in more detail on the magnetization pattern around a vacancy at small finite @xmath63 .
[ cormag](c ) shows the variation of magnetization @xmath157 due to the presence of the vacancy .
indeed , the spatial variations of this quantity presents similarities with those of correlations in the vacancy - free system : @xmath157 is positive everywhere and decreases with increasing distance to the vacancy .
moreover , the spatial oscillations of correlation function appear in the magnetization pattern . indeed ,
if one compares values of @xmath158 at relative position either @xmath159 or @xmath160 from the vacancy [ see values @xmath161 and @xmath162 respectively in fig .
[ cormag](c ) ] , their respective ratio is comparable to that of corresponding correlations in the vacancy - free system .
while @xmath163 correlations scale as @xmath164 , the magnetization enhancement @xmath158 on further neighbors of the vacancy grows approximately linearly with @xmath63 ; in contrast , for @xmath165 , it decreases from @xmath153 at @xmath69 towards small negative values in the polarized phase ( see @xmath166 in fig .
[ difmz ] ) .
although the response observed here is obtained for a finite cluster size , we may expect that in a large system , the spin at site @xmath147 acquires a large magnetization due to the combined effect of the neighboring vacancy and the small field , and in turn polarizes its other neighbors in a way that the @xmath157 reproduces the above discussed correlation function pattern ( fig .
[ cormag](a ) ) .
+ if one considers the effect of a vacancy on the magnetic response of the system in a field @xmath63 , a relevant quantity to consider is @xmath167-(n-1)m_0(h_z)$ ] , with @xmath59 the magnetization _ per site _ of the vacancy - free system , and @xmath168 $ ] is the total magnetization of the system with @xmath169 spins plus a vacant site .
this quantity , shown in fig .
[ difmz ] , vanishes in the infinite field limit ( where all spins are polarized ) , but also at small field in the tl . in finite systems ,
the previously discussed finiteness of @xmath153 goes along with a deviation from the expected logarithmic behavior of @xmath170 ; yet , we notice on large enough systems ( @xmath22 ) a non - linearity of the magnetization curve in the small field range @xmath171 , the slope @xmath172 decreasing moderately with increasing small @xmath63 .
but when @xmath63 becomes close to the value @xmath115 of the liquid / polarized transition for the undoped system , the slope @xmath173 increases strongly , so that @xmath174 peaks near @xmath175 .
the exact position of the peak depends slightly on the cluster size , and the peak amplitude grows with increasing system size ( roughly linearly : compare amplitudes for @xmath24 and @xmath22 on fig .
[ difmz ] ) .
the existence of this peak can be understood in the following way : not only does the presence of a vacancy cause an increase of the magnetization in the liquid phase , but it shifts slightly the liquid / polarized transition towards smaller fields @xmath176 ( the difference between both fields being likely proportional to the vacancy concentration in a diluted limit ) , and the peak signals a field range @xmath177 .
this behavior can loosely remind of spin polarons formed around mobile holes in a paramagnetic background , so that the temperature of transition to the ferromagnetic ( fm ) phase shifts up with increasing hole density@xcite . here , instead of temperature
the parameter driving the transition is the magnetic field , and moreover a vacancy is here not mobile , but its effect is as well to polarize its surroundings more and more when approaching the fm phase ; consequently , at small finite vacancy concentration , the transition field should be slightly smaller than in the vacancy - free case .
spin - vacancy induced magnetization @xmath174 as a function of magnetic field @xmath63 , for different @xmath60 and cluster sizes .
shown also is a partial contribution ( open circles ) of the spin at site @xmath147 , being the @xmath15-neighbor of the vacancy , for @xmath31 and @xmath22.,width=294 ]    an analogy between the spin correlations in the vacancy - free system and the magnetization pattern around a vacancy can also be noticed in the case where the kitaev hamiltonian is perturbed by small heisenberg interactions [ see fig .
[ cormag](b ) and fig .
[ cormag](d ) ] .
note that in fig .
[ cormag](d ) , we also included an infinitesimal field @xmath63 , the effect of which is only to break the kramers degeneracy related to the odd number of spins .
there , both the correlations and the vacancy - enhanced magnetizations @xmath158 decrease rapidly with distance to the vacancy .
moreover the sign of these quantities alternates consistently with the pattern of @xmath15-oriented stripes shown in inset of fig .
[ slrstripes ] , and their decrease is much faster when going perpendicularly to stripes than along stripes : for the magnetizations @xmath178 ( or equivalently @xmath158 , since @xmath179 ) , one can again compare in fig .
[ cormag](d ) the absolute values at positions @xmath160 ( @xmath180 ) and @xmath159 ( @xmath181 ) from the vacancy .
the staggered @xmath15-bond magnetization , which measures the @xmath15-stripe order parameter , m^z(_z)=_r e^i _ z r ^z_r , a+^z_r , b , is shown in fig .
[ slrstripes ] .
it also evidences that , in the liquid phase , the short - range stripe pattern is induced around the vacancy by small heisenberg interactions .
starting from the @xmath31 limit , the magnetizations at further neighbors increase ( in absolute values ) significantly with increasing amplitude of heisenberg interactions , i.e. with decreasing @xmath60 .
we expect that the extensive quantity @xmath182 diverges in the tl for @xmath183 .
moreover , we note that the vacancy - induced magnetization @xmath174 , for fixed small field @xmath63 , is reduced by the presence of heisenberg interactions .
it is significantly smaller for @xmath113 than for @xmath31 , as the tendency to stripe order around the vacancy , at small field , competes with the aligning effect of the field .
consequently , the peak in @xmath184 shifts slightly towards higher fields with decreasing @xmath60 : the effect of heisenberg interactions is to increase the position @xmath185 of transition to the polarized phase ( as for @xmath115 in the vacancy - free case , see fig .  [ diag ] ) .
we briefly comment here on the fact that in the polarized phase ( and unlike in small field ) , the vacancy - induced magnetization change @xmath166 at the @xmath15-neighbor is negative .
this can be understood with the perturbative treatment of transverse interactions described in section [ phas ] ; but here , to evaluate the quantum correction to the magnetization at a given site , we consider the off - diagonal ( perturbing ) couplings only on neighboring bonds  see ref .   for a similar approach . within this approach ,
the magnetization at this site is @xmath186 ; and the predicted value for @xmath166 is of the same order as seen in fig .
[ difmz ] and it decreases in a similar way as @xmath63 increases , sufficiently far away from the transition field to either the liquid or the stripe phase .
qualitatively , one can interpret this by the fact that at high field , fluctuations from a fully - polarized state are effectively stronger at the vacancy s @xmath15-neighbor than elsewhere , since the absence of ferromagnetic interaction on the missing @xmath15-bond makes it easier to flip this spin than other spins .
with a similar reasoning one can explain that , in the same field range , the magnetization at other nearest neighbors of the impurity is enhanced compared to @xmath59 .
if , instead of a field along an easy axis , a field in the @xmath29 $ ] direction is applied , one can assume@xcite that the magnetic response is approximately a linear combination of responses to each field component : thus each neighbor of a vacancy acquires an anomalously large magnetization - the neighbor spin at site @xmath84 , that would be connected to the missing spin by a @xmath187-bond , has a magnetization @xmath188 , and in turn brings a smaller polarization , also parallel to @xmath189 , to other spins in its vicinity - mainly on neighboring sites which are close to the spatial axis in the prolongation of the missing @xmath187-bond . hence
, three domains should coexist in the system , each with a different magnetization direction .
note that a vector @xmath190 $ ] can be defined and used to characterize the polarization of spins around the vacancy .
qualitatively , @xmath191 should point in a direction closer to the easy axis @xmath187 for which @xmath189 is the largest than to other easy axes ; but it can have three non - zero components simultaneously .      in the stripe - ordered phase , among the 6 allowed stripe patterns ( 3 possible stripe directions times 2 spin orientations ) , a vacancy in combination with an infinitesimal field @xmath192 selects the @xmath15-stripe pattern for which @xmath193 up to small quantum corrections .
this selection is due to the necessary imbalance between the number of spins respectively aligned and antialigned with the field ( see inset of fig .  [ slrstripes ] ) .
this state is stable versus increasing @xmath63 , until the field reaches a critical value ( or anisotropy field @xmath194 ) , where canted stripe patterns  those with stripe directions along @xmath32 or @xmath33  become more favorable .
this effect is similar to the local rotation of the staggered magnetization close to a vacancy in @xmath6 antiferromagnets @xcite . here
, the main difference is that there is only a finite number of stripe patterns which can not be continuously connected to each other without going out of the ground state manifold ; consequently , the system has to choose between allowed stripe patterns , tilted or not by the field , and @xmath191 is forced to point exactly along one easy axis @xmath187 , corresponding to ( one of ) the most favorable pattern(s ) . to understand the existence of a finite anisotropy field , one has to compare the energy of the @xmath15-stripe pattern ( _ _ untiltable__@xcite in a @xmath66$]-field ) to that of the _ tiltable _ @xmath32- and @xmath33-stripe patterns .
the former , thanks to its total magnetization @xmath195 , acquires locally an energy gain @xmath196 with respect to the zero - field value but it has zero susceptibility ( @xmath197 ) thus no contribution @xmath198 to its energy ; in contrary , the field lowers the energy of the latter patterns by an extensive quantity @xmath199 .
this comparison gives @xmath200 for @xmath201 , not far from the position of the jump in @xmath202 in fig .
[ slrstripes ] .    in the vicinity of the stripe - liquid transition
, the magnetization pattern induced by a small @xmath66$]-field ( @xmath203 ) around a vacancy can be of special interest : assuming that vacancies are diluted enough such that their mutual interactions can be neglected , the patterns expected on both sides of the transition are clearly distinct from each other . in the stripe phase ,
all neighboring spins of the vacancy have a large , almost field - independent magnetization , with two of them aligned with the field and the third one antialigned with @xmath57 .
in contrast , in the liquid phase only one nearest neighbor spin is significantly polarized ( with a strongly field - dependent polarization ) and aligned with the field ; and the two others are weakly polarized opposite to the field . thus , if one would have a local probe to determine the ( site - resolved ) magnetization in the neighborhood of a vacancy  this probe could be @xmath204 nmr for instance  the measured pattern could indicate on which side of the transition the system is .
we have also seen that a vacancy can block the spin canting mechanism in the stripe phase , under a small field in ( or close to ) the @xmath66 $ ] axis .
this phenomenon is absent for a field in the @xmath29 $ ] direction , where all stripe patterns remain equivalent and continuously tiltable by the field .
consequently , in a finite system , the situation concerning the susceptibility anisotropy at zero - field : @xmath205}$ ] , is now opposite to the previously discussed case without vacancies where we found @xmath90}<\chi^z$ ] .
concerning now the thermodynamic limit with a finite concentration @xmath206 of vacancies , if one assumes that @xmath206 is large enough to influence the magnetic response at low temperatures but small enough so that vacancies behave independently from each other , a similar effect could account for the experimentally observed @xmath207 .
however , given that in a small field , within the kitaev model , the magnetization enhancement caused by two vacancies residing on the same sublattice is much larger than twice the magnetization enhancement caused by a single vacancy @xcite , one should take care about cooperative effects between the vacancies .
( color online ) the @xmath15-stripe order parameter @xmath208 computed on the @xmath22 cluster , as a function of @xmath60 at different small fields @xmath63 .
right - inset : the magnetization pattern @xmath209 around the vacancy for @xmath201 and @xmath210 .
left - inset : total magnetization @xmath170 as a function of a field @xmath63 applied in the liquid ( @xmath113 ) and in the stripe ( @xmath201 ) phases , with a magnetization jump at @xmath211.,title="fig:",width=302 ] +
we turn now to the case where several vacancies are present in the system ; we first consider their effect on the magnetization curves ; then , we will focus on their effective interaction in the stripe phase and show how a vacancy pair can select an ordered pattern .
it is instructive to compare the effect of vacancies or pairs of vacancies on magnetization curves , for given field direction and system size . to do this comparison we choose a @xmath66$]-field and consider the periodic @xmath22 cluster , but discuss as well implications for the thermodynamic limit .
magnetization curves @xmath202 without or with vacancies , in different phases of the kitaev - heisenberg model , are shown in fig .
[ magnsl ] .
_ in the nel phase _
, although the total magnetization @xmath212 is not a conserved quantity in presence of small kitaev interactions ( @xmath213 ) , the magnetization curve @xmath202 on finite clusters displays a succession of plateaux ( almost flat for @xmath214 ) , as seen in fig .
[ magnsl](a ) .
they are separated by jumps @xmath215 corresponding to a flip of one spin @xmath2 . at fixed @xmath60 the step - like increase of the magnetization per site
@xmath72 disappears only in the tl , the steps becoming smaller and closer to each other as @xmath80 increases .
moreover the steps disappear , at fixed cluster size , as @xmath60 gets closer to @xmath216 .    within this phase ,
one is tempted to define an effective spin @xmath2 , @xmath217 , associated to the vacancy , such that @xmath218 corresponds to @xmath217 pointing along @xmath219 ; with two vacancies , their effective @xmath220 s add up if vacancies are on the same sublattice ( e.g. , n.n.n .
case ) and cancel each other if on opposite sublattices ( n.n .
case ) . from the behavior of
@xmath221 one obtains similar conclusions concerning effective couplings between other components of @xmath222 .
one could therefore think that an effective interaction between two vacancy spins , if it can be defined , should be @xmath6-invariant , and possibly of heisenberg type @xmath223 ; yet , the gapless nature of the undoped system and the continuous behaviour of magnetization curves in the tl make it unclear how to formulate the effect of vacancies in terms of effective spins .
_ in the liquid phase _ , in contrast , the magnetization does not show any plateau , but vacancies cause peculiar features at small fields : we have already seen that , compared to the zero - vacancy curve , the slope of @xmath202 at small field is significantly increased by the presence of one vacancy in the system .
if two vacancies are present , their effect depends strongly on their respective sublattices ( see fig .  [ magnsl](c ) ) : if they are on opposite sublattices the @xmath202 curve is , as in the vacancy - free case , linear for small @xmath63 [ and @xmath224 is only slightly modified ] ; but if they are on the same sublattice , this slope is greatly increased  for a large range of @xmath63 the vacancy - induced increase in @xmath170 is more than twice that of the single - impurity case . again , even though we can not identify a logarithmic behavior @xmath225 as derived in ref .   for @xmath31 , we see that vacancies , as soon as they are not equally distributed on both sublattices , act like _ partially - unbound moments _ causing non - linearities in the magnetization properties at low field .
_ in the stripe phase _ , as in the liquid phase , kitaev interactions are strong enough to prohibit the existence of plateaux in the magnetization curves @xmath202 , except at small fields in presence of vacancies : with one vacancy , the previously discussed @xmath193 plateau is again suggestive of an effective spin @xmath2 of the vacancy ; and with two vacancies on the same sublattice , the @xmath226 plateau seen on fig .
[ magnsl](b ) ( pair p1 ) extends up to a field @xmath227 , that is more than twice the anisotropy field @xmath194 for a single vacancy .
this indicates a cooperative effect of both vacancies , which stabilize the @xmath15-stripe pattern more than if they behaved independently from each other . to describe this cooperative behavior
, we will , similarly as in section [ sosf ] , analyze the effect of vacancies on the energies of the various stripe patterns .
( color online ) magnetization curves @xmath202 on the cluster @xmath22 with 0 , 1 , or 2 vacancies . the labels @xmath228 correspond to relative positions of 2 vacancies as shown in panel ( c ) .
panels ( a),(b ) , and ( c ) correspond to @xmath229 , respectively.,title="fig:",width=306 ] +   ( color online ) magnetization curves @xmath202 on the cluster @xmath22 with 0 , 1 , or 2 vacancies .
the labels @xmath228 correspond to relative positions of 2 vacancies as shown in panel ( c ) .
panels ( a),(b ) , and ( c ) correspond to @xmath229 , respectively.,title="fig:",width=294 ] +     ( color online ) structure factors ( a ) @xmath230 and ( b ) @xmath67 at @xmath231 calculated for various respective positions of 2 vacancies @xmath232 ,
@xmath233 , @xmath234 , indicated in the insets . also shown
is the structure factor of the vacancy - free system ( open circles ) , for which @xmath235 .
all results are obtained on the @xmath22 cluster.,title="fig:",width=272 ] +   ( color online ) structure factors ( a ) @xmath230 and ( b ) @xmath67 at @xmath231 calculated for various respective positions of 2 vacancies @xmath232 , @xmath233 , @xmath234 , indicated in the insets .
also shown is the structure factor of the vacancy - free system ( open circles ) , for which @xmath235 .
all results are obtained on the @xmath22 cluster.,title="fig:",width=272 ]      we focus now on the stripe - ordered phase , and will explore how at @xmath236 , a pair of defects can , depending on their relative position , favor or disfavor one of the three stripe directions . in terms of symmetries , this can be understood as a breaking of @xmath47 symmetry related to the spatial orientation of stripes , while the @xmath45 symmetry ( @xmath237 ) remains unbroken .
remarkably , as shown in fig .
[ mstripy_2 ] , the selected stripe pattern for a given vacancy pair depends on the coupling constant @xmath60 . + this effect is easy to understand for nearest - neighbor defects on a @xmath238bond ( configuration @xmath232 : see fig.[mstripy_2 ] ) ; classically , the @xmath238stripe pattern has an energy lower by e^z_1 = 2 - 4than the one of other stripe patterns .
indeed , in the @xmath15-pattern only one fm bond - here a @xmath15-bond , with a classical contribution to the energy @xmath239 - is lost due to vacancies , instead of two for other patterns ( e.g. two @xmath32-bonds , in a @xmath32-stripe pattern ) , and for @xmath240 fm interactions contribute more than af ones ( whose classical contribution to the energy is @xmath241 ) . for @xmath242 the situation
is reversed , which explains that @xmath230 becomes small and @xmath243 .    for second neighbors ( case @xmath233 on fig .
[ mstripy_2 ] ) , stripe patterns are classically degenerate , but this degeneracy is lifted by quantum fluctuations .
we evaluate in second order perturbation theory the contribution to the energy of various bonds in which at least one site is a neighbor of the vacancy ; actually only af bonds give a contribution in this approximation .
for instance fluctuations on a bond belonging to the same hexagon as the two impurities yield a contribution @xmath244 $ ] to the energy of the @xmath15-stripe state . collecting contributions from all relevant bonds ,
the energy of the @xmath15-stripe state ( relative to that of other stripe orientations ) is : e^z_2=2(1 - 2)^2 ( + - ) , which is positive for @xmath60 in the vicinity of @xmath2 , and changes sign at @xmath245 .
this explains the drop of the @xmath15-stripe structure factor ( dotted curve on fig .
[ mstripy_2](a ) ) , from large values for @xmath246 $ ] , to almost zero for smaller @xmath60 .
the smaller values of @xmath247 in the vicinity of @xmath248 are related to the smaller energy scale separating the different stripe patterns : the above expression , with its prefactor @xmath249 , indicates that the 3 patterns become degenerate at @xmath248 .    a similar energy comparison between different patterns for third neighbor vacancy pairs ( configuration @xmath234 : see fig .  [ mstripy_2 ] )
predicts a change of stripe orientation at @xmath250 between @xmath238stripes ( @xmath251 ) and @xmath32-stripes ( @xmath252 ) , following from the perturbative result : e^z_3 = , whereas this orientational change is seen on @xmath22 at a slightly smaller value @xmath253 .
notice that here ( and more generally for pairs of vacancies at further neighbor sites which are not related to each other by any reflection symmetry of the honeycomb lattice ) the three stripe orientations are all non - equivalent to each other .
finally we note that at zero - field , the conclusions above hold also for pairs of neighbor vacancies with other relative positions : for instance with next - nearest neighbor vacancies separated by a @xmath32- and a @xmath15-bond and @xmath254 , according to the previous reasoning the favored pattern is a @xmath33-stripe pattern .    in presence of a small magnetic field , the @xmath45 symmetry characterizing each stripe orientation ( e.g. , @xmath15-stripes ) is broken .
moreover , similar effects occur as in the previously discussed single vacancy case : if the most favorable pattern in an infinitesimal @xmath66 $ ] field ( e.g. , @xmath15-stripe pattern around a vacancy pair in configuration @xmath232 , for @xmath240 ) is untiltable for this field direction , as the field amplitude @xmath63 gets larger than a critical value @xmath255 this pattern becomes less favorable than other patterns which are tiltable . in the limit of large systems @xmath256 ,
a classical estimate of @xmath255 gives a result @xmath257 which , for @xmath201 and @xmath22 , has a value @xmath258 comparable to the observed position of the magnetization jump .    the relevant quantity for the discussion of either low temperature properties or magnetic ordering at small vacancy concentrations is the energy scale @xmath259 , namely the energy which favors a particular stripe pattern in the vicinity of a vacancy pair  @xmath260 standing for the number of bonds separating both vacancies .
far away from the stripe / liquid transition ( say , for @xmath261 ) @xmath259 decreases rapidly with @xmath260 ( e.g. , @xmath262 for @xmath263 ) .
closer to @xmath35 , larger quantum fluctuations make the stripe - orientation selection by @xmath264 neighbor vacancy pairs as efficient as by nearest neighbor vacancy pairs .
several consequences can result from this : ( i ) in the simpler case with only two vacancies , the energy lowering of some stripe patterns for @xmath265 also stabilizes these patterns with respect to the spin liquid , which leads to a slight shift of the stripe / liquid transition towards larger @xmath60  possibly @xmath266 to @xmath267 depending on the relative position of vacancies , as can be seen in fig .  [ mstripy_2](a ) .
( ii ) more importantly , if one focuses on systems with substantial vacancy concentrations , say , @xmath268 , neighboring vacancy pairs as described in fig .
[ mstripy_2 ] are abundant enough to select locally stripe orientations that will frustrate each other at the global scale .
( iii ) as a result , the system could display a glassy behavior , especially in the vicinity of the stripe - liquid transition .
the resulting glassy phase would consist of stripe - ordered nanoscale domains , with a vacancy pair in each domain , selecting its stripe orientation . in order to perform a modeling of such a system and to reproduce a possible spin - glass behavior
, one should obviously need to consider systems of ( at least ) hundreds of ir atoms , and to take into account the energy cost of domain walls between regions of different stripe orientations .
moreover , obviously not all vacancies are grouped into distinct first , second or third nearest neighbor vacancy pairs ; in particular the effect of unpaired vacancies , even at zero external field , could be subtle , although one can imagine single vacancies pinning walls between domains of different stripe orientation .
we have considered here several important aspects of the kitaev - heisenberg model , which , despite its simple structure , displays in its @xmath1 phase diagram three distinct phases , two of which are antiferromagnetically ordered while the third one is a spin liquid .
we employed numerical methods ( exact diagonalization ) for our investigation of the properties of the model in a magnetic field and the effect of vacancies on the magnetic properties .
the numerical methods were complemented by classical and perturbative analysis particularly in the frame of the stripe phase .
we showed how the ground state phases evolve in a magnetic field oriented along the @xmath66 $ ] direction ( one of the local easy axes of the kitaev - heisenberg hamiltonian ) .
we found that the three phases are robust to small fields .
next , we analyzed the influence of the field direction .
the magnetic response is qualitatively distinct in the different phases . in the stripe phase ,
the magnetization can be described by a spin canting mechanism , which also provides an understanding of the anisotropy in the numerically obtained susceptibilities . in the liquid phase
, we compared the effect of small magnetic fields and of small heisenberg interactions , each considered as a perturbation of the kitaev hamiltonian .
such perturbations are shown to induce specific patterns in the spin - spin correlations beyond nearest neighbors , while at larger fields a transition into the spin - polarized phase occurs .
we determined the phase diagram of the model in a @xmath66$]-field , with good agreement between numerical results and analytical estimates .
for the latter analysis , we took into account both the spin canting and the quantum fluctuations in the ordered and polarized phases .
this phase diagram has similarities with its counterpart for a @xmath29 $ ] field orientation @xcite , however , a major difference between them is that for a @xmath66$]-field we do not see any quantum critical point in the vicinity of the stripe - liquid phase transition .
we presume that this is related to the absence of a field - induced gap in the low - energy spectrum of the liquid phase for a @xmath66$]-field@xcite .
subsequently , we focused on the effect of non - magnetic vacancies on the spin correlations and the magnetic response of the model . for a single vacancy in the liquid phase
, we found that the spatially anisotropic magnetization pattern around the vacancy is related to the anisotropy of spin correlations in the vacancy - free system .
up to finite - size effects inherent to our approach , these results are also in agreement with predictions of a non - linear response to a small @xmath66$]-field , whose effect is mostly seen on one specific neighbor of the vacancy . in the stripe phase ,
a vacancy coupled to a small [ 001]-field breaks both the @xmath45 ( time reversal ) and the @xmath47 ( orientational ) symmetry characterizing the ground state manifold ; a single stripe pattern is selected , which is not continuously tiltable by the field in contrast to the vacancy - free situation .
having in mind the description of systems at finite vacancy concentration , we considered the interplay of two vacancies sitting either on nearest or further neighbor sites .
again the magnetic response to an easy - axis field depends strongly on the nature of the phase , but it also depends sensitively on the relative position of vacancies . in the liquid phase a non - linear behavior , much stronger than for the single - impurity case , is seen if the two vacancies are on the same sublattice , while the response remains linear if they reside on different sublattices . in the stripe phase , vacancy pairs select a specific stripe orientation even at zero field .
this selection mechanism , depending on the value of coupling parameter @xmath60 , can be well understood by a perturbative approach evaluating locally the effect of quantum fluctuations .
the response to a small field depends on whether the selected stripe pattern is tiltable or not for this field direction@xcite .
although the strongest selection effect is caused by nearest neighbor vacancy pairs for a wide range of @xmath60 , further neighbor pairs have to be considered as well in the vicinity of the stripe - liquid transition .
therefore we conclude that the intrinsic randomness of vacancies , in combination with the selection mechanism induced by vacancy pairs , may lead to spin - glass behavior in the stripe regime . there , a vacancy concentration of a few percents should be enough for the system to show a glassy behavior , with a tendency to the formation of stripe - ordered nanodomains around vacancy pairs .    finally , we comment on the possible relevance of this work to layered iridates , which motivated the derivation of the kitaev - heisenberg model .
the experimental data ( magnetic susceptibility , heat capacity@xcite ) on na@xmath3iro@xmath4 indicate antiferromagnetism at low temperature and a substantial magnetic anisotropy . yet
, the angular dependence of the measured susceptibility is qualitatively opposite to that found for a vacancy - free model .
we then showed that non - magnetic vacancies in the dilute limit , i.e. , where interactions between vacancies are irrelevant , can reverse locally this tendency by blocking the spin - canting for some field directions .
the recent x - ray magnetic scattering data @xcite ruled out the nel state , implying that the kitaev interaction in na@xmath3iro@xmath4 may dominate over a simple heisenberg coupling .
in addition to the stripe phase intrinsic to the kitaev - heisenberg model , the so - called  zig - zag  spin order is also consistent with the data @xcite .
an element in favor of a description of layered iridates a@xmath3iro@xmath4 ( a@xmath269li , na ) by the kitaev - heisenberg model is the comparison of the low temperature susceptibilities of these compounds@xcite . using our result @xmath138 for the stripe phase
, we found a relation between @xmath60 values in these two compounds , which indicates that the li compound , as the na - based one , is most probably in the stripe phase but much closer to the transition towards the spin liquid .
this conclusion is in agreement with previous work@xcite estimating the relative @xmath60 values from a comparison of the paramagnetic curie temperatures .
the application of uniaxial pressure on a single crystal of li@xmath3iro@xmath4 might even drive the system into the liquid phase .    to describe the experimental situation in layered iridates on a quantitative level
, other factors such as lattice distortions have to be considered as well .
in particular , the lattice distortions may bring about substantial anisotropy in the ground state @xmath270-factors and in the van vleck contribution to the magnetic susceptibility , originating from transitions to higher lying spin - orbital quartet split by non - cubic ( tetragonal and/or trigonal ) crystal fields .
furthermore , these distortions may lead to spin interactions not included in the kitaev - heisenberg model . for intermediate strength of spin - orbit coupling like in iridates
, these effects might be essential for a quantitative description of magnetic properties .
moreover , our study suggests that a glassy behavior could occur in the vicinity of the spin liquid phase ; in this context , considering the present model at larger scales , with a finite density of vacancies , could help to explore the effect of temperature .                                                                    a quantitative comparison of our calculated susceptibilities to those measured in ref .
would include also the effect of temperature , especially as the measured anisotropy was also large in the paramagnetic phase .