Search

We use the bond fluctuation model (BFM) to study the pore-blockade times of a translocating polymer of length N in two dimensions, in the absence of external forces on the polymer (i.e., unbiased translocation) and hydrodynamic interactions (i.e., the polymer is a Rouse polymer), through a narrow pore. Earlier studies using the BFM concluded that the pore-blockade time scales with polymer length as τd∼Nβ, with β=1+2ν, whereas some recent studies using different polymer models produce results consistent with β=2+ν, originally predicted by us. Here ν is the Flory exponent of the polymer; ν=0.75 in 2D. In this paper we show that for the BFM if the simulations are extended to longer polymers, the purported scaling τd∼N1+2ν ceases to hold. We characterize the finite-size effects, and study the mobility of individual monomers in the BFM. In particular, we find that in the BFM, in the vicinity of the pore the individual monomeric mobilities are heavily suppressed in the direction perpendicular to the membrane. After a modification of the BFM which counters this suppression (but possibly introduces other artifacts in the dynamics), the apparent exponent β increases significantly. Our conclusion is that BFM simulations do not rule out our theoretical prediction for unbiased translocation, namely, β=2+ν. [ABSTRACT FROM AUTHOR]

In this paper, we construct an effective model for the dynamics of an excluded-volume chain under confinement by extending the formalism of Rouse modes. We make specific predictions about the behavior of the modes for a single polymer confined to a tube. The results are tested against Monte Carlo simulations using the bond-fluctuation algorithm which uses a lattice representation of the polymer chain with excluded-volume effects. [ABSTRACT FROM AUTHOR]

This paper is aimed to develop a Brownian dynamics simulation method for strongly confined semiflexible polymers where numerical simulation plays an indispensable role in complementing theory and experiments. A wormlike chain under strong confinement is modeled as a string of virtual spherical beads connected by inextensible rods with length varying according to the confinement intensity of the chain measured by the Odijk deflection length. The model takes hydrodynamic interactions into account. The geometrical constraints associated with the inextensible rods are realized by the so-called linear constraint solver. The model parameters are studied by quantitatively comparing the simulated properties of a double-stranded DNA chain with available experimental data and theoretical predictions. [ABSTRACT FROM AUTHOR]

This work presents an immersed boundary method that allows fast Brownian dynamics simulation of solutions of polymer chains and other Brownian objects in complex geometries with fluctuating hydrodynamics. The approach is based on the general geometry Ewald-like method, which solves the Stokes equation with distributed regularized point forces in O(N) or O(NlogN) operations, where N is the number of point forces in the system. Time-integration is performed using a midpoint algorithm and Chebyshev polynomial approximation proposed by Fixman. This approach is applied to the dynamics of a genomic DNA molecule driven by flow through a nanofluidic slit with an array of nanopits on one wall of the slit. The dynamics of the DNA molecule was studied as a function of the Péclet number and chain length (the base case being λ-DNA). The transport characteristics of the hopping dynamics in this device differ at low and high Péclet number, and for long DNA, relative to the pit size, the dynamics is governed by the segments residing in the pit. By comparing with results that neglect them, hydrodynamic interactions are shown to play an important quantitative role in the hopping dynamics. [ABSTRACT FROM AUTHOR]

The dynamics of a collapsing polymer under a temperature quench in dilute solution is investigated by dissipative particles dynamics. Hydrodynamic interactions and many-body interaction are preserved naturally by incorporating explicit solvent particles in this approach. Our simulation suggests a four-stage collapse pathway: localized clusters formation, cluster coarsening in situ, coarsening involving global backbone conformation change into a crumpled globule, and compaction of the globule. For all the quench depths and chain lengths used in our study, collapse proceeds without the chain getting trapped in a metastable 'sausage' configuration, as reported in some earlier studies. We obtain the time scales for each of the first three stages, as well as its scaling with the quench depths ξ and chain lengths N. The total collapse time scales as τc ∼ ξ-0.46 ± 0.04N0.98 ± 0.09, with the quench depth and degree of polymerization. [ABSTRACT FROM AUTHOR]

The influence of hydrodynamic screening near a surface on the dynamics of a single semiflexible polymer is studied by means of Brownian dynamics simulations and hydrodynamic mean field theory. The polymer motion is characterized in terms of the mean squared displacements of the end-monomers, the end-to-end vector, and the scalar end-to-end distance. In order to control hydrodynamic screening effects, the polymer is confined to a plane at a fixed separation from the wall. When gradually decreasing this separation, a crossover from Zimm-type towards Rouse (free-draining) polymer dynamics is induced. However, this crossover is rather slow and the free-draining limit is not completely reached-substantial deviations from Rouse-like dynamics are registered in both simulations and theory-even at distances of the polymer from the wall on the order of the monomer size. Remarkably, the effect of surface-induced screening of hydrodynamic interactions sensitively depends on the type of dynamic observable considered. For vectorial quantities such as the end-to-end vector, hydrodynamic interactions are important and therefore surface screening effects are sizeable. For a scalar quantity such as the end-to-end distance, on the other hand, hydrodynamic interactions are less important, but a pronounced dependence of dynamic scaling exponents on the persistence length to contour length ratio becomes noticeable. Our findings are discussed against the background of single-molecule experiments on f-actin [L. Le Goff et al., Phys. Rev. Lett. 89, 258101 (2002)]. [ABSTRACT FROM AUTHOR]

We consider the dynamics of charged polymers in free solution in a slit geometry under the influence of an electrical field, applied at an angle to the plane parallel walls of the confinement. The simulations are carried out using the Brownian dynamics method with explicit counterions and implicit hydrodynamics. The hydrodynamic interactions between all the particles and the plane parallel walls are taken into account using a diffusion matrix which depends on slit geometry and the actual polyelectrolyte-solute conformations. We observe a selective transport of the charged polymers, as a function of the degree of polymerization and slit height. [ABSTRACT FROM AUTHOR]

Using a lattice-based Monte Carlo code for simulating self-avoiding flexible polymers in three dimensions in the absence of explicit hydrodynamics, we study their Rouse modes. For self-avoiding polymers, the Rouse modes are not expected to be statistically independent; nevertheless, we demonstrate that numerically these modes maintain a high degree of statistical independence. Based on high-precision simulation data we put forward an approximate analytical expression for the mode amplitude correlation functions for long polymers. From this, we derive analytically and confirm numerically several scaling properties for self-avoiding flexible polymers, such as (i) the real-space end-to-end distance, (ii) the end-to-end vector correlation function, (iii) the correlation function of the small spatial vector connecting two nearby monomers at the middle of a polymer, and (iv) the anomalous dynamics of the middle monomer. Importantly, expanding on our recent work on the theory of polymer translocation, we also demonstrate that the anomalous dynamics of the middle monomer can be obtained from the forces it experiences, by the use of the fluctuation-dissipation theorem. [ABSTRACT FROM AUTHOR]

Within Kirkwood theory, we study the translational diffusion coefficient of a single polymer chain in dilute solution, and focus on the small difference between the short-time Kirkwood value D[SUPK] and the asymptotic long-time value D. We calculate this correction term by highly accurate large-scale Brownian dynamics simulations, and show that it is in perfect agreement with the rigorous variational result D