arxiv: v3 [cond-mat.stat-mech] 16 Aug 2012

Transcription

1 The Fortuin-Kasteleyn and Damage Sreading transitions in Random bond Ising lattices P. H. Lundow 1 and I. A. Cambell 2 1 Deartment of Theoretical Physics, Kungliga Tekniska högskolan, SE Stockholm, Sweden 2 Laboratoire Charles Coulomb, Université Montellier II, Montellier, France arxiv: v3 [cond-mat.stat-mech] 16 Aug 2012 The Fortuin-Kasteleyn and heat-bath damage sreading temeratures T FK() and T ds () are studied on random bond Ising models of dimension two to five and as functions of the ferromagnetic interaction robability ; the conjecture that T ds () T FK() is tested. It follows from a statement by Nishimori that in any such system exact coordinates can be given for the intersection oint between the Fortuin-Kasteleyn T FK() transition line and the Nishimori line, [ NL,FK,T NL,FK]. There are no finite size corrections for this intersection oint. In dimension three, at the intersection concentration [ NL,FK] the damage sreading T ds () is found to be equal to T FK() to within 0.1%. For the other dimensions however T ds () is observed to be systematically a few ercent lower than T FK(). PACS numbers: Lk, q, Cn, Cx INTRODUCTION The hysics of Ising sin glasses (ISGs) in the ordered regime below the freezing temerature T g has been intensively investigated for decades; the aramagnetic regime above T g has attracted less attention. However, in addition to the standard ordering (or freezing) transitions other significant critical temeratures within the aramagnetic regime can be defined oerationally and estimated numerically with high recision. We have studied random-bond Ising models (RBIM) in dimensions two, three, four and five having ferromagnetic near neighbor interactions with robability and antiferromagnetic interactions with robability 1 over the whole range of in the aramagnetic regime. The Hamiltonian is H = ij J ijs i S j where the sum is taken over all nearest neighbour bonds ij. For the discrete ±J interaction distribution the bond values are thus chosen according to P(J ij ) = δ(j ij J)+(1 )δ(j ij +J) (1) We will set J = 1 and use when convenient inverse temeraturesβ = 1/T. The ure Isingmodel is recoveredfor = 1 and the ure antiferromagnetic model for = 0, while the standard strong disorder bimodal sin glass case corresonds to = 1/2. In addition to the standard ferromagnetic and sin glass ordering transitions, on the RBIM hase diagram in the aramagnetic regime other hysically significant lines can be oerationally defined as functions of. In addition to the Griffiths line (T gr () = T c ) [1, 2] and the Nishimori line [3] there is the Fortuin-Kasteleyn (FK) transition line T FK () [4 9] and the heat-bath damage sreading transition line T ds () [10 13]. Early large scale relaxation measurements on the ISG in dimension three were interreted [14] in terms of a dynamic transition at the Griffiths temerature (i.e. T G () = T c ), with a qualitative change in relaxation behavior of the autocorrelation function from non-exonential to exonential. It can alternatively be considered that T FK or T ds definesadynamictransition. AtT FK thefkclustersizediverges so the standard cluster fliing algorithms [15, 16] break down. At T ds, the time scale for the couling from the ast equilibration criterion D(t) 0 diverges at large L [17, 18] so this criterion similarly breaks down. A riori it seems lausible that the two dynamic breakdowns should occur at similar, if erhas not identical, temeratures as has been conjectured by a number of authors [7, 19, 20]. In the following discussion it is shown that there are universal analytic relations for the intersection oint between T FK () and the Nishimori line [3]. With these exact values in hand together with the resent accurate numerical damage sreading data a critical test can be made of the conjectured rule T ds () T FK (). This is shown to be a good aroximation, almost exact in dimension three and accurate to a few ercent in dimensions two, four and five. However the equality between the temeratures is not a general rule. TRANSITION DEFINITIONS Nishimori [3] has shown for the RBIM that, due to extra symmetries of the roblem, a number of quantities may be comuted exactly when the equality (1 NL )/ NL = ex( 2β NL ) (2) holds; this condition defines the Nishimori line (NL), a line traversing the entire β() hase diagram including both aramagnetic and ferromagnetic regimes. In articular on the NL the internal energy of the system er bond (or edge) is U NL = [2 NL 1] = tanh(β NL ) (3)

2 2 This energy er bond has exactly the value that indeendent bonds would have at the same temerature. We will see below that indeed on the NL the ositions of the satisfied bonds are uncorrelated [21, 22] and that bond ositions remain very close to random over a wide strong disorder regime around = 1/2. In other words anywhere on the NL the satisfied bonds are uncorrelated on average, so a fortiori the FK active bonds are distributed at random. The rule for the Fortuin-Kasteleyn(FK) transition line [4 9], corresonding initially to the ure ferromagnets, is to first select for some articular equilibrium configuration the entire set of satisfied bonds ij where J ij S i S j is ositive. These bonds are then decimated at random leaving a fraction [1 ex( 2/T)] of active satisfied bonds. The FK stochastic transition at T FK occurs when the set of active bonds ercolates through the lattice. From the way in which the bonds are laid down, in the general case this is correlated bond ercolation, as oosed to the standard uncorrelated bond ercolation with random bond occuation. Remarkably, for ure ferromagnets it can be roved that the FK transition β FK ( = 1)coincidesexactlywiththeCurietemerature β c ( = 1) [4]. The Swendsen-Wang and Wolff cluster algorithms [15, 16] seed u equilibration dramatically as long as the FK clusters remain of finite size. The same oerational definition can also be used elsewhere in the hase diagram; in articular T FK ( = 1/2) was estimated numerically for strong disorder RBIM on square, cubic, and triangular lattices [6 9]. It turned out thatt FK ( = 1/2)wasineachcasemuchhigherthanthe standard sin glass ordering temerature T g ( = 1/2). Even in the two-dimensional case where the glass temerature is zero, T FK ( = 1/2) is similar to but lower than the Griffiths temerature T gr = T c ( = 1) [8, 9]. The hysical basis for the ercolation transition temerature was exlained in terms of FK drolets for the case of the fully frustrated lattice and general Potts q in Refs. [23, 24]. Imaoka et al. [25] estimated T FK () numerically over the entire range of for the square and triangular RBIM lattices. BelowwewillintroduceafurthertemeratureT FKr () closely related to T FK (). By construction the fraction of satisfied bonds n s (,T) in equilibrium is related to the equilibriuminternalenergyerbond U(,T)through U(,T) = 1 2n s (,T) so the fraction of all bonds which are FK active bonds is just [19] P a (,T) = (1 U(,T))[1 ex( 2/T)]/2. (4) U(,T) and hence P a (,T) can be readily measured numerically to high recision. It was noted [19] that if the satisfied bond ositions are assumed to be uncorrelated, the FK transition would occur at a temerature such that P a (,T FK ) P c where P c is the standard random bond ercolation concentration for the lattice. We will refer to the temerature where the condition P a (,T FKr ) = P c (5) is satisfied as T FKr () (r standing for random). This conjecture led to estimates for T FKr ( = 1/2) on square and on cubic lattices in higher dimensions [19] which were in good agreement with the numerical estimates for T FK ( = 1/2) available at the time. It can be noticed that for the ure Ising square lattice ferromagnet at criticality U(β c ) = 1/2 1/2 = while tanh(β c ) = which is very different because the ositions of the satisfied bonds are strongly correlated. However the critical FK concentration ofactive bonds is P a ( = 1,T c ) = 1/2which is accidentally equal to P c = 1/2 for this lattice. For the ure Ising ferromagnet on the triangular lattice at criticality on the other hand P a ( = 1,T c ) = which is not quite equal to P c = for this lattice. Finally, a heat-bath damage sreading transition, T ds () can be defined [10 13]. Heat-bath udate rules are alied at fixed temerature to two initially nonidentical sin configurations A(0) and B(0) (which are not necessarily equilibrated) of a given samle (meaning A and B have exactly the same sets of interactions J ij ); the random number used in each subsequent single sin udate ste is the same for both configurations. At high enough temeratures, on annealing under this rocedure for sufficient time t, A(t) and B(t) will become identical; below the damage temerature T ds the damage (D(t) the Hamming distance between A(t) and B(t) divided by L d ) will stabilize for long times at a temerature deendent non-zero value. For the ure ferromagnet case T ds ( = 1) is equal to the Curie temerature T c ( = 1) [10, 11]. For the strong disorder sin glass, just as T FK ( = 1/2) is much higher than the glass temerature T g ( = 1/2), T ds ( = 1/2) is also much higher than T g ( = 1/2) [10 13]. The damage sreading transition for a given couling algorithm, in articular heat bath (HB), can be described as a regular to chaotic dynamic transition from the viewoint of the couling from the ast aroach. When the damage falls to zero after a sufficient anneal time it is a guarantee that the system has been strictly equilibrated [18, 26, 27]. In systems with frustration this guarantee breaks down when T < T ds (), so erfect equilibration in this sense for the regime near the critical temerature cannot be achieved. It should be noted that T ds deends on the udating rotocol, with for instance Glauber udating giving very different results from HB udating. It turns out that even within the HB rotocol the recise value obtained for T ds changes slightly deending on whether sequential or random udating is used. The results reorted here are for random udating. With sequential udating the observed damage sreading temerature is of the order of 1% lower.

3 3 Early comarisons of numerical data from different grous indicated that T ds ( = 1/2) T FK ( = 1/2) for dimensions two and three [7, 19], and on this basis it has been conjectured by a number of authors [7, 19, 20] that T ds () T FK () is a general rule defining a joint dynamic transition temerature above which relaxation is exonential in the long time limit and below which relaxation is chaotic. It should be noted that in ractice it is very hard to identify such a transition directly from autocorrelation function decay q(t) data. Recently Yamaguchi[20] focused attention on behavior on the NL. He rovided conjectures suggesting that the intersection of the T FK () line and the NL would occur when and NL,FK = (1+P c )/2 (6) β NL,FK = ln[(1+p c )/(1 P c )]/2 (7) where again P c is the random ercolation concentration for the lattice. It turns out that because of the analytic value for the energy U NL () on the NL, Eq. (3), the Yamaguchi condition is strictly equivalent to the equality T NL,FK = T NL,FKr, meaning that this condition holds if the FK active bonds are distributed at random. The random bond conjecture for β FK, Eq. (4) of Ref. [19], is identical to Eq. (3.3) of Ref. [20]. In fact Nishimori [21] years earlier had made the strict statement : We have also roved indeendence of the local internal energy of different bonds, which indicates that the system effectively slits into uncorrelated sets of bonds on the [NL] in the hase diagram. In other words anywhere on the NL the satisfied bonds are uncorrelated on average, so a fortiori the FK active bonds are distributed at random. Hence at the NL-FK intersection concentration NL,FK, we have T FK T FKr (), meaning that the conditions for this intersection conjectured byyamaguchi, Eq.(6) and Eq.(7), holdexactlyfor any RBIM lattice if the FK transition is a random active bond ercolation transition. Remarkably, a further consequence of Nishimori s statement is that this equality should hold for the mean over many samles indeendently of the samle size L. The numerical results below show that for dimension two, and resumably in other dimensions also, in addition to this identity at NL,FK, for a wide range of around = 1/2, T FK () T FKr () remains a very good aroximation, as was conjectured for = 1/2 in Ref. [19]. The active bonds are thus essentially uncorrelated at β FK () in this wide high disorder regime. It should be noted however that as a general rule the ositions of active bonds are correlated so T FK () is not equal to T FKr (). Yamaguchi [20], following Refs. [7, 19] for the case of = 1/2, also conjectured that T ds (NL) = T NL,FK at the intersection. The resent measurements show that in dimension two T ds () is in fact systematically lower than T FK () by a few % over the entire range of, excet at and near the ure ferromagnet limit = 1. At NL,FK for cubic models in dimensions three to five the conjecture holds to within about 0.1% in dimension three, to 2% in dimension four and to 4% in dimension five. The general relation T ds () T FK () is therefore simly a reasonably good aroximation. In dimensions two, three and four the resent data show that the damage sreading transition temerature like the FK transition temerature is very insensitive to for a wide range of around = 1/2. FK NUMERICAL RESULTS IN DIMENSION TWO The two dimensional lattices studied were the square lattice and the triangular lattice. On both lattices the standard hase diagrams as functions of are well established. The ure ferromagnetic Curie temeratures T c ( = 1) are exactly T c (sq) = 2/ln[2 1/2 1] = and T c (tri) = 4/ln3 = resectively. When is lowered from = 1, the Curie temerature T c () drosgraduallyuntil acritical oint isreached at r on the Nishimori line where the Curie temerature tends suddenly to zero with weak re-entrant behavior [3, 28 32]. Between the re-entrant regime and = 1/2 the system can be considered a sin glass but with no finite temerature ordering transition. The square lattice hase diagram for < 1/2 is the exact mirror image of the hase diagram for > 1/2 with antiferromagnetic order taking the lace of ferromagnetic order as the lattice is biartite. For the triangular lattice on the other hand there is no finite temerature order below = 1/2 as the fully frustrated antiferromagnetic limit at = 0 is aroached. The internal energy er bond U(,T,L) at equilibrium was estimated numerically. For small triangular lattices, 3 L 11, the energy as a function of temerature for each secific samle was calculated exactly (u to numerical recision) using transfer matrices. Once a samle is icked we slit the lattice into L arts, each art corresonding to a transfer matrix A i. After choosing a numerical value of T we evaluate the matrices and comute the trace of their roduct to obtain Z(,T,L) = tr(a 1 A 2,A L ). This standard aroach is described in great detail for the Ising case in section 3 of [33] and is of course easily adated to the sin glass case. We thus comuted lnz(,t,l) for = 0,0.02,0.04,...,0.50 and = 0.51,...,0.98,0.99,1 and some values of T. The number of samles ranged between for L = 3,4,5 and then down to only 1024 for L = 11. Taking the average lnz/ β then rovides us with U(,T,L). Values at (,T) outside

4 4 the comuted data grid were obtained through 3rd order interolation. For larger L (for the square and triangular lattice) the Monte Carlo data were collected after equilibration using standard Metroolis udating. At each (,T), and for each samle, we collected at least a few million measurements of U, and of course even more for moderate L. For the square lattices we used L = 16,32,64 on 64 different samles with 0.5 < 1 in stes of For the triangular lattice we used L = 16,32 on 32 different samles but with 0 < < 1 in stes of (slightly denser at high and low ). Some 50 values of T were used for both lattices. Again, intermediate values of U were obtained by interolating the values at the (,T) grid oints. By interolation of oints derived from Eq. (4) and the consistency condition Eq. (5) the temerature T FKr () where the fraction of FK active bonds P a (,T) is equal to P c can be estimated to high recision. In Figures 1 and 2 the T FKr () values are comared to the directly measured T FK () values [8, 25] for the square and triangular lattices resectively. The damage sreading temeratures to be discussed later are also shown. The error in the T FKr -estimates is found by first solving Equation (4) for each individual samle (using the interolation we mentioned above) which gives us a samle-to-samle standard deviation from which we obtain the standard error of the average-samle T FKr. Once the standard error is established we are free to lot a smoothed version (from fitting a olynomial of high degree) of T FKr versus. Needless to say we have not estimated the error in U, and hence T FKr, for each individual samle, but we simly assume this error should be reflected in the samle-to-samle variation. In dimension two, P c = 1/2 for the square lattice and P c = 2sin(π/18) = for the triangular lattice. The exact NL-FK intersection values from Eqs. (6) and (7) are then : NL,FK = 3/4, β NL,FK = ln(3)/2 = for the square lattice, and NL,FK = , β NL,FK = for the triangular lattice. The estimated T FK () and T FKr () lines run directly through these exact NL-FK intersection oints as they should, Figures 1 and 2. The measured equilibrium energy er bond is slightly higher than the random bond energy tanh(β) for β > β FKr and slightly lower for β < β FKr. To within the high numerical recision the estimated β FKr ( NL,FK,L) at the NL-FK intersection oint is indeendent of L down to L = 3, Figures 3 and 4. The observed absence of finite scaling corrections at the NL also follows from Nishimori s general statement quoted above [21]. For -values to the left and right of NL,FK there are weak finite size effects of oosite signs. This absence of mean finite size scaling deviations arises because T FKr deends only on the energy. Other arameters can still show finite size scaling deviations[34]. The large size limits for T FK () and T FKr () remain equal to each other T FIG. 1: (Color online) The square lattice transitions, temeratures T as functions of ferromagnetic interaction concentrations. The ink dashed line on the right indicates the ferromagnetic Curie temeratures, data taken from Ref. [32]. The green solid traversal line is the exact Nishimori line (NL). Red circles : FK transitions T FK, Ref. [25]. Blue squares : random active bond ercolation line T FKr. Black diamonds : heat bath damage sreading transition T ds. The errors on T FK and T ds are of the size of the oints. Large black square : exact intersection oint T NL,FK. to within the numerical recision of the T FK () oints [25] over the strong disorder range of extending from = 1/2to NL,FK inthe squarelattice, andfrom 0.3 to NL,FK in the triangular lattice. Over these ranges of the satisfied bonds at T FK are very close to being uncorrelated; for instance at = 1/2 on the square lattice U(,β FKr )/tanh(β FKr ) = which remains near to the uncorrelated value of 1. The further randomness introduced by the FK decimation is sufficient to render the active bond ositions essentially uncorrelated. The T FK () [25](or T FKr () which has been measured here with higher recision) is rather insensitive to within this range, so that the measured T FK (1/2) is similar to the exact T NL,FK on both lattices. For both lattices in the more strongly ferromagnetic ranges > NL,FK, and in the range < 0.3 near to the fully frustrated = 0 limit in the triangular lattice, T FK () > T FKr (). At these concentrations the satisfied bonds are significantly correlated at T FK () and the FK decimation does not sufficiently comensate the correlations so as to roduce randomness among the active bonds.

5 5 T T FKr () FIG. 2: (Color online) The triangular lattice transitions, temeratures T as functions of ferromagnetic interaction concentrations. The ink dashed line on the right indicates the ferromagnetic Curie temeratures. The green solid traversal line is the exact Nishimori line (NL). The red dashed vertical line in the centre incates = 0.5. Red circles : FK transitions T FK, Ref. [25]. Blue squares : random active bond ercolation line T FKr. Black diamonds : heat bath damage sreading transition T ds. The errors on T FK and T ds are of the size of the oints. Large black square : exact intersection oint T NL,FK. T FKr () FIG. 3: (Color online) The size effect on triangular lattice temeratures T FKr at which the fraction of FK sites is equal to the random site ercolation concentration for the lattice P c, Eq. (5). Temeratures T FKr as functions of ferromagnetic interaction concentrations for lattice sizes L from 3 to 11. L = 3,4,5,11 black, red, green, blue; from to to bottom on the right hand side and from bottom to to on the left hand side. The black square is the exact T NL,FK intersection Eqs. (6) and (7). FIG. 4: (Color online) As for Fig. 3, a closeu of the region around the T NL,FK intersection (square). L = 3,4,5,11 (black square, red circle, green triangle, blue diamond; from to to bottom on the right hand side and from bottom to to on the left hand side). The error bars are the standard error when measuring T FKr after solving Eqs. (4) and (5) and corresond to samles for L = 3,4,5 and 1024 for L = 11. The curves are 10th degree olynomials fitted to the whole range from = 0 to = 1. DAMAGE SPREADING NUMERICAL RESULTS IN DIMENSION TWO For dimension two the critical temeratures T ds () belowwhichthelongtimelargelheatbathdamagesreading D(,T,t) tends to a non-zero value are also shown in Figures 1 and 2. The values given are obtained from extraolating data at increasing sizes, Figures 5 and 6, to infinite size. On the basis of data for square and cubic lattices at = 1/2 it was conjectured earlier that T FK () T ds () [7, 19, 20] as is the case in the ure ferromagnet = 1 limit. However the resent data demonstrate that while the curves for T ds () and T FK () lie near together and are of very similar shae, in articular both being almost indeendent of for the range near = 1/2, it is clear that T ds () < T FK () excet when tends to the ferromagnetic limit = 1. So unfortunately it is not ossible to define a unique joint dynamic transition temerature. We will see that this is true also in higher dimensions. The resent accurate critical value for the strong disorder square lattice T ds (1/2) = 1.69(2) is in good agreement with the earlier value T ds (1/2) = 1.70 [13] and is very close to the regular to chaotic dynamic transition temerature estimated for the same lattice in Ref. [18] from the divergence of the couling time with L 2, where the data indicate a transition at T ds (1/2) 1.72.

6 D(,L) D(,L) FIG. 5: (Color online) The damage sreading D(β,L) for the triangular lattice at = The lattice sizes are L = 48,64,96,128,256,384,512,1024 (black square, red circle, green triangle, blue inverted triangle, cyan diamond, yellow left triangle, ink right triangle, brown star. Right to left in the region where the curves searate). Shown error bars are the standard error. The exact FK transition inverse temerature is β = so clearly distinct from the limiting damage sreading inverse temerature. DIMENSIONS THREE, FOUR AND FIVE It is numerically much more demanding to estimate T FK () recisely through direct measurements (as in [25]), in articular allowing for finite size effects, than to estimate T ds () to the same level of accuracy. For cubic lattices in dimensions three, four and five the random bond critical concentrations P c though not exact have been estimated to very high recision, P c (3) = (5), P c (4) = (10) and P c (5) = (3) resectively [35, 36]. Using the exact NL-FK intersection oint exressions, Eqs. (6) and (7), one thus has NL,FK = and β NL,FK = in dimension three, NL,FK = and β NL,FK = in dimension four, and NL,FK = and β NL,FK = in dimension five. Hence as NL,FK and β FK ( NL,FK ) can be taken as known almost exactly, it is sufficient to estimate β ds ( NL,FK ) numerically (allowing carefully for finite size corrections) to obtain an accurate estimate of theratiot ds ( NL,FK )/T FK ( NL,FK )foreachdimension. The equilibrium damage D(L, β) was measured for given L as a function of β, and the results were extraolated to obtain an estimate of the β ds value at which D(,β) falls to zero, Figures 5 and 6. There is a clear enveloecurve forall Lin eachcase, with the dataoints leaving the curve later and later as L increases. There are various ways to extraolate to infinite L. One effi- FIG. 6: (Color online) The cubic lattice equilibrium damage sreading D(β,L) as a function of size and inverse temerature for = Sizes : L = 12,16,24,32,48,64 (black square, red circle, blue triangle, green inverted triangle, cyan diamond, ink star ; right to left in the region where the curves searate). Shown error bars are the standard error. The extraolated intersection with the β axis gives the infinite size critical β ds. In this case β ds = (15) can hardly be distinguished from the exact FK transition inverse temerature β FK = , red arrow. cient method is to fit the enveloe oints assuming that near β ds the behavior follows D(β) = A(β β ds ) B and adjusting B to obtain a straight line. We have also used a simle scaling formula β ds (L) = β ds + CL λ to verify. Of course, a certain statistical uncertainty will enter deending on which L to include in the fitting rocess which we take into account in the final error estimate. The infinite L damage sreading temeratures T ds ( NL,FK ) were estimated for the central values of NL,FK in the three dimensions : T ds ( NL,FK ) = 3.932(2) in dimension three, T ds ( NL,FK ) = 6.057(10) in dimension four, and T ds ( NL,FK ) = 8.13(1) in dimension five. TheobservedratiosT FK ( NL,FK )/T ds ( NL,FK )are equal to (5), 1.022(2) and 1.036(10) resectively. IndimensionthreeT ds ( NL,FK )isindistinguishablefrom T FK ( NL,FK ), while in dimensions four and five the values aear tantalizingly close but not identical. We have no exlanation for the striking similarity between the two temeratures for the articular case of dimension three. Adding a little more detail for this articular case, the individual β ds (L) were estimated to an accuracy of ± for L = 48,64, ± for L = 32, ± for L = 16,20,24 and ±0.001 for L = 12. Using the exonent λ = 2.05 gives a rojected β ds that deends only to a very small degree on which L are included in the fit (though we always include L = 48,64). We receive a median value of β ds = and a standard deviation of

7 giving us T ds = 3.932(2). Similar methods were used to estimate T ds for the other dimensions. As in dimension two, the curve for T ds () is almost flat over a wide range of around = 1/2 for dimensions three and four, Figures 7 and 8 (the case of dimension five was not studied). These data are for fixed L (L = 32 and L = 16 resectively) so the T ds () values are not quite equivalent to the values quoted for T ds at NL,FK. Judging from the behavior observed in dimension two, it seems very lausible to assume that the ratio T FK ()/T ds () is always ractically indeendent of in the strong disorder regime 1/2 < < NL,FK. Itcanbenotedthatnumericaldataforthetimedeendence of the autocorrelation function q(t) at all temeratures in two dimensional fully frustrated systems have been interreted in terms of exonential relaxation with logarithmic corrections due to vortex-vortex interactions [37]. The results indicate that in these systems (which have well established FK and damage transitions) there is no dynamic critical temerature in the Ogielski sense [14]. T ds () FIG. 7: The damage sreading critical temerature T ds () for a 3d cubic lattice size L = 32. The error bars corresond mainly to the residual uncertainty in the extraolation to the intersection with the β axis. DAMAGE CLUSTERS It is of interest to examine in some detail the damage clusters in dimension two where they can be readily visualized. Keeing β > β ds and letting t (or at least allowing for equilibration) how does the damage actually sread? Note that we define a cluster as a connected comonent in the lattice induced by the damaged sites, i.e. a maximal set of damaged sites such that there is a lattice ath(made u of horizontal and vertical stes) between each air. The sum of the individual cluster sizes is thus the damage without normalization, i.e. L d D(β). We have collected data on cluster sizes for d = 2,3,4,5. For d = 2 we have only used the square lattice, not the triangular lattice. What we see in a square lattice snashot at β near β ds and at any arbitrary fixed time is not, as in the case for ercolation(including FK ercolation) or random grahs near c, the formation of a giant cluster, but rather a large number of very small clusters. At subsequent times the clusters evolve and flutter through the lattice. An examle is shown in Fig. 9 The measurements take lace as above for the damage sreading; starting from random (infinite temerature) sin configurations, at each time ste we udate L d randomly selected sites, then search for all clusters and their sizes. On a measurement we collect the number of clusters, the average cluster size, the size of the largest cluster and the size of a randomly selected cluster. As before, we collect 10 6 measurements from each samle after discarding between and time stes (deending on L) to allow for equilibration of the damage sreading. We then average over the time stes to get the average T ds () FIG. 8: The damage sreading critical temerature T ds () for a 4d cubic lattice size L = 16. The error bars corresond mainly to the residual uncertainty in the extraolation to the intersection with the β axis. for a articular samle. The data are then averaged over the samles; we haveonly used eightsamlesin eachcase for the cluster measurements (for the damage seading we used between 8 and 128 samles deending on lattice size). In all cases we have set = NL,FK. First we discuss the exected number of clusters, which wedenoten c (β,l). In Figure10welotn c (β,l)/l 3 versus β for the simle cubic lattice. Equivalent behavior is found also for d = 2,4,5. Note that for the square

8 number of clusters FIG. 9: An illustrative snashot of an instantaneous damage site configuration for an L = 128 square lattice at = 0.75 and temerature T = 1.61, just below the damage sreading transition. FIG. 10: (Color online) The number of clusters n c(β,l) on the 3d cubic lattice normalized by L 3 at = NL,FK. Sizes L = 12,16,24,32,48 shown as black squares, red circles, blue triangles, green inverted triangles, ink left triangles. The errors are smaller than the size of the oints. lattice all lattice sizes agree on a global maximum robability located at β max = 0.656(1), where n c (β max,l) (1)L 2. For d = 3 we obtain β max = 0.313(1) and n c (β max,l) (1)L 3, for d = 4 we get β max = 0.211(1) and n c (β max,l) (1)L 4, while d = 5 gave β max = 0.177(1) and n c (β max,l) (1)L 5. Thus the damage is always distributed over O(L d ) clusters. With so many clusters the average cluster at each time ste must be rather small. We measure the number of damaged sites DL d and the number of clusters n c, giving us the mean cluster size at each time ste. The normalized time average (and then the samle average, though the samle variation is very small), is then s(β,l) = DL d /n c. In Figure 11 we lot s(β,l) versus β, again for d = 3. The curve clearly suggests a ositive right limit at β ds. To estimate this limit we fit a simle exression a 0 +a 1 ex(a 2 β) to the oints resulting in s(β ds ) = lim β β + lim L s(β,l) = 3.46(1). We ds estimated the right limits corresonding to d = 2,3,4,5 to be resectively 9.00(2), 3.46(1), 2.76(1) and 1.95(1) sites at β ds. One could alternatively define s as the average damage divided by the average number of clusters. This is not strictly the same as our resent definition (mean ratio versus ratio of the means) but the difference is of course vanishingly small here. For examle, the maximum in Figure 10 at β max = gives n c L 3 and the number of damaged sites, see Figure 6, is 0.139(1)L 3. Hence the average cluster size is mean cluster size FIG. 11: (Color online) The mean cluster size s(β,l) on the 3d cubic lattice at = NL,FK. Sizes L = 12,16,24,32,48 shown as black squares, red circles, blue triangles, green inverted triangles, ink left triangles. The errors are smaller than the size of the oints. roughly 6.1 which matches Figure 11. We end this section with a final remark concerning the distribution of cluster sizes. At each time ste we here ick one cluster uniformly at random and measure its size. In Figure 12 we show a set of size distributions (or density functions) as logpr vs logs for a range of β

9 9 ln((n s )) ln(n s ) FIG. 12: (Color online) The distribution of cluster sizes on the 3d cubic lattice for L = 12 at inverse temeratures β = 0.295, 0.335, 0.370, (black squares, red circles, green triangles, and blue inverted triangles). The errors are of the size of the oints excet for very small where they can be judged by the scatter (note the logarithmic scales). for d = 3,L = 12. They all have a high value at size 1 (many isolated sites) and then dro quickly. They can be exected to behave like this since the mean size is between 4 and 70 in this temerature range. However, for β > β max (or thereabout) the distributions show a second, and rather wide, maximum located in the neighborhoodofthetotal damaged(β,l). This canbe understood as follows. Suose for the sake of argument that the damaged sites are distributed at random in sace. (This is only aroximate as correlations between damage sites should be allowed for). Then when D(β) exceeds the site ercolation concentration there will exist a single giant ercolating cluster of damaged sites together with residual small clusters. As D(β) increases further the ercolating cluster will contain almost all the damaged sites. On this criterion the eak should aear in the distribution at D(β) 0.59, 0.31, 0.20 and 0.14 in dimensions 2,3,4 and 5 resectively which gives an indication in rough agreement with the data. Indeed in the caseofthe squarethe eak neveraearswhich is consistent with the fact that D(β) never aroaches For the dimensions where the eak does aear the global arameters, in articular D(β), increase smoothly with β and show no sign of any critical behavior as the giant cluster forms. CONCLUSION The exact values of the coordinates of the intersection oint where the Fortuin-Kasteleyn transition line crosses the Nishimori line can be derived for an RBIM from the analytic condition that satisfied bonds are uncorrelatedonthe NL [21]. On anylattice this leadsto the exact exressions NL,FK = (1 + P c )/2 and T NL,FK = 2/ln[(1+P c )/(1 P c )] [19, 20] where P c is the standard random bond ercolation concentration for the articular lattice. For lattices in dimension two (and robably in higher dimensions also) the uncorrelated bond condition remains a very good aroximation at the FK transition temerature over a wide strong-disorder region sanning = 1/2. In ure ferromagnets T FK ( = 1) = T ds ( = 1) = T c ( = 1). The conjectured equivalence for the RBIM between the FK transition temerature and the heat bath damage sreading temerature, T FK () T ds () [7, 19, 20], searating an exonential from a chaotic dynamic regime has been tested at = NL,FK on simle cubic lattices in dimensions two, three, four and five. It holds to within 0.1% in dimension three, to within 2% in dimension four, and to 3% in dimension five. The equivalence aears always to be a good aroximation. We have no exlanation to roose for the quasi-equality in the case of dimension three. For the square and triangle lattices the difference is larger : 5% and 16% resectively. The FK transition in the strong-disorder regime close to = 1/2 can thus be said to be well understood. However the basic hysical condition determining the the damage sreading transition temerature, which lays an imortant role in limiting erfect equilibration in RBIMs, and its roximity to the FK transition, remains unclear. ACKNOWLEDGEMENTS We are very grateful for enlightening remarks by H. Nishimori, C. Yamaguchi, R. Ziff and W. Krauth. The comutations were erformed on resources rovided by the Swedish National Infrastructure for Comuting (SNIC) at High Performance Comuting Center North (HPC2N). [1] R.B. Griffiths, Phys. Rev. Lett. 23, 17 (1969). [2] Y. Matsuda, H. Nishimori and K. Hukushima, J. Phys. A: Math. Theor. 41, (2008). [3] H. Nishimori, Prog. Theor. Phys. 66, 1169 (1981). [4] C. M. Fortuin and P. W. Kasteleyn, Physica 57, 536 (1972). [5] A. Coniglio and W. Klein, J. Phys. A 13, 2775 (1980).

10 10 [6] A. Coniglio, F. di Liberto, G. Monroy, and F. Peruggi, Phys. Rev. B 44, (1991) [7] L. de Arcangelis, A. Coniglio, and F. Peruggi, Eurohys. Lett. 14, 515 (1991). [8] V. Cataudella, Physica A 183, 249 (1991). [9] M. Zhang, C.Z. Yang, Eurohys. Lett. 22 (1993) 505. [10] B. Derrida and G. Weisbuch, Eurohys. Lett. 4, 657 (1987). [11] B. Derrida, Phys. Re. 184, 207 (1989). [12] L. de Arcangelis, A. Coniglio, and H. Herrmann, Eurohys. Lett. 9, 749 (1989). [13] I. A. Cambell and L. de Arcangelis, Physica 178, 29 (1991). [14] A. T. Ogielski, Phys. Rev. B 32, 7384 (1985). [15] R. H. Swendsen and J.-S. Wang, Phys. Rev. Lett. 58, 86 (1987). [16] U. Wolff, Phys. Rev. Lett. 62, 361 (1989). [17] J. Pro and D. Wilson, Random Struct. Algorithms, 9, 223 (1996). [18] E.P. Bernard, C. Chanal, and W. Krauth, Eurohysics Lett. 92, (2010). [19] I.A. Cambell and L. Bernardi, Phys. Rev. B 50, (1994). [20] C. Yamaguchi, arxiv: [21] H. Nishimori, Prog. Th. Phys. 76, 305 (1986). [22] H. Nishimori, Statistical Physics of Sin Glasses and Information Processing: An Introduction (Oxford University Press, 2001),.51. [23] S. Prakash, A. Coniglio, and H.E. Stanley, Phys. Rev. E 49, 2742 [24] V. Cataudella, A. Coniglio, L. de Arcangelis, and F. di Liberto, Physica A 192, 167 (1993). [25] H. Imaoka, H.Ikeda, and Y.Kasai, PhysicaA246 (1997) 18 [26] C. Chanal and W. Krauth, Phys. Rev. Lett. 100, (2008). [27] C. Chanal and W. Krauth, Phys. Rev. E 81, (2010). [28] F.D. Nobre, Phys. Rev. E 64, (2001). [29] C. Wang et al., Ann. Phys. 303, 31 (2003). [30] C. Amoruso and A. K. Hartmann, Phys. Rev. B 70, (2004). [31] F. Parisen Toldin et al., J. Stat. Phys. 135, 1039 (2009). [32] C.K. Thomas and H.G. Katzgraber, Phys. Rev. E 84, (R) (2011) [33] R. Häggkvist and P.H. Lundow, J. Stat. Phys. 108,429 (2002). [34] M. Hasenbusch, F.P. Toldin, A. Pelissetto, and E. Vicari, Phys. Rev. E 77, (2008). [35] C.D. Lorenz and R.M. Ziff, Phys. Rev. E 57, 230 (1998). [36] S.M. Dammer and H. Hinrichsen, J. Stat. Mech: Theory Ex. 7, (2004). [37] J.-C. Walter and C. Chatelain, arxiv: