Effective medium approach for stiff polymer networks with flexible cross-links

Recent experiments have demonstrated that the nonlinear elasticity of in vitro networks of the biopolymer actin is dramatically altered in the presence of a flexible cross-linker such as the abundant cytoskeletal protein filamin. The basic principles of such networks remain poorly understood. Here we describe an effective medium theory of flexibly cross-linked stiff polymer networks. We argue that the response of the cross-links can be fully attributed to entropic stiffening, while softening due to domain unfolding can be ignored. The network is modeled as a collection of randomly oriented rods connected by flexible cross-links to an elastic continuum. This effective medium is treated in a linear elastic limit as well as in a more general framework, in which the medium self-consistently represents the nonlinear network behavior. This model predicts that the nonlinear elastic response sets in at strains proportional to cross-linker length and inversely proportional to filament length. Furthermore, we find that the differential modulus scales linearly with the stress in the stiffening regime. These results are in excellent agreement with bulk rheology data.Comment: 12 pages, 8 figure

The bend stiffness of S-DNA

We formulate and solve a two-state model for the elasticity of nicked, double-stranded DNA that borrows features from both the Worm Like Chain and the Bragg--Zimm model. Our model is computationally simple, and gives an excellent fit to recent experimental data through the entire overstretching transition. The fit gives the first value for the bending stiffness of the overstretched state as about 10 nm*kbt, a value quite different from either B-form or single-stranded DNA.Comment: 7 pages, 1 figur

Nonlinear elasticity of composite networks of stiff biopolymers with flexible linkers

Motivated by recent experiments showing nonlinear elasticity of in vitro networks of the biopolymer actin cross-linked with filamin, we present an effective medium theory of flexibly cross-linked stiff polymer networks. We model such networks by randomly oriented elastic rods connected by flexible connectors to a surrounding elastic continuum, which self-consistently represents the behavior of the rest of the network. This model yields a crossover from a linear elastic regime to a highly nonlinear elastic regime that stiffens in a way quantitatively consistent with experiment.Comment: 4 pages, 3 figure

Reply to "Comment on 'Theory of high-force DNA stretching and overstretching'"

In his Comment to an earlier paper [Phys. Rev. E 67, 051906 (2003)] Lam points out an error in Eq. (20) of the original paper. Here we show that use of the corrected expression produces results very similar to those presented in our original paper, so our qualitative conclusions are unchanged

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.

Evaluation of SIR-A space radar for geologic interpretation: United States, Panama, Colombia, and New Guinea

Comparisons between LANDSAT MSS imagery, and aircraft and space radar imagery from different geologic environments in the United States, Panama, Colombia, and New Guinea demonstrate the interdependence of radar system geometry and terrain configuration for optimum retrieval of geologic information. Illustrations suggest that in the case of space radars (SIR-A in particular), the ability to acquire multiple look-angle/look-direction radar images of a given area is more valuable for landform mapping than further improvements in spatial resolution. Radar look-angle is concluded to be one of the most important system parameters of a space radar designed to be used for geologic reconnaissance mapping. The optimum set of system parameters must be determined for imaging different classes of landform features and tailoring the look-angle to local topography

Semiflexible Filamentous Composites

Inspired by the ubiquity of composite filamentous networks in nature we investigate models of biopolymer networks that consist of interconnected floppy and stiff filaments. Numerical simulations carried out in three dimensions allow us to explore the microscopic partitioning of stresses and strains between the stiff and floppy fractions c_s and c_f, and reveal a non-trivial relationship between the mechanical behavior and the relative fraction of stiff polymer: when there are few stiff polymers, non-percolated stiff ``inclusions`` are protected from large deformations by an encompassing floppy matrix, while at higher fractions of stiff material the stiff network is independently percolated and dominates the mechanical response.Comment: Phys. Rev. Lett, to appear (4 pages, 2 figures

A Microstructure-based Graph Neural Network for Accelerating Multiscale Simulations

Simulating the mechanical response of advanced materials can be done more accurately using concurrent multiscale models than with single-scale simulations. However, the computational costs stand in the way of the practical application of this approach. The costs originate from microscale Finite Element (FE) models that must be solved at every macroscopic integration point. A plethora of surrogate modeling strategies attempt to alleviate this cost by learning to predict macroscopic stresses from macroscopic strains, completely replacing the microscale models. In this work, we introduce an alternative surrogate modeling strategy that allows for keeping the multiscale nature of the problem, allowing it to be used interchangeably with an FE solver for any time step. Our surrogate provides all microscopic quantities, which are then homogenized to obtain macroscopic quantities of interest. We achieve this for an elasto-plastic material by predicting full-field microscopic strains using a graph neural network (GNN) while retaining the microscopic constitutive material model to obtain the stresses. This hybrid data-physics graph-based approach avoids the high dimensionality originating from predicting full-field responses while allowing non-locality to arise. By training the GNN on a variety of meshes, it learns to generalize to unseen meshes, allowing a single model to be used for a range of microstructures. The embedded microscopic constitutive model in the GNN implicitly tracks history-dependent variables and leads to improved accuracy. We demonstrate for several challenging scenarios that the surrogate can predict complex macroscopic stress-strain paths. As the computation time of our method scales favorably with the number of elements in the microstructure compared to the FE method, our method can significantly accelerate FE2 simulations

Theory of High-Force DNA Stretching and Overstretching

Single molecule experiments on single- and double stranded DNA have sparked a renewed interest in the force-extension of polymers. The extensible Freely Jointed Chain (FJC) model is frequently invoked to explain the observed behavior of single-stranded DNA. We demonstrate that this model does not satisfactorily describe recent high-force stretching data. We instead propose a model (the Discrete Persistent Chain, or ``DPC'') that borrows features from both the FJC and the Wormlike Chain, and show that it resembles the data more closely. We find that most of the high-force behavior previously attributed to stretch elasticity is really a feature of the corrected entropic elasticity; the true stretch compliance of single-stranded DNA is several times smaller than that found by previous authors. Next we elaborate our model to allow coexistence of two conformational states of DNA, each with its own stretch and bend elastic constants. Our model is computationally simple, and gives an excellent fit through the entire overstretching transition of nicked, double-stranded DNA. The fit gives the first values for the elastic constants of the stretched state. In particular we find the effective bend stiffness for DNA in this state to be about 10 nm*kbt, a value quite different from either B-form or single-stranded DNAComment: 33 pages, 11 figures. High-quality figures available upon reques