No Arabic abstract
We study the relaxation dynamics of a coarse-grained polymer chain at different degrees of stretching by both analytical means and numerical simulations. The macromolecule is modelled as a string of beads, connected by anharmonic springs, subject to a tensile force applied at the end monomer of the chain while the other end is fixed at the origin of coordinates. The impact of bond non-linearity on the relaxation dynamics of the polymer at different degrees of stretching is treated analytically within the Gaussian self-consistent approach (GSC) and then compared to simulation results derived from two different methods: Monte-Carlo (MC) and Molecular Dynamics (MD). At low and medium degrees of chain elongation we find good agreement between GSC predictions and the Monte-Carlo simulations. However, for strongly stretched chains the MD method, which takes into account inertial effects, reveals two important aspects of the nonlinear interaction between monomers: (i) a coupling and energy transfer between the damped, oscillatory normal modes of the chain, and (ii) the appearance of non-vanishing contributions of a continuum of frequencies around the characteristic modes in the power spectrum of the normal mode correlation functions.
Using analytical techniques and Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a narrow channel of width $R$ embedded in two dimensions, driven by a force proportional to the number of monomers in the channel. Such a setup mimics typical experimental situations in nano/micro-fluidics. During the the translocation process if the monomers in the channel can sufficiently quickly assume steady state motion, we observe the scaling $tausim N/F$ of the translocation time $tau$ with the driving force $F$ per bead and the number $N$ of monomers per chain. With smaller channel width $R$, steady state motion cannot be achieved, effecting a non-universal dependence of $tau$ on $N$ and $F$. From the simulations we also deduce the waiting time distributions under various conditions for the single segment passage through the channel entrance. For different chain lengths but the same driving force, the curves of the waiting time as a function of the translocation coordinate $s$ feature a maximum located at identical $s_{mathrm{max}}$, while with increasing the driving force or the channel width the value of $s_{mathrm{max}}$ decreases.
We investigate the crystallization mechanism of a single, flexible homopolymer chain with short range attractions. For a sufficiently narrow attractive well, the system undergoes a first-order like freezing transition from an expanded disordered coil to a compact crystalline state. Based on a maximum likelihood analysis of committor values computed for configurations obtained by Wang-Landau sampling, we construct a non-linear string reaction coordinate for the coil-to-crystal transition. In contrast to a linear reaction coordinate, the string reaction coordinate captures the effect of different degrees of freedom controlling different stages of the transition. Our analysis indicates that a combination of the energy and the global crystallinity parameter Q6 provide the most accurate measure for the progress of the transition. While the crystallinity parameter Q6 is most relevant in the initial stages of the crystallization, the later stages are dominated by a decrease in the potential energy.
We investigate the ejection dynamics of a ring polymer out of a cylindrical nanochannel using both theoretical analysis and three dimensional Langevin dynamics simulations. The ejection dynamics for ring polymers shows two regimes like for linear polymers, depending on the relative length of the chain compared with the channel. For long chains with length $N$ larger than the critical chain length $N_{c}$, at which the chain just fully occupies the nanochannel, the ejection for ring polymers is faster compared with linear chains of identical length due to a larger entropic pulling force; while for short chains ($N<N_c$), it takes longer time for ring polymers to eject out of the channel due to a longer distance to be diffused to reach the exit of the channel before experiencing the entropic pulling force. These results can help understand many biological processes, such as bacterial chromosome segregation.
We study the melting of a double stranded DNA in the presence of stretching forces, via 3D Monte-Carlo simulations, exactly solvable models and heuristic arguments. The resulting force-temperature phase diagram is dramatically different for the cases where the force is applied to only one strand or to both. Different assumptions on the monomer size of single and double stranded DNA lead to opposite conclusions as to whether DNA melts or not as it overstretches.
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.