No Arabic abstract
Semiflexible polymers in concentrated lyotropic solution are studied within a bead-spring model by molecular dynamics simulations, focusing on the emergence of a smectic A phase and its properties. We systematically vary the density of the monomeric units for several contour lengths that are taken smaller than the chain persistence length. The difficulties concerning the equilibration of such systems and the choice of appropriate ensemble (constant volume versus constant pressure, where all three linear dimensions of the simulation box can fluctuate independently) are carefully discussed. Using HOOMD-blue on graphics processing units, systems containing more than a million monomeric units are accessible, making it possible to distinguish the order of the phase transitions that occur. While in this model the nematic-smectic transition is continuous, the transition from the smectic phase to a related crystalline structure with true three-dimensional long-range order is clearly of first order. Further, both orientational and positional correlations of monomeric units are studied as well as the order parameters characterizing the nematic, smectic A, and crystalline phases. The analogy between smectic order and one-dimensional harmonic crystals with respect to the behavior of the structure factor is also explored. Finally, the results are put in perspective with pertinent theoretical predictions and possible experiments.
Using Molecular Dynamics simulations, we study the force-induced detachment of a coarse-grained model polymer chain from an adhesive substrate. One of the chain ends is thereby pulled at constant speed off the attractive substrate and the resulting saw-tooth profile of the measured mean force $< f >$ vs height $D$ of the end-segment over the plane is analyzed for a broad variety of parameters. It is shown that the observed characteristic oscillations in the $< f >$-$D$ profile depend on the bending and not on the torsional stiffness of the detached chains. Allowing for the presence of hydrodynamic interactions (HI) in a setup with explicit solvent and DPD-thermostat, rather than the case of Langevin thermostat, one finds that HI have little effect on the $< f >$-$D$ profile. Also the change of substrate affinity with respect to the solvent from solvophilic to solvophobic is found to play negligible role in the desorption process. In contrast, a changing ratio $epsilon_s^A / epsilon_s^B$ of the binding energies of $A$- and $B$-segments in the detachment of an $AB$-copolymer from adhesive surface strongly changes the $< f >$-$D$ profile whereby the $B$-spikes vanish when $epsilon_s^A / epsilon_s^B < 0.15$. Eventually, performing an atomistic simulation of a (bio)-polymer {it polyglycine}, we demonstrate that the simulation results, derived from our coarse-grained model, comply favorably with those from the all-atom simulation.
We study the dynamical properties of semiflexible polymers with a recently introduced bead-spring model. We focus on double-stranded DNA. The two parameters of the model, $T^*$ and $ u$, are chosen to match its experimental force-extension curve. The bead-spring Hamiltonian is approximated in the first order by the Hessian that is quadratic in the bead positions. The eigenmodels of the Hessian provide the longitudinal (stretching) and transverse (bending) eigenmodes of the polymer, and the corresponding eigenvalues match well with the established phenomenology of semiflexible polymers. Using the longitudinal and transverse eigenmodes, we obtain analytical expressions of (i) the autocorrelation function of the end-to-end vector, (ii) the autocorrelation function of a bond (i.e., a spring, or a tangent) vector at the middle of the chain, and (iii) the mean-square displacement of a tagged bead in the middle of the chain, as sum over the contributions from the modes. We also perform simulations with the full dynamics of the model. The simulations yield numerical values of the correlation functions (i-iii) that agree very well with the analytical expressions for the linearized dynamics. We also study the mean-square displacement of the longitudinal component of the end-to-end vector that showcases strong nonlinear effects in the polymer dynamics, and we identify at least an effective $t^{7/8}$ power-law regime in its time-dependence. Nevertheless, in comparison to the full mean-square displacement of the end-to-end vector the nonlinear effects remain small at all times --- it is in this sense we state that our results demonstrate that the linearized dynamics suffices for dsDNA fragments that are shorter than or comparable to the persistence length. Our results are consistent with those of the wormlike chain (WLC) model, the commonly used descriptive tool of semiflexible polymers.
The cytoskeleton of eukaryotic cells provides mechanical support and governs intracellular transport. These functions rely on the complex mechanical properties of networks of semiflexible protein filaments. Recent theoretical interest has focused on mesoscopic properties of such networks and especially on the effect of local, non-affine bending deformations on mechanics. Here, we study the impact of local network deformations on the scale-dependent mobility of probe particles in entangled networks of semiflexible actin filaments by high-bandwidth microrheology. We find that micron-sized particles in these networks experience two opposing non-continuum elastic effects: entropic depletion reduces the effective network rigidity, while local non-affine deformations of the network substantially enhance the rigidity at low frequencies. We show that a simple model of lateral bending of filaments embedded in a viscoelastic background leads to a scaling regime for the apparent elastic modulus G(omega) sim omega^{9/16}, closely matching the experiments. These results provide quantitative evidence for how different a semiflexible polymer network can feel for small objects, and they demonstrate how non-affine bending deformations can be dominant for the mobility of vesicles and organelles in the cell.
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.
We develop a theory of Smectic A - Smectic C phase transition with anomalously weak smectic layer contraction. We construct a phenomenological description of this transition by generalizing the Chen-Lubensky model. Using a mean-field molecular model, we demonstrate that a relatively simple interaction potential suffices to describe the transition. The theoretical results are in excellent agreement with experimental data.