Generic rigidity percolation in two dimensions

Rigidity transitions in anisotropic networks: a crossover scaling analysis

We study how the rigidity transition in a triangular lattice changes as a function of anisotropy by preferentially filling bonds on the lattice in one direction. We discover that the onset of rigidity in anisotropic spring networks on a regular triangular lattice arises in at least two steps, reminiscent of the two-step melting transition in two dimensional crystals. In particular, our simulations demonstrate that the percolation of stress-supporting bonds happens at different critical volume fractions along different directions. By examining each independent component of the elasticity tensor, we determine universal exponents and develop universal scaling functions to analyze isotropic rigidity percolation as a multicritical point. Our crossover scaling approach is applicable to anisotropic biological materials (e.g. cellular cytoskeletons, extracellular networks of tissues like tendons), and extensions to this analysis are important for the strain stiffening of these materials.

Short-range correlation of stress chains near solid-to-liquid transition in active monolayers

Using a three-dimensional model of cell monolayers, we study the spatial organization of active stress chains as the monolayer transitions from a solid to a liquid state. The critical exponents that characterize this transition map the isotropic stress percolation onto the two-dimensional random percolation universality class, suggesting short-range stress correlations near this transition. This mapping is achieved via two distinct, independent pathways: (i) cell-cell adhesion and (ii) active traction forces. We unify our findings by linking the nature of this transition to high-stress fluctuations, distinctly linked to each pathway. The results elevate the importance of the transmission of mechanical information in dense active matter and provide a new context for understanding the non-equilibrium statistical physics of phase transition in active systems.

Effects of coordination and stiffness scale separation in disordered elastic networks

Many fibrous materials are modeled as elastic networks featuring a substantial separation between the stiffness scales that characterize different microscopic deformation modes of the network's constituents. This scale separation has been shown to give rise to emergent complexity in these systems' linear and nonlinear mechanical response. Here we study numerically a simple model featuring said stiffness scale separation in two-dimensions and show that its mechanical response is governed by the competition between the characteristic stiffness of collective nonphononic soft modes of the stiff subsystem, and the characteristic stiffness of the soft interactions. We present and rationalize the behavior of the shear modulus of our complex networks across the unjamming transition at which the stiff subsystem alone loses its macroscopic mechanical rigidity. We further establish a relation in the soft-interaction-dominated regime between the shear modulus, the characteristic frequency of nonphononic vibrational modes, and the mesoscopic correlation length that marks the crossover from a disorder-dominated response to local mechanical perturbations in the near field, to a linear, continuumlike response in the far field. The effects of spatial dimension on the observed scaling behavior are discussed, in addition to the interplay between stiffness scales in strain-stiffened networks, which is relevant to understanding the nonlinear mechanics of non-Brownian fibrous biomatter.

Properties of packings and dispersions of superellipse sector particles

Superellipse sector particles (SeSPs) are segments of superelliptical curves that form a tunable set of hard-particle shapes for granular and colloidal systems. SeSPs allow for continuous parametrization of corner sharpness, aspect ratio, and particle curvature; rods, circles, rectangles, and staples are examples of shapes SeSPs can model. We compare three computational processes: pair-wise Monte Carlo simulations that explore particle-particle geometric constraints, Monte Carlo simulations that reveal how these geometric constraints play out over dispersions of many particles, and Molecular Dynamics simulations that form random loose and close packings. We investigate the dependence of critical random loose and close packing fractions on particle parameters, finding that both values increase with opening aperture and decrease with increasing corner sharpness. The identified packing fractions are compared with the mean-field prediction of the random contact model; we find deviations from the model's prediction due to correlations between particle orientations. The complex interaction of spatial proximity and orientational alignment is also explored with a generalized spatioorientational distribution area (SODA) plot, which shows how higher density packings are achieved through particles assuming a small number of preferred configurations that depend sensitively on particle shape and system preparation.

Free volume theory explains the unusual behavior of viscosity in a non-confluent tissue during morphogenesis

A recent experiment on zebrafish blastoderm morphogenesis showed that the viscosity (η) of a non-confluent embryonic tissue grows sharply until a critical cell packing fraction (ϕS). The increase in η up to ϕS is similar to the behavior observed in several glass-forming materials, which suggests that the cell dynamics is sluggish or glass-like. Surprisingly, η is a constant above ϕS. To determine the mechanism of this unusual dependence of η on ϕ, we performed extensive simulations using an agent-based model of a dense non-confluent two-dimensional tissue. We show that polydispersity in the cell size, and the propensity of the cells to deform, results in the saturation of the available free area per cell beyond a critical packing fraction. Saturation in the free space not only explains the viscosity plateau above ϕS but also provides a relationship between equilibrium geometrical packing to the dramatic increase in the relaxation dynamics.

Rigidity percolation in a random tensegrity via analytic graph theory

Functional structures from across the engineered and biological world combine rigid elements such as bones and columns with flexible ones such as cables, fibers, and membranes. These structures are known loosely as tensegrities, since these cable-like elements have the highly nonlinear property of supporting only extensile tension. Marginally rigid systems are of particular interest because the number of structural constraints permits both flexible deformation and the support of external loads. We present a model system in which tensegrity elements are added at random to a regular backbone. This system can be solved analytically via a directed graph theory, revealing a mechanical critical point generalizing that of Maxwell. We show that even the addition of a few cable-like elements fundamentally modifies the nature of this transition point, as well as the later transition to a fully rigid structure. Moreover, the tensegrity network displays a collective avalanche behavior, in which the addition of a single cable leads to the elimination of multiple floppy modes, a phenomenon that becomes dominant at the transition point. These phenomena have implications for systems with nonlinear mechanical constraints, from biopolymer networks to soft robots to jammed packings to origami sheets.

Mixing in Langmuir Monolayers: Perfluorotetradecanoic Acid and a Gemini Surfactant without a Linker

A recently reported anionic gemini surfactant, a member of the so-called "gemini without a linker" family, has recently been reported to form closely packed crystalline monolayers at the air-water interface. In this work, the impact on monolayer properties of the compound, C18-0-C18, that result from its mixing with a benchmark perfluorinated surfactant, perfluorotetradecanoic acid (PF), is explored. The films exhibit nonideal mixing, as determined by surface pressure-area (π-A) isotherms and surface potential measurements, and phase-separation between the two components was observed by the direct visualization of the monolayers, and grazing-incident X-ray diffraction at the air-water interface. The pure and mixed films follow similar trends in the order of C18-0-C18 < PF < χPF = 0.50 mixed films for both their extent of hysteresis and their stability at the air-water interface. Further, crystallographic data for the mixed film emerge as a simple combination of distinct diffraction patterns characteristic of both the individual components, consistent with the other findings reported here and thus clarify the intermolecular behavior of the binary mixture at the surface.

Discontinuous rigidity transition associated with shear jamming in granular simulations

We investigate the rigidity transition associated with shear jamming in frictionless, as well as frictional, disk packings in the quasi-static regime and at low shear rates. For frictionless disks, the transition under quasi-static shear is discontinuous, with an instantaneous emergence of a system spanning rigid clusters at the jamming transition. For frictional systems, the transition appears continuous for finite shear rates, but becomes sharper for lower shear rates. In the quasi-static limit, it is discontinuous as in the frictionless case. Thus, our results show that the rigidity transition associated with shear jamming is discontinuous, as demonstrated in the past for isotropic jamming of frictionless particles, and therefore a unifying feature of the jamming transition in general.

Heterogeneous Dynamics of Sheared Particle-Laden Fluid Interfaces with Janus Particle Doping

The formation of particle clusters can substantially modify the dynamics and mechanical properties of suspensions in both two and three dimensions. While it has been well established that large network-spanning clusters increase the rigidity of particle systems, it is still unclear how the presence of localized nonpercolating clusters affects the dynamics and mechanical properties of particle suspensions. Here, we introduce self-assembled localized particle clusters at a fluid-fluid interface by mixing a fraction of Janus particles in a monolayer of homogeneous colloids. Each Janus particle binds to a few nearby homogeneous colloids, resulting in numerous small clusters uniformly distributed across the interface. Using a custom magnetic rod interfacial stress rheometer, we apply linear oscillatory shear to the particle-laden fluid interface. By analyzing the local affine deformation of particles from optical microscopy, we show that particles in localized clusters experience substantially lower shear-induced stretching than their neighbors outside clusters. We hypothesize that such heterogeneous dynamics induced by particle clusters increase the effective surface coverage of particles, which in turn enhances the shear moduli of the interface, as confirmed by direct interfacial rheological measurements. Our study illustrates the microscopic dynamics of small clusters in a shear flow and reveals their profound effects on the macroscopic rheology of particle-laden fluid interfaces. Our findings open an avenue for designing interfacial materials with improved mechanical properties via the control of formation of localized particle clusters.

Evidences of Conformal Invariance in 2D Rigidity Percolation

The rigidity transition occurs when, as the density of microscopic components is increased, a disordered medium becomes able to transmit and ensure macroscopic mechanical stability, owing to the appearance of a space-spanning rigid connected component, or cluster. As a second-order phase transition it exhibits a scale invariant critical point, at which the rigid clusters are random fractals. We show, using numerical analysis, that these clusters are also conformally invariant, and we use conformal field theory to predict the form of universal finite-size effects. Furthermore, although connectivity and rigidity percolation are usually thought to be of fundamentally different natures, we provide evidence of unexpected similarities between the statistical properties of their random clusters at criticality. Our work opens a new research avenue through the application of the powerful 2D conformal field theory tools to understand the critical behavior of a wide range of physical and biological materials exhibiting such a mechanical transition.

Explosive Percolation Obeys Standard Finite-Size Scaling in an Event-Based Ensemble

Explosive percolation in the Achlioptas process, which has attracted much research attention, is known to exhibit a rich variety of critical phenomena that are anomalous from the perspective of continuous phase transitions. Hereby, we show that, in an event-based ensemble, the critical behaviors in explosive percolation are rather clean and obey the standard finite-size scaling theory, except for the large fluctuation of pseudo-critical points. In the fluctuation window, multiple fractal structures emerge and the values can be derived from a crossover scaling theory. Further, their mixing effects account well for the previously observed anomalous phenomena. Making use of the clean scaling in the event-based ensemble, we determine with a high precision the critical points and exponents for a number of bond-insertion rules and clarify ambiguities about their universalities. Our findings hold true for any spatial dimensions.

Explosive rigidity percolation in kirigami

Controlling the connectivity and rigidity of kirigami, i.e. the process of cutting paper to deploy it into an articulated system, is critical in the manifestations of kirigami in art, science and technology, as it provides the resulting metamaterial with a range of mechanical and geometric properties. Here, we combine deterministic and stochastic approaches for the control of rigidity in kirigami using the power of k choices, an approach borrowed from the statistical mechanics of explosive percolation transitions. We show that several methods for rigidifying a kirigami system by incrementally changing either the connectivity or the rigidity of individual components allow us to control the nature of the explosive transition by a choice of selection rules. Our results suggest simple lessons for the design of mechanical metamaterials.

Localized growth and remodelling drives spongy mesophyll morphogenesis

The spongy mesophyll is a complex, porous tissue found in plant leaves that enables carbon capture and provides mechanical stability. Unlike many other biological tissues, which remain confluent throughout development, the spongy mesophyll must develop from an initially confluent tissue into a tortuous network of cells with a large proportion of intercellular airspace. How the airspace in the spongy mesophyll develops while the tissue remains mechanically stable is unknown. Here, we use computer simulations of deformable polygons to develop a purely mechanical model for the development of the spongy mesophyll tissue. By stipulating that cell wall growth and remodelling occurs only near void space, our computational model is able to recapitulate spongy mesophyll development observed in Arabidopsis thaliana leaves. We find that robust generation of pore space in the spongy mesophyll requires a balance of cell growth, adhesion, stiffness and tissue pressure to ensure cell networks become porous yet maintain mechanical stability. The success of this mechanical model of morphogenesis suggests that simple physical principles can coordinate and drive the development of complex plant tissues like the spongy mesophyll.

Scaling to zero of compressive modulus in disordered isostatic cubic networks

Networks with as many mechanical constraints as degrees of freedom and no redundant constraints are minimally rigid or isostatic. Isostatic networks are relevant in the study of network glasses, soft matter, and sphere packings. Because of being at the verge of mechanical collapse, they have anomalous elastic and dynamical properties not found in the more commonly occurring hyperstatic networks. In particular, while hyperstatic networks are only slightly affected by geometric disorder, the elastic properties of isostatic networks are dramatically altered by it. In this paper, we show how disorder and system size strongly affect the ability of isostatic networks to sustain a compressive load. We develop an analytic method to calculate the bulk compressive modulus B for various boundary conditions as a function of disorder strength and system size. For simplicity, we consider square and cubic lattices with L^{d} sites, each having d mechanical degrees of freedom, and dL^{d} rotatable springs in the presence of hot-solid disorder of magnitude ε. Additionally, ∼L^{θ} sites may be fixed, thus introducing a nonextensive number of redundancies, either in the bulk or on the boundaries of the system. In all cases, B is analytically and numerically shown to decay as L^{-μ} with μ_{large}=d-θ for large disorder and μ_{small}=max{(d-θ-1),0} for small disorder. Furthermore B(L,ε)L^{μ_{small}}=g(λ) with λ=L^{(μ_{large}-μ_{small})}ε^{2} a scaling variable such that λ<<1 is small disorder and λ>1 is large disorder. The faster decay to zero of B in the large disorder regime results from a broad distribution of spring tensions, including tensions of both signs in equal proportions, which is remarkable since the system is under a purely compressive load. Notably, the bulk modulus is discontinuous at ε=0, a consequence of the fact that the regular network sits at an unstable degenerate configuration.

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.

Rigidity percolation uncovers a structural basis for embryonic tissue phase transitions

Embryo morphogenesis is impacted by dynamic changes in tissue material properties, which have been proposed to occur via processes akin to phase transitions (PTs). Here, we show that rigidity percolation provides a simple and robust theoretical framework to predict material/structural PTs of embryonic tissues from local cell connectivity. By using percolation theory, combined with directly monitoring dynamic changes in tissue rheology and cell contact mechanics, we demonstrate that the zebrafish blastoderm undergoes a genuine rigidity PT, brought about by a small reduction in adhesion-dependent cell connectivity below a critical value. We quantitatively predict and experimentally verify hallmarks of PTs, including power-law exponents and associated discontinuities of macroscopic observables. Finally, we show that this uniform PT depends on blastoderm cells undergoing meta-synchronous divisions causing random and, consequently, uniform changes in cell connectivity. Collectively, our theoretical and experimental findings reveal the structural basis of material PTs in an organismal context.

Spongelike Rigid Structures in Frictional Granular Packings

We show how rigidity emerges in experiments on sheared two-dimensional frictional granular materials by using generalizations of two methods for identifying rigid structures. Both approaches, the force-based dynamical matrix and the topology-based rigidity percolation, agree with each other and identify similar rigid structures. As the system becomes jammed, at a critical contact number z_{c}=2.4±0.1, a rigid backbone interspersed with floppy, particle-filled holes of a broad range of sizes emerges, creating a spongelike morphology. While the pressure within rigid structures always exceeds the pressure outside the rigid structures, they are not identified with the force chains of shear jamming. These findings highlight the need to focus on mechanical stability arising through arch structures and hinges at the mesoscale.

Correlated rigidity percolation in fractal lattices

Rigidity percolation (RP) is the emergence of mechanical stability in networks. Motivated by the experimentally observed fractal nature of materials like colloidal gels and disordered fiber networks, we study RP in a fractal network where intrinsic correlations in particle positions is controlled by the fractal iteration. Specifically, we calculate the critical packing fractions of site-diluted lattices of Sierpiński gaskets (SG's) with varying degrees of fractal iteration. Our results suggest that although the correlation length exponent and fractal dimension of the RP of these lattices are identical to that of the regular triangular lattice, the critical volume fraction is dramatically lower due to the fractal nature of the network. Furthermore, we develop a simplified model for an SG lattice based on the fragility analysis of a single SG. This simplified model provides an upper bound for the critical packing fractions of the full fractal lattice, and this upper bound is strictly obeyed by the disorder averaged RP threshold of the fractal lattices. Our results characterize rigidity in ultralow-density fractal networks.

Two-Scale Scenario of Rigidity Percolation of Sticky Particles

In the presence of attraction, the jamming transition of packings of frictionless particles corresponds to the rigidity percolation. When the range of attraction is long, the distribution of the size of rigid clusters, P(s), is continuous and shows a power-law decay. For systems with short-range attractions, however, P(s) appears discontinuous. There is a power-law decay for small cluster sizes, followed by a low probability gap and a peak near the system size. We find that this appearing "discontinuity" does not mean that the transition is discontinuous. In fact, it signifies the coexistence of two distinct length scales, associated with the largest cluster and smaller ones, respectively. The comparison between the largest and second largest clusters indicates that their growth rates with system size are rather different. However, both cluster sizes tend to diverge in the large system size limit, suggesting that the jamming transition of systems with short-range attractions is still continuous. In the framework of the two-scale scenario, we also derive a generalized hyperscaling relation. With robust evidence, our work challenges the former single-scale view of the rigidity percolation.

Nonequilibrium dynamics of isostatic spring networks

Marginally stable systems exhibit rich critical mechanical behavior. Such isostatic assemblies can be actively driven, but it is unclear how their critical nature affects their nonequilibrium dynamics. Here, we study the influence of isostaticity on the nonequilibrium dynamics of active spring networks. In our model, heterogeneously distributed white or colored, motorlike noise drives the system into a nonequilibrium steady state. We quantify the nonequilibrium dynamics of pairs of network nodes by the characteristic cycling frequency ω-an experimentally accessible measure of the circulation of the associated phase space currents. The distribution of these cycling frequencies exhibits critical scaling, which we approximately capture by a mean-field theory. Finally, we show that the scaling behavior of ω with distance is controlled by a diverging length scale. Overall, we provide a theoretical approach to elucidate the role of marginality in active disordered systems.

Corner-Sharing Tetrahedra for Modeling Micro-structure

This paper introduces Corner-Sharing Tetrahedra (CoSTs), a minimalist, constraint-graph representation of micro-structure. CoSTs have built-in structural guarantees, such as connectivity and minimal rigidity. CoSTs form a space, fully accessible via local operations, that is rich enough to design regular or irregular micro-structure at multiple scales within curved objects. All operations are based on efficient local graph manipulation, which also enables efficient analysis and adjustment of static physical properties. Geometric and material detail, parametric or solid splines, can be added locally, on-demand, for example, for printing.

Mechanical Heterogeneity in Tissues Promotes Rigidity and Controls Cellular Invasion

We study the influence of cell-level mechanical heterogeneity in epithelial tissues using a vertex-based model. Heterogeneity is introduced into the cell shape index (p_{0}) that tunes the stiffness at a single-cell level. The addition of heterogeneity can always enhance the mechanical rigidity of the epithelial layer by increasing its shear modulus, hence making it more rigid. There is an excellent scaling collapse of our data as a function of a single scaling variable f_{r}, which accounts for the overall fraction of rigid cells. We identify a universal threshold f_{r}^{*} that demarcates fluid versus solid tissues. Furthermore, this rigidity onset is far below the contact percolation threshold of rigid cells. These results give rise to a separation of rigidity and contact percolation processes that leads to distinct types of solid states. We also investigate the influence of heterogeneity on tumor invasion dynamics. There is an overall impedance of invasion as the tissue becomes more rigid. Invasion can also occur in an intermediate heterogeneous solid state that is characterized by significant spatial-temporal intermittency.

Correlated Rigidity Percolation and Colloidal Gels

Rigidity percolation (RP) occurs when mechanical stability emerges in disordered networks as constraints or components are added. Here we discuss RP with structural correlations, an effect ignored in classical theories albeit relevant to many liquid-to-amorphous-solid transitions, such as colloidal gelation, which are due to attractive interactions and aggregation. Using a lattice model, we show that structural correlations shift RP to lower volume fractions. Through molecular dynamics simulations, we show that increasing attraction in colloidal gelation increases structural correlation and thus lowers the RP transition, agreeing with experiments. Hence, the emergence of rigidity at colloidal gelation can be understood as a RP transition, but occurs at volume fractions far below values predicted by the classical RP, due to attractive interactions which induce structural correlation.

Excess floppy modes and multibranched mechanisms in metamaterials with symmetries

Floppy modes-deformations that cost zero energy-are central to the mechanics of a wide class of systems. For disordered systems, such as random networks and particle packings, it is well-understood how the number of floppy modes is controlled by the topology of the connections. Here we uncover that symmetric geometries, present in, e.g., mechanical metamaterials, can feature an unlimited number of excess floppy modes that are absent in generic geometries, and in addition can support floppy modes that are multibranched. We study the number Δ of excess floppy modes by comparing generic and symmetric geometries with identical topologies, and show that Δ is extensive, peaks at intermediate connection densities, and exhibits mean-field scaling. We then develop an approximate yet accurate cluster counting algorithm that captures these findings. Finally, we leverage our insights to design metamaterials with multiple folding mechanisms.

Rigidity percolation and geometric information in floppy origami

Origami structures with a large number of excess folds are capable of storing distinguishable geometric states that are energetically equivalent. As the number of excess folds is reduced, the system has fewer equivalent states and can eventually become rigid. We quantify this transition from a floppy to a rigid state as a function of the presence of folding constraints in a classic origami tessellation, Miura-ori. We show that in a fully triangulated Miura-ori that is maximally floppy, adding constraints via the elimination of diagonal folds in the quads decreases the number of degrees of freedom in the system, first linearly and then nonlinearly. In the nonlinear regime, mechanical cooperativity sets in via a redundancy in the assignment of constraints, and the degrees of freedom depend on constraint density in a scale-invariant manner. A percolation transition in the redundancy in the constraints as a function of constraint density suggests how excess folds in an origami structure can be used to store geometric information in a scale-invariant way.

Normal stress anisotropy and marginal stability in athermal elastic networks

Hydrogels of semiflexible biopolymers such as collagen have been shown to contract axially under shear strain, in contrast to the axial dilation observed for most elastic materials. Recent work has shown that this behavior can be understood in terms of the porous, two-component nature and consequent time-dependent compressibility of hydrogels. The apparent normal stress measured by a torsional rheometer reflects only the tensile contribution of the axial component σzz on long (compressible) timescales, crossing over to the first normal stress difference, N1 = σxx - σzz at short (incompressible) times. While the behavior of N1 is well understood for isotropic viscoelastic materials undergoing affine shear deformation, biopolymer networks are often anisotropic and deform nonaffinely. Here, we numerically study the normal stresses that arise under shear in subisostatic, athermal semiflexible polymer networks. We show that such systems exhibit strong deviations from affine behavior and that these anomalies are controlled by a rigidity transition as a function of strain.

RIGID GRAPH COMPRESSION: MOTIF-BASED RIGIDITY ANALYSIS FOR DISORDERED FIBER NETWORKS

Using particle-scale models to accurately describe property enhancements and phase transitions in macroscopic behavior is a major engineering challenge in composite materials science. To address some of these challenges, we use the graph theoretic property of rigidity to model mechanical reinforcement in composites with stiff rod-like particles. We develop an efficient algorithmic approach called rigid graph compression (RGC) to describe the transition from floppy to rigid in disordered fiber networks ("rod-hinge systems"), which form the reinforcing phase in many composite systems. To establish RGC on a firm theoretical foundation, we adapt rigidity matroid theory to identify primitive topological network motifs that serve as rules for composing interacting rigid particles into larger rigid components. This approach is computationally efficient and stable, because RGC requires only topological information about rod interactions (encoded by a sparse unweighted network) rather than geometrical details such as rod locations or pairwise distances (as required in rigidity matroid theory). We conduct numerical experiments on simulated two-dimensional rod-hinge systems to demonstrate that RGC closely approximates the rigidity percolation threshold for such systems, through comparison with the pebble game algorithm (which is exact in two dimensions). Importantly, whereas the pebble game is derived from Laman's condition and is only valid in two dimensions, the RGC approach naturally extends to higher dimensions.

Entropy favors heterogeneous structures of networks near the rigidity threshold

The dynamical properties and mechanical functions of amorphous materials are governed by their microscopic structures, particularly the elasticity of the interaction networks, which is generally complicated by structural heterogeneity. This ubiquitous heterogeneous nature of amorphous materials is intriguingly attributed to a complex role of entropy. Here, we show in disordered networks that the vibrational entropy increases by creating phase-separated structures when the interaction connectivity is close to the onset of network rigidity. The stress energy, which conversely penalizes the heterogeneity, finally dominates a smaller vicinity of the rigidity threshold at the glass transition and creates a homogeneous intermediate phase. This picture of structures changing between homogeneous and heterogeneous phases by varying connectivity provides an interpretation of the transitions observed in chalcogenide glasses.

Normal Stresses, Contraction, and Stiffening in Sheared Elastic Networks

When elastic solids are sheared, a nonlinear effect named after Poynting gives rise to normal stresses or changes in volume. We provide a novel relation between the Poynting effect and the microscopic Grüneisen parameter, which quantifies how stretching shifts vibrational modes. By applying this relation to random spring networks, a minimal model for, e.g., biopolymer gels and solid foams, we find that networks contract or develop tension because they vibrate faster when stretched. The amplitude of the Poynting effect is sensitive to the network's linear elastic moduli, which can be tuned via its preparation protocol and connectivity. Finally, we show that the Poynting effect can be used to predict the finite strain scale where the material stiffens under shear.

Geometry and the onset of rigidity in a disordered network

Disordered spring networks that are undercoordinated may abruptly rigidify when sufficient strain is applied. Since the deformation in response to applied strain does not change the generic quantifiers of network architecture, the number of nodes and the number of bonds between them, this rigidity transition must have a geometric origin. Naive, degree-of-freedom-based mechanical analyses such as the Maxwell-Calladine count or the pebble game algorithm overlook such geometric rigidity transitions and offer no means of predicting or characterizing them. We apply tools that were developed for the topological analysis of zero modes and states of self-stress on regular lattices to two-dimensional random spring networks and demonstrate that the onset of rigidity, at a finite simple shear strain γ^{★}, coincides with the appearance of a single state of self-stress, accompanied by a single floppy mode. The process conserves the topologically invariant difference between the number of zero modes and the number of states of self-stress but imparts a finite shear modulus to the spring network. Beyond the critical shear, the network acquires a highly anisotropic elastic modulus, resisting further deformation most strongly in the direction of the rigidifying shear. We confirm previously reported critical scaling of the corresponding differential shear modulus. In the subcritical regime, a singular value decomposition of the network's compatibility matrix foreshadows the onset of rigidity by way of a continuously vanishing singular value corresponding to the nascent state of self-stress.

Cation-induced monolayer collapse at lower surface pressure follows specific headgroup percolation

A Langmuir monolayer can be considered as a two-dimensional (2D) sheet at higher surface pressure which structurally deform with mechanical compression depending upon the elastic nature of the monolayer. The deformed structures formed after a certain elastic limit are called collapsed structures. To explore monolayer collapses at lower surface pressure and to see the effect of ions on such monolayer collapses, out-of-plane structures and in-plane morphologies of stearic acid Langmuir monolayers have been studied both at lower (≈6.8) and higher (≈9.5) subphase pH in the presence of Mg^{2+},Ca^{2+},Zn^{2+},Cd^{2+}, and Ba^{2+} ions. At lower subphase pH and in the presence of all cations, the stearic acid monolayer remains as a monolayer before collapse, which generally takes place at higher surface pressure (π_{c}>50mN/m). However, at higher subphase pH, structural changes of stearic acid monolayers occur at relatively lower surface pressure depending upon the specific dissolved ions. Among the same group elements of Mg^{2+},Ca^{2+}, and Ba^{2+}, only for Ba^{2+} ions does monolayer to multilayer transition take place from a much lower surface pressure of the monolayer, remaining, however, as a monolayer for Mg^{2+} and Ca^{2+} ions. For another same group elements of Zn^{2+} and Cd^{2+} ions, a less covered bilayer structure forms on top of the monolayer structure at lower surface pressure, which is evidenced from both x-ray reflectometry and atomic force microscopy. Fourier transform infrared spectroscopy confirms the presence of two coexisting conformations formed by the two different metal-headgroup coordinations and the monolayer to trilayer or multilayer transformation takes place when the coverage ratio of the two molecular conformations changes from the critical value (p_{c}) of ≈0.66. Such ion-specific monolayer collapses are correlated with the 2D lattice percolation model.

Critical behaviour in the nonlinear elastic response of hydrogels

In this paper we study the elastic response of synthetic hydrogels to an applied shear stress. The hydrogels studied here have previously been shown to mimic the behaviour of biopolymer networks when they are sufficiently far above the gel point. We show that near the gel point they exhibit an elastic response that is consistent with the predicted critical behaviour of networks near or below the isostatic point of marginal stability. This point separates rigid and floppy states, distinguished by the presence or absence of finite linear elastic moduli. Recent theoretical work has also focused on the response of such networks to finite or large deformations, both near and below the isostatic point. Despite this interest, experimental evidence for the existence of criticality in such networks has been lacking. Using computer simulations, we identify critical signatures in the mechanical response of sub-isostatic networks as a function of applied shear stress. We also present experimental evidence consistent with these predictions. Furthermore, our results show the existence of two distinct critical regimes, one of which arises from the nonlinear stretch response of semi-flexible polymers.

The loss of tension in an infinite membrane with holes distributed according to a Poisson law

What is the effect of punching holes at random in an infinite tensed membrane? When will the membrane still support tension? This problem was introduced by Connelly in connection with applications of rigidity theory to natural sciences. The answer clearly depends on the shapes and the distribution of the holes. We briefly outline a mathematical theory of tension based on graph rigidity theory and introduce a probabilistic model for this problem. We show that if the centers of the holes are distributed in ℝ2 according to a Poisson law with density λ &gt; 0, and the shapes are i.i.d. and independent of the locations of their centers, the tension is lost on all of ℝ2 for any λ. After introducing a certain step-by-step dynamic for the loss of tension, we establish the existence of a nonrandom N = N(λ) such that N steps are almost surely enough for the loss of tension. Also, we prove that N(λ) &gt; 2 almost surely for sufficiently small λ. The processes described in the paper are related to bootstrap and rigidity percolation.

The loss of tension in an infinite membrane with holes distributed according to a Poisson law

What is the effect of punching holes at random in an infinite tensed membrane? When will the membrane still support tension? This problem was introduced by Connelly in connection with applications of rigidity theory to natural sciences. The answer clearly depends on the shapes and the distribution of the holes. We briefly outline a mathematical theory of tension based on graph rigidity theory and introduce a probabilistic model for this problem. We show that if the centers of the holes are distributed in ℝ2 according to a Poisson law with density λ > 0, and the shapes are i.i.d. and independent of the locations of their centers, the tension is lost on all of ℝ2 for any λ. After introducing a certain step-by-step dynamic for the loss of tension, we establish the existence of a nonrandom N = N(λ) such that N steps are almost surely enough for the loss of tension. Also, we prove that N(λ) > 2 almost surely for sufficiently small λ. The processes described in the paper are related to bootstrap and rigidity percolation.

Topological similarity of random cell complexes and applications

Although random cell complexes occur throughout the physical sciences, there does not appear to be a standard way to quantify their statistical similarities and differences. The various proposals in the literature are usually motivated by the analysis of particular physical systems and do not necessarily apply to general situations. The central concepts in this paper-the swatch and the cloth-provide a description of the local topology of a cell complex that is general (any physical system that can be represented as a cell complex is admissible) and complete (any statistical question about the local topology can be answered from the cloth). Furthermore, this approach allows a distance to be defined that measures the similarity of the local topology of two cell complexes. The distance is used to identify a steady state of a model grain boundary network, quantify the approach to this steady state, and show that the steady state is independent of the initial conditions. The same distance is then employed to show that the long-term properties in simulations of a specific model of a dislocation network do not depend on the implementation of dislocation intersections.

Interlocking-induced stiffness in stochastically microcracked materials beyond the transport percolation threshold

We study the mechanical behavior of two-dimensional, stochastically microcracked continua in the range of crack densities close to, and above, the transport percolation threshold. We show that these materials retain stiffness up to crack densities much larger than the transport percolation threshold due to topological interlocking of sample subdomains. Even with a linear constitutive law for the continuum, the mechanical behavior becomes nonlinear in the range of crack densities bounded by the transport and stiffness percolation thresholds. The effect is due to the fractal nature of the fragmentation process and is not linked to the roughness of individual cracks.

Viscoelastic properties of marginal networks in a solvent

Polymer networks at the margins of mechanical stability are known to be highly sensitive to applied forces and fields and to exhibit an anomalously large resistance to deformation. In this paper, we study the effects of hydrodynamic interactions on the behavior of marginal networks using a hybrid molecular dynamics and multiparticle collision dynamics simulation technique. We examine how the filament and solvent properties affect the response of marginal networks to shear. We find that the stiffening of the network shows a stronger dependence on the shear frequency when hydrodynamic interactions are present than when they are not. The network shear modulus scales as G'∼ω(α(c)), with a critical stiffening exponent α(c) that can be controlled by varying the relative concentrations of the network and the solvent. Our results show that this arises due to the solvent aiding the relaxation of the network and suppressing the network nonaffinity, with the system deforming more affinely when hydrodynamic interactions are maximized.

Finite-temperature mechanical instability in disordered lattices

Mechanical instability takes different forms in various ordered and disordered systems and little is known about how thermal fluctuations affect different classes of mechanical instabilities. We develop an analytic theory involving renormalization of rigidity and coherent potential approximation that can be used to understand finite-temperature mechanical stabilities in various disordered systems. We use this theory to study two disordered lattices: a randomly diluted triangular lattice and a randomly braced square lattice. These two lattices belong to two different universality classes as they approach mechanical instability at T=0. We show that thermal fluctuations stabilize both lattices. In particular, the triangular lattice displays a critical regime in which the shear modulus scales as G∼T(1/2), whereas the square lattice shows G∼T(2/3). We discuss generic scaling laws for finite-T mechanical instabilities and relate them to experimental systems.

Anchored boundary conditions for locally isostatic networks

Finite pieces of locally isostatic networks have a large number of floppy modes because of missing constraints at the surface. Here we show that by imposing suitable boundary conditions at the surface the network can be rendered effectively isostatic. We refer to these as anchored boundary conditions. An important example is formed by a two-dimensional network of corner sharing triangles, which is the focus of this paper. Another way of rendering such networks isostatic is by adding an external wire along which all unpinned vertices can slide (sliding boundary conditions). This approach also allows for the incorporation of boundaries associated with internal holes and complex sample geometries, which are illustrated with examples. The recent synthesis of bilayers of vitreous silica has provided impetus for this work. Experimental results from the imaging of finite pieces at the atomic level need such boundary conditions, if the observed structure is to be computer refined so that the interior atoms have the perception of being in an infinite isostatic environment.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Phonons and elasticity in critically coordinated lattices

Much of our understanding of vibrational excitations and elasticity is based upon analysis of frames consisting of sites connected by bonds occupied by central-force springs, the stability of which depends on the average number of neighbors per site z. When z  <  zc  ≈  2d, where d is the spatial dimension, frames are unstable with respect to internal deformations. This pedagogical review focuses on the properties of frames with z at or near zc, which model systems like randomly packed spheres near jamming and network glasses. Using an index theorem, N0  -NS  =  dN  -NB relating the number of sites, N, and number of bonds, NB, to the number, N0, of modes of zero energy and the number, NS, of states of self stress, in which springs can be under positive or negative tension while forces on sites remain zero, it explores the properties of periodic square, kagome, and related lattices for which z  =  zc and the relation between states of self stress and zero modes in periodic lattices to the surface zero modes of finite free lattices (with free boundary conditions). It shows how modifications to the periodic kagome lattice can eliminate all but trivial translational zero modes and create topologically distinct classes, analogous to those of topological insulators, with protected zero modes at free boundaries and at interfaces between different topological classes.

Rigidity loss in disordered systems: three scenarios

We reveal significant qualitative differences in the rigidity transition of three types of disordered network materials: randomly diluted spring networks, jammed sphere packings, and stress-relieved networks that are diluted using a protocol that avoids the appearance of floppy regions. The marginal state of jammed and stress-relieved networks are globally isostatic, while marginal randomly diluted networks show both overconstrained and underconstrained regions. When a single bond is added to or removed from these isostatic systems, jammed networks become globally overconstrained or floppy, whereas the effect on stress-relieved networks is more local and limited. These differences are also reflected in the linear elastic properties and point to the highly effective and unusual role of global self-organization in jammed sphere packings.

Rigidity percolation by next-nearest-neighbor bonds on generic and regular isostatic lattices

We study rigidity percolation transitions in two-dimensional central-force isostatic lattices, including the square and the kagome lattices, as next-nearest-neighbor bonds ("braces") are randomly added to the system. In particular, we focus on the differences between regular lattices, which are perfectly periodic, and generic lattices with the same topology of bonds but whose sites are at random positions in space. We find that the regular square and kagome lattices exhibit a rigidity percolation transition when the number of braces is ∼LlnL, where L is the linear size of the lattice. This transition exhibits features of both first-order and second-order transitions: The whole lattice becomes rigid at the transition, and a diverging length scale also exists. In contrast, we find that the rigidity percolation transition in the generic lattices occur when the number of braces is very close to the number obtained from Maxwell's law for floppy modes, which is ∼L. The transition in generic lattices is a very sharp first-order-like transition, at which the addition of one brace connects all small rigid regions in the bulk of the lattice, leaving only floppy modes on the edge. We characterize these transitions using numerical simulations and develop analytic theories capturing each transition. Our results relate to other interesting problems, including jamming and bootstrap percolation.

Mechanics of anisotropic spring networks

We construct and analyze a model for a disordered linear spring network with anisotropy. The modeling is motivated by, for example, granular systems, nematic elastomers, and ultimately cytoskeletal networks exhibiting some underlying anisotropy. The model consists of a triangular lattice with two different bond occupation probabilities, p(x) and p(y), for the linear springs. We develop an effective medium theory (EMT) to describe the network elasticity as a function of p(x) and p(y). We find that the onset of rigidity in the EMT agrees with Maxwell constraint counting. We also find beyond linear behavior in the shear and bulk modulus as a function of occupation probability in the rigid phase for small strains, which differs from the isotropic case. We compare our EMT with numerical simulations to find rather good agreement. Finally, we discuss the implications of extending the reach of effective medium theory as well as draw connections with prior work on both anisotropic and isotropic spring networks.

Evolution of covalent networks under cooling: contrasting the rigidity window and jamming scenarios

We study the evolution of structural disorder under cooling in supercooled liquids, focusing on covalent networks. We introduce a model for the energy of networks that incorporates weak noncovalent interactions. We show that at low temperature these interactions considerably affect the network topology near the rigidity transition that occurs as the coordination increases. As a result, this transition becomes mean field and does not present a line of critical points previously argued for, the "rigidity window." Vibrational modes are then not fractons but instead are similar to the anomalous modes observed in packings of particles near jamming. These results suggest an alternative interpretation for the intermediate phase observed in chalcogenides.