Langevin dynamics simulations of polymer translocation through nanopores

We investigate the dynamics of polymer translocation through a nanopore using two-dimensional Langevin dynamics simulations. In the absence of an external driving force, we consider a polymer which is initially placed in the middle of the pore and study the escape time tau(e) required for the polymer to completely exit the pore on either side. The distribution of the escape times is wide and has a long tail. We find that tau(e) scales with the chain length N as tau(e) approximately N(1+2nu), where nu is the Flory exponent. For driven translocation, we concentrate on the influence of the friction coefficient xi, the driving force E, and the length of the chain N on the translocation time tau, which is defined as the time duration between the first monomer entering the pore and the last monomer leaving the pore. For strong driving forces, the distribution of translocation times is symmetric and narrow without a long tail and tau approximately E(-1). The influence of xi depends on the ratio between the driving and frictional forces. For intermediate xi, we find a crossover scaling for tau with N from tau approximately N(2nu) for relatively short chains to tau approximately N(1+nu) for longer chains. However, for higher xi, only tau approximately N(1+nu) is observed even for short chains, and there is no crossover behavior. This result can be explained by the fact that increasing xi increases the Rouse relaxation time of the chain, in which case even relatively short chains have no time to relax during translocation. Our results are in good agreement with previous simulations based on the fluctuating bond lattice model of polymers at intermediate friction values, but reveal additional features of dependency on friction.     

Polymer chains in a soft nanotube: a Monte Carlo study

We study the equilibrium properties of flexible polymer chains confined in a soft tube by means of extensive Monte Carlo simulations. The tube wall is that of a single sheet six-coordinated self-avoiding tethered membrane. Our study assumes that there is no adsorption of the chain on the wall. By varying the length N of the polymer and the tube diameter D we examine the variation of the polymer gyration radius Rg and diffusion coefficient Ddiff in soft and rigid tubes of identical diameter and compare them to scaling theory predictions. We find that the swollen region of the soft tube surrounding the chain exhibits a cigarlike cylindrical shape for sufficiently narrow tubes with D相似文献   

Polymer translocation through a nanopore under an applied external field

We investigate the dynamics of polymer translocation through a nanopore under an externally applied field using the two-dimensional fluctuating bond model with single-segment Monte Carlo moves. We concentrate on the influence of the field strength E, length of the chain N, and length of the pore L on forced translocation. As our main result, we find a crossover scaling for the translocation time tau with the chain length from tau approximately N2nu for relatively short polymers to tau approximately N1+nu for longer chains, where nu is the Flory exponent. We demonstrate that this crossover is due to the change in the dependence of the translocation velocity v on the chain length. For relatively short chains v approximately N-nu, which crosses over to v approximately N(-1) for long polymers. The reason for this is that with increasing N there is a high density of segments near the exit of the pore, which slows down the translocation process due to slow relaxation of the chain. For the case of a long nanopore for which R parallel, the radius of gyration Rg along the pore, is smaller than the pore length, we find no clear scaling of the translocation time with the chain length. For large N, however, the asymptotic scaling tau approximately N1+nu is recovered. In this regime, tau is almost independent of L. We have previously found that for a polymer, which is initially placed in the middle of the pore, there is a minimum in the escape time for R parallel approximately L. We show here that this minimum persists for weak fields E such that EL is less than some critical value, but vanishes for large values of EL.     

Polymer translocation through a nanopore: a two-dimensional Monte Carlo study

We investigate the problem of polymer translocation through a nanopore in the absence of an external driving force. To this end, we use the two-dimensional fluctuating bond model with single-segment Monte Carlo moves. To overcome the entropic barrier without artificial restrictions, we consider a polymer which is initially placed in the middle of the pore and study the escape time tau required for the polymer to completely exit the pore on either end. We find numerically that tau scales with the chain length N as tau approximately N(1+2nu), where nu is the Flory exponent. This is the same scaling as predicted for the translocation time of a polymer which passes through the nanopore in one direction only. We examine the interplay between the pore length L and the radius of gyration R(g). For LR(g), we find tau approximately N. In addition, we numerically find the scaling function describing crossover between short and long pores. We also show that tau has a minimum as a function of L for longer chains when the radius of gyration along the pore direction R( parallel) approximately L. Finally, we demonstrate that the stiffness of the polymer does not change the scaling behavior of translocation dynamics for single-segment dynamics.     

The effect of hydrodynamic interactions on the dynamics of DNA translocation through pores

In this work, we investigate the effect of hydrodynamic interactions on the dynamics of DNA translocation through micropores. We simulate DNA as a bead-spring chain and use a lattice Boltzmann method to simulate the flow field that arises from the motion of the molecule. We investigate the free-draining entrance of DNA to the pore by diffusion and find that, consistent with experiments, molecules have a higher probability of entering the pore from one end. We then consider the electric-field driven translocation of 21-210 microm DNA with and without hydrodynamic interactions. Consistent with experiments, we study translocation events that are much shorter than the relaxation time of DNA. We find that the effect of hydrodynamic interactions on this process is to cause different regions of a molecule, other than the ones pulled by voltage or chain connectivity into the pore, to move toward the pore. We quantify this effect and show that it is smaller than the difference in the translocation dynamics of chains that arises from different initial configurations of the molecules. A power-law scaling of translocation time with chain length is observed, with exponents of 1.28+/-0.03 and 1.31+/-0.03 in simulations with and without hydrodynamic interactions, respectively. Our results are in good agreement with recent translocation experiments conducted in small pores and show that, for the regime considered in this work, hydrodynamic interactions play a minor role in the relation of the translocation time to chain length. For fast translocation processes, the effect of hydrodynamic interactions is local and the main factor determining the dynamics of DNA is the initial configuration of the molecules.     

Polymer translocation in solid-state nanopores: dependence of scaling behavior on pore dimensions and applied voltage

We investigate unforced and forced translocation of a Rouse polymer (in the absence of hydrodynamic interactions) through a silicon nitride nanopore by three-dimensional Langevin dynamics simulations, as a function of pore dimensions and applied voltage. Our nanopore model consists of an atomistically detailed nanopore constructed using the crystal structure of β-Si(3)N(4). We also use realistic parameters in our simulation models rather than traditional dimensionless quantities. When the polymer length is much larger than the pore length, we find the translocation time versus chain length scales as τ ～ N(2+ν) for the unforced case and as τ ～ N((1+2ν)/(1+ν)) for the forced case. Our results agree with theoretical predictions which indicate that memory effects and tension on the polymer chain play an important role during the translocation process. We also find that the scaling exponents are highly dependent on the applied voltage (force). When the length of the polymer is on the order of the length of the pore, we do not find a continuous scaling law, but rather scaling exponents that increase as the length of the polymer increases. Finally, we investigate the scaling behavior of translocation time versus applied voltage for different polymer and pore lengths. For long pores, we obtain the theoretical scaling law of τ ～ 1/V(α), where α ? 1 for all voltages and polymer lengths. For short pores, we find that α decreases for very large voltages and/or small polymer lengths, indicating that the value of α = 1 is not universal. The results of our simulations are discussed in the context of experimental measurements made under different conditions and with differing pore geometries.     

Dynamical scaling of single chains on adsorbing substrates: diffusion processes

We study the dynamics of tethered chains of length N on adsorbing surfaces, considering the dilute case; for this we use the bond fluctuation model and scaling concepts. In particular, we focus on the mean-square displacement of single monomers and of the center of mass of the chains. The characteristic time tau of the fluctuations of a free chain in a good solvent grows as tau approximately N(a), where the coefficient a obeys a=2nu+1. We show that the same coefficient also holds at the critical point of adsorption. At intermediate time scales single monomers show subdiffusive behavior; this concurs with the behavior calculated from scaling arguments based on the dynamical exponent a. In the adsorbed state tau(perpendicular), the time scale for the relaxation in the direction perpendicular to the surface, becomes independent of N; tau(perpendicular) is then the relaxation time of an adsorption blob. In the direction parallel to the surface the motion is similar to that of a two-dimensional chain and is controlled by a time scale given by tau(parallel) approximately N(2nu(2)+1)L(-2Delta(nu/nu)), where nu(2) is the Flory exponent in two dimensions, nu is the Flory exponent in three dimensions, and Deltanu=nu(2)-nu. For the motion parallel to the surface we find dynamical scaling over a range of about four decades in time.     

Dissipative particle dynamics simulations of polymer melts. I. Building potential of mean force for polyethylene and cis-polybutadiene

We present the derivation of coarse-grained force fields for two types of polymers, polyethylene (PE), and cis-polybutadiene (cis-PB), using the concept of potential of mean force. Coarse-grained force fields were obtained from microscopic simulations for several coarse-graining levels, i.e., different number of monomers lambda per mesoscopic unit called "bead." These force fields are then used in dissipative particle dynamics (DPD) simulations to study structural and dynamical properties of polymer melts of PE and cis-PB. The radial distribution functions g(R), the end-to-end distance R0, the end-to-end vector relaxation time tau, and the chain center of mass self-diffusion D(CM), are computed for different chain lengths at different coarse-graining factor lambda. Scaling laws typical of the Rouse regime are obtained for both polymers for chain lengths ranging from 6 to 50 beads. It is found that the end-to-end distance R0 obtained from DPD simulations agree well with values obtained from both microscopic simulations and experiments. The dependence of the friction coefficient used in DPD simulations versus the coarse-graining level is discussed in view of the overall scaling of the dynamical properties.     

Nondriven Polymer Translocation Through a Nanopore:Scaling for Translocation Time with Chain Length

We investigated the dynamics of the passage for a polymer chain through a nanopore in the absence of any external driving force with Weeks-Chandler-Andersen potential in two-dimensional simulations, in particular, focused our attention on the scaling law of the mean translocation time. We found that the effect of hydrodynamic interactions is the major factor in determining the scaling exponents with increasing pore size. The scaling close to N1+2v was observed when the hydrodynamic interactions were screened in the cases of small pore sizes, while the scaling close to N3v was obtained when the hydrodynamic interactions were present in the cases of large pore sizes.     

Universal properties of a single polymer chain in slit: Scaling versus molecular dynamics simulations

We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic aspects. We confirm a number of well known scaling predictions and analyze their range of validity by means of comprehensive molecular dynamics simulations using a coarse-grained bead-spring model of a flexible polymer chain. The normal and parallel components of the average end-to-end distance, mean radius of gyration and their distributions, the density profile, the force exerted on the slit walls, and the local bond orientation characteristics are obtained in slits of width D=4/10 (in units of the bead diameter) and for chain lengths N=50/300. We demonstrate that a wide range of static chain properties in normal direction can be described quantitatively by analytic model-independent expressions in perfect agreement with computer experiment. In particular, the observed profile of confinement-induced bond orientation is shown to closely match theory predictions. The anisotropy of confinement is found to be manifested most dramatically in the dynamic behavior of the polymer chain. We examine the relation between characteristic times for translational diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and normal relaxation times are in good agreement with our observations. A novel feature is the observed coupling of normal and lateral modes with two vastly different relaxation times. We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain length.     

Polymer chain dynamics at interfaces: role of boundary conditions at solid interface

Using classical molecular dynamics simulations, we study the dynamical properties of a single polymer chain dissolved in an explicit solvent and strongly adsorbed at solid-liquid interface. To circumvent a serious challenge posed by finite size effects due to long-range hydrodynamic effects, we developed a correction procedure that substantially limits the finite size effects. Concurrently, we provide an analysis of distinctly different size effects in the directions transverse and normal to the interface. We find that on analytically smooth interfaces, corresponding to the slip boundary condition, the motions of the polymer chain and the surrounding solvent are hydrodynamically coupled. This leads to the chain diffusion coefficient D scaling with the chain degree of polymerization N as D approximately N(-3/4), consistent with the Zimm dynamics for strongly adsorbed chains. Introduction of transverse forces at the interface results in loss of correlation between the motion of the polymer chain and the solvent. Consequently, D approximately N(-1), which is a characteristic of the Rouse dynamics.     

Kinetics of loop formation in polymer chains

We investigate the kinetics of loop formation in ideal flexible polymer chains (the Rouse model), and polymers in good and poor solvents. We show for the Rouse model, using a modification of the theory of Szabo, Schulten, and Schulten, that the time scale for cyclization is tau(c) approximately tau(0)N(2) (where tau(0) is a microscopic time scale and N is the number of monomers), provided the coupling between the relaxation dynamics of the end-to-end vector and the looping dynamics is taken into account. The resulting analytic expression fits the simulation results accurately when a, the capture radius for contact formation, exceeds b, the average distance between two connected beads. Simulations also show that when a < b, tau(c) approximately N(alpha)(tau), where 1.5 < alpha(tau) < or = 2 in the range 7 < N < 200 used in the simulations. By using a diffusion coefficient that is dependent on the length scales a and b (with a < b), which captures the two-stage mechanism by which looping occurs when a < b, we obtain an analytic expression for tauc that fits the simulation results well. The kinetics of contact formation between the ends of the chain are profoundly effected when interactions between monomers are taken into account. Remarkably, for N < 100, the values of tau(c) decrease by more than 2 orders of magnitude when the solvent quality changes from good to poor. Fits of the simulation data for tau(c) to a power law in N (tau(c) approximately N(alpha)(tau)) show that alpha(tau) varies from about 2.4 in a good solvent to about 1.0 in poor solvents. The effective exponent alpha(tau) decreases as the strength of the attractive monomer-monomer interactions increases. Loop formation in poor solvents, in which the polymer adopts dense, compact globular conformations, occurs by a reptation-like mechanism of the ends of the chain. The time for contact formation between beads that are interior to the chain in good solvents changes nonmonotonically as the loop length varies. In contrast, the variation in interior loop closure time is monotonic in poor solvents. The implications of our results for contact formation in polypeptide chains, RNA, and single-stranded DNA are briefly outlined.     

The rate constant of polymer reversal inside a pore

Translocation of biopolymers through pores is implicated in many biological phenomena. Confinement within a pore often breaks ergodicity on experimental and/or biological time scales by creating large entropic barriers to conformational rearrangements of the chain. Here, we study one example of such hindered rearrangement, in which the chain reverses its direction inside a long pore. Our goal is twofold. First, we study the dependence of the time scale of polymer reversal on the pore size and on the polymer length. Second, we examine the ability of simple one-dimensional theories to quantitatively describe a transition in a system with a complex energy landscape by comparing them with the exact rate constant obtained using brute-force simulations and the forward flux sampling method. We find that one-dimensional transition state theory (TST) using the polymer extension along the pore axis as the reaction coordinate adequately accounts for the exponentially strong dependence of the reversal rate constant on the pore radius r and the polymer length N, while the transmission factor, i.e., the ratio of the exact rate and the TST approximation, has a much weaker power law r and N dependence. We have further attempted to estimate the transmission factor from Kramer's theory, which assumes the reaction coordinate dynamics to be governed by a Langevin equation. However, such an approximation was found to be inadequate. Finally, we examine the scaling behavior of the reversal rate constant with N and r and show that finite size effects are important even for chains with N up to several hundreds.     

Dynamic aspects of mercury porosimetry: a lattice model study

Grand canonical Monte Carlo simulations using both Glauber dynamics and Kawasaki dynamics have been carried out for a recently developed lattice model of a nonwetting fluid confined in a porous material. The calculations are aimed at investigating the molecular scale mechanisms leading to mercury retention encountered during mercury porosimetry experiments. We first describe a set of simulations on slit and ink-bottle pores. We have studied the influence of the pore width parameter on the intrusion/extrusion curve shapes and investigated the corresponding mechanisms. Entrapment appears during Kawasaki dynamics simulations of extrusion performed on ink-bottle pores when the system is studied for short relaxation times. We then consider the more realistic and complex case of a Vycor glass building on recent work on the dynamics of adsorption of wetting fluids (Woo, H. J.; Monson, P. A. Phys. Rev. E 2003, 67, 041207). Our results suggest that mercury entrapment is caused by a decrease in the rate of mass transfer associated with the fragmentation of the liquid during extrusion.     

Hopping time of a hard disk fluid in a narrow channel

We use Monte Carlo (MC) and molecular dynamics (MD) methods to study the self-diffusion of hard disk fluids, confined within a narrow channel. The channels have a pore radius of Rp, above the passing limit of hard disk diameter (sigma(hd)). We focus on the average time (tau(hop)) needed for a hard disk to hop past a nearest neighbor in the longitudinal direction. This parameter plays a key role in a recent theory of the crossover from single-file diffusion to the bulk limit. For narrow channels near the hopping threshold (Rp=1 in units of sigma(hd)), both MC and MD results for tau(hop) diverge as approximately (Rp-1)(-2). Our results indicate that the scaling law exponent does not appear to be dependent on the differences between the two dynamics. This exponent is consistent with the prediction of an approximate transition state theory.     

Polymer translocation through a nanopore induced by adsorption: Monte Carlo simulation of a coarse-grained model

Dynamic Monte Carlo simulation of a bead-spring model of flexible macromolecules threading through a very narrow pore in a very thin rigid membrane are presented, assuming at the cis side of the membrane a purely repulsive monomer-wall interaction, while the trans side is attractive. Two choices of monomer-wall attraction epsilon are considered, one choice is slightly below and the other slightly above the "mushroom to pancake" adsorption threshold epsilon(c) for an infinitely long chain. Studying chain lengths N=32, 64, 128, and 256 and varying the number of monomers N(trans) (time t=0) that have already passed the pore when the simulation started, over a wide range, we find for epsilonepsilon(c) a finite number N(trans)(t=0) suffices that the translocation probability is close to unity. In the case epsilonepsilon(c), we find that the translocation time scales as tau proportional, variant N(1.65+/-0.08). We suggest a tentative scaling explanation for this result. Also the distribution of translocation times is obtained and discussed.     

Static and dynamic properties of tethered chains at adsorbing surfaces: a Monte Carlo study

We present extensive Monte Carlo simulations of tethered chains of length N on adsorbing surfaces, considering the dilute case in good solvents, and analyze our results using scaling arguments. We focus on the mean number M of chain contacts with the adsorbing wall, on the chain's extension (the radius of gyration) perpendicular and parallel to the adsorbing surface, on the probability distribution of the free end and on the density profile for all monomers. At the critical adsorption strength epsilon(c) one has M(c) approximately N(phi), and we find (using the above results) as best candidate phi to equal 0.59. However, slight changes in the estimation of epsilon(c) lead to large deviations in the resulting phi; this might be a possible reason for the difference in the phi values reported in the literature. We also investigate the dynamical scaling behavior at epsilon(c), by focusing on the end-to-end correlation function and on the correlation function of monomers adsorbed at the wall. We find that at epsilon(c) the dynamic scaling exponent a (which describes the relaxation time of the chain as a function of N) is the same as that of free chains. Furthermore, we find that for tethered chains the modes perpendicular to the surface relax quicker than those parallel to it, which may be seen as a splitting in the relaxation spectrum.     

Numerical and theoretical study on the mechanism of biopolymer translocation process through a nano-pore

We conducted a numerical study on the translocation of a biopolymer from the cis side to the trans side of a membrane through a synthetic nano-pore driven by an external electric field in the presence of hydrodynamic interactions (HIs). The motion of the polymer is simulated by 3D Langevin dynamics technique using a worm-like chain model of N identical beads, while HI between the polymer and fluid are incorporated by the lattice Boltzmann equation. The translocation process is induced by electrophoretic force, which sequentially straightens out the folds of the initial random configuration of the polymer chain on the cis side. Our simulation results on translocation time and velocity are in good quantitative agreement with the corresponding experimental ones when the surface charge on the nano-pore and the HI effect are considered explicitly. We found that the translocation velocity of each bead inside the nano-pore mainly depends upon the length of the straightened portion of the polymer in forced motion near the pore. We confirmed this by a theoretical formula. After performing simulations with different pore lengths, we observed that translocation velocity mainly depends upon the applied potential difference rather than upon the electric field inside the nano-pore.     

Excluded volume effects on the intrachain reaction kinetics

On the basis of the recently developed optimized Rouse-Zimm theory of chain polymers with excluded volume interactions, we calculate the long-time first-order rate constant k(1) for end-to-end cyclization of linear chain polymers. We first find that the optimized Rouse-Zimm theory provides the longest chain relaxation times tau(1) of excluded volume chains that are in excellent agreement with the available Brownian dynamics simulation results. In the free-draining limit, the cyclization rate is diffusion-controlled and k(1) is inversely proportional to tau(1), and the k(1) values calculated using the Wilemski-Fixman rate theory are in good agreement with Brownian dynamics simulation results. However, when hydrodynamic interactions are included, noticeable deviations are found. The main sources of errors are fluctuating hydrodynamic interaction and correlation hole effects as well as the non-Markovian reaction dynamic effect. The physical natures of these factors are discussed, and estimates for the magnitudes of required corrections are given. When the corrections are included, the present theory allows the prediction of accurate k(1) values for the cyclization of finite-length chains in good solvents as well as the correct scaling exponent in the long-chain limit.     

Computer simulation of linkage of two ring chains

We performed off-lattice Monte Carlo simulations of links of two model ring chains with chain length N up to 32,768 in the theta solution or amorphous bulk state by using a random walk model (Model I), and molecular dynamics simulations of two model ring chains in solution with excluded volume interaction (Model II) to investigate topological effects on the geometry of link and ring conformation. In the case of Model I, the mean squared linking number, its distribution, and the size of two chains with fixed linking number are investigated. Our simulation results confirm the previous theoretical prediction that the mean squared linking number decays as pe(-qs(2)) with the distance of centers of chain mass s, where p and q are found to be chain length dependent and q asymptotically approaches to 0.75 as chain length increases. The linking number distribution of two chains has a universal form for long chains, but our simulation results clearly show that the distribution function deviates from the Gaussian distribution, a fact not predicted by any previous theoretical work. A scaling prediction is proposed to predict the link size, and is checked for our simulations for the Model II. The simulation results confirmed the scaling prediction of the blob picture that the link with linking number m occupies a compact volume of m blobs, and the size of the link is asymptotic to R(L) ≈ bN(ν)m(1/3-ν), where N is the chain length, and v is the Flory exponent of polymer in solutions.