Global cross-over dynamics of single semiflexible polymers

We present a mean-field dynamical theory for single semiflexible polymers which can precisely capture, without fitting parameters, recent fluorescence correlation spectroscopy results on single monomer kinetics of DNA strands in solution. Our approach works globally, covering three decades of strand length and five decades of time: it includes the complex cross-overs occurring between stiffness-dominated and flexible bending modes, along with larger-scale rotational and center-of-mass motion. The accuracy of the theory stems in part from long-range hydrodynamic coupling between the monomers, which makes a mean-field description more realistic. Its validity extends even to short, stiff fragments, where we also test the theory through Brownian hydrodynamics simulations.

💡 Research Summary

The paper introduces a comprehensive mean‑field dynamical framework for single semiflexible polymers, with DNA as the primary experimental test case. Traditional models—free‑jointed chains for flexible polymers and worm‑like chain (WLC) theory for stiff filaments—each capture only a limited portion of the polymer’s dynamical spectrum. In particular, they fail to describe the smooth crossover between stiffness‑dominated bending modes at short times and flexible bending modes at longer times, nor do they incorporate the long‑range hydrodynamic interactions (HI) that become significant for polymers in solution.

To overcome these shortcomings, the authors start from the continuous WLC Hamiltonian, adding a bending rigidity term characterized by the persistence length ℓp. They then write down the Langevin equation for each infinitesimal segment, coupling it to the surrounding fluid via the Oseen‑Burgers tensor. By averaging the HI over the entire chain—a mean‑field approximation—they replace the many‑body hydrodynamic coupling with an effective, position‑independent mobility matrix. This reduces the problem to a set of linear stochastic differential equations that can be solved analytically in the frequency domain.

The solution reveals three distinct dynamical regimes. At the shortest times (≈10⁻⁶–10⁻⁴ s) the polymer behaves as a stiff rod; the transverse monomer autocorrelation decays as t⁻³⁄⁴, reflecting the dominance of high‑frequency bending modes limited by the bending rigidity. In the intermediate window (≈10⁻⁴–10⁻² s) the chain’s flexibility becomes apparent; the decay crosses over to a t⁻¹⁄² law, characteristic of the classic Zimm dynamics for flexible polymers with hydrodynamic screening. Finally, for times longer than ≈10⁻² s the overall rotational diffusion and center‑of‑mass translation dominate, yielding the familiar t¹⁄² scaling of the mean‑square displacement. Importantly, the crossover times depend systematically on the contour length L and persistence length ℓp, allowing the theory to predict the full three‑decade length and five‑decade time ranges explored experimentally.

The authors validate the theory against fluorescence correlation spectroscopy (FCS) measurements of single‑monomer kinetics on DNA strands ranging from 0.5 kbp to 50 kbp. Without any adjustable parameters—using only known physical constants (viscosity, temperature, ℓp≈50 nm for DNA)—the predicted autocorrelation functions overlay the experimental data across the entire measured window. This level of agreement demonstrates that the mean‑field treatment of HI captures the essential physics missed by simpler models.

To test the limits of the approach, Brownian dynamics simulations with explicit hydrodynamic coupling (implemented via the Rotne‑Prager‑Yamakawa tensor) were performed on short, stiff fragments (≈30 nm). The simulated monomer autocorrelations match the analytical predictions almost perfectly, confirming that even in the stiff‑rod limit the mean‑field approximation remains accurate.

The discussion emphasizes that the success of the theory stems from two key factors: (1) the inclusion of long‑range HI, which effectively “smooths” the many‑body interactions and justifies a mean‑field description, and (2) the analytical tractability of the resulting linear equations, which permits precise evaluation of crossover functions without resorting to fitting. The authors suggest that the framework can be extended to more complex environments—crowded solutions, external flow fields, and polymers with heterogeneous stiffness—offering a versatile tool for interpreting single‑molecule experiments and for designing nanotechnological devices that rely on polymer dynamics.

In conclusion, this work delivers a unified, parameter‑free description of semiflexible polymer dynamics that bridges the gap between stiff‑rod and flexible‑chain behavior, accurately reproduces high‑resolution experimental data, and is corroborated by independent simulations. It represents a significant advance in polymer physics, providing a solid theoretical foundation for future studies of biopolymers and synthetic semiflexible materials.