Semiflexible polymer dynamics with a bead-spring model

Abstract

We study the dynamical properties of semiflexible polymers with a recently introduced bead-spring model. We focus on double-stranded DNA (dsDNA). The two parameters of the model, T* and ν, are chosen to match its experimental force-extension curve. In comparison to its groundstate value, the bead-spring Hamiltonian is approximated in the first order by the Hessian that is quadratic in the bead positions. The eigenmodes 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. At the Hessian approximation of the Hamiltonian, the polymer dynamics is linear. Using the longitudinal and transverse eigenmodes, for the linearized problem, 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 the sum over the contributions from the modes—the so-called ‘mode sums’. We also perform simulations with the full dynamics of the model. The simulations yield numerical values of the correlations functions (i–iii) that agree very well with the analytical expressions for the linearized dynamics. This does not however mean that the nonlinearities are not present. In fact, 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 t7/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.