No Arabic abstract
The Kremer-Grest (KG) model is a standard for studying generic polymer properties. Here we have equilibrated KG melts up to and beyond $200$ entanglements per chain for varying chain stiffness. We present methods for estimating the Kuhn length corrected for incompressibility effects, for estimating the entanglement length corrected for chain stiffness, for estimating bead frictions and Kuhn times taking into account entanglement effects. These are the key parameters for enabling quantitative, accurate, and parameter free comparisons between theory, experiment and simulations of KG polymer models with varying stiffness. We demonstrate this for the mean-square monomer displacements in moderately to highly entangled melts as well as for the shear relaxation modulus for unentangled melts, which are found to be in excellent agreement with the predictions from standard theories of polymer dynamics.
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 Kremer-Grest (KG) bead-spring model is a near standard in Molecular Dynamic simulations of generic polymer properties. It owes its popularity to its computational efficiency, rather than its ability to represent specific polymer species and conditions. Here we investigate how to adapt the model to match the universal properties of a wide range of chemical polymers species. For this purpose we vary a single parameter originally introduced by Faller and Muller-Plathe, the chain stiffness. Examples include polystyrene, polyethylene, polypropylene, cis-polyisoprene, polydimethylsiloxane, polyethyleneoxide and styrene-butadiene rubber. We do this by matching the number of Kuhn segments per chain and the number of Kuhn segments per cubic Kuhn volume for the polymer species and for the Kremer-Grest model. We also derive mapping relations for converting KG model units back to physical units, in particular we obtain the entanglement time for the KG model as function of stiffness allowing for a time mapping. To test these relations, we generate large equilibrated well entangled polymer melts, and measure the entanglement moduli using a static primitive-path analysis of the entangled melt structure as well as by simulations of step-strain deformation of the model melts. The obtained moduli for our model polymer melts are in good agreement with the experimentally expected moduli.
We simulate a relaxation process of non-brownian particles in a sheared viscous medium; the small shear strain is initially applied to a system, which then undergoes relaxation. The relaxation time and the correlation length are estimated as functions of density, which algebraically diverge at the jamming density. This implies that the relaxation time can be scaled by the correlation length using the dynamic critical exponent, which is estimated as 4.6(2). It is also found that shear stress undergoes power-law decay at the jamming density, which is reminiscent of critical slowing down.
Self-avoiding polymers in two-dimensional ($d=2$) melts are known to adopt compact configurations of typical size $R(N) sim N^{1/d}$ with $N$ being the chain length. Using molecular dynamics simulations we show that the irregular shapes of these chains are characterized by a perimeter length $L(N) sim R(N)^{dpm}$ of fractal dimension $dpm = d-Theta_2 =5/4$ with $Theta_2=3/4$ being a well-known contact exponent. Due to the self-similar structure of the chains, compactness and perimeter fractality repeat for subchains of all arc-lengths $s$ down to a few monomers. The Kratky representation of the intramolecular form factor $F(q)$ reveals a strong non-monotonous behavior with $q^2F(q) sim 1/(qN^{1/d})^{Theta_2}$ in the intermediate regime of the wavevector $q$. Measuring the scattering of labeled subchains %($s F(q) sim L(s)$) the form factor may allow to test our predictions in real experiments.