Search

The mechanical properties of biological materials are spatially heterogeneous. Typical tissues are made up of a spanning fibrous extracellular matrix in which various inclusions, such as living cells, are embedded. While the influence of inclusions on the stiffness of common elastic materials such as rubber has been studied for decades and can be understood in terms of the volume fraction and shape of inclusions, the same is not true for disordered filamentous and fibrous networks. Recent work has shown that, in isolation, such networks exhibit unusual viscoelastic behavior indicative of an underlying mechanical phase transition controlled by network connectivity and strain. How this behavior is modified when inclusions are present is unclear. Here, we present a theoretical and computational study of the influence of rigid inclusions on the mechanics of disordered elastic networks near the connectivity-controlled central force rigidity transition. Combining scaling theory and coarse-grained simulations, we predict and confirm an anomalously strong dependence of the composite stiffness on inclusion volume fraction, beyond that seen in ordinary composites. This stiffening exceeds the well-established volume fraction-dependent stiffening expected in conventional composites, e.g., as an elastic analogue of the classic volume fraction dependent increase in the viscosity of liquids first identified by Einstein. We show that this enhancement is a consequence of the interplay between inter-particle spacing and an emergent correlation length, leading to an effective finite-size scaling imposed by the presence of inclusions. We outline the expected scaling of the shear modulus and strain fluctuations with the inclusion volume fraction and network connectivity, confirm these predictions in simulations, and discuss potential experimental tests and implications for our predictions in real systems.

Stretching an elastic material along one axis typically induces contraction along the transverse axes, a phenomenon known as the Poisson effect. From these strains, one can compute the specific volume, which generally either increases or, in the incompressible limit, remains constant as the material is stretched. However, in networks of semiflexible or stiff polymers, which are typically highly compressible yet stiffen significantly when stretched, one instead sees a significant reduction in specific volume under finite strains. This volume reduction is accompanied by increasing alignment of filaments along the strain axis and a nonlinear elastic response, with stiffening of the apparent Young's modulus. For semiflexible networks, in which entropic bending elasticity governs the linear elastic regime, the nonlinear Poisson effect is caused by the nonlinear force-extension relationship of the constituent filaments, which produces a highly asymmetric response of the constituent polymers to stretching and compression. The details of this relationship depend on the geometric and elastic properties of the underlying filaments, which can vary greatly in experimental systems. Here, we provide a comprehensive characterization of the nonlinear Poisson effect in an affine network model and explore the influence of filament properties on essential features of the macroscopic response, including strain-driven alignment and volume reduction.

In active materials, uncoordinated internal stresses lead to emergent long-range flows. An understanding of how the behavior of active materials depends on mesoscopic (hydrodynamic) parameters is developing, but there remains a gap in knowledge concerning how hydrodynamic parameters depend on the properties of microscopic elements. In this work, we combine experiments and multiscale modeling to relate the structure and dynamics of active nematics composed of biopolymer filaments and molecular motors to their microscopic properties, in particular motor processivity, speed, and valency. We show that crosslinking of filaments by both motors and passive crosslinkers not only augments the contributions to nematic elasticity from excluded volume effects but dominates them. By altering motor kinetics we show that a competition between motor speed and crosslinking results in a nonmonotonic dependence of nematic flow on motor speed. By modulating passive filament crosslinking we show that energy transfer into nematic flow is in large part dictated by crosslinking. Thus motor proteins both generate activity and contribute to nematic elasticity. Our results provide new insights for rationally engineering active materials.

Fibrous networks such as collagen are common in biological systems. Recent theoretical and experimental efforts have shed light on the mechanics of single component networks. Most real biopolymer networks, however, are composites made of elements with different rigidity. For instance, the extracellular matrix in mammalian tissues consists of stiff collagen fibers in a background matrix of flexible polymers such as hyaluronic acid (HA). The interplay between different biopolymer components in such composite networks remains unclear. In this work, we use 2D coarse-grained models to study the nonlinear strain-stiffening behavior of composites. We introduce a local volume constraint to model the incompressibility of HA. We also perform rheology experiments on composites of collagen with HA. We demonstrate both theoretically and experimentally that the linear shear modulus of composite networks can be increased by approximately an order of magnitude above the corresponding moduli of the pure components. Our model shows that this synergistic effect can be understood in terms of the local incompressibility of HA, which acts to suppress density fluctuations of the collagen matrix with which it is entangled.

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 highly strain-sensitive mechanical properties. In this work, we perform simulations of three-dimensional networks produced by the irreversible formation of crosslinks 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 crosslinker availability and detail the associated connectivity dependence of both linear elasticity and nonlinear strain stiffening behavior, drawing comparisons with prior experimental measurements of the crosslinker concentration-dependent elasticity of actin gels.

Flocking in $d=2$ is a genuine non-equilibrium phenomenon for which irreversibility is an essential ingredient. We study a class of minimal flocking models whose only source of irreversibility is self-propulsion and use the entropy production rate (EPR) to quantify the departure from equilibrium across their phase diagrams. The EPR is maximal in the vicinity of the order-disorder transition, where reshuffling of the interaction network is fast. We show that signatures of irreversibility come in the form of asymmetries in the steady state distribution of the flock's microstates. They occur as consequences of the time reversal symmetry breaking in the considered self-propelled systems, independently of the interaction details. In the case of metric pairwise forces, they reduce to local asymmetries in the distribution of pairs of particles. This study suggests a possible use of pair asymmetries both to quantify the departure from equilibrium and to learn relevant information about aligning interaction potentials from data., Comment: 8 pages + Appendix; 6 figures

Networks and dense suspensions frequently reside near a boundary between soft (or fluid-like) and rigid (or solid-like) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.