No Arabic abstract
We present a theoretical framework for the linear and nonlinear visco-elastic properties of reversibly crosslinked networks of semiflexible polymers. In contrast to affine models where network strain couples to the polymer end-to-end distance, in our model strain rather serves to locally distort the network structure. This induces bending modes in the polymer filaments, the properties of wich are slaved to the surrounding network structure. Specifically, we investigate the frequency-dependent linear rheology, in particular in combination with crosslink binding/unbinding processes. We also develop schematic extensions to describe the nonlinear response during creep measurements as well as during constant-strainrate ramps.
We have developed a new technique to measure viscoelasticity in soft materials such as polymer solutions, by monitoring thermal fluctuations of embedded probe particles using laser interferometry in a microscope. Interferometry allows us to obtain power spectra of fluctuating beads from 0.1 Hz to 20 kHz, and with sub-nanometer spatial resolution. Using linear response theory, we determined the frequency-dependent loss and storage shear moduli up to frequencies on the order of a kHz. Our technique measures local values of the viscoelastic response, without actively straining the system, and is especially suited to soft biopolymer networks. We studied semiflexible F-actin solutions and, as a control, flexible polyacrylamide (PAAm) gels, the latter close to their gelation threshold. With small particles, we could probe the transition from macroscopic viscoelasticity to more complex microscopic dynamics. In the macroscopic limit we find shear moduli at 0.1 Hz of G=0.11 +/- 0.03 Pa and 0.17 +/- 0.07 Pa for 1 and 2 mg/ml actin solutions, close to the onset of the elastic plateau, and scaling behavior consistent with G(omega) as omega^(3/4) at higher frequencies. For polyacrylamide we measured plateau moduli of 2.0, 24, 100 and 280 Pa for crosslinked gels of 2, 2.5, 3 and 5% concentration (weight/volume) respectively, in agreement to within a factor of two with values obtained from conventional rheology. We also found evidence for scaling of G(omega) as omega^(1/2), consistent with the predictions of the Rouse model for flexible polymers.
We present a method to generate realistic, three-dimensional networks of crosslinked semiflexible polymers. The free energy of these networks is obtained from the force-extension characteristics of the individual polymers and their persistent directionality through the crosslinks. A Monte Carlo scheme is employed to obtain isotropic, homogeneous networks that minimize the free energy, and for which all of the relevant parameters can be varied: the persistence length, the contour length as well as the crosslinking length may be chosen at will. We also provide an initial survey of the mechanical properties of our networks subjected to shear strains, showing them to display the expected non-linear stiffening behavior. Also, a key role for non-affinity and its relation to order in the network is uncovered.
Reversible crosslinking is a design paradigm for polymeric materials, wherein they are microscopically reinforced with chemical species that form transient crosslinks between the polymer chains. Besides the potential for self-healing, recent experimental work suggests that freely diffusing reversible crosslinks in polymer networks, such as gels, can enhance the toughness of the material without substantial change in elasticity. This presents the opportunity for making highly elastic materials that can be strained to a large extent before rupturing. Here, we employ Gaussian chain theory, molecular simulation, and polymer self-consistent field theory for networks to construct an equilibrium picture for how reversible crosslinks can toughen a polymer network without affecting its linear elasticity. Maximisation of polymer entropy drives the reversible crosslinks to bind preferentially near the permanent crosslinks in the network, leading to local molecular reinforcement without significant alteration of the network topology. In equilibrium conditions, permanent crosslinks share effectively the load with neighbouring reversible crosslinks, forming multi-functional crosslink points. The network is thereby globally toughened, while the linear elasticity is left largely unaltered. Practical guidelines are proposed to optimise this design in experiment, along with a discussion of key kinetic and timescale considerations.
Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydrodynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymers fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5%$ of several key observables, the model allows for single Langevin time steps of $approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.
Motivated by the structure of networks of cross-linked cytoskeletal biopolymers, we study the orientationally ordered phases in two-dimensional networks of randomly cross-linked semiflexible polymers. We consider permanent cross-links which prescribe a finite angle and treat them as quenched disorder in a semi-microscopic replica field theory. Starting from a fluid of un-cross-linked polymers and small polymer clusters (sol) and increasing the cross-link density, a continuous gelation transition occurs. In the resulting gel, the semiflexible chains either display long range orientational order or are frozen in random directions depending on the value of the crossing angle, the crosslink concentration and the stiffness of the polymers. A crossing angle $thetasim 2pi/M$ leads to long range $M$-fold orientational order, e.g., hexatic or tetratic for $theta=60^{circ}$ or $90^{circ}$, respectively. The transition is discontinuous and the critical cross-link density depends on the bending stiffness of the polymers and the cross-link geometry: the higher the stiffness and the lower $M$, the lower the critical number of cross-links. In between the sol and the long range ordered state, we always observe a gel which is a statistically isotropic amorphous solid (SIAS) with random positional and random orientational localization of the participating polymers.