Normal stress anisotropy and marginal stability in athermal elastic networks

Molecular dynamics simulations of microscopic structural transition and macroscopic mechanical properties of magnetic gels

Magnetic gels with embedded micro-/nano-sized magnetic particles in cross-linked polymer networks can be actuated by external magnetic fields, with changes in their internal microscopic structures and macroscopic mechanical properties. We investigate the responses of such magnetic gels to an external magnetic field, by means of coarse-grained molecular dynamics simulations. We find that the dynamics of magnetic particles are determined by the interplay of magnetic dipole-dipole interactions, polymer elasticity, and thermal fluctuations. The corresponding microscopic structures formed by the magnetic particles, such as elongated chains, can be controlled by the external magnetic field. Furthermore, the magnetic gels can exhibit reinforced macroscopic mechanical properties, where the elastic modulus increases algebraically with the magnetic moments of the particles in the form of ∝(m-mc)2 when magnetic chains are formed. This simulation work can not only serve as a tool for studying the microscopic and the macroscopic responses of the magnetic gels, but also facilitate future fabrications and practical controls of magnetic composites with desired physical properties.

Effective medium theory for mechanical phase transitions of fiber networks

Networks of stiff fibers govern the elasticity of biological structures such as the extracellular matrix of collagen. These networks are known to stiffen nonlinearly under shear or extensional strain. Recently, it has been shown that such stiffening is governed by a strain-controlled athermal but critical phase transition, from a floppy phase below the critical strain to a rigid phase above the critical strain. While this phase transition has been extensively studied numerically and experimentally, a complete analytical theory for this transition remains elusive. Here, we present an effective medium theory (EMT) for this mechanical phase transition of fiber networks. We extend a previous EMT appropriate for linear elasticity to incorporate nonlinear effects via an anharmonic Hamiltonian. The mean-field predictions of this theory, including the critical exponents, scaling relations and non-affine fluctuations qualitatively agree with previous experimental and numerical results.

Cross-linker Mobility Governs Fracture Behavior of Catch-Bonded Networks

While most chemical bonds weaken under the action of mechanical force (called slip bond behavior), nature has developed bonds that do the opposite: their lifetime increases as force is applied. While such catch bonds have been studied quite extensively at the single molecule level and in adhesive contacts, recent work has shown that they are also abundantly present as crosslinkers in the actin cytoskeleton. However, their role and the mechanism by which they operate in these networks have remained unclear. Here, we present computer simulations that show how polymer networks crosslinked with either slip or catch bonds respond to mechanical stress. Our results reveal that catch bonding may be required to protect dynamic networks against fracture, in particular for mobile linkers that can diffuse freely after unbinding. While mobile slip bonds lead to networks that are very weak at high stresses, mobile catch bonds accumulate in high stress regions and thereby stabilize cracks, leading to a more ductile fracture behavior. This allows cells to combine structural adaptivity at low stresses with mechanical stability at high stresses.

Structural Features and Nonlinear Rheology of Self-Assembled Networks of Cross-Linked Semiflexible Polymers

Disordered networks of semiflexible filaments are common support structures in biology. Familiar examples include fibrous matrices in blood clots, bacterial biofilms, and essential components of cells and tissues of plants, animals, and fungi. Despite the ubiquity of these networks in biomaterials, we have only a limited understanding of the relationship between their structural features and their highly strain-sensitive mechanical properties. In this work, we perform simulations of three-dimensional networks produced by the irreversible formation of cross-links between linker-decorated semiflexible filaments. We characterize the structure of networks formed by a simple diffusion-dependent assembly process and measure their associated steady-state rheological features at finite temperature over a range of applied prestrains that encompass the strain-stiffening transition. We quantify the dependence of network connectivity on cross-linker availability and detail the associated connectivity dependence of both linear elasticity and nonlinear strain-stiffening behavior, drawing comparisons with prior experimental measurements of the cross-linker concentration-dependent elasticity of actin gels.

Mechanics of fiber networks under a bulk strain

Biopolymer networks are common in biological systems from the cytoskeleton of individual cells to collagen in the extracellular matrix. The mechanics of these systems under applied strain can be explained in some cases by a phase transition from soft to rigid states. For collagen networks, it has been shown that this transition is critical in nature and it is predicted to exhibit diverging fluctuations near a critical strain that depends on the network's connectivity and structure. Whereas prior work focused mostly on shear deformation that is more accessible experimentally, here we study the mechanics of such networks under an applied bulk or isotropic extension. We confirm that the bulk modulus of subisostatic fiber networks exhibits similar critical behavior as a function of bulk strain. We find different nonmean-field exponents for bulk as opposed to shear. We also confirm a similar hyperscaling relation to what was previously found for shear.

Stretchy and disordered: Toward understanding fracture in soft network materials via mesoscopic computer simulations

Soft network materials exist in numerous forms ranging from polymer networks, such as elastomers, to fiber networks, such as collagen. In addition, in colloidal gels, an underlying network structure can be identified, and several metamaterials and textiles can be considered network materials as well. Many of these materials share a highly disordered microstructure and can undergo large deformations before damage becomes visible at the macroscopic level. Despite their widespread presence, we still lack a clear picture of how the network structure controls the fracture processes of these soft materials. In this Perspective, we will focus on progress and open questions concerning fracture at the mesoscopic scale, in which the network architecture is clearly resolved, but neither the material-specific atomistic features nor the macroscopic sample geometries are considered. We will describe concepts regarding the network elastic response that have been established in recent years and turn out to be pre-requisites to understand the fracture response. We will mostly consider simulation studies, where the influence of specific network features on the material mechanics can be cleanly assessed. Rather than focusing on specific systems, we will discuss future challenges that should be addressed to gain new fundamental insights that would be relevant across several examples of soft network materials.

Shear-induced phase transition and critical exponents in three-dimensional fiber networks

When subject to applied strain, fiber networks exhibit nonlinear elastic stiffening. Recent theory and experiments have shown that this phenomenon is controlled by an underlying mechanical phase transition that is critical in nature. Growing simulation evidence points to non-mean-field behavior for this transition and a hyperscaling relation has been proposed to relate the corresponding critical exponents. Here, we report simulations on two distinct network structures in three dimensions. By performing a finite-size scaling analysis, we test hyperscaling and identify various critical exponents. From the apparent validity of hyperscaling, as well as the non-mean-field exponents we observe, our results suggest that the upper critical dimension for the strain-controlled phase transition is above three, in contrast to the jamming transition that represents another athermal, mechanical phase transition.

Strain Stiffening and Negative Normal Force of Agarose Hydrogels

Inspired by the specific strain stiffening and negative normal force phenomena in several biological networks, herein, we show strain stiffening and negative normal force in agarose hydrogels. We use both pre-strain and strain amplitude sweep protocols in dynamic rheological measurements where the gel slip was suppressed by the in situ gelation in the cross-hatched parallel plate rheometer geometry. Within the stiffening region, we show the scaling relation for the differential modulus K ∝ σ1, where σ is stress. The strain at the onset of stiffening is almost constant throughout the concentration range. The gels show negative apparent normal stress difference when sheared as a result of the gel contraction. The pore size of the hydrogel is large enough to allow water to move with respect to the network to balance the pressure difference caused by the hoop stress. The rheological analysis together with scanning electron microscopy suggests that the agarose gels can be described using subisostatic athermal network models where the connectivity dictates the stiffening behavior. Therefore, the simple agarose gels appear to capture several of the viscoelastic properties, which were previously thought to be characteristic to biological protein macromolecules.

Finite size effects in critical fiber networks

Fibrous networks such as collagen are common in physiological systems. One important function of these networks is to provide mechanical stability for cells and tissues. At physiological levels of connectivity, such networks would be mechanically unstable with only central-force interactions. While networks can be stabilized by bending interactions, it has also been shown that they exhibit a critical transition from floppy to rigid as a function of applied strain. Beyond a certain strain threshold, it is predicted that underconstrained networks with only central-force interactions exhibit a discontinuity in the shear modulus. We study the finite-size scaling behavior of this transition and identify both the mechanical discontinuity and critical exponents in the thermodynamic limit. We find both non-mean-field behavior and evidence for a hyperscaling relation for the critical exponents, for which the network stiffness is analogous to the heat capacity for thermal phase transitions. Further evidence for this is also found in the self-averaging properties of fiber networks.

Nonlinear Poisson Effect Governed by a Mechanical Critical Transition

Under extensional strain, fiber networks can exhibit an anomalously large and nonlinear Poisson effect accompanied by a dramatic transverse contraction and volume reduction for applied strains as small as a few percent. We demonstrate that this phenomenon is controlled by a collective mechanical phase transition that occurs at a critical uniaxial strain that depends on network connectivity. This transition is punctuated by an anomalous peak in the apparent Poisson's ratio and other critical signatures such as diverging nonaffine strain fluctuations.

Athermal Fracture of Elastic Networks: How Rigidity Challenges the Unavoidable Size-Induced Brittleness

By performing extensive simulations with unprecedentedly large system sizes, we unveil how rigidity influences the fracture of disordered materials. We observe the largest damage in networks with connectivity close to the isostatic point and when the rupture thresholds are small. However, irrespective of network and spring properties, a more brittle fracture is observed upon increasing system size. Differently from most of the fracture descriptors, the maximum stress drop, a proxy for brittleness, displays a universal nonmonotonic dependence on system size. Based on this uncommon trend it is possible to identify the characteristic system size L^{*} at which brittleness kicks in. The more the disorder in network connectivity or in spring thresholds, the larger L^{*}. Finally, we speculate how this size-induced brittleness is influenced by thermal fluctuations.

Force chains in cell-cell mechanical communication

Force chains (FCs) are a key determinant of the micromechanical properties and behaviour of heterogeneous materials, such as granular systems. However, less is known about FCs in fibrous materials, such as the networks composing the extracellular matrix (ECM) of biological systems. Using a finite-element computational model, we simulated the contraction of a single cell and two nearby cells embedded in two-dimensional fibrous elastic networks and analysed the tensile FCs that developed in the ECM. The role of ECM nonlinear elasticity on FC formation was evaluated by considering linear and nonlinear, i.e. exhibiting 'buckling' and/or 'strain-stiffening', stress-strain curves. The effect of the degree of cell contraction and network coordination value was assessed. We found that nonlinear elasticity of the ECM fibres influenced the structure of the FCs, facilitating the transition towards more distinct chains that were less branched and more radially oriented than the chains formed in linear elastic networks. When two neighbouring cells contract, a larger number of FCs bridged between the cells in nonlinear networks, and these chains had a larger effective rigidity than the chains that did not reach a neighbouring cell. These results suggest that FCs function as a route for mechanical communication between distant cells and highlight the contribution of ECM fibre nonlinear elasticity to the formation of FCs.

Strain-Stiffening of Agarose Gels

Strain-stiffening is one of the characteristic properties of biological hydrogels and extracellular matrices, where the stiffness increases upon increased deformation. Whereas strain-stiffening is ubiquitous in protein-based materials, it has been less observed for polysaccharide and synthetic polymer gels. Here we show that agarose, that is, a common linear polysaccharide, forms helical fibrillar bundles upon cooling from aqueous solution. The hydrogels with these semiflexible fibrils show pronounced strain-stiffening. However, to reveal strain-stiffening, suppressing wall slippage turned as untrivial. Upon exploring different sample preparation techniques and rheological architectures, the cross-hatched parallel plate geometries and in situ gelation in the rheometer successfully prevented the slippage and resolved the strain-stiffening behavior. Combining with microscopy, we conclude that strain-stiffening is due to the semiflexible nature of the agarose fibrils and their geometrical connectivity, which is below the central-force isostatic critical connectivity. The biocompatibility and the observed strain-stiffening suggest the potential of agarose hydrogels in biomedical applications.

A minimal-length approach unifies rigidity in underconstrained materials

We present an approach to understand geometric-incompatibility-induced rigidity in underconstrained materials, including subisostatic 2D spring networks and 2D and 3D vertex models for dense biological tissues. We show that in all these models a geometric criterion, represented by a minimal length [Formula: see text], determines the onset of prestresses and rigidity. This allows us to predict not only the correct scalings for the elastic material properties, but also the precise magnitudes for bulk modulus and shear modulus discontinuities at the rigidity transition as well as the magnitude of the Poynting effect. We also predict from first principles that the ratio of the excess shear modulus to the shear stress should be inversely proportional to the critical strain with a prefactor of 3. We propose that this factor of 3 is a general hallmark of geometrically induced rigidity in underconstrained materials and could be used to distinguish this effect from nonlinear mechanics of single components in experiments. Finally, our results may lay important foundations for ways to estimate [Formula: see text] from measurements of local geometric structure and thus help develop methods to characterize large-scale mechanical properties from imaging data.