arxiv: v1 [cond-mat.stat-mech] 21 Mar 2018

Transcription

1 APS Elastic backbone defines a new transition in the ercolation model Cesar I. N. Samaio Filho 1, José S. Andrade Jr. 1,2, Hans J. Herrmann 1,2, André A. Moreira 1 1 Deartamento de Física, Universidade Federal do Ceará, Fortaleza, Ceará, Brasil 2 Comutational Physics for Engineering Materials, IfB, ETH Zurich, Schafmattstrasse, 8093 Zurich, Switzerland arxiv: v1 [cond-mat.stat-mech] 21 Mar 2018 The elastic backbone is the set of all shortest aths. We found a new hase transition at eb above the classical ercolation threshold at which the elastic backbone becomes dense. At this transition in 2d its fractal dimension is 1.750±0.003, and one obtains a novel set of critical exonents β eb = 0.50 ± 0.02, = 1.97 ± 0.05, and ν eb = 2.00 ± 0.02 fulfilling consistent critical scaling laws. Interestingly, however, the hyerscaling relation is violated. Using Binder s cumulant, we determine, with high recision, the critical robabilities eb for the triangular and tilted square lattice for site and bond ercolation. This transition describes a sudden rigidification as a function of density when stretching a damaged tissue. PACS numbers:.0.ah,.0.al, q, Da Desite being one of the simlest and most studied models, classical ercolation [1 ] still bears yet uncovered surrises. It is well known that at the ercolation threshold ( c ) the shortest ath between two oosite sides of the system is fractal, with a fractal dimension that is only known numerically to be d s = ±0.000 [5 8] in two dimension and increases with dimension, until becoming two at and above the critical dimension d =. In fact the shortest ath is not unique: several shortest aths can exist simultaneously and the set of all shortest aths has been called elastic backbone in the ast [9], because it is the subset of the backbone that, when elongated, would give the first contribution to a restoring force. The elastic backbone indeed determines the first resistance that is felt, when stretching damaged [10 12] or biological tissues [13 18] and thus has exerimental relevance. It has been numerically established that at the ercolation threshold c its fractal dimension is indistinguishable from the one of the shortest ath [9]. Here we reort on the discovery that, above the classical ercolation threshold c, there exists another critical robability eb > c at which the elastic backbone becomes dense. While its dimension is known to be d s at c, our results in two dimensions show that it becomes unity between c and eb, it is 1.750±0.003 at eb, and is equal to two above. At eb we reveal critical scaling laws and a set of exonents, however, violating hyerscaling. We simulated two-dimensional ercolation configurations at occuation robability for systems of size with eriodic boundary conditions in horizontal direction and oen boundaries at to and bottom. Using the burning algorithm, we identifieor each configuration the elastic backbone, i.e. the set of all shortest aths [9]. We considered site and bond ercolation on the triangular lattice and on two tyes of square lattices, namely, tilted and non-tilted. In Fig. 1 we see the elastic back- Corresondence to: cesar@fisica.ufc.br (a) (b) (c) FIG. 1: Shown in red are tyical elastic backbones obtained for site ercolation on a tilted square lattice of linear size = 512 and calculated with: (a) = ( c < < eb ), (b) = ( = eb ), and (c) = ( > eb ). Occuied sites that are in the sanning cluster, but do not belong to the elastic backbone, are shown in blue. Other sites are reresented in white. bones for site ercolation on a tilted square lattice for three different robabilities below, at and above eb. We define the quantity M eb as the mass of the elastic backbone, i.e., the number of sites that belong to the elastic backbone, and its density as m () = M eb /N, where N is the total number of sites of the lattice. In Fig. 2 we lot m () against > c for site ercolation

2 2 m () = 25 = 512 = 102 = 208 = 09 m () β eb /ν eb FIG. 2: Density m () of the elastic backbone as a function of the occuation robability for site ercolation on the triangular lattice for different system sizes. Inset: determination of the threshold using finite-size scaling. on the triangular lattice anind that there exists a clear hase transition at a value eb > c at which the elastic backbone becomes dense, with m acting as order arameter of this transition. Using finite size scaling, eb can be determined more recisely as shown in the inset, yielding also an estimate for the exonent β eb /ν eb = 0.25 ± The threshold eb can be determined even more recisely using Binder s cumulant defined as m eb U () = (1) 3 m 2 eb 2 In Fig. 3 we show the analysis of Binder s cumulant for site ercolation on the triangular lattice, where we obtain eb = ± Moreover, by alying the same analysis, we finor site ercolation eb = 1 on the normal square lattice and eb = ± on the tilted square lattice. For bond ercolation, we obtain eb = ± on the tilted square lattice and eb = 0.50 ± on the triangular lattice. At = 1, for the square lattice, starting from a given node in the border, there is just one shortest ath of size reaching the other side. In this way, the square lattice leaves the dense hase for any fraction of nodes removed. For the triangular and tilted lattices, on the other hand, there are 2 different shortest aths (of size ) reaching the other side. This much larger number of shortest aths increases the mass of the elastic backbone, and therefore the triangular and tilted lattices remain dense, even after some sites have been removed. It should be noted that a similar deendence on the lattice is observed in other models. In the case of Rigidity Percolation [19 22], a transition on the triangular lattice is observed at finite, while rigidity is only attained at = 1 for the square lattice [23 25]. The fractal dimension of the elastic backbone is ob- ln(1-1.5u (q)) = 25 = 512 = 102 = 208 = 09 eb () FIG. 3: Binder s cumulant U () for different sizes as a function of the occuation robability for site ercolation on the triangular lattice. Inset: The deendence of eb () on 1 for site ercolation on the triangular lattice. eb () (circles) is obtainerom the crossing oint of U and U 2 for each air of successive values. Here we consider =, 128, 25, 512, 102, 208, 09, and Extraolating through the data oints (blue line) to the thermodynamic limit, we obtain eb ( ) = ± (red dashed line). tained directly from the double-logarithmic lot of the mass of the elastic backbone M eb against the linear size of the considered lattice, at the threshold eb. We find that the elastic backbone is fractal with a fractal dimension 0 ± for all studied models, while for c < < eb it is one-dimensional, as shown in Fig.. In this case, although the sanning cluster is dense, the existence of various holes revents the effective coalescence of the shortest aths resent in the system, therefore leading to an elastic backbone which is a fractal. At the threshold eb, the exonent β/ν is then obtainerom the relation β/ν = d, resulting, for all studied models, in β/ν = ± Moreover, the resonse function of the order arameter, i.e. what would corresond in magnetic systems to the suscetibility, χ = N( meb 2 meb 2 ), diverges at the critical threshold eb with an exonent /ν eb = 0 ± 0.02 for all models considered, as shown in Fig. 5 for site ercolation on the tilted square and triangular lattices. Finally, we also erform a full finite-size scaling analysis for m () anor χ () of the form, m () = β eb/ν eb m(ε 1/ν eb ), (2) χ () = /ν eb χ(ε 1/ν eb ), (3) where ε = ( eb ) is the distance from the critical threshold. The exonents β eb /ν eb, /ν eb, and ν eb are, resectively, associated with the decay of the order arameter,

3 3 M eb d = d = χ ( c ) /ν eb = /ν eb = FIG. : ogarithmic lot of the mass M eb of the elastic backbone as a function of the lattice size for site ercolation at eb on the tilted square lattice (orange circles) and on the triangular lattice (blue squares). The same is shown for the value of = 0.3, between c and eb, on the tilted square lattice for site ercolation (red stars) anor bond ercolation (green triangles), with ranging from to 138 sites. At = eb, the least-squares fit to the data of a ower law, M eb, gives the exonent 00 ± for the tilted square lattice (black line) and 00 ± for the triangular lattice (blue line). At = 0.3 on the tilted square lattice, the least-squares fit to the data of a ower law, M eb d, gives the exonent d = 000 ± for site ercolation (red line) and d = 000 ± for bond ercolation (green line). In all cases, the errors are smaller than the symbols. the divergence of the suscetibility, and the finite-size effects. As shown in Fig., for the secific case of site ercolation on the triangular lattice, we obtain excellent data collase for values β eb = 0.50 ± 0.03, = 1.97 ± 0.05 and ν eb = 2.00±0.0. We note that hyerscaling relation 2β eb + = dν eb is violated [2, 27]. Similar data collase with the same exonents have been founor the other considered lattices. We also studied the elastic backbone transitions on a lattice model that mixes the features of normal square and triangular lattices, by adding to the normal square lattice with robability q some additional diagonals going from to-left to bottom-right. In this way, for q = 0, we obtain a normal square lattice and, for q = 1, a triangular lattice. For every value of q > 0, we found a value of c (q) < 1 at which the elastic backbone becomes dense. In Fig. 7a we show the finite-size scaling of the mass M eb of the elastic backbone for site ercolation at eb (q), for q = 0.0 and q = The results show that the fractal dimension is within error bars the same for both cases. Moreover, we calculated these fractal dimensions at c (q) for different values of q > 0 anound that they do not deend on q. This underlines, on one FIG. 5: ogarithmic lot of the suscetibility χ of the elastic backbone as a function of the lattice size for site ercolation at eb, with ranging from to 138 sites. The leastsquares fit to the data of a ower law, χ( eb ) γ/ν, gives the exonent γ/ν = 0 ± 0.01 for the tilted square lattice (orange circles) and γ/ν = 0±0.02 for the triangular lattice (blue squares). hand, the universality of d eb and questions, on the other hand, any relation to Rigidity Percolation [19 22], which does exhibit a transition at an intermediate value of q. In Fig. 7b, we show the hase diagram between ercolation robability and the density of diagonals q, where three hases are identified. In the non-ercolating hase, the sanning cluster is absent and consequently the elastic backbone does not ercolate. The non-ercolating hase is bounded by the curve c (q) which reresents the classical ercolation thresholds. The non-dense hase is characterized by a linear scaling with of the elastic backbone and it is bounded by the curves c (q) on the bottom and eb (q) on the to, which defines the critical line along which the order arameter m () vanishes and the elastic backbone is fractal. In the dense hase, the dimension of the elastic backbone is equal to the dimension of the lattice considered. Concluding, we discovered that the mass of the elastic backbone serves as order arameter for a new transition within the connected hase of classical ercolation, exhibiting a new set of critical exonents β eb 1/2, 2, ν eb 2, and 7/ in two dimensions. Interestingly, however, hyerscaling is violated, being this to our knowledge the first examle for a violation of hyerscaling in a urely geometrical model. It would be also interesting to investigate higher dimensions and try to formulate a mean-field aroximation. Similar transitions for elastic backbones could be exected in models with tunable disorder [28 3]. Our findings have direct consequences to the stretching of random fibrous materials like biological tissues: when the first restoring force is felt, the resistance will grow

4 m () β eb /ν eb a) = 25 = 512 = 102 = 208 = 09 β eb = ε 1/ν eb M () q = 0.0 q = 0.80 (a) b) 0.8 Dense Phase χ () - /ν eb χ max ε 1/ν eb = 1.97 FIG. : (a) Results from the finite-size scaling analysis of the order arameter m obtaineor site ercolation on the triangular lattice. The growing and decaying curves corresond to values above and below eb, resectively. The dashed line is the least-squares fit to data in the scaling region above eb for Eq. 2. (b) The same as in (a), but for the suscetibility χ. Blue symbols for < eb and red symbols for > eb. Here, the dashed line has sloe Inset of (b): og-log lot of the maximum of χ as a function of system size. very gently with dislacement below the threshold eb, while above eb the system will then to be instead very stiff. It is thus a transition between two very different stress-strain relations for a damaged tissue. Furthermore, there exists an interesting similarity between the elastic backbone ercolation, as introduced here, and the hase transition associated to Rigidity Percolation, since the stressed backbone [37 39] is a fractal whose dimension ( = 1.78 ± 0.02) is very close to the fractal dimension of the elastic backbone at eb. For non-tilted square lattices, in both cases, the transition is shifted to = 1 [23 25]. Nevertheless the two transitions describe different henomena, since they are in general Non-dense Phase Non-ercolating Phase (b) q FIG. 7: (a) ogarithmic lot of the mass M eb of the elastic backbone as a function of the lattice size for site ercolation at eb (q) for two different values of q. The least-squares fit to the data of a ower law, M eb ( eb (q)), gives the exonent 0 ± for q = 0.0 (orange circles) and 0 ± for q = 0.80 (red squares). The results show that the fractal dimension is within error bars the same for both cases. (b) Phase diagram between ercolation robability and the density of diagonals q. Three hases can be identified. The non-ercolating hase, where the sanning cluster is absent, bounded on to by the curve c(q) (filled circles) which reresents the classical ercolation thresholds. The non-dense hase is bounded by the curve c(q) on the bottom and the curve eb (q) on the to (stars) that is the critical line along which the order arameter m () vanishes and the elastic backbone is fractal. Finally, in the dense hase on the to the elastic backbone becomes dense. located at different thresholds values and have a different set of critical exonents. Acknowledgments We thank the Brazilian agencies CNPq, CAPES, FUN- CAP, the National Institute of Science and Technology

5 5 for Comlex Systems, and the Euroean Research Council (ERC) Advanced Grant FlowCCS for financial suort. [1] S. R. Broadbent and J. M. Hammersley, in Math. Proc. Cambridge (Cambridge University Press, 1957), vol. 53, [2] S. Kirkatrick, Rev. Mod. Phys. 5, 57 (1973). [3] D. Stauffer and A. Aharony, Tailor & Francis, ondon (1985). [] M. Sahimi, Alications of Percolation Theory (Taylor & Francis Grou, ondon, UK, 199). [5] H. J. Herrmann and H. E. Stanley, J. Phys. A-Math. Gen. 21, 829 (1988). [] P. Grassberger, J. Phys. A-Math. Gen. 25, 575 (1992). [7] N. V. Dokholyan, Y. ee, S. V. Buldyrev, S. Havlin, P. R. King, and H. E. Stanley, J. Stat. Phys. 93, 03 (1998). [8] M. E. J. Newman and R. M. Ziff, Phys. Rev. ett. 85, 10 (2000). [9] H. J. Herrmann, D. C. Hong, and H. E. Stanley, J. Phys. A-Math. Gen. 17, 21 (198). [10] D. J. Jacobs and M. F. Thore, Phys. Rev. ett. 75, 051 (1995). [11] V. K. de Souza and P. Harrowell, Proc. Natl. Acad. Sci. USA 10, 1513 (2009). [12]. Zhang, D. Z. Rocklin,. M. Sander, and X. Mao, Phys. Rev. Materials 1, (2017). [13] J. Bates, G. Maksym, D. Navajas, and B. Suki, Ann. Biomed. Eng. 22, 7 (199). [1] H. Yuan, E. P. Ingenito, and B. Suki, J. Al. Physiol. 83, 120 (1997). [15] J. W. Baish and R. K. Jain, Cancer Res. 0, 383 (2000). [1] B. Suki, Am. J. Resir. Crit. Care Med. 1, 133 (2002). [17] M. C. Ritter, R. Jesudason, A. Majumdar, D. Stamenović, J. A. Buczek-Thomas, P. J. Stone, M. A. Nugent, and B. Suki, Proc. Natl. Acad. Sci. USA 10, 1081 (2009). [18] F. E. ennon, G. C. Cianci, N. A. Ciriani, T. A. Hensing, H. J. Zhang, C.-T. Chen, S. D. Murgu, E. E. Vokes, M. W. Vannier, and R. Salgia, Nat. Rev. Clin. Oncol. 12, (2015). [19] Y. Kantor and I. Webman, Phys. Rev. ett. 52, 1891 (198). [20] W. Bresser, P. Boolchand, and P. Suranyi, Phys. Rev. ett. 5, 293 (198). [21] C. Moukarzel, P. M. Duxbury, and P. eath, Phys. Rev. ett. 78, 180 (1997). [22] W. G. Ellenbroek, V. F. Hagh, A. Kumar, M. F. Thore, and M. Van Hecke, Phys. Rev. ett. 11, (2015). [23] C. F. Moukarzel, Phys. Rev. ett. 81, 13 (1998). [2] A. Souslov, A. J. iu, and T. C. ubensky, Phys. Rev. ett. 103, (2009). [25] X. Mao, N. Xu, and T. C. ubensky, Phys. Rev. ett. 10, (2010). [2] K. Binder and J.-S. Wang, J. Stat. Phys. 55, 87 (1989). [27] R.. C. Vink, T. Fischer, and K. Binder, Phys. Rev. E 82, (2010). [28] M. Porto, S. Havlin, S. Schwarzer, and A. Bunde, Phys. Rev. ett. 79, 00 (1997). [29]. A. Braunstein, S. V. Buldyrev, R. Cohen, S. Havlin, and H. E. Stanley, Phys. Rev. ett. 91, (2003). [30] J. S. Andrade, E. A. Oliveira, A. A. Moreira, and H. J. Herrmann, Phys. Rev. ett. 103, (2009). [31] S. Pradhan, A. Hansen, and B. K. Chakrabarti, Rev. Mod. Phys. 82, 99 (2010). [32] E. Fehr, D. Kadau, J. S. Andrade, and H. J. Herrmann, Phys. Rev. ett. 10, (2011). [33] A. A. Moreira, C.. N. Oliveira, A. Hansen, N. A. M. Araújo, H. J. Herrmann, and J. S. Andrade Jr, Phys. Rev. ett. 109, (2012). [3] A. Shekhawat, S. Zaeri, and J. P. Sethna, Phys. Rev. ett. 110, (2013). [35] C.. N. Oliveira, P. A. Morais, A. A. Moreira, and J. S. Andrade, Phys. Rev. ett. 112, (201). [3] C. I. N. Samaio Filho, A. A. Moreira, N. A. M. Araújo, J. S. Andrade Jr, and H. J. Herrmann, Phys. Rev. ett. 117, (201). [37] C. Moukarzel and P. M. Duxbury, Phys. Rev. ett. 75, 055 (1995). [38] D. J. Jacobs and M. F. Thore, Phys. Rev. ett. 80, 551 (1998). [39] C. Moukarzel and P. M. Duxbury, Phys. Rev. E 59, 21 (1999).