low - dimensional quantum heisenberg antiferromagnets are known to exhibit non - trivial properties due to quantum fluctuations .
a most striking example comes from spin ladders : interpolating between the chain geometry and the square lattice geometry , they have provided an interesting new playground to explore the physics of low - dimensional strongly correlated electron systems @xcite . specifically , it has been shown that af heisenberg spin-@xmath1 ladders with an even number of legs are quantum spin liquids with purely short - range spin correlations .
their spin correlations decay exponentially due to a finite singlet - triplet gap in the energy spectrum .
in contrast , heisenberg ladders with an odd number of legs display properties similar to those of single chains at low temperature , namely gapless spin excitations , and power - law spin correlations .
this difference between even - leg and odd - leg ladders was predicted by theory and confirmed by experiments in a variety of systems @xcite .    in this paper
we focus on disordered spin ladders . for the case of doping with static non - magnetic impurities ,
one distinguishes between site and bond disorder .
site dilution occurs upon doping the magnetic ions with non - magnetic ones ( as for instance substitution of cu@xmath2 with zn@xmath2 ) .
this leads to the formation of local moments ( lms ) close to the dopant site @xcite .
effective residual interactions between the lms lead to long - range antiferromagnetic ordering@xcite , experimentally detected in systems such as @xmath3 ( ref .  ) and bi(cu@xmath4zn@xmath5po@xmath6 ( ref .  ) .
bond disorder occurs instead when the dopant ions replace the ions which act as bridges between the magnetic ions .
for instance , bond disorder can be introduced in ipa-@xmath7 by a partial substitution of non - magnetic @xmath8 for the likewise non - magnetic @xmath9 , affecting the bond angles in the @xmath10-halogen - halogen-@xmath10 superexchange pathways @xcite .
the reduced strength of magnetic interactions on the affected bonds leads to bose - glass behavior when a strong magnetic field is applied .    from the theory side ,
a number of studies have addressed the effects of site and bond disorder on heisenberg ladders . ref .   examined the correlation length of randomly site - diluted spin-1/2 heisenberg 2-leg ladders at weak and intermediate interchain couplings , showing an apparent divergence of the spin correlations at low temperatures , due to the presence of the impurities .
a related divergence of the staggered susceptibility has been reported in ref .  .
this behavior has been related to that of the random - exchange heisenberg model ( rehm ) @xcite , describing the effective random interactions between lms @xcite . in a more recent study ,
three of us have investigated bond - diluted 2-leg ladders @xcite and observed enhancement of spin correlations due to bond dilution . however , a general understanding of the mechanisms for enhancement of correlations in site- and bond - diluted heisenberg ladders is still lacking .
it is also unclear whether such an enhancement of spin correlations generally exists for even - leg ladders with a number of legs greater than two .    in this paper
, we present a study of site- and bond - diluted spin-1/2 af heisenberg ladders with n=2 , 4 and 6 legs . making use of quantum monte carlo simulations
, we can address the correlation length of the system down to extremely small temperatures at which the asymptotic @xmath11 behavior sets in .
our main findings can be summarized as follows :    1 .   for site - diluted even - leg ladders ,
we observe a strong enhancement of correlations upon doping up to @xmath12 legs . in the specific case of 2-leg ladders ,
we find that the system realizes the physics of the rehm with a _ discrete _ distribution of the effective couplings , leading to a _ non - universal _ behavior at low temperatures ; hence the known predictions for the universal regime of the rehm with a continuous distribution of couplings are expected not to apply to realistic models of doped spin ladders .
2 .   for bond - diluted ladders
, we find that correlations are _ suppressed _ by a low level of dilution , due to dilution - induced dimerization .
hence the system is first driven to a gapless phase with short - range correlations .
the correlation length becomes logarithmically divergent for vanishing temperatures only beyond a critical dilution .
we study the heisenberg model on @xmath13-leg ladders @xmath14 where @xmath15 is the quantum spin operator at site @xmath16 along the @xmath17-th leg .
the first term in the hamiltonian describes interactions along the legs with coupling @xmath18 , whereas the second term represents the rung interaction with coupling @xmath19 .
the random variables @xmath20 , @xmath21 express the site or bond dilution : in the case of site dilution , @xmath22 and @xmath23 , where @xmath24 takes value 1 if the site @xmath25 is occupied ( with probability @xmath26 ) or 0 if it is empty ( with probability @xmath27 ) . in the case of bond dilution ,
@xmath28 take values 0 with probability @xmath29 if the bond is not occupied , and take values 1 otherwise .
the simulations are performed using the stochastic series expansion ( sse ) quantum monte carlo ( qmc ) method based on the directed loop algorithm @xcite .
periodic boundary conditions ( pbc ) along the leg direction are used . to access the regime of very low temperatures , at which the asymptotic low-@xmath30 behavior of the correlation length sets in , we have made use a @xmath31-doubling scheme @xcite , allowing us to efficiently access inverse temperatures up to @xmath32 .
the main objective of this paper is to study the correlation length along the leg direction , which is calculated via the disorder - averaged second - moment estimator : @xcite : @xmath33_{\rm av}}{[s(\pi+2\pi / l,\pi)]_{\rm av}}-1 } .
\label{e.2ndmoment}\ ] ] here @xmath34 is the time - averaged structure factor @xmath35 with @xmath36 corresponding to the total number of sites .
@xmath37_{\rm av}$ ] denotes the disorder average , which is performed over 300 - 600 disorder realizations .
antiferromagnetic even - leg ladders without doping display a rung - singlet ground state @xcite with a finite gap @xmath38 to triplet excitations . in this state ,
the @xmath13 spins on the same rung preferentially form a singlet state , and therefore effectively decouple from the rest of the ladder .
this leads to exponentially decaying correlations in the direction of the legs , characterized by a finite correlation length @xmath39 .    at low enough concentration of dopants , the main effect of site dilution on an even - leg ladder
is that of turning the number of spins on a rung from even to odd : in this situation , the state on the rung turns from a singlet into a _ doublet _ , which corresponds to an effective @xmath0 localized moment ( lm ) .
this @xmath0 moment remains exponentially localized ( over a characteristic length @xmath39 ) close to the impurity site @xcite , but its finite overlap with other lms leads to an effective interaction which generically decreases exponentially with the distance . in the case of dominant rung interactions , @xmath40 , the interaction between lms is appropriately described within second - order perturbation theory ( in @xmath41 ) as resulting from the exchange of virtual massive triplets between two lms @xcite . in the case of a 2-leg ladder , in which we can identify the location of a lm with a spin site next to a missing rung partner ( dangling spin ) , the effective coupling between lms has the form of a su(2)-invariant heisenberg interaction , and at large separation between lms , @xmath42 , the coupling strength takes the form @xcite @xmath43 where the staggering factor takes value @xmath44 if the two lms belong to the same sublattice and @xmath45 otherwise .
eq .   holds strictly speaking only in the case in which @xmath46 , but a similar exponential decay of the effective lm coupling ( _ without _ the square - root denominator ) has been observed numerically in the case of much stronger leg interaction @xcite , as _ e.g. _ for @xmath47 , for which @xmath48
. moreover the above expression is numerically found to account for the decay of the effective couplings already for moderate distances , @xmath49 .
hence , towards an effective model of the site - diluted ladders , we will assume in the following that the exponential decay of lm couplings with decay rate @xmath39 remains valid even when perturbation theory is no longer applicable  which is the case of interest . in this paper
, we deliberately choose @xmath50 to have a sizable correlation length @xmath39 , for reasons which will become clear in the following .
our assumption then implies that non - perturbative effects only affect the prefactor to the exponential in eq .  .
in addition , we will assume that the exponential decay present in eq .
sets in for distances between lms of the order of @xmath39 .    in the case of @xmath13-leg ladders with @xmath51
a detailed theory of effective lm couplings is not available to our knowledge , but for widely spaced vacancies , the extended structure of the localized doublet becomes irrelevant , and for @xmath52 a perturbation approach analogous to that of 2-leg ladders should be applicable , leading to effective staggering couplings between lms with exponential dependence on the distance .
( color online ) . probability distribution ( @xmath53 ) and cumulative distribution ( @xmath54 ) of the effective couplings @xmath55 between nearest - neighbor localized moments in site - diluted 2-leg ladders , with parameters @xmath56 , @xmath57 ( upper panels ) , and @xmath58 , @xmath59 ( lower panels ) .
the red dashed lines show the probability distributions in the continuum approximation according to eq .
.,title="fig:",width=359 ]   ( color online ) .
probability distribution ( @xmath53 ) and cumulative distribution ( @xmath54 ) of the effective couplings @xmath55 between nearest - neighbor localized moments in site - diluted 2-leg ladders , with parameters @xmath56 , @xmath57 ( upper panels ) , and @xmath58 , @xmath59 ( lower panels ) .
the red dashed lines show the probability distributions in the continuum approximation according to eq .
.,title="fig:",width=359 ]    in the case of low doping , @xmath60 , the probability distribution of having two nearest - neighboring lms at a distance @xmath61 along the leg direction in a 2-leg ladder is given by @xmath62 with a corresponding average @xmath63 .
as found numerically in ref .  , the effective coupling has the form @xmath64 for @xmath65 . as an approximation
, one may assume that the same exponential behavior survives for @xmath66 .
then the probability distribution for the strength of the couplings between nearest neighboring lms reads @xmath67.\end{aligned}\ ] ] this result suggests the possibility of extracting analytically the probability distribution for @xmath68 , as done in ref .  .
to this end , one may further take the continuum approximation , namely , approximate the summation in eq .   by an integral : @xmath69~. \label{e.pdpj}\end{aligned}\ ] ] this allows then to write @xmath70 where @xmath71 .
hence , within the continuum approximation the n.n .
couplings obey a simple power - law distribution.@xcite the average absolute value of the coupling between neighboring lms takes then the form @xcite : @xmath72_{\rm av } = \frac{1-\gamma}{2-\gamma } ~j_0 = \frac{2z\xi_0}{1 + 2z\xi_0 } ~j_0 .
\label{e.jav}\ ] ]    however , a critical analysis of the above derivation shows that the continuum approximation is problematic .
in fact , it requires that the distances @xmath61 giving a significant contribution to the sum in eq .
be @xmath73 , which is in contradiction with the fact that @xmath74 decreases exponentially with @xmath61 ; the characteristic decay length is @xmath75 , which for @xmath76 takes values @xmath77 .
[ f.pjeff-numerical ] shows the distribution @xmath78 determined numerically according to eq .  , _
i.e. _ , by sampling the _ discrete _ distribution lengths @xmath74 , for @xmath79 , @xmath57 ( corresponding to @xmath80 , and giving @xmath81 ) , and @xmath82 , @xmath59 ( corresponding to @xmath83 , and giving @xmath84 )  these parameter sets will be relevant for our study of correlations in the following .
we notice that the distribution shown in fig .
[ f.pjeff-numerical ] is only quantitatively correct when @xmath85 , namely for @xmath86 , while it is only an approximation otherwise .
it is clear that for @xmath85 , the distribution of n.n .
couplings obtained numerically deviates strongly from the prediction of eq .  .
over most of its support , the distribution of n.n .
couplings has a fundamentally _ discrete _ structure , due to the fact that the largest values of @xmath68 , which also take the largest probabilities , are associated with short separations @xmath61 among n.n .
these probabilities _ grow _ as a power law of the absolute value of the coupling , opposed to what predicted by eq .
for @xmath87 .
obviously a coarse graining on the @xmath88-peaks of the exact distribution would gradually reconstruct the shape of eq .  , but given the finite separation between the @xmath88-peaks even in the thermodynamic limit , this coarse graining is not justified over a large range of @xmath68 .
in fact the @xmath88-peaks become dense only in the limit @xmath89 , which corresponds to larger and larger separations among the lms , and hence to lesser and lesser probabilities . only in that limit
the continuum approximation appears legitimate ( in fact numerically we recover the analytical prediction for @xmath90 due to the finite width of our bins , @xmath91 ) .    in conclusion , for @xmath92 and @xmath93
, @xmath78 has a remarkable structure : it assigns the largest probability to the range in which @xmath68 is a _ discrete _ variable , preventing a straightforward approximation of the distribution with a continuous function .
this will have significant consequences on the analysis of doped ladder systems in relation with models of randomly coupled spins .
( color online ) .
correlation length @xmath94 of site - diluted ladders with @xmath95 and @xmath96 .
upper panel : comparison between 2-leg , 4-leg and 6-leg ladders at inverse temperature @xmath97 .
lower panel : correlation length of the 2-leg ladder at different temperatures.,width=283 ]      in both cases of @xmath98 and @xmath51 the unfrustrated effective couplings between lms can lead to the formation of a _
_ ground state with algebraically decaying correlations between the sites , and hence an infinite correlation length .
when looking at a fixed finite temperature @xmath30 the correlation length is necessarily finite for one - dimensional su(2 ) invariant systems with short - range interactions , but it is expected to increase because the average coupling between lms increases with @xmath27 , as in eq .  ,
so that the effective temperature @xmath99_{\rm av}$ ] decreases .    despite the increase of correlations due to disorder , the range of correlations is fatally bounded by the fact that any finite amount of site dilution will break an infinite @xmath13-leg ladders into finite segments , of characteristic length @xmath100 ( corresponding to the inverse linear density of rungs which are fully removed by doping )
when considering doped ladders at finite temperatures no geometric bound on the correlation length ( neither coming from finite - size effects nor from ladder fragmentation ) is present as long as the correlation length satisfies the condition @xmath101 it is clear that such a condition can only be satisfied when the doping is sufficiently weak .
hence , for any fixed temperature , upon increasing the doping concentration the correlation length @xmath102 will cross over from a growing behavior at low doping @xmath103 to a decreasing behavior for @xmath104 , hence going through a maximum for a given optimal doping value @xmath105 , which is quite stable at sufficiently low temperatures .
this is indeed seen in fig .
[ f.sitedil](a ) , in which we observe @xcite that , for @xmath106 , @xmath107 for @xmath98 , 4 , and 6 .    for any finite value of doping and at arbitrarily low temperatures , the correlation length is upper bounded by the average length of the ladder segments , @xmath108 . yet , as shown in fig .
[ f.sitedil](b ) , at low enough doping the upper bound is reached extremely slowly in temperature - in fact , for the 2-leg ladder at @xmath109 and @xmath47 , the correlation length is still not saturated to its finite @xmath110 value for @xmath111 .
it is easy to understand that this occurs because the average effective coupling @xmath112_{\rm av}$ ] decreases as @xmath27 decreases , leading to exceedingly weak correlations even at extremely low temperatures .
this means that for sufficiently low doping , the condition @xmath113 is actually satisfied over all numerically accessible low temperatures ( and quite possibly also over all experimentally accessible temperatures ) , so that the segmented nature of the doped ladders become in fact _
irrelevant_. in the following we will focus our discussion on this regime .      as discussed in section [ s.pjeff ] , from the statistical point of view the system of interacting lms randomly distributed on the bipartite 2-leg ladder can be approximated by a system of @xmath0 spins interacting with randomly distributed couplings , following an intrinsically _ discrete _ distribution . in general terms , a doped 2-leg ladder represents a physical realization of the so - called random - exchange heisenberg model ( rehm ) @xcite , but with a special structure of the coupling distribution .
( color online ) .
temperature dependence of the correlation length of site - diluted 2-leg ladders with @xmath114 , and power - law fits at low temperature.,width=283 ]    numerical real - space renormalization group studies @xcite show that the rehm colorred for a _ continuous _ distribution of couplings @xmath115 has two regimes :    * _ universal regime _ :
if the distribution is not singular , or has a power - law singularity @xmath116 with @xmath117 ( where @xmath118 ) , the low - temperature thermodynamics of the rehm is dominated by a unique fixed point ; the correlation length diverges as @xmath119 with a universal scaling exponent , @xmath120 with @xmath121 .
this result has been confirmed by quantum monte carlo @xcite for a non - singular distribution ( square box , @xmath122 ) .
* _ non - universal regime _ : for a more strongly singular distribution @xmath123 the renormalization group flows to a non - universal fixed point , which depends on the initial coupling distribution .
studies on doped coupled spin - peierls chains using real space renormalization group and sse qmc in this regime also confirm that the exponent @xmath124 depends on the doping concentration@xcite .    as discussed in sec .
[ s.pjeff ] , the distribution of couplings between lms in a doped 2-leg ladder is intrinsically discrete , approximating a continuous distribution only for vanishing couplings .
it is therefore interesting to ask whether the doped 2-leg ladder reproduces the known physics of the rehm model with a continuous distribution of couplings .
if a continuous approximation to the coupling distribution is in order , as in eq .  , the system should exhibit a _
universal rehm behavior _ for @xmath125 , and a non - universal behavior for @xmath126 .
in other words , ladders with sufficiently strong doping @xmath27 or sufficiently large correlation length @xmath39 in the undoped limit should exhibit universal scaling properties at low temperature . to our knowledge
the only numerical study of correlations in site - diluted ladders is ref .  , which , making use of quantum monte carlo , investigates ladders with @xmath127 and dopings @xmath128 , 0.04 and 0.1 .
given that @xmath129 in this case @xcite , for @xmath130 and 0.1 one obtains @xmath131 and @xmath132 respectively , which appear to be both safely on the universal side of the rehm with continuous couplings .
indeed a fit to the low - temperature behavior of the correlation length to the scaling law @xmath133 gives an estimate of @xmath134 ( independent of @xmath27 ) which appears consistent with the predictions of refs .  .
we have performed qmc simulations to extensively study the low - temperature behavior of the correlations in site - diluted 2-leg ladders .
our qmc results , however , do not support the conclusion in ref .   about a universal rehm behavior of diluted 2-leg ladders .
we investigate doped ladders at similar concentrations with respect to those of ref .  , down to significantly lower temperatures - ref .
stops at @xmath135 , while we descend down to @xmath136 . a collection of our results is shown in fig .  [ f.xit ] .
in qualitative agreement with what expected from the theory of the rehm model , we also observe the clear onset of a low - temperature power law scaling of the type @xmath137 .
yet our central observation is that this power - law scaling exhibits a _
non - universal , doping dependent _ exponent @xmath138 .
for all the parameter sets we considered , a fit of the low - temperature correlation length provides values of @xmath124 which are systematically below the one exhibited by the universal regime of the rehm with continuous couplings .    the data shown in fig .
[ f.xit ] refer to 2-leg ladders with @xmath139 and @xmath140 and @xmath141 , with corresponding @xmath142 and @xmath143 respectively ( here we use @xmath59 , which we estimate independently with qmc ) ; and 2-leg ladders with @xmath127 ( as in ref .  ) and @xmath144 , @xmath145 , and @xmath146 , corresponding to @xmath147 , @xmath148 , and @xmath149 respectively .
although all the values of @xmath150 fall below the critical value @xmath151 , as shown in fig .
[ f.xit ] , in fact , all the corresponding data sets for the correlation length display a _ non - universal _ power - law scaling .
in particular , we observe that data sets with close values of the parameter @xmath150 ( _ e.g. _ @xmath139 , @xmath140 with @xmath152 , and @xmath127 , @xmath153 , with @xmath154 ) show a clearly different scaling .
we also critically revisit the parameter sets explored in ref .  .
a comparison to our data , as done in fig .
[ f.greven ] , shows that in the data of ref .
the significant scattering due to the statistical uncertainty is masking the correct asymptotic scaling regime at low temperature .    a comprehensive summary of our results on the low - temperature scaling of the correlation length is provided in fig .
[ ph.diagram ] , where we report the values of the @xmath124 exponent obtained for different doping values and @xmath155 ratios .
the values of doping concentrations and ratios were chosen such that the average distance between lms @xmath156 is at least larger than @xmath39 .
the lms , therefore , are not overlapping on average and can be described by an effective rehm with couplings obeying eq .  .
furthermore , for all the points shown in fig .
[ ph.diagram ] , @xmath117 , we do not observe a universal @xmath124 value .
indeed , it is found that @xmath124 depends not only on the doping concentration but also on the ratio of couplings @xmath155 . and
an @xmath124 exponent close to that of the universal rehm regime is only observed in a narrow region around @xmath157 and @xmath76 .
( color online ) .
values of the @xmath124 exponent extracted from the low - temperature behavior of site - diluted 2-leg ladders.,width=283 ]      the above results lead us to conclude that doped 2-leg ladders do not realize in general the universal physics of the rehm with a continuous coupling distribution .
in fact , the coupling distribution for the lms in doped 2-leg ladders has a fundamentally discrete structure , approximating a continuous one only for vanishing couplings .
the weight assigned by the distribution to strong discrete couplings is dominant over that of weak , quasi - continuous couplings in the case of sizable doping @xmath27 ( suppressing the probability of large separation between the lms ) and of large correlation length @xmath158 ( increasing the strength of the couplings between lms , the distance being fixed ) . nonetheless , this is the regime in which the parameter @xmath159 is small , and in which one would expect the universal rehm physics to be manifested in the continuum approximation . on the other hand ,
the quasi - continuous part of the distribution ( corresponding to very small couplings ) acquires a bigger weight for small @xmath27 and @xmath39 , which determine a small average coupling according to eq .   and hence increase the probability of the small values of the coupling . in this regime @xmath160 is closer to 1 , which leads to too strong a singularity in the probability distribution for the universal regime of the rehm with continuous couplings to be manifested .    as a consequence ,
we generally observe a low - temperature power - law scaling of the correlation length with doping - dependent exponents .
as far as the numerical results are concerned , this scaling appears to be the asymptotic one for @xmath11 , given that it sets in for temperatures well below the average lm coupling .
on the other hand , according to a real - space renormalization group ( rsrg ) approach @xcite , the system at a lower temperature is increasingly sensitive to the part of the distribution related to the weaker couplings ; given that this part is the one which best approximates a continuous distribution , one could naturally expect that at very low temperatures the physics of the rehm with continuous couplings be reproduced by the doped 2-leg ladders . this argument would then suggest that , when lowering the temperature , the system will sooner or later attain the universal regime of the rehm with continuous couplings if @xmath117 .
nonetheless one can argue that in a doped 2-leg ladder the rsrg flow has to be necessarily stopped at a finite length , corresponding roughly to the average segment length @xmath108 . in practice
this imposes an upper bound to the correlation length and a lower bound to the temperature through the condition eq .  ; these bounds might prevent one from attaining the low - temperature regime at which universal rehm physics would manifest itself .
finally , our study leaves one open question concerning the role of discreteness of the initial coupling distribution for the physics of the rehm . the numerical rg study of ref .   only addresses initial continuous distributions ; a systematic rg study of the flow starting from realistic distributions stemming from the doped - ladder physics would be highly desirable .
in the following , we consider the doping- and temperature dependence of correlations in bond - diluted ladders , which show a marked difference with respect to site dilution . this is due to a fundamental geometric aspect which distinguishes site and bond dilution : diluting a site leaves a single unpaired spin , giving rise to a lm , while eliminating a bond leaves always _ two _ unpaired spins , located on different sublattices .
the corresponding lms are therefore interacting with an effective antiferromagnetic coupling , mediated by the the shortest path of bonds connecting them . in the following
, we will focus for simplicity on the case @xmath96 .
[ f.bonddil-fixt ] shows the evolution of the correlation length at a fixed , low temperature @xmath161 . a striking difference with respect to the case of site dilution , shown in fig .
[ f.bonddil-fixt ] , is that the correlation length does not increase immediately with doping : it remains essentially constant for @xmath98 , while it even _ decreases _ for @xmath51 .    in the case of a 2-leg ladder
, we can understand qualitatively the behavior of correlations by considering that bond dilution has two fundamentally different effects ( sketched in fig .
[ f.bonddil-cartoon](b ) ) :    * if a rung bond is diluted , two lms are liberated but still interact antiferromagnetically , so that they can screen each other to form a rung singlet at a lower energy @xmath162 ( corresponding to the effective coupling between lms which are rung neighbors ) . this screening is very effective at low dilution because in that case the @xmath162 interaction is largely dominant over the interaction between two lms belonging to different rungs , due to the large spacing of two diluted rung bonds  this spacing is @xmath163 ( given that only one bond in three is a rung bond ) ; * if a leg bond is diluted , the rung - singlet state on the two adjacent rungs is further reinforced because its connectivity to the other rungs is lowered .
( color online ) .
correlation length of bond - diluted 2- , 4- , and 6-leg ladders with @xmath164 , @xmath96 , and temperature @xmath165,width=283 ]    the dilution of rung bonds has therefore the effect of lowering locally the gap over the ground state , whereas the dilution of leg bonds has the effect of increasing it @xcite .
the fact that there are twice as many leg bonds as rung bonds suggests that the second effect dominates over the first .
therefore the spin gap of the undoped ladder is locally preserved or even enhanced , leading to short - range correlations even in presence of bond dilution .
nonetheless , standard consideration on rare - event physics lead to conclude that even a very weak bond dilution leads the system to a _
gapless , griffiths - like phase_. in fact , the gap can close locally if rare regions appear in which only rung bonds are diluted  the limiting case being that of the formation of local strands made of two uncoupled @xmath0 chains .
these exponentially rare regions lead to the closing of the spin gap in the thermodynamic limit ; yet global correlations remain short ranged , due to the localized nature of the rare , locally gapless regions .
this griffiths phase is reminiscent of what has been observed by some of us in the heisenberg model on the inhomogeneously bond - diluted square lattice @xcite . in that case
, the predominance of correlation - suppressing dilution was guaranteed by the inhomogeneity of doping probabilities ( favoring the appearance of ladder - like or dimer - like structures ) .
( color online ) .
sketch of the bond dilution effects on 2-leg ladders . a ) pure ladder ; b ) ladder at low dilution : enhancement of rung singlets by leg - bond dilution , and low - energy singlets formed between lms after rung - bond dilution ; c ) ladder at stronger dilution : couplings between lms belonging to different diluted rungs.,width=302 ]     ( color online ) .
correlation length for 2-leg and 4-leg ladders with dilution of ( center-)rung bonds .
the average length of ladders segments @xmath108 is obtained upon homogeneous dilution .
, width=283 ]    for a sufficiently large bond dilution @xmath166 ( @xmath167 ) , the correlation length starts increasing with increasing @xmath27 .
this implies necessarily that the liberated lms belonging to different rungs begin to correlate with each other ( fig .
[ f.bonddil-cartoon](c ) , see also refs . ) .
such correlations appear because the spacing between diluted rung bonds decreases , and hence the couplings between lms belonging to different rungs become of the same order of @xmath162 .
we have isolated the effect of enhancement of correlations through rung - bond dilution by considering uniquely this form of dilution , namely by taking @xmath168 with probability @xmath27 , while @xmath169 with probability 1 .
the resulting correlations @xmath170 ( r = rung ) as a function of @xmath27 for 2- and 4-leg ladders are shown in fig .
[ f.rungdilution ] , and they are seen to increase quite fast with dilution ( following the approximate form @xmath171 for sufficiently low dilution ) . yet in the original system the correlations induced by rung - bond dilution are upper - bounded by the characteristic length of the ladder segments @xmath108 , which we have estimated by generalizing an efficient algorithm recently developed for homogeneous percolation @xcite .
it is found that at finite @xmath27 , @xmath172 .
it follows the same asymptotic behavior as from the naive estimate @xmath173 , but with a renormalized factor @xmath174 .
hence the enhancement of correlations due to rung - bond dilution competes with this geometrical restriction on correlations for sufficiently strong doping .
quite remarkably , we observe that the optimal doping for the enhancement of correlations in the homogeneously doped 2-leg ladder coincides with the value of dilution ( @xmath175 ) at which the correlation length of the rung - only diluted ladder equals the average length of ladder segments @xmath108
. therefore we can conclude that the non - monotonic behavior of the correlation length at intermediate dilution values is completely determined by the competition between the correlation enhancement through rung - bond dilution and the correlation suppression due to ladder fragmentation .
( color online ) .
( a ) correlation length of bond - diluted 2-leg ladders with @xmath95 and @xmath176 at various temperatures .
( b ) temperature dependence of the correlation length of above system .
the curves with @xmath177 are omitted due to @xmath178 is limited by the geometric length .
( c ) and ( d ) fitting parameters of correlation length @xmath179 , as a function of bond dilution.,width=321 ]      as discussed in the previous subsection , for @xmath180 correlations start to increase as a function of bond dilution .
[ f.bonddilution.t](a ) shows the doping dependence of the correlation length of 2-leg ladders for higher temperatures than those considered in fig .
[ f.bonddil-fixt ] : here we observe that the correlation length does not appear to evolve significantly with temperature for doping values @xmath181 , whereas it becomes significantly dependent on temperature for @xmath166 .
moreover , the low - temperature dependence of the correlation length , see fig .
[ f.bonddilution.t](b ) , indicates that in this regime the correlation length grows logarithmically with decreasing temperature as @xmath182     ( color online ) .
( a ) temperature dependence of correlation lengths for four cases of rung - bond and random dilution .
, width=283 ]    the dependence of the coefficients @xmath183 and @xmath184 on bond dilution is shown in fig .
[ f.bonddilution.t](c ) and ( d ) .
while @xmath183 decreases with increasing bond dilution , @xmath184 increases .
in particular @xmath184 appears to vanish for @xmath181 , and it seems to suggest that a phase transition occurs at finite @xmath180 from a phase with correlations converging to a finite value at low @xmath30 to a phase with correlations diverging logarithmically with decreasing temperature .
obviously the divergent behavior of correlation length only persists up to lengths of the order @xmath108 , corresponding to the characteristic length of ladder segments , but the exceedingly slow growth of @xmath178 with decreasing @xmath30 keeps it safely below @xmath108 for all the temperatures and doping values we explored .    to gain a deeper understanding on the above results , we study the temperature dependence of rung - bond dilution and random dilution separately .
we consider rung - bond dilution with concentrations @xmath185 such that the fraction of diluted rung bonds is the same as for the cases of random dilution that we have investigated .
this allows us to quantitatively ascertain the separate effect of leg - bond and rung - bond dilution . in presence of rung - bond dilution only , the correlation length is seen to converge to a finite value as @xmath186 , even at relatively high concentration .
two examples with concentration @xmath187 and @xmath188 are shown in fig .
[ para.z ] ( data for higher concentration , not shown here , exhibit a similar behavior ) .
this implies that the system with rung - bond dilution only is gapless ( because of the unbounded size of rare regions with _
e.g. _ diluted adjacent rung bonds ) and it has a finite correlation length , namely it is in a griffiths - like phase . adding leg - bond dilution ( with concentrations @xmath189 and @xmath190 respectively ) , we recover a randomly bond - diluted ladder ( with concentrations @xmath191 and @xmath192 respectively ) . for large enough dilution
the correlation is found to grow logarithmically with decreasing temperature , as in eq .  , which seems atypical for a griffiths - like phase .
on the other hand , we observe that , for dilution @xmath185 , the correlation length for the rung - bond diluted ladder with concentration @xmath193 is an upper bound to that of the randomly diluted ladder with concentration @xmath27 , a fact which is somewhat intuitive , given that rung - bond dilution is the main mechanism leading to the enhancement of correlations .
therefore it is a priori unclear whether the logarithmic growth of the randomly diluted ladder persists to even lower temperatures than the ones explored here , because this would eventually lead the correlation length of the randomly diluted ladder to exceed that of the rung - bond diluted one .
we can then conclude that the correlation length results are _ a priori _ consistent with two different scenarios for the randomly diluted 2-leg ladders . in a first scenario
the system transitions from a griffiths - like phase with a finite correlations for @xmath194 to a new phase with logarithmically diverging correlations , which is reminiscent of the behavior of a system controlled by an infinite - randomness fixed point ( irfp )  see the discussion below . in a second scenario
the correlation length exhibits different temperature dependences at intermediate temperatures for increasing bond dilution . yet
the correlation length of the randomly diluted ladders converges always to a finite value , because it is expected to be upper - bounded by that of rung - bond diluted ladders , which is seen to converge to a finite value . in this second scenario
the system is therefore in a griffiths - like phase for all the values of bond dilution we explored . in the following
we find that the analysis of the temperature dependence of the uniform susceptibility helps clarifying which scenario is the most appropriate .
[ xu.z](a ) shows the uniform susceptibility of a 2-leg ladder for various values of bond dilution . in the undoped system
, the susceptibility of 2-leg ladders vanishes exponentially at low temperatures , in agreement with well - known previous results , and showing the gapped nature of the spectrum .
in contrast , the susceptibility of the doped system diverges following a curie law for low temperatures .
furthermore , fig .
[ xu.z](b ) shows that the curie coefficient increases as a function of bond dilution concentration following a power law , @xmath195 .
one might expect that curie paramagnetism comes trivially from spins which have remained isolated after bond dilution . yet
the concentration of such spins would scale as @xmath196 with dilution , given that one needs to dilute at least 3 bonds to decouple one spin from the rest of the ladder .
therefore the paramagnetism observed in the system has a collective nature .
in fact , two contributions add to the one of free spins :     ( color online ) .
( a ) low - temperature uniform susceptibility of bond - diluted ladders with @xmath197 and @xmath198 .
( b ) the curie coefficient increases as a power - law function of the bond dilution concentration @xmath195.,width=302 ]    * odd - numbered clusters , obtained after fragmentation of the ladder into two ( or more ) pieces , have a doublet ground state , behaving as a collective spin 1/2 . in our simulations on systems with an even total number of sites , three contiguous diluted bonds lead to an isolated site plus an odd - numbered cluster , so that the latter clusters have a probability @xmath199 .
alternatively odd - numbered clusters can be obtained by cutting the ladder at two distinct locations into two odd - numbered extended clusters ( bigger than one single site ) .
yet one can easily see that this would require at least 5 diluted bonds at specific locations , giving rise to a probability @xmath200 , and consequently to a very small contribution to the uniform susceptibility .
therefore the contribution to collective paramagnetism can not come from such clusters .
* clusters with an even number of spins , but a non - equal number of spins on a and b sublattices , can also contribute to a diverging susceptibility .
they represent therefore a further candidate for collective paramagnetism .
( color online ) .
( a ) uniform susceptibility of even - numbered spin clusters of randomly bond - diluted 2-leg ladders
. the solid lines are ( a ) logarithmic fits of the form @xmath201 and ( b ) power - law fits @xmath202.,width=302 ]    to isolate the second contribution , we calculate the uniform susceptibility coming from _ even - numbered _ spin clusters only , @xmath203 . this quantity is shown in fig .
[ xuc.t ] .
we find that even - numbered uniform susceptibility at low temperatures can be equally well described by two fitting laws : a logarithmically - corrected curie behavior : @xmath204 and a power - law corrected curie behavior : @xmath205 in both cases we can conclude that the contribution to the uniform susceptibility coming from even - numbered clusters diverges more weakly than for a curie law , and therefore that our system exhibits a highly non - trivial collective paramagnetism .
the fitting coefficients @xmath206 , and @xmath207 for both cases are shown in fig .
[ fit.para1 ] .
we observe that all these fitting coefficients evolve smoothly with @xmath27 , apparently contradicting the picture of a possible phase transition suggested by the behavior of the correlation length .
in particular the non - universal power - law dependence of eq .   is typical of griffiths behavior , as observed _
e.g. _ in anisotropically bond - diluted square lattices @xcite .
therefore the uniform susceptibility seems to favor the scenario for which the system is in a griffiths phase for all values of the bond dilution we explored . for completeness we should also mention a third , less likely scenario ,
for which the correlation length is diverging logarithmically as @xmath11 even for infinitesimal doping ( but the prefactor @xmath184 of the logarithmic divergence becomes lower than our resolution ) .
therefore the system would be controlled by an irfp for all doping values explored here ; this scenario is somewhat consistent with the fact that a logarithmically corrected curie law is consistent with the susceptibility at all doping values considered in this section , and that it could be another sign of irfp behavior .
( color online ) .
fitting coefficients as a function of bond dilution .
left panels are fitting coefficients of the logarithmic function @xmath201 .
right panels are fitting coefficients of the power - law function @xmath202.,width=302 ]      as sketched in fig .
[ f.bonddil-cartoon ] , the rung lms in a bond - diluted ladder realize a ladder with weaker random couplings which are at once rung , leg and diagonal ones but without frustration .
the rung couplings are randomized by the variety of different local environments which mediate the effective interaction between the rung lms ; while the leg couplings and the diagonal couplings are obviously randomized ( in a correlated way ) by the positional disorder of the diluted rung bonds .    for ladders with @xmath47 we can not easily rely on perturbation theory to extract accurate expressions for the effective couplings , and
an explicit calculation of the disorder statistics goes beyond the purpose of this paper . in any instance
it is clear that the rung couplings obey a fundamentally different probability distribution with respect to the leg and the diagonal couplings : the rung couplings are always antiferromagnetic and weakly randomized , while the leg and diagonal couplings always take opposite signs and can be either ferromagnetic or antiferromagnetic - they are hence strongly disordered . to the best of our knowledge , such a peculiar model of a random ladder has not been investigated before .    in a critical quantum spin system
the relationship between energy scales ( @xmath208 ) and length scales ( @xmath209 at finite temperature ) is governed by the so - called dynamical critical exponent @xmath210 , @xmath211 .
the logarithmic dependence of the correlation length on temperature , eq .
, implies that @xmath212 : as mentioned before , this behavior , together with the logarithmically corrected curie law in the susceptibility , could be suggestive of the fact that the bond - diluted ladder at sufficiently strong disorder is governed by an irfp , similarly to what happens to _ e.g. _ random antiferromagnetic @xmath0 chains @xcite .
nonetheless , it is not the same irfp as in random antiferromagnetic chains , given that the correlation length and the susceptibility diverge differently ( as @xmath213 and as @xmath214 respectively ) in the latter .
existing rg studies of disordered ladder systems focus primarily on the case of antiferromagnetic ladders , either with nearest neighbor couplings only or with frustrated ones @xcite .
the occurrence of extended parameter regions with irfp is only observed in the case of frustrated ladders .
therefore it is quite remarkable to observe an irfp - like behavior in our unfrustrated system for a large interval of doping values .
on the other hand , the absence of a documented irfp for unfrustrated random ladders suggest that the logarithmic divergence observed in our system might not be the true asymptotic behavior for @xmath11 . nonetheless , a definite conclusion on such an asymptotic behavior would require a precise determination of the distribution of the effective couplings between lms , which goes beyond the scopes of this paper .
[ f.bonddil-fixt ] shows that the dependence of the correlation length on bond dilution in 4-leg and 6-leg ladders exhibits a qualitatively different behavior compared to 2-leg ladders . indeed for 4-leg and 6-leg ladders
@xmath178 is initially _ decreasing _ with increasing bond dilution . in the case
@xmath176 the correlation length reaches a minimum for @xmath215 , and then starts growing with doping , up to a maximum attained for @xmath175 . in the following
we will analyze this phenomenon in details .
( color online ) .
bond - dilution effects on 4-leg ladders : ( 1 ) rung - singlet enhancement upon leg - bond dilution ; ( 2 ) center - rung dilution leading to the formation of two 2-spin rung singlets ; ( 3 ) outer - rung dilution with formation of a low - energy singlets between a lms and a 3-spin doublet.,width=302 ]    similarly to what done for 2-leg ladders , we can distinguish among different effects of doping depending on which bond is diluted ; we will focus here on the case of a 4-leg ladder ( which is sketched in fig .
[ f.bonddil-cartoon-4leg ] ) :    1 .
the dilution of a _ leg bond _ leads to the the enhancement of the 4-spin rung singlet on the two adjacent rungs ; 2 .   the dilution of a _
center rung bond _ leads to the formation of two 2-spin rung singlets on the intact outer rungs ; 3 .
the dilution of an _ outer rung bond _ separates an outer lm from three rung spins , which tend to form locally a spin - doublet .
this doublet should have then tendency to form again a low - energy singlet with the lm through the effective antiferromagnetic coupling which binds them .
hence we observe that bond dilution of type 1 ) and 2 ) have the tendency to enhance locally the spin gap of the ladder , while dilution of type 3 ) has the tendency to suppress it .
the case of 6-leg ladders is analogous , with the only difference that dilution of the central rung bond leads to the formation of two rung doublets .
therefore we will hereafter focus on the 4-leg case for simplicity .
we have disentangled the competing effects of dilution on 4-leg ladders by considering only one at a time .
first of all we consider the effect of enhancement of correlations due to outer rung - bond dilution only , whose results are shown in fig .
[ f.rungdilution ] .
analogously to what has been seen for rung - bond dilution in the 2-leg ladders , outer - rung - bond dilution in 4-leg ladders leads to a strong enhancement of correlations ( with a similar exponential scaling @xmath216 , or = outer rung ) . in homogeneously diluted ladders
this effect of enhancement is fundamentally limited by the finiteness of ladder segments ( due to leg - bond dilution ) : indeed a numerical estimate of the average length of ladder segments , @xmath108 , shows that @xmath108 crosses the correlation length @xmath217 for a value of @xmath27 ( @xmath218 ) which nicely corresponds to the optimal doping in the homogeneously bond - diluted ladder .
hence , similarly to 2-leg ladders , we can conclude that the optimal doping is completely dictated by the competing effect of outer - rung - bond dilution and ladder fragmentation .
( color online ) . correlation length with selective bond dilution .
bond dilution is distinguished as : ( a ) rung - bond and leg - bond dilution in 2-leg ladders ; ( b ) center - rung bond , outer - rung bond and leg bond dilution in 4-leg ladders .
the dashed line is the average correlation length from above selectively bond - diluted mechanisms , width=302 ]    finally , we study the effect of selective dilution on the correlation length . in fig .
[ f.selective-dilution](a ) , we investigate the effect of rung - bond and leg - bond dilution in 2-leg ladders . in fig .
[ f.selective-dilution](b ) , we distinguish dilution of outer - rung bonds , center - rung bonds , and leg bonds in 4-leg ladders . as anticipated , latter forms of dilution in above cases lead to a significant suppression of correlations , again with approximate exponential scaling , @xmath219 ( l = leg , cr = center rung ) .
we observe that @xmath220 ( @xmath221 for 2-led ladder ) , so that the simultaneous suppression and enhancement of correlations due to dilution of different types of bonds are in strong competition .
the fact that at low doping the suppression of correlation wins is simply due to combinatorics : in 2-leg ladders ( a ) , rung bonds only hold a fraction @xmath222 while leg bonds contain @xmath223 of all bonds ; in 4-leg ladders ( b ) , outer - rung bonds only represent a fraction @xmath224 of all the bonds , whereas leg bonds and center - rung bonds are respectively @xmath225 and @xmath226 .
hence correlation - enhancing bond dilution has a global lower probability than correlation - suppressing one . in the limit of very low dilution ( @xmath227 ) we can easily imagine to separate a ladder into correlated regions ( of characteristic length @xmath39 ) which are statistically independent from each other . upon dilution , each correlated region
will only be affected by a single missing bond on average , with a probability @xmath228 for the @xmath229-th bond type ( @xmath229 = r , l in 2-leg ladders ; and @xmath229 = or , l , cr in 4-leg ladders ) . therefore the local correlations will be enhanced or suppressed in the same way as for the case of selective bond doping of type @xmath229 .
for @xmath227 it is then tempting to write the correlation length of the homogeneously diluted ladder as a spatial average of different local correlation lengths . assuming that the ratio of structure factors entering in eq .
can be written as @xmath230_{\rm av}}{[s(\pi+2\pi / l,\pi)]_{\rm av } } \approx \sum_i f_i \frac{[s(\pi,\pi)]_{{\rm av},i}}{[s(\pi+2\pi / l,\pi)]_{{\rm av},i}}\ ] ] where @xmath37_{{\rm av},i}$ ] denotes averaging over disorder realizations with dilution of bond of @xmath229-th type only ( and with dilution concentration @xmath231 ) , we obtain that @xmath232^{1/2}~.\ ] ] this average correlation length is presented by the dashed lines in fig .
[ f.selective-dilution ] .
it clearly shows a quantitative agreement with random bond dilution results at small dilution .
in this paper we have discussed the effects of site and bond dilution on the low - temperature correlations of even - leg @xmath0 ladders with antiferromagnetic heisenberg interactions .
site dilution is found to prune rung singlets and thus create localized moments which interact through unfrustrated , distance dependent couplings .
the hamiltonian describing the effective interaction between these moments is a random exchange heisenberg model ( rehm ) , with a gapless spectrum and power - law diverging correlations as temperature decreases to zero .
we find that the distribution of the effective couplings has an intrinsically discrete structure , which prevents the system from realizing the universal regime of the rehm with continuous couplings .
further studies would hence be desired to clarify the role of a discrete distribution of effective couplings in determining the fixed point which governs the low temperature physics of this system .
bond dilution , on the other hand , can either enhance or suppress locally the spin gap in even - leg ladders , depending on the location of the bond which is diluted . as a result of these competing effects , weak bond dilution is not enhancing the correlations of pure ladders .
in fact , correlations can even be suppressed by disorder , especially in 4-leg and 6-leg ladders .
the resulting short - range correlated phase has a gapless , griffiths nature , due to the appearance of exponentially rare , but locally gapless regions . beyond a critical concentration , the correlation - enhancing effect of bond dilution leads to a phase in 2-leg ladders , with a logarithmic scaling of the correlation length with decreasing temperature .
this behavior may suggests that the low - temperature behavior of the system is governed by an infinite randomness fixed point over an extended range of parameters .
nonetheless , based on the investigation of the uniform susceptibility of the system , and by comparison with existing results on disordered ladders , we argue that this might not be the true asymptotic behavior for @xmath11 , and that the system is more likely to remain in a griffiths phase for all values of bond dilution .    in order to connect our results to experiments on site - diluted and bond - disordered ladders ,
a fundamental question to address is the role of finite inter - ladder couplings , which are unavoidable in real materials .
if @xmath233 is their energy scale , such couplings drive the system towards a three - dimensional magnetically ordered phase below a temperature @xmath234 , which can be estimated via a mean - field approach as @xmath235 , where @xmath178 is the correlation length of the uncoupled ladders . for a gapless critical phase , as that induced by site dilution ( at any concentration ) and by bond dilution ( beyond a critical concentration ) , @xmath102 diverges for @xmath11 , so that the system orders at finite temperature . on the other hand
, the divergence of @xmath178 with decreasing temperature appears to be relatively weak ( power law with a small power @xmath236 for site dilution , or logarithmic for bond dilution ) .
hence a difference of several orders of magnitude between the intra - ladder and the interladder couplings , which is common to various ladder materials , allows the asymptotic low - temperature behavior of @xmath178 for the uncoupled doped ladders to be manifested at temperatures which lie well above the critical temperature @xmath234 for magnetic ordering . for bond - diluted ladders , moreover
, short - range correlations proper of the pure ladder persist up to a critical doping , so that three - dimensional ordering in coupled ladders should be absent at weak dilution .
increasing the bond dilution leads to an increase of correlations at low temperature , so that the system is expected to undergo a disorder - induced quantum phase transition from a griffiths phase to a magnetically ordered phase ( at @xmath110 ) , but with exponentially small transition temperatures , @xmath237 , which might be hardly detectable in experiments . as a result the temperature scaling of correlations investigated here for uncoupled ladders
is expected to be relevant to the behavior of real ladder materials with non - magnetic doping .
we would like to thank n. bray - ali , s. garnerone , and b. normand for valuable discussions .
we also acknowledge financial support by the department of energy , grant number de - fg03 - 01er45908 .
t. r. acknowledges the hospitality of the mpipks - dresden , where part of this work was completed .
r. y. was supported by nsf grant no .
dmr-1006985 and robert a. welch foundation grant no .
the numerical computations were carried out on the university of southern california high - performance supercomputer cluster .
200 e. dagotto , t. m. rice , science * 271 * , 618 ( 1996 ) . m. sigrist and a. furusaki , j. phys .
. jpn . * 65 * , 2385 ( 1996 ) .
sandvik , e. dagotto , and d. j. scalapino , phys . rev .
b * 56 * , 11701 ( 1997 ) . g.b .
martins , m. laukamp , j. riera , and e. dagotto , phys .
* 78 * , 3563 ( 1997 ) .
s. wessel , b. normand , m. sigrist , and s. haas , phys .
* 86 * , 1086 ( 2001 ) . n. laflorencie and d. poilblanc , phys .
lett . * 90 * , 157202 ( 2003 ) . c. yasuda , s. todo , m. matsumoto , and h. takayama , phys
b * 64 * , 092405 ( 2001 ) .
m. imada and y. iino , j. phys .
66 * , 568 ( 1997 ) .
h. weber and m. vojta , eur .
j. b * 53 * , 185 ( 2006 ) .
m. azuma , y. fujishiro , m. takano , m. nohara and h. takagi , phys . rev .
b. * 55 * , r8658 ( 1997 ) .
j. bobroff , n. laflorencie , l. k. alexander , a. v. mahajan , b. koteswararao , and p. mendels , phys .
lett . * 103 * , 047201 ( 2009 ) .
h. manaka , a. v. kolomiets , and t. goto , phys .
* 101 * 077204 ( 2008 ) .
h. manaka , h. a. katori , o. v. kolomiets , and t. goto , phys .
b * 79 * 092401 ( 2009 ) .
tao hong , a. zheludev , h. manaka , and l .- p .
regnault , phys .
b * 81 * 060410 ( 2010 )
. m. greven and r. j. birgeneau , phys .
lett . * 81 * , 1945 ( 1998 ) .
t. miyazaki , m. troyer , m. ogata , k. ueda , and d. yoshioka , j. phys .
* 66 * , 2580 ( 1997 ) .
a. furusaki , m. sigrist , p. a. lee , k. tanaka , and n. nagaosa , phys .
lett . * 73 * , 2622 ( 1994 ) .
e. westerberg , a. furusaki , m. sigrist , and p. a. lee , phys .
lett . * 75 * , 4302 ( 1995 ) .
n. nagaosa , a. furusaki , m. sigrist , and h. fukuyama , j. phys .
. jpn . * 67 * , 3724 ( 1996 ) .
r. yu , t. roscilde , and s. haas , phys .
lett . * 94 * , 197204 ( 2005 ) .
r. yu , t. roscilde , and s. haas , phys .
b * 73 * , 064406 ( 2006 ) . c. yasuda , s. todo , and h takayama , j. phys . soc . jpn . * 75 * , 124704 ( 2006 ) .
o. f. syljusen and a. w. sandvik , phys .
e * 66 * , 046701 ( 2002 ) .
a. w. sandvik , phys .
b * 66 * 024418 ( 2002 ) . f. cooper , b. freedman , and d. preston , nucl .
b * 210 * , 210 ( 1982 ) .
mikeska , u. neugebauer , and u. schollwck , phys .
b * 55 * , 2955 ( 1997 ) .
e. westerberg , a. furusaki , m. sigrist , and p. a. lee , phys .
b * 55 * , 12578 ( 1997 ) .
b. frischmuth , m. sigrist , b. ammon , and m. troyer , phys .
b * 60 * , 3388 ( 1999 ) .
n. laflorencie , d. poilblanc , and m. sigrist , phys .
b. * 71 * , 212403 ( 2005 ) .
the notion of local gap is associated to the local susceptibility at site @xmath229 , @xmath238 . assuming that @xmath239 one obtains that @xmath240 , where @xmath241 is the local gap .
d. stauffer , `` introduction to percolation theory '' , taylor and francis , london 1985 .
m. e. j. newman and r. m. ziff , phys .
85 * , 4104 ( 2000 ) ; m. e. j. newman and r. m. ziff , phys .
e * 64 * , 016706 ( 2001 ) .
d. s. fisher , phys .
b * 50 * , 3799 ( 1994 ) .
e. yusuf and k. yang , phys .
b * 65 * , 224428 ( 2002 ) .
r. mlin , y .- c .
lin , p. lajk , h. rieger , and f. igli , phys .
b * 65 * , 104415 ( 2002 ) . j. a. hoyos and e. miranda , phys . rev .
b * 69 * , 214411 ( 2004 ) . f. igli and c. monthus , phys .
rep . * 412 * , 277 ( 2005 ) .
it is to be said that , at the low temperature to which the figure refers ( @xmath242 ) , the correlation length on the doped 4-leg and 6-leg ladders grows to values which are no longer satisfying the condition @xmath243 , and hence it is limited by the system size ( @xmath164 ) explored there .