arxiv: v1 [cond-mat.dis-nn] 6 Feb 2008

Transcription

1 Reentrant and orward hase Diagrams of the nisotroic Three-Dimensional Ising Sin Glass an Güven,. Nihat Berker, Michael Hinczewski, and Hidetoshi Nishimori Deartment of hysics, Koç University, Sarıyer 5, Istanbul, Turkey, Deartment of hysics, Massachusetts Institute of Technology, ambridge, Massachusetts 9, U.S.., eza Gürsey Research Institute, TÜBITK - Boshorus University, Çengelköy, Istanbul, Turkey, and Deartment of hysics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, Tokyo 5-55, Jaan arxiv:.77v [cond-mat.dis-nn] eb The satially uniaxially anisotroic d = Ising sin glass is solved exactly on a hierarchical lattice. ive different ordered hases, namely ferromagnetic, columnar, layered, antiferromagnetic, and singlass hases, are found in the global hase diagram. The sin-glass hase is more extensive when randomness is introduced within the lanes than when it is introduced in lines along one direction. hase diagram cross-sections, with no Nishimori symmetry, with Nishimori symmetry lines, or entirely imbedded into Nishimori symmetry, are studied. The boundary between the ferromagnetic and sin-glass hases can be either reentrant or forward, that is either receding from or enetrating into the sin-glass hase, as temerature is lowered. However, this boundary is always reentrant when the multicritical oint terminating it is on the Nishimori symmetry line. S numbers: 75..Nr,..aq, 5.7.h, 5..c I. INTRODUTION The Ising sin glass [] yields a hase diagram with a distinctively comlex ordered hase, in d =. wide accumulation of methods and results has occurred for this system. Most remarkably, in site of its high satial dimension and comlex ordering behavior, exact or recise information is being obtained for this system.[,,, 5,, 7,, 9,,, ] Thus, in the hase diagram in terms of temerature and concentration of antiferromagnetic bonds, the occurrence of the Nishimori symmetry line has been deduced [, ] and the accurate location of the multicritical oint has been redicted [, ]. urthermore, in systems with the Nishimori symmetry, it has been shown that the ferromagnetic hase cannot extend to antiferromagnetic bond concentrations beyond that of the multicritical oint.[, ] The two remaining otions being a straight line or a reentrance situation, subsequent works [, ] on hierarchical lattices have shown that for these systems, the singlass hase diagram is reentrant, namely that below the multicritical oint, the ferromagnetic hase recedes from the sin-glass hase as temerature is lowered. Exact results recently have also been extended to otts sin glasses.[5] satially uniaxially anisotroic d = system is studied in this work, to our knowledge the first study of quenched randomness and frustration in a satially anisotroic higher-dimensional system. In fact, both anisotroy and quenched randomness have acquired increased relevance from high-temerature suerconductivity results.[, 7] Our calculation is exact for a hierarchical lattice and aroximate for a cubic lattice. We find a rich hase diagram (e.g., ig.) with five different ordered hases, namely with ferromagnetic, antiferromagnetic, layered, columnar, and sin-glass order. The sin-glass hase is more extensive when randomness is introduced within the lanes than when it is introduced in lines along one direction. The global hase diagram includes cross-sections with no Nishimori symmetry, cross-sections with Nishimori symmetry lines, and a cross-section entirely imbedded within Nishimori symmetry. Thus, the multicritical oint between the sin-glass, ferromagnetic, and aramagnetic hases, reviously found to occur on the Nishimori symmetry line, is also found here at oints with no Nishimori symmetry, but renormalizes to a fixed distribution of interaction robabilities that obeys Nishimori symmetry. Nevertheless, we find that the boundary between the ferromagnetic and sin-glass hases can be either reentrant or forward, that is either receding from or enetrating into the sin-glass hase, as temerature is lowered. When the multicritical oint is not on the Nishimori symmetry line, the ferromagnetic-sin glass boundary can be reentrant or forward. However, when the multicritical oint is on the Nishimori symmetry line, this boundary is always reentrant [, ], consistently with the rigorous result [, ]. II. UNIXIY NISOTROI SIN GSS The uniaxially anisotroic Ising sin-glass system has the Hamiltonian βh = Kijs u i s j, () u ij u where s i = ± at each site i, ij u denotes a sum over nearest-neighbor airs of sites along the z direction (u = z) or in the xy lane (u = xy), and the bond strengths K u ij are equal to Ku > with robability u and K u with robability u, resectively corresonding to ferromagnetic and antiferromagnetic interaction. When imbedded into a cubic lattice, the Hamiltonian () yields a uniaxially anisotroic d = system.

2 = =5. =. =.5 =. = = =. =.9 =. = = =. =.7 =. = IG. : (olor on-line) To line: onstant-temerature cross-sections of the global hase diagram for / =.5, as a function of and, which are the concentrations of antiferromagnetic xy and z bonds, resectively. t low temeratures (high ), the central sin-glass () hase searates the corner ferromagnetic (), columnar (), antiferromagnetic (), and layered () hases. The diagrams are twofold symmetric along each axis, but not fourfold symmetric, due to the difference between longitudinal ( = ) and transverse ( = ) sin glasses. s temerature increases, the aramagnetic () hase aears at the central oint, first reaches the transverse sin-glass system and eliminates the sin-glass hase, then reaches the longitudinal sin-glass system and eliminates the sin-glass hase. In the latter system, the sin-glass and aramagnetic hases simultaneously occur for a very narrow range of temeratures, as also seen in the inset in the lower left anel of ig.. Middle line: onstant-temerature cross-sections for =. Bottom line: onstant-temerature cross-sections for / =. ll hase transitions in this figure are second order. Hierarchical lattices are d-dimensional lattices yielding exact renormalization-grou solutions to comlex statistical mechanics roblems. These lattices are constructed by the reeated self-imbedding of a grah into a bond [, 9, ]. The shortest ath between the external vertices of the grah gives the length rescaling factor b and the number of bonds in the grah gives the volume rescaling factor b d, from which the dimension d is determined. Hierarchical lattices have been used to study a wide variety of roblems, including chaotic rescaling [, ], sin-glass [], random-field [], Schrödinger equation [], lattice-vibration [5], dynamic scaling [], randomresistor network [7], aeriodic magnet [], comlex hase diagram [9], directed-ath [], heteroolymer [], directed-olymer [], and, most recently, scale-free and small-world network [,, 5,, 7,, 9, ] systems, etc. More recently, hierarchical lattices have been created [] for the study of satially anisotroic systems. The mutual reeated self-imbedding of two aroriately chosen grahs, with differentiated interac-

3 The quenched robability distribution ( K) of the relacing bond is calculated by the convolution ( K) = dk I dk II I (K I ) II (K II )δ( K R(K I, K II )), () IG. : onstruction of the uniaxially anisotroic d= hierarchical model. Two grahs are mutually and reeatedly selfimbedded. Note that for =, =, and =, the system reduces resectively to the d =, isotroic d = and d = systems. tions, yields a uniaxially anisotroic system, whereas a higher number of grahs is needed to achieve higher satial anisotroy.[] These hierarchical systems must reduce to isotroy and/or lower satial dimensions when corresonding interactions are set equal to each other or to zero, as illustrated in ig.. n anisotroic hierarchical lattice has already been used to obtain the hase diagram of the uniaxially anisotroic d = tj model of electronic conduction.[] When imbedded into the hierarchical lattice of ig., the Hamiltonian () yields a uniaxially anisotroic d = sin-glass system that is exactly soluble. III. EXT RENORMIZTION-GROU SOUTION: OWS O THE QUENHED DISTRIBUTIONS O THE NISOTROI SIN-GSS INTERTIONS The renormalization-grou solution roceeds in the direction oosite to the construction of a hierarchical model. Each grah is relaced by a renormalized bond via summation over the sins on the internal sites of the grah. This is achieved by a combination of two tyes of stes: the relacement, by a single bond K ij, of two bonds that are either in arallel, referred to as bondmoving: or in series, referred to as decimation: K ij = K I ij + KII ij, () K ik = [ ] ln cosh(kij + K jk ). () cosh(k ij K jk ) where R(K I, K II ) is the right-hand side of Eq.() or (), K I and K II are the interactions entering the right-hand side of either of these equations, with quenched robability distributions I (K I ) and II (K II ).[, ] ccordingly, the renormalization of xy is obtained as follows, following the uer ig. in the direction oosite to the arrow: (i) from the bond-moving of xy with itself, obtaining ; (ii) from the bond-moving of z with itself, obtaining ; (iii) from the decimation of and, obtaining ; (iv) from the decimation of and, obtaining ; (v) from the decimation of xy and xy, obtaining 5 ; (vi) from the decimation of and, obtaining ; (vii) from the decimation of and, obtaining 7 ; (viii) from the decimation of 5 and xy, obtaining ; (ix) from the bond-moving of and 7, obtaining 9 ; (x) from the bond-moving of 7 and 7, obtaining ; (xi) from the bond-moving of 9 and, obtaining (xii) finally, from the bond-moving of and, obtaining the renormalized quenched distribution xy. Thus, in each renormalization-grou ste, the renormalized distribution xy is obtained from the convolutions of 7 unrenormalized distributions xy and z. The renormalized distribution z is similarly obtained from the convolutions of 7 unrenormalized distributions xy and z, but with a different sequencing dictated by the lower ig.. The renormalization-grou transformations of the quenched robability distributions xy and z, given in the receding aragrah, are imlemented numerically, resulting in a distribution of interaction-strength values and a robability associated with each value, namely a histogram. Thus, the initial ±K u doubledelta distribution functions, described after Eq.(), are of course not conserved under the scale coarsening of the renormalization-grou transformation. The number of histograms increases after each convolution. When a maximum number of histograms, set by us, is reached, a binning rocedure is alied [, ]: Before each convolution, the range of interaction values is divided into bins, searately for ositive and negative interactions. The interactions falling into the same bin are combined according to their relative robabilities. The convolution then restores the set maximum number of histograms. In this work, we have used the maximum number of 9, for histograms for each distribution xy and z.

4 / / xy = = = =.5 IG. : (olor on-line) Temerature-concentration hase diagrams for isotroically mixed (uer left), transverse (uer right), longitudinal (lower left), and =.5 sin-glass systems. In all cases, / =.5. The uer left and right hase diagrams are seen to be, resectively, reentrant and forward, namely with a ferromagnetic hase that, resectively, recedes from or roceeds towards the sin-glass hase as temerature is lowered, as clearly seen in the insets. There are no oints obeying Nishimori symmetry in the hase diagrams of this figure. Note the remarkably narrow sin-glass hase, reaching zero-temerature, in the longitudinal singlass system, as also seen in the inset. ll hase transitions in this figure are second order. hase < + > < Kxy > < Kz + > < K z > erro + + ntiferro olumnar + ayered + Sin Glass + + ara TBE I: Sinks of the renormalization-grou flows in the different hases. These sinks are characterized here in terms of the average ositive and negative interactions of their limiting quenched robability distribution. / / / / = = = =.... IG. : (olor on-line) Zero-temerature hase diagrams of the longitudinal (left column) and transverse (right column) sin-glass systems. With the aroriate reversal in variables, the transverse and longitudinal sin-glass hase diagrams are seen here to be qualitatively similar, but quantitatively different. The sin-glass hase is more extensive in the transverse case. ll hase transitions in this figure are second order. transitions.[, ] We find six different hases for this system, with corresonding sinks characterized in Table I in terms of the average ositive and negative interactions of the limiting distribution. These hases are the ferromagnetic, antiferromagnetic, layered, columnar, sin-glass ordered hases and the disordered aramagnetic hase. In the layered hase, the sins are mutually aligned in each xy lane; these lanes of mutually aligned sins form an antiferromagnetic attern along the z direction. In the columnar hase, the sins are mutually aligned along the z direction; these lines of mutually aligned sins form an antiferromagnetic attern along the xy directions. Both of these hases are thus distinct from the antiferromagnetic hase, which is antiferromagnetic in all three directions. There is a single sin-glass hase, extending to anisotroic systems. / / IV. HSE DIGRMS ND IXED DISTRIBUTIONS We have obtained the global hase diagram of the uniaxially anisotroic d = sin-glass system in terms of the original interactions and robabilities (,,, ). In each thermodynamic hase, quenched robability distributions flow, under reeated renormalization-grou transformations, to a limiting behavior (sink) characteristic of that thermodynamic hase. hase boundary oints flow to their own characteristic (unstable) fixed distributions, shown below. nalysis at these unstable fixed distributions yields the order of the hase. hase Diagrams with no Nishimori Symmetry ross-sections of the global hase diagram are given in igs.,, and. ll hase transitions in these figures are second order. ig. shows constant-temerature cross-sections of the global hase diagram as a function of and. t low temeratures (high ), the central sin-glass () hase searates the corner ferromagnetic (), columnar (), antiferromagnetic (), and layered () hases. The diagrams are twofold symmetric along each axis, but not fourfold symmetric, due to the difference between transverse ( = ) and longitudinal ( = ) sin glasses. s temerature in-

5 5 creases, the aramagnetic () hase aears at the central oint, first reaches the transverse sin-glass system and eliminates the sin-glass hase, then reaches the longitudinal sin-glass system and eliminates the sin-glass hase. ig. shows temerature-concentration hase diagrams for isotroically mixed, transverse, longitudinal, and =.5 sin-glass systems. The uer left and right hase diagrams are seen to be, resectively, reentrant and forward, namely with a ferromagnetic hase that, resectively, recedes from or roceeds towards the sin-glass hase as temerature is lowered, as clearly seen in the insets. The Nishimori symmetry (see below) is obeyed only at four isolated ordinary oints in each cross-section in ig. and is not obeyed at any oint in the hase diagrams in igs.,, so that the forward behavior is not excluded by the rigorous results [, ]. remarkably narrow sin-glass hase, reaching zerotemerature, occurs in the longitudinal sin-glass system. Zero-temerature hase diagrams are shown in ig. for the longitudinal (left column) and transverse (right column) sin-glass systems. With the aroriate reversal in variables, the longitudinal and transverse singlass hase diagrams are seen in this figure to be qualitatively similar, but quantitatively different. The singlass hase is more extensive in the transverse case. This can be understood from the more extensive intermixing of the ferromagnetic and antiferromagnetic bonds. B. Temerature-oncentration hase Diagrams with Nishimori Symmetry urved ines The Nishimori symmetry condition [, ] for isotroic systems, = e ±K, (5) generalizes, for uniaxially anisotroic sin-glass systems, to = e ±Kxy and = e ±Kz. () or Nishimori symmetry to obtain, both equations have to be satisfied, but the signs in the exonents can be chosen indeendently. The Nishimori condition, in its general form u ( K u ) u (K u ) = e ±Ku (7) for each histogram air of each distribution, is invariant (closed) under our renormalization-grou transformation. If one of the two conditions in Eq.() is fixed, hase diagram cross-sections are obtained, in which Nishimori / / / xy IG. 5: (olor on-line) hase diagrams with Nishimori symmetry lines (dashed) for different anisotroy arameters: The ratio / is, and.5 from to to bottom. In the left column, satisfies the Nishimori condition. In the right column, satisfies the Nishimori condition. ll hase transitions in this figure are second order. symmetry holds along a line. Thus, throughout the three hase diagrams on the left in ig.5, the condition on (, ) is fixed. The condition on (, ), and therefore Nishimori symmetry, is satisfied along the dashed lines on the left in ig.5. In these temerature versus concentration hase diagrams, it is seen that the multicritical oints between the ferromagnetic, sin-glass, and aramagnetic hases lie on the Nishimori symmetry line. urthermore, it has been roven [, ] that a forward hase diagram cannot occur below such a multicritical oint that is on the symmetry line. On the left in ig.5, this is indeed the case, with reentrant hase diagrams, as also seen in isotroic sin glasses [, ]. Recall that in Sec.IV, multicritical oints, between the same hases as here, that do not lie on Nishimori symmetry occur with both reentrant and forward hase diagrams. However, the latter non-symmetric multicritical oints flow, under renormalization-grou transformations, to the (doubly unstable) fixed distribution of the symmetric multicritical oints, therefore being in the same universality class and having the same critical exonents. In the three hase diagrams on the right of ig.5, the condition on (, ) is fixed. In these concentrationconcentration hase diagrams, the multicritical oints between the ordered (ferromagnetic, antiferromagnetic, layered, or columnar), sin-glass, and aramagnetic hases again lie on the Nishimori symmetry lines

6 ... =.7 =.7 Nishimori Surface....9 K/<K> (a) K (b) IG. : (olor on-line) The Nishimori condition for is held throughout the leftmost figure and for throughout the center figure. The comlimentary Nishimori condition, for and resectively, is held along the dashed straight lines, which intersect the ordered (,,, or )-singlassaramagnetic multicritical oints. In the rightmost figure both conditions are satisfied throughout the figure. In this figure, the hase boundaries around the aramagnetic hase are actually lines of the multicritical oints where the aramagnetic, ordered (,,, or ), and sin-glass (not seen in this cross-section) hases meet. In the side figures, firstorder boundaries (dotted) occur between the ferromagnetic and layered hases, and between the antiferromagnetic and columnar hases, terminating at d = critical oints. ll other hase transitions (full lines) in this figure are second order.. oncentration-oncentration hase Diagrams with Nishimori Symmetry Straight ines In the hase diagrams in ig.5, the ratio / is held constant. On the left and center of ig., again the condition in Eq.() on one interaction is fixed and the other interaction strength is held constant. Thus, the Nishimori symmetry lines becomes straight lines. The multicritical oints between the ordered (ferromagnetic, antiferromagnetic, layered, or columnar), sin-glass, and aramagnetic hases again lie on the Nishimori symmetry lines. In the left hase diagram, due to the enforced Nishimori symmetry condition, = along the line =.5 and the system reduces to d =. long this line, first-order transitions between ferromagnetic and layered hases and between antiferromagnetic and columnar hases terminate at d = critical oints. rom.5, d = second-order boundaries between each ordered hase and the aramagnetic hase terminate on the d = critical oints. In the center hase diagram, due to the enforced Nishimori symmetry condition, = along the line =.5 and the system reduces to d =. ccordingly, the system is disordered (aramagnetic) along the entire length of this line. D. The hase Diagram Entirely Imbedded in Nishimori Symmetry In the rightmost ig., both conditions of Eq.() are satisfied throughout the figure. With two symmetry constraints, this is a unique surface in the global hase diagram of our model. The system reduces to d = and... (d) K/< K > (c) K... IG. 7: (olor on-line) ixed distributions, with circles and crosses showing one renormalization-grou transformation and thereby by their exact suerosition attesting to the fixed nature of the distributions. The distributions have been binned for exhibition uroses. (a) or the ferromagneticsinglass hase boundary, a runaway to infinite couling; (b) for the aramagnetic-singlass hase boundary. Both of these fixed distributions are satially isotroic, attracting isotroic and anisotroic boundaries, and do not obey Nishimori symmetry. (c) or the ferromagnetic-singlassaramagnetic multicritical oint. This fixed distribution is satially isotroic and obeys Nishimori symmetry. This fixed distribution attracts the isotroic multicritical oint, which obeys Nishimori symmetry, and anisotroic multicritical oints, which obey and do not obey Nishimori symmetry. The fixed distributions for the antiferromagnetic-singlass, columnar-singlass, layered-singlass hase boundaries and for the antiferromagnetic-singlass-aramagnetic, columnarsinglass-aramagnetic, layered-singlass-aramagnetic multicritical oints are as shown here in (a) and (c) resectively, but with the aroriate and/or reflections. (d) ixed distribution for the sin-glass hase. This hase sink is an isotroic runaway, attracting both satially isotroic and anisotroic sin-glass hase oints, and does not obey Nishimori symmetry. d =, as exlained above, for =.5 and =.5 resectively. The hase boundaries around the aramagnetic hases are actually lines of the multicritical oints where the aramagnetic, ordered (ferromagnetic, layered, antiferromagnetic, or columnar), and sin-glass (not seen in this cross-section) hases meet. No sin-glass hase occurs within the Nishimorisymmetric subsace. The hase transitions seen in the rightmost ig., namely ordered-singlassaramagnetic multicritical and ferromagnetic-layered, antiferromagnetic-layered first-order transitions, are the only hase transitions of the system that occur under Nishimori symmetry.

7 7 E. ixed Distributions The fixed distributions underinning the hase diagrams of this system are given in ig.7. The fixed distributions for the ferromagnetic-singlass boundary, aramagnetic-singlass boundary, and the ferromagnetic-singlass-aramagnetic multicritical oints are satially isotroic, but attract both satially isotroic and anisotroic hase transitions. The fixed distribution for the ferromagnetic-singlassaramagnetic multicritical oints obeys Nishimori symmetry, but attracts multicritical oints that obey and do not obey Nishimori symmetry. In the latter cases, as seen above, both reentrant and forward hase diagrams occur. The fixed distributions for the antiferromagneticsinglass, columnar-singlass, layered-singlass hase boundaries and for the antiferromagnetic-singlassaramagnetic, columnar-singlass-aramagnetic, layered-singlass-aramagnetic multicritical oints are as shown in ig.7 (a) and (c) resectively, but with the aroriate and/or reflections. V. ONUSION The exact solution of the satially uniaxially anisotroic sin glass on a d = hierarchical lattice yields new hase diagrams. In view of the semiquantitative agreement between satially isotroic sin-glass results on cubic and hierarchical lattices [], it would certainly be worthwhile to investigate on cubic lattices the new henomena found in the resent study. urthermore, the exact study of sin glasses on fully anisotroic d = hierarchical lattices [] may yield even more new hase transition henomena. cknowledgments This research was suorted by the Scientific and Technological Research ouncil (TÜBİTK) and by the cademy of Sciences of Turkey. [] H. Nishimori, Statistical hysics of Sin Glasses and Information rocessing (Oxford University ress, ). [] H. Nishimori, J. hys., 7 (9). [] H. Nishimori, rog. Theor. hys., 9 (9). [] H. Nishimori, J. hys. Soc. Jaan 55, 5 (9). [5] Y. Ozeki and H. Nishimori, J. hys., 99 (99). [] H. Nishimori, hysica 5, (99). [7] H. Nishimori and K. Nemoto, J. hys. Soc. Jaan 7, 9 (). [] J.-M. Maillard, K. Nemoto, and H. Nishimori, J. hys. : Math. Gen., 9799 (). [9] K. Takeda and H. Nishimori, Nucl. hys. B, 77 (). [] K. Takeda, T. Sasamoto, and H. Nishimori, J. hys. : Math. Gen., 75 (5). [] M. Hinczewski and.n. Berker, hys. Rev. B 7, (5). [] H. Nishimori, J. Stat. hys., 977 (7). [] G. Migliorini and.n. Berker, hys. Rev. B 57, (99). [].D. Nobre, hys. Rev. E, (). [5] M. Ohzeki, J. hys. Soc. Jaan 7, (7). [] M. Hinczewski and.n. Berker, Eur. hys. J. B 5, (). [7] M. Hinczewski and.n. Berker, arxiv:cond-mat/77v [cond-mat.str-el]. [].N. Berker and S. Ostlund, J. hys., 9 (979). [9] R.B. Griffiths and M. Kaufman, hys. Rev. B, 5R (9). [] M. Kaufman and R.B. Griffiths, hys. Rev. B, (9). [] S.R. McKay,.N. Berker, and S. Kirkatrick, hys. Rev. ett., 77 (9). []. Monthus and T. Garel, arxiv:7.7v [condmat.dis-nn]. []. alicov,.n. Berker, and S.R. McKay, hys. Rev. B 5, (995). [] E. Domany, S. lexander, D. Bensimon, and.. Kadanoff, hys. Rev. B, (9). [5] J.-M. anglois,.-m.s. Tremblay, and B.W. Southern, hys. Rev. B, (9). [] R.B. Stinchcombe and.. Maggs, J. hys. 9, 99 (9). [7] R.. ngulo and E. Medina, J. Stat. hys. 75, 5 (99). [] T..S. Haddad, S.T.R. inho, and S.R. Salinas, hys. Rev. E, (). [9] J.-X. e and Z.R. Yang, hys. Rev. E 9, 7 (). [] R.. da Silveira and J.-. Bouchaud, hys. Rev. ett. 9, 59 (). [].-H. Tang and H. haté, arxiv:cond-mat/75v [cond-mat.stat-mech]. []. Monthus and T. Garel, J. Stat. Mech.: Theory and Exeriment, (). [] M. Hinczewski and.n. Berker, hys. Rev. E 7, (). [] M. Hinczewski, hys. Rev E 75, (7). [5] Z.Z. Zhang,.. Rong, and S.G. Zhou, hysica 77, 9 (7). [] Z.Z. Zhang, S.G. Zhou, and T. Zou, Eur. hys. J. B 5, 59 (7). [7] Z.Z. Zhang, Z.G. Zhou, and.. hen, Eur. hys. J. B 5, 7 (7). [] H.D. Rozenfeld, S. Havlin, and D. ben-vraham, New J. hys. 9, 75 (7). [9] H.D. Rozenfeld and D. ben-vraham, hys. Rev. E 75, (7). [] E. Khajeh, S.N. Dorogovtsev, and J... Mendes, hys.

8 Rev. E 75, (7). []. Erbaş,. Tuncer, B. Yücesoy, and.n. Berker, hys. Rev. E 7, 9 (5).