Soft modes and elasticity of nearly isostatic lattices: randomness and dissipation

Testing the heterogeneous-elasticity theory for low-energy excitations in structural glasses

Understanding the statistical mechanics of low-energy excitations in structural glasses has been the focus of extensive research efforts in the past decades due to their key roles in determining the low-temperature mechanical and transport properties of these intrinsically nonequilibrium materials. While it is established that glasses feature low-energy nonphononic excitations that follow a non-Debye vibrational density of states, we currently lack a well-founded theory of these fundamental objects and their vibrational spectra. A recent theory-that builds on the so-called heterogeneous-elasticity theory (HET) and its extensions-provides explicit predictions for the scaling of the low-frequency tail of the nonphononic spectrum of glasses, the localization properties of the vibrational modes that populate this tail, and its connections to glass formation history and to the form of the distribution of weak microscopic (interatomic) stresses. Here, we employ computer models of structural glasses to quantitatively test these predictions. Our findings do not support the HET's predictions regarding the nature and statistics of low-energy excitations in glasses, highlighting the need for additional theoretical developments.

Achieving adjustable elasticity with non-affine to affine transition

For various engineering and industrial applications it is desirable to realize mechanical systems with broadly adjustable elasticity to respond flexibly to the external environment. Here we discover a topology-correlated transition between affine and non-affine regimes in elasticity in both two- and three-dimensional packing-derived networks. Based on this transition, we numerically design and experimentally realize multifunctional systems with adjustable elasticity. Within one system, we achieve solid-like affine response, liquid-like non-affine response and a continuous tunability in between. Moreover, the system also exhibits a broadly tunable Poisson's ratio from positive to negative values, which is of practical interest for energy absorption and for fracture-resistant materials. Our study reveals a fundamental connection between elasticity and network topology, and demonstrates its practical potential for designing mechanical systems and metamaterials.

Directed force propagation in semiflexible networks

We consider the propagation of tension along specific filaments of a semiflexible filament network in response to the application of a point force using a combination of numerical simulations and analytic theory. We find the distribution of force within the network is highly heterogeneous, with a small number of fibers supporting a significant fraction of the applied load over distances of multiple mesh sizes surrounding the point of force application. We suggest that these structures may be thought of as tensile force chains, whose structure we explore via simulation. We develop self-consistent calculations of the point-force response function and introduce a transfer matrix approach to explore the decay of tension (into bending) energy and the branching of tensile force chains in the network.

Mechanics of Metric Frustration in Contorted Filament Bundles: From Local Symmetry to Columnar Elasticity

Bundles of filaments are subject to geometric frustration: certain deformations (e.g., bending while twisted) require longitudinal variations in spacing between filaments. While bundles are common-from protein fibers to yarns-the mechanical consequences of longitudinal frustration are unknown. We derive a geometrically nonlinear formalism for bundle mechanics, using a gaugelike symmetry under reptations along filament backbones. We relate force balance to orientational geometry and assess the elastic cost of frustration in twisted-toroidal bundles.

Revealing viscoelastic bending relaxation dynamics of isolated semiflexible colloidal polymers

The viscoelastic properties of filaments and biopolymers play a crucial role in soft and biological materials from biopolymer networks to novel synthetic metamaterials. Colloidal particles with specific valency allow mimicking polymers and more complex molecular structures at the colloidal scale, offering direct observation of their internal degrees of freedom. Here, we elucidate the time-dependent viscoelastic response in the bending of isolated semi-flexible colloidal polymers, assembled from dipatch colloidal particles by reversible critical Casimir forces. By tuning the patch-patch interaction strength, we adjust the polymers' viscoelastic properties, and follow spontaneous bending modes and their relaxation directly on the particle level. We find that the elastic response is well described by that of a semiflexible rod with persistence length of order 1000 μm, tunable by the critical Casimir interaction strength. We identify the viscous relaxation on longer timescales to be due to internal friction, leading to a wavelength-independent relaxation time similar to single biopolymers, but in the colloidal case arising from the contact mechanics of the bonded patches. These tunable mechanical properties of assembled colloidal filaments open the door to "colloidal architectures", rationally designed (network) structures with desired topology and mechanical properties.

Elasticity of colloidal gels: structural heterogeneity, floppy modes, and rigidity

Rheological measurements of model colloidal gels reveal that large variations in the shear moduli as colloidal volume-fraction changes are not reflected by simple structural parameters such as the coordination number, which remains almost a constant. We resolve this apparent contradiction by conducting a normal-mode analysis of experimentally measured bond networks of gels of colloidal particles with short-ranged attraction. We find that structural heterogeneity of the gels, which leads to floppy modes and a nonaffine-affine crossover as frequency increases, evolves as a function of the volume fraction and is key to understanding the frequency-dependent elasticity. Without any free parameters, we achieve good qualitative agreement with the measured mechanical response. Furthermore, we achieve universal collapse of the shear moduli through a phenomenological spring-dashpot model that accounts for the interplay between fluid viscosity, particle dissipation, and contributions from the affine and non-affine network deformation.

Continuum Theory for Topological Edge Soft Modes

Topological edge zero modes and states of self stress have been intensively studied in discrete lattices at the Maxwell point, offering robust properties concerning surface and interface stiffness and stress focusing. In this Letter, we present a topological elasticity theory for general continuous media where a gauge-invariant bulk topological index independent of microscopic details is defined. This index directly predicts the number of zero modes on edges at long length scales, and it naturally extends to media that deviate from the Maxwell point, depicting how topological zero modes turn into topological soft modes.

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.

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.

Autonomous Actuation of Zero Modes in Mechanical Networks Far from Equilibrium

A zero mode, or floppy mode, is a nontrivial coupling of mechanical components yielding a degree of freedom with no resistance to deformation. Engineered zero modes have the potential to act as microscopic motors or memory devices, but this requires an internal actuation mechanism that can overcome unwanted fluctuations in other modes and the dissipation inherent in real systems. In this Letter, we show theoretically and experimentally that complex zero modes in mechanical networks can be selectively mobilized by nonequilibrium activity. We find that a correlated active bath actuates an infinitesimal zero mode while simultaneously suppressing fluctuations in higher modes compared to thermal fluctuations, which we experimentally mimic by high frequency shaking of a physical network. Furthermore, self-propulsive dynamics spontaneously mobilize finite mechanisms as exemplified by a self-propelled topological soliton. Nonequilibrium activity thus enables autonomous actuation of coordinated mechanisms engineered through network topology.

Elastic Backbone Defines a New Transition in the Percolation Model

The elastic backbone is the set of all shortest paths. We found a new phase transition at p_{eb} above the classical percolation 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 exponents β_{eb}=0.50±0.02, γ_{eb}=1.97±0.05, and ν_{eb}=2.00±0.02, fulfilling consistent critical scaling laws. Interestingly, however, the hyperscaling relation is violated. Using Binder's cumulant, we determine, with high precision, the critical probabilities p_{eb} for the triangular and tilted square lattice for site and bond percolation. This transition describes a sudden rigidification as a function of density when stretching a damaged tissue.

Topological Edge Floppy Modes in Disordered Fiber Networks

Disordered fiber networks are ubiquitous in a broad range of natural (e.g., cytoskeleton) and manmade (e.g., aerogels) materials. In this Letter, we discuss the emergence of topological floppy edge modes in two-dimensional fiber networks as a result of deformation or active driving. It is known that a network of straight fibers exhibits bulk floppy modes which only bend the fibers without stretching them. We find that, interestingly, with a perturbation in geometry, these bulk modes evolve into edge modes. We introduce a topological index for these edge modes and discuss their implications in biology.

Unjamming in models with analytic pairwise potentials

Canonical models for studying the unjamming scenario in systems of soft repulsive particles assume pairwise potentials with a sharp cutoff in the interaction range. The sharp cutoff renders the potential nonanalytic but makes it possible to describe many properties of the solid in terms of the coordination number z, which has an unambiguous definition in these cases. Pairwise potentials without a sharp cutoff in the interaction range have not been studied in this context, but should in fact be considered to understand the relevance of the unjamming phenomenology in systems where such a cutoff is not present. In this work we explore two systems with such interactions: an inverse power law and an exponentially decaying pairwise potential, with the control parameters being the exponent (of the inverse power law) for the former and the number density for the latter. Both systems are shown to exhibit the characteristic features of the unjamming transition, among which are the vanishing of the shear-to-bulk modulus ratio and the emergence of an excess of low-frequency vibrational modes. We establish a relation between the pressure-to-bulk modulus ratio and the distance to unjamming in each of our model systems. This allows us to predict the dependence of other key observables on the distance to unjamming. Our results provide the means for a quantitative estimation of the proximity of generic glass-forming models to the unjamming transition in the absence of a clear-cut definition of the coordination number and highlight the general irrelevance of nonaffine contributions to the bulk modulus.

Transformable topological mechanical metamaterials

Mechanical metamaterials are engineered materials whose structures give them novel mechanical properties, including negative Poisson's ratios, negative compressibilities and phononic bandgaps. Of particular interest are systems near the point of mechanical instability, which recently have been shown to distribute force and motion in robust ways determined by a nontrivial topological state. Here we discuss the classification of and propose a design principle for mechanical metamaterials that can be easily and reversibly transformed between states with dramatically different mechanical and acoustic properties via a soft strain. Remarkably, despite the low energetic cost of this transition, quantities such as the edge stiffness and speed of sound can change by orders of magnitude. We show that the existence and form of a soft deformation directly determines floppy edge modes and phonon dispersion. Finally, we generalize the soft strain to generate domain structures that allow further tuning of the material.

Here Rocklin et al. propose a design principle using operations that cost little energy and realize mechanical metamaterials that can be easily and reversibly transformed between states with different mechanical and acoustic properties.

Signal transmissibility in marginal granular materials

We examine the 'transmissibility' of a simulated two-dimensional pack of frictionless disks formed by confining dilute disks in a shrinking, periodic box to the point of mechanical stability. Two opposite boundaries are then removed, thus allowing a set of free motions. Small free displacements on one boundary then induce proportional displacements on the opposite boundary. Transmissibility is the ability to distinguish different perturbations by their distant responses. We assess transmissibility by successively identifying free orthonormal modes of motion that have the smallest distant responses. The last modes to be identified in this 'pessimistic' basis are the most transmissive. The transmitted amplitudes of these most transmissive modes fall off exponentially with mode number. Similar exponential falloff is seen in a simple elastic medium, though the responsible modes differ greatly in structure in the two systems. Thus the marginal pack's transmissibility is qualitatively similar to that of a simple elastic medium. We compare our results with recent findings based on the projection of the space of free motion onto interior sites.

The role of rigidity in controlling material failure

We investigate how material rigidity acts as a key control parameter for the failure of solids under stress. In both experiments and simulations, we demonstrate that material failure can be continuously tuned by varying the underlying rigidity of the material while holding the amount of disorder constant. As the rigidity transition is approached, failure due to the application of uniaxial stress evolves from brittle cracking to system-spanning diffuse breaking. This evolution in failure behavior can be parameterized by the width of the crack. As a system becomes more and more floppy, this crack width increases until it saturates at the system size. Thus, the spatial extent of the failure zone can be used as a direct probe for material rigidity.

Zero-frequency and slow elastic modes in phononic monolayer granular membranes

We theoretically study the dispersion properties of elastic waves in hexagonal and honeycomb monolayer granular membranes with either out-of-plane or in-plane particle motion. The particles interact predominantly via normal and transverse contact rigidities. When rotational degrees of freedom are taken into account, the bending and torsional rigidities of the intergrain contacts can control some of the phononic modes. The existence of zero-frequency modes, zero-group-velocity modes and their transformation into slow propagating phononic modes due to weak bending and torsional intergrain interactions are investigated. We also study the formation and manipulation of Dirac cones and multiple degenerated modes. This could motivate variety of potential applications in elastic waves control by manipulating the contact rigidities in granular phononic crystals.

Nonlinear elasticity of disordered fiber networks

Disordered biopolymer gels have striking mechanical properties including strong nonlinearities. In the case of athermal gels (such as collagen-I) the nonlinearity has long been associated with a crossover from a bending dominated to a stretching dominated regime of elasticity. The physics of this crossover is related to the existence of a central-force isostatic point and to the fact that for most gels the bending modulus is small. This crossover induces scaling behavior for the elastic moduli. In particular, for linear elasticity such a scaling law has been demonstrated [Broedersz et al. Nat. Phys., 2011 7, 983]. In this work we generalize the scaling to the nonlinear regime with a two-parameter scaling law involving three critical exponents. We test the scaling law numerically for two disordered lattice models, and find a good scaling collapse for the shear modulus in both the linear and nonlinear regimes. We compute all the critical exponents for the two lattice models and discuss the applicability of our results to real systems.

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.

Finite-temperature buckling of an extensible rod

Thermal fluctuations can play an important role in the buckling of elastic objects at small scales, such as polymers or nanotubes. In this paper, we study the finite-temperature buckling transition of an extensible rod by analyzing fluctuation corrections to the elasticity of the rod. We find that, in both two and three dimensions, thermal fluctuations delay the buckling transition, and near the transition, there is a critical regime in which fluctuations are prominent and make a contribution to the effective force that is of order √T. We verify our theoretical prediction of the phase diagram with Monte Carlo simulations.

Lattice mechanics of origami tessellations

Origami-based design holds promise for developing materials whose mechanical properties are tuned by crease patterns introduced to thin sheets. Although there have been heuristic developments in constructing patterns with desirable qualities, the bridge between origami and physics has yet to be fully developed. To truly consider origami structures as a class of materials, methods akin to solid mechanics need to be developed to understand their long-wavelength behavior. We introduce here a lattice theory for examining the mechanics of origami tessellations in terms of the topology of their crease pattern and the relationship between the folds at each vertex. This formulation provides a general method for associating mechanical properties with periodic folded structures and allows for a concrete connection between more conventional materials and the mechanical metamaterials constructed using origami-based design.

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 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.

Force distribution affects vibrational properties in hard-sphere glasses

We theoretically and numerically study the elastic properties of hard-sphere glasses and provide a real-space description of their mechanical stability. In contrast to repulsive particles at zero temperature, we argue that the presence of certain pairs of particles interacting with a small force f soften elastic properties. This softening affects the exponents characterizing elasticity at high pressure, leading to experimentally testable predictions. Denoting P(f) ~ f(θ(e)), the force distribution of such pairs and ϕ(c) the packing fraction at which pressure diverges, we predict that (i) the density of states has a low-frequency peak at a scale ω*, rising up to it as D(ω) ~ ω(2+a), and decaying above ω* as D(ω) ~ ω(-a) where a = (1 - θ(e))/(3 + θ(e)) and ω is the frequency, (ii) shear modulus and mean-squared displacement are inversely proportional with ⟨δR²⟩ ~ 1/μ ~ (ϕ(c) - ϕ)(κ), where κ = 2 - 2/(3 + θ(e)), and (iii) continuum elasticity breaks down on a scale ℓ(c) ~ 1/√(δz) ~ (ϕ(c) - ϕ)(-b), where b = (1 + θ(e))/(6 + 2θ(e)) and δz = z - 2d, where z is the coordination and d the spatial dimension. We numerically test (i) and provide data supporting that θ(e) ≈ 0.41 in our bidisperse system, independently of system preparation in two and three dimensions, leading to κ ≈ 1.41, a ≈ 0.17, and b ≈ 0.21. Our results for the mean-square displacement are consistent with a recent exact replica computation for d = ∞, whereas some observations differ, as rationalized by the present approach.

Penrose tilings as jammed solids

Penrose tilings form lattices, exhibiting fivefold symmetry and isotropic elasticity, with inhomogeneous coordination much like that of the force networks in jammed systems. Under periodic boundary conditions, their average coordination is exactly four. We study the elastic and vibrational properties of rational approximants to these lattices as a function of unit-cell size N(S) and find that they have of order sqrt[N(S)] zero modes and states of self-stress and yet all their elastic moduli vanish. In their generic form, obtained by randomizing site positions, their elastic and vibrational properties are similar to those of particulate systems at jamming with a nonzero bulk modulus, vanishing shear modulus, and a flat density of states.

Effects of coordination and pressure on sound attenuation, boson peak and elasticity in amorphous solids

Connectedness and applied stress strongly affect elasticity in solids. In various amorphous materials, mechanical stability can be lost either by reducing connectedness or by increasing pressure. We present an effective medium theory of elasticity that extends previous approaches by incorporating the effect of compression, of amplitude e, allowing one to describe quantitative features of sound propagation, transport, the boson peak, and elastic moduli near the elastic instability occurring at a compression ec. The theory disentangles several frequencies characterizing the vibrational spectrum: the onset frequency where strongly-scattered modes appear in the vibrational spectrum, the pressure-independent frequency ω* where the density of states displays a plateau, the boson peak frequency ωBP found to scale as , and the Ioffe-Regel frequency ωIR where scattering length and wavelength become equal. We predict that sound attenuation crosses over from ω(4) to ω(2) behaviour at ω0, consistent with observations in glasses. We predict that a frequency-dependent length scale ls(ω) and speed of sound ν(ω) characterize vibrational modes, and could be extracted from scattering data. One key result is the prediction of a flat diffusivity above ω0, in agreement with previously unexplained observations. We find that the shear modulus does not vanish at the elastic instability, but drops by a factor of 2. We check our predictions in packings of soft particles and study the case of covalent networks and silica, for which we predict ωIR ≈ ωBP. Overall, our approach unifies sound attenuation, transport and length scales entering elasticity in a single framework where disorder is not the main parameter controlling the boson peak, in agreement with observations. This framework leads to a phase diagram where various glasses can be placed, connecting microscopic structure to vibrational properties.

Rheology of semiflexible bundle networks with transient linkers

We present a theoretical and computational analysis of the rheology of networks made up of bundles of semiflexible filaments bound by transient cross-linkers. Such systems are ubiquitous in the cytoskeleton and can be formed in vitro using filamentous actin and various cross-linkers. We find that their high-frequency rheology is characterized by a scaling behavior that is quite distinct from that of networks of the well-studied single semiflexible filaments. This regime can be understood theoretically in terms of a length-scale-dependent bending modulus for bundles. Next, we observe new dissipative dynamics associated with the shear-induced disruption of the network at intermediate frequencies. Finally, at low frequencies, we encounter a region of non-Newtonian rheology characterized by power-law scaling. This regime is dominated by bundle dissolution and large-scale rearrangements of the network driven by equilibrium thermal fluctuations.

Physics in ordered and disordered colloidal matter composed of poly(N-isopropylacrylamide) microgel particles

This review collects and describes experiments that employ colloidal suspensions to probe physics in ordered and disordered solids and related complex fluids. The unifying feature of this body of work is its clever usage of poly(N-isopropylacrylamide) (PNIPAM) microgel particles. These temperature-sensitive colloidal particles provide experimenters with a 'knob' for in situ control of particle size, particle interaction and particle packing fraction that, in turn, influence the structural and dynamical behavior of the complex fluids and solids. A brief summary of PNIPAM particle synthesis and properties is given, followed by a synopsis of current activity in the field. The latter discussion describes a variety of soft matter investigations including those that explore formation and melting of crystals and clusters, and those that probe structure, rearrangement and rheology of disordered (jammed/glassy) and partially ordered matter. The review, therefore, provides a snapshot of a broad range of physics phenomenology which benefits from the unique properties of responsive microgel particles.

Shear shocks in fragile networks

A minimal model for studying the mechanical properties of amorphous solids is a disordered network of point masses connected by unbreakable springs. At a critical value of its mean connectivity, such a network becomes fragile: it undergoes a rigidity transition signaled by a vanishing shear modulus and transverse sound speed. We investigate analytically and numerically the linear and nonlinear visco-elastic response of these fragile solids by probing how shear fronts propagate through them. Our approach, which we tentatively label shear front rheology, provides an alternative route to standard oscillatory rheology. In the linear regime, we observe at late times a diffusive broadening of the fronts controlled by an effective shear viscosity that diverges at the critical point. No matter how small the microscopic coefficient of dissipation, strongly disordered networks behave as if they were overdamped because energy is irreversibly leaked into diverging nonaffine fluctuations. Close to the transition, the regime of linear response becomes vanishingly small: the tiniest shear strains generate strongly nonlinear shear shock waves qualitatively different from their compressional counterparts in granular media. The inherent nonlinearities trigger an energy cascade from low to high frequency components that keep the network away from attaining the quasi-static limit. This mechanism, reminiscent of acoustic turbulence, causes a superdiffusive broadening of the shock width.

Convolution theory for dynamic systems: a bottom-up approach to the viscoelasticity of polymeric networks

Biological materials such as the cytoskeleton are characterized by remarkable viscoelastic properties and therefore represent the subject of numerous micro- and macrorheological experimental studies. By generalizing the previously introduced dynamic convolution theory (DCT) to two dimensions, we devise a bottom-up approach for the viscoelastic properties of extended, crosslinked semiflexible polymer networks. Brownian dynamics (BD) simulations serve to determine the dynamic linear self- and cross-response properties of isolated semiflexible polymers to externally applied forces and torques; these response functions are used as input to the DCT. For a given network topology, the frequency-dependent response of the network subject to a given external force/torque distribution is calculated via the DCT allowing to resolve both micro- and macrorheological properties of the networks. A mapping on continuum viscoelastic theory yields the corresponding viscoelastic bulk moduli. Special attention is drawn to the flexibility of crosslinkers, which couple angular degrees of freedom at the network nodes and which are found to sensitively affect the resulting rheological properties of the polymeric meshwork.

Nanoindentation of 35 virus capsids in a molecular model: relating mechanical properties to structure

A coarse-grained model is used to study the mechanical response of 35 virus capsids of symmetries T = 1, T = 2, T = 3, pseudo T = 3, T = 4, and T = 7. The model is based on the native structure of the proteins that constitute the capsids and is described in terms of the C atoms associated with each amino acid. The number of these atoms ranges between 8 460 (for SPMV – satellite panicum mosaic virus) and 135 780 (for NBV – nudaureli virus). Nanoindentation by a broad AFM tip is modeled as compression between two planes: either both flat or one flat and one curved. Plots of the compressive force versus plate separation show a variety of behaviors, but in each case there is an elastic region which extends to a characteristic force . Crossing results in a drop in the force and irreversible damage. Across the 35 capsids studied, both and the elastic stiffness are observed to vary by a factor of 20. The changes in mechanical properties do not correlate simply with virus size or symmetry. There is a strong connection to the mean coordination number , defined as the mean number of interactions to neighboring amino acids. The Young's modulus for thin shell capsids rises roughly quadratically with , where 6 is the minimum coordination for elastic stability in three dimensions.

Effective-medium theory of a filamentous triangular lattice

We present an effective-medium theory that includes bending as well as stretching forces, and we use it to calculate the mechanical response of a diluted filamentous triangular lattice. In this lattice, bonds are central-force springs, and there are bending forces between neighboring bonds on the same filament. We investigate the diluted lattice in which each bond is present with a probability p. We find a rigidity threshold p(b) which has the same value for all positive bending rigidity and a crossover characterizing bending, stretching, and bend-stretch coupled elastic regimes controlled by the central-force rigidity percolation point at p(CF)=/~2/3 of the lattice when fiber bending rigidity vanishes.

Actively stressed marginal networks

We study the effects of motor-generated stresses in disordered three-dimensional fiber networks using a combination of a mean-field theory, scaling analysis, and a computational model. We find that motor activity controls the elasticity in an anomalous fashion close to the point of marginal stability by coupling to critical network fluctuations. We also show that motor stresses can stabilize initially floppy networks, extending the range of critical behavior to a broad regime of network connectivities below the marginal point. Away from this regime, or at high stress, motors give rise to a linear increase in stiffness with stress. Finally, we demonstrate that our results are captured by a simple, constitutive scaling relation highlighting the important role of nonaffine strain fluctuations as a susceptibility to motor stress.

Surface phonons, elastic response, and conformal invariance in twisted kagome lattices

Model lattices consisting of balls connected by central-force springs provide much of our understanding of mechanical response and phonon structure of real materials. Their stability depends critically on their coordination number z. d-dimensional lattices with z = 2d are at the threshold of mechanical stability and are isostatic. Lattices with z < 2d exhibit zero-frequency "floppy" modes that provide avenues for lattice collapse. The physics of systems as diverse as architectural structures, network glasses, randomly packed spheres, and biopolymer networks is strongly influenced by a nearby isostatic lattice. We explore elasticity and phonons of a special class of two-dimensional isostatic lattices constructed by distorting the kagome lattice. We show that the phonon structure of these lattices, characterized by vanishing bulk moduli and thus negative Poisson ratios (equivalently, auxetic elasticity), depends sensitively on boundary conditions and on the nature of the kagome distortions. We construct lattices that under free boundary conditions exhibit surface floppy modes only or a combination of both surface and bulk floppy modes; and we show that bulk floppy modes present under free boundary conditions are also present under periodic boundary conditions but that surface modes are not. In the long-wavelength limit, the elastic theory of all these lattices is a conformally invariant field theory with holographic properties (characteristics of the bulk are encoded on the sample boundary), and the surface waves are Rayleigh waves. We discuss our results in relation to recent work on jammed systems. Our results highlight the importance of network architecture in determining floppy-mode structure.

A unified framework for non-brownian suspension flows and soft amorphous solids

While the rheology of non-brownian suspensions in the dilute regime is well understood, their behavior in the dense limit remains mystifying. As the packing fraction of particles increases, particle motion becomes more collective, leading to a growing length scale and scaling properties in the rheology as the material approaches the jamming transition. There is no accepted microscopic description of this phenomenon. However, in recent years it has been understood that the elasticity of simple amorphous solids is governed by a critical point, the unjamming transition where the pressure vanishes, and where elastic properties display scaling and a diverging length scale. The correspondence between these two transitions is at present unclear. Here we show that for a simple model of dense flow, which we argue captures the essential physics near the jamming threshold, a formal analogy can be made between the rheology of the flow and the elasticity of simple networks. This analogy leads to a new conceptual framework to relate microscopic structure to rheology. It enables us to define and compute numerically normal modes and a density of states. We find striking similarities between the density of states in flow, and that of amorphous solids near unjamming: both display a plateau above some frequency scale ω(∗) ∼ |z(c) - z|, where z is the coordination of the network of particle in contact, z(c) = 2D where D is the spatial dimension. However, a spectacular difference appears: the density of states in flow displays a single mode at another frequency scale ω(min) ≪ ω(∗) governing the divergence of the viscosity.

Tuning jammed frictionless disk packings from isostatic to hyperstatic

We perform extensive computational studies of two-dimensional static bidisperse disk packings using two distinct packing-generation protocols. The first involves thermally quenching equilibrated liquid configurations to zero temperature over a range of thermal quench rates r and initial packing fractions followed by compression and decompression in small steps to reach packing fractions φ(J) at jamming onset. For the second, we seed the system with initial configurations that promote micro- and macrophase-separated packings followed by compression and decompression to φ(J). Using these protocols, we generate more than 10(4) static packings over a wide range of packing fraction, contact number, and compositional and positional order. We find that disordered, isostatic packings exist over a finite range of packing fractions in the large-system limit. In agreement with previous calculations, the most dilute mechanically stable packings with φ min ≈ 0.84 are obtained for r > r*, where r* is the rate above which φ(J) is insensitive to rate. We further compare the structural and mechanical properties of isostatic versus hyperstatic packings. The structural characterizations include the contact number, several order parameters, and mixing ratios of the large and small particles. We find that the isostatic packings are positionally and compositionally disordered (with only small changes in a number of order parameters), whereas bond-orientational and compositional order increase strongly with contact number for hyperstatic packings. In addition, we calculate the static shear modulus and normal mode frequencies (in the harmonic approximation) of the static packings to understand the extent to which the mechanical properties of disordered, isostatic packings differ from partially ordered packings. We find that the mechanical properties of the packings change continuously as the contact number increases from isostatic to hyperstatic.

Coherent potential approximation of random nearly isostatic kagome lattice

The kagome lattice has coordination number 4, and it is mechanically isostatic when nearest-neighbor sites are connected by central-force springs. A lattice of N sites has O(√N) zero-frequency floppy modes that convert to finite-frequency anomalous modes when next-nearest-neighbor (NNN) springs are added. We use the coherent potential approximation to study the mode structure and mechanical properties of the kagome lattice in which NNN springs with spring constant κ are added with probability P=Δz/4, where Δz=z-4 and z is the average coordination number. The effective medium static NNN spring constant κ(m) scales as P(2) for P≪κ and as P for P≫κ, yielding a frequency scale ω*~Δz and a length scale l*~(Δz)(-1). To a very good approximation at small nonzero frequency, κ(m)(P,ω)/κ(m)(P,0) is a scaling function of ω/ω*. The Ioffe-Regel limit beyond which plane-wave states become ill-defined is reached at a frequency of order ω*.