Site percolation on central-force elastic networks

Elastic moduli of body-centered cubic lattice near rigidity percolation threshold: Finite-size effects and evidence for first-order phase transition

Extensive numerical simulations of rigidity percolation with only central forces in large three-dimensional lattices have indicated that many of their topological properties undergo a first-order phase transition at the rigidity percolation threshold p_{ce}. In contrast with such properties, past numerical calculations of the elastic moduli of the same lattices had provided evidence for a second-order phase transition. In this paper we present the results of extensive simulation of rigidity percolation in large body-centered cubic (bcc) lattices, and show that as the linear size L of the lattice increases, the elastic moduli close to p_{ce} decrease in a stepwise, discontinuous manner, a feature that is absent in lattices with L<30. The number and size of such steps increase with L. As p_{ce} is approached, long-range, nondecaying orientational correlations are built up, giving rise to compact, nonfractal clusters. As a result, we find that the backbone of the lattice at p_{ce} is compact with a fractal dimension D_{bb}≈3. The absence of fractal, scale-invariant clusters, the hallmark of second-order phase transitions, together with the stairwise behavior of the elastic moduli, provide strong evidence that, at least in bcc lattices, many of the topological properties of rigidity percolation as well as its elastic moduli may undergo a first-order phase transition at p_{ce}. In relatively small lattices, however, the boundary effects interfere with the nonlocal nature of the rigidity percolation. As a result, only when such effects diminish in large lattices does the true nature of the phase transition emerge.

Relationship between vibrations and dynamical heterogeneity in a model glass former: extended soft modes but local relaxation

We study the relation between short-time vibrational modes and long-time relaxational dynamics in a kinetically constrained lattice gas with harmonic interactions between neighbouring particles. We find a correlation between the location of the low- (high-) frequency vibrational modes and regions of high (low) propensity for motion. This is similar to what was observed in continuous force systems, but our interpretation is different: in our case relaxation is due to localised excitations which propagate through the system; these localised excitations act as background disorder for the elastic network, giving rise to anomalous vibrational modes. Our results provide an example whereby a correlation between spatially extended low-frequency modes and high-propensity regions does not imply that relaxational dynamics originates in extended soft modes but rather belies their common origin. We consider other measures of elastic heterogeneity, such as non-affine displacement fields and mode localisation lengths, and discuss implications of our results to interpretations of dynamic heterogeneity more generally.

Self-healing materials: a review

The ability of materials to self-heal from mechanical and thermally induced damage is explored in this paper and has significance in the field of fracture and fatigue. The history and evolution of several self-repair systems is examined including nano-beam healing elements, passive self-healing, autonomic self-healing and ballistic self-repair. Self-healing mechanisms utilized in the design of these unusual materials draw much information from the related field of polymer-polymer interfaces and crack healing. The relationship of material damage to material healing is examined in a manner to provide an understanding of the kinetics and damage reversal processes necessary to impart self-healing characteristics. In self-healing systems, there are transitions from hard-to-soft matter in ballistic impact and solvent bonding and conversely, soft-to-hard matter transitions in high rate yielding materials and shear-thickening fluids. These transitions are examined in terms of a new theory of the glass transition and yielding, viz., the twinkling fractal theory of the hard-to-soft matter transition. Success in the design of self-healing materials has important consequences for material safety, product performance and enhanced fatigue lifetime.

Algorithms for three-dimensional rigidity analysis and a first-order percolation transition

A fast computer algorithm, the pebble game, has been used successfully to analyze the rigidity of two-dimensional (2D) elastic networks, as well as of a special class of 3D networks, the bond-bending networks, and enabled significant progress in studies of rigidity percolation on such networks. Application of the pebble game approach to general 3D networks has been hindered by the fact that the underlying mathematical theory is, strictly speaking, invalid in this case. We construct an approximate pebble game algorithm for general 3D networks, as well as a slower but exact algorithm, the relaxation algorithm, that we use for testing the new pebble game. Based on the results of these tests and additional considerations, we argue that in the particular case of randomly diluted central-force networks on bcc and fcc lattices, the pebble game is essentially exact. Using the pebble game, we observe an extremely sharp jump in the largest rigid cluster size in bond-diluted central-force networks in 3D, with the percolating cluster appearing and taking up most of the network after a single bond addition. This strongly suggests a first-order rigidity percolation transition, which is in contrast to the second-order transitions found previously for the 2D central-force and 3D bond-bending networks. While a first order rigidity transition has been observed previously for Bethe lattices and networks with "chemical order," here it is in a regular randomly diluted network. In the case of site dilution, the transition is also first order for bcc lattices, but results for fcc lattices suggest a second-order transition. Even in bond-diluted lattices, while the transition appears massively first order in the order parameter (the percolating cluster size), it is continuous in the elastic moduli. This, and the apparent nonuniversality, make this phase transition highly unusual.

Rigidity percolation in a field

Rigidity percolation with g degrees of freedom per site is analyzed on randomly diluted Erdös-Renyi graphs, with average connectivity gamma, in the presence of a field h. In the (gamma,h) plane, the rigid and flexible phases are separated by a line of first-order transitions whose location is determined exactly. This line ends at a critical point with classical critical exponents. Analytic expressions are given for the densities n(F) of uncanceled degrees of freedom and gamma(r) of redundant bonds. Upon crossing the coexistence line, gamma(r) and n(F) are continuous, although their first derivatives are discontinuous. We extend, for the case of nonzero field, a recently proposed hypothesis, namely, that the density of uncanceled degrees of freedom is a "free energy" for rigidity percolation. Analytic expressions are obtained for the energy, entropy, and specific heat. Some analogies with a liquid-vapor transition are discussed. Particularizing to zero field, we find that the existence of a (g+1) core is a necessary condition for rigidity percolation with g degrees of freedom. At the transition point gamma(c), Maxwell counting of degrees of freedom is exact on the rigid cluster and on the (g+1) rigid core, i.e., the average coordination of these subgraphs is exactly 2g, although gamma(c), the average coordination of the whole system, is smaller than 2g. gamma(c) is found to converge to 2g for large g, i.e., in this limit Maxwell counting is exact globally as well.

Elasticity of randomly diluted central force networks under tension

We study the rigidity of two-dimensional site-diluted central force triangular networks under tension. We calculate the shear modulus micro directly and fit it with a power law of the form mu approximately (p-p*)(f), where p is the concentration of sites, p* its critical value, and f the critical exponent. We find that the critical behavior of mu is quite sensitive to tension. As the tension is increased there is at first a sharp drop in the values of both p* and f, followed by a slower decrease towards the values of the diluted Gaussian spring network (or random resistor network). We find that the size of the critical region is also sensitive to tension. The tension-free system has a narrower critical regime with the power law failing for p>0.8. In contrast, a small tension is sufficient to extend the power law to near p=1. The physical basis for these behaviors is discussed.

Directed rigidity and bootstrap percolation in 1+1 dimensions

We study directed rigidity percolation (equivalent to directed bootstrap percolation) on three different lattices: square, triangular, and augmented triangular. The first two of these display a first-order transition at p=1, while the augmented triangular lattice shows a continuous transition at a nontrivial p(c). On the augmented triangular lattice we find, by extensive numerical simulation, that the the directed rigidity percolation transition belongs to the same universality class as the directed percolation. The same conclusion is reached by studying its surface critical behavior, i.e., the spreading of rigidity from finite clusters close to a nonrigid wall. Near the discontinuous transition at p=1 on the triangular lattice, we are able to calculate the finite-size behavior of the density of rigid sites analytically. Our results are confirmed by numerical simulation.