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.     

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.     

Kinetics of polymer translocation through a pore

We theoretically study kinetics of a polymer threading through a pore embedded in a flat membrane. We numerically solve three coupled kinetic equations for the number n(1) of polymer segments in one side of the membrane and expansion factors of the polymer chain in each side of the membrane. We find the time evolution n(1) proportional to t(1/(1+nu)) at late stages and the translocation time tau(t) is scaled as tau(t) proportional to 1+nu) for large number n of the polymer segments, where nu is the effective size exponent of the radius of gyration of the polymer. When the polymer is translocated into a region with a good solvent condition (nu=3/5), we obtain n(1) proportional to t(5/8) and tau(t) proportional to n(8/5).     

Structure, dynamics, and phase transitions of tethered membranes: a Monte Carlo simulation study

A coarse-grained model of a self-avoiding tethered membrane with hexagonal coordination, embedded in three-dimensional space, is studied by means of extensive Monte Carlo computer simulations. The simulations are performed at various temperatures for membranes with linear size 5< or =L< or =50. We find that the membrane undergoes several folding transitions from a high-temperature flat phase to multiple-folded structure as the temperature is steadily decreased. Using a suitable order parameter and finite size scaling analysis, these phase transitions are shown to be of first order. The equilibrium shape of the membranes is analyzed by calculating the eigenvalues lambda(max) (2)> or =lambda(med) (2)> or =lambda(min) (2) of the inertia tensor. We present a systematic finite size scaling analysis of the radius of gyration and the eigenvalues of the inertia tensor at different phases of the observed folding transitions. In the high-temperature flat phase, the radius of gyration R(g) grows with the linear size of the membrane L as R(g) proportional to L(nu), where the exponent nu is approximately equal to 1.0. The eigenvalues of the inertia tensor scale as lambda(max) proportional to lambda(med) proportional to L(nu) and lambda(min) proportional to L(nu(min) ), whereby the roughness exponent nu(min) is approximately equal to 0.7. We also find that the time tau(R) of a self-avoiding membrane to diffuse a distance R(g) scales as tau(R) proportional to L(2nu+2), which is in good agreement with the theoretical predictions.     

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.     

Monte Carlo simulation of dendrimers in variable solvent quality

We study via lattice Monte Carlo simulation and Flory theory the properties of g=1-6 dendrimers in variable solvent quality. For all the generations studied, we find that the radius of gyration R(g) collapses significantly (factor of 2) going from athermal to extreme poor solvent conditions, indicating that varying solvent quality is an effective means of controlling dendrimer size. We also find that in athermal, theta, and extreme poor solvent conditions, the radius of gyration of dendrimers scales with the total number of monomers roughly as R(g) approximately N(1/3). However, a more careful analysis shows that in athermal and theta solvents, there is, in fact, a small but systematic deviation of R(g) from R(g) approximately N(1/3) scaling and the simulation data is described better by the Flory theory prediction of R(g) approximately N(1/5)[(g+1)m](2/5) in athermal solvents and R(g) approximately N(1/4)[(g+1)m](1/4) in theta solvents. We also find for our simulation data that stronger deviations from constant density scaling are possible, with scaling behavior as shallow as R(g) approximately N(0.26) possible for solvent conditions in between theta and the completely collapsed state. It is evident therefore that dendrimers do not obey (or even approximately obey) R(g) approximately N(1/3) scaling under all solvent conditions. Under all solvent conditions, we find that the intramolecular density is dense corelike (i.e., the density maximum is in the interior of the dendrimer) and terminal groups are delocalized throughout the dendrimer.     

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相似文献   

A Monte Carlo algorithm to study polymer translocation through nanopores. II. Scaling laws

In the first paper of this series, we developed a new one-dimensional Monte Carlo approach for the study of flexible chains that are translocating through a small channel. We also presented a numerical scheme that can be used to obtain exact values for both the escape times and the escape probabilities given an initial pore-polymer configuration. We now present and discuss the fundamental scaling behaviors predicted by this Monte Carlo method. Our most important result is the fact that, in the presence of an external bias E, we observe a change in the scaling law for the translocation time tau as function of the polymer length N: In the general expression tau approximately N(beta)E, the exponent changes from beta=1 for moderately long chains to beta=1+nu or beta=2nu for very large values of N (for Rouse and Zimm dynamics, respectively). We also observe an increase in the effective diffusion coefficient due to the presence of entropic pulling on unbiased polymer chains.     

Dynamics of polymer translocation into a circular nanocontainer through a nanopore

Using Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a circular nanocontainer through a nanopore under a driving force F. We observe that the translocation probability initially increases and then saturates with increasing F, independent of φ, which is the average density of the whole chain in the nanocontainer. The translocation time distribution undergoes a transition from a Gaussian distribution to an asymmetric distribution with increasing φ. Moreover, we find a nonuniversal scaling exponent of the translocation time as chain length, depending on φ and F. These results are interpreted by the conformation of the translocated chain in the nanocontainer and the time of an individual segment passing through the pore during translocation.     

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.     

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.     

Unexpected relaxation dynamics of a self-avoiding polymer in cylindrical confinement

We report extensive simulations of the relaxation dynamics of a self-avoiding polymer confined inside a cylindrical pore. In particular, we concentrate on examining how confinement influences the scaling behavior of the global relaxation time of the chain, tau, with the chain length N and pore diameter D. An earlier scaling analysis based on the de Gennes blob picture led to tau approximately N(2)D(13). Our numerical effort that combines molecular dynamics and Monte Carlo simulations, however, consistently produces different tau results for N up to 2000. We argue that the previous scaling prediction is only asymptotically valid in the limit N"D(53)"1, which is currently inaccessible to computer simulations and, more interestingly, is also difficult to reach in experiments. Our results are thus relevant for the interpretation of recent experiments with DNA in nano- and microchannels.     

Unforced translocation of a polymer chain through a nanopore: the solvent effect

The authors have performed the Langevin dynamics simulation to investigate the unforced polymer translocation through a narrow nanopore in an impermeable membrane. The effects of solvent quality controlled by the attraction strength lambda of the Lennard-Jones cosine potential between polymer beads and beads on two sides of the membrane on the translocation processes are extensively examined. For polymer translocation under the same solvent quality on both sides of the membrane, the two-dimensional and three-dimensional simulations confirm the scaling law of tautrans approximately N1+2upsilon for the translocation in the good solvent, where tautrans is the translocation time, N is the chain length, and upsilon is the Flory exponent. For the three-dimensional polymer translocation under different solvent qualities on two sides of the membrane, the translocation efficiency may be notably improved. The scaling law between tautrans and N varies from tautrans approximately N1+2upsilon to tautrans approximately N with the increase of the difference of solvent qualities, and the crossover occurs at the theta temperature point, where a scaling law of tautrans approximately N1.27 is found. The simulation results here also show that the translocation time changes from a wide and asymmetric distribution with a long tail to a narrow and symmetric distribution with the increase of the difference of the solvent qualities.     

Shear thinning in dilute polymer solutions

We use bead-spring models for a polymer coupled to a solvent described by multiparticle collision dynamics to investigate shear thinning effects in dilute polymer solutions. First, we consider the polymer motion and configuration in a shear flow. For flexible polymer models we find a sharp increase in the polymer radius of gyration and the fluctuations in the radius of gyration at a Weissenberg number approximately 1. We then consider the polymer viscosity and the effect of solvent quality, excluded volume, hydrodynamic coupling between the beads, and finite extensibility of the polymer bonds. We conclude that the excluded volume effect is the major cause of shear thinning in polymer solutions. Comparing the behavior of semiflexible chains, we find that the fluctuations in the radius of gyration are suppressed when compared to the flexible case. The shear thinning is greater and, as the rigidity is increased, the viscosity measurements tend to those for a multibead rod.     

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.     

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.     

Equilibrium conformational dynamics of a polymer in a solvent

Molecular dynamics simulations were used to study the conformational dynamics of a bead-spring model polymer in an explicit solvent under good solvent conditions. The dynamics of the polymer chain were investigated using an analysis of the time autocorrelation functions of the Rouse coordinates of the polymer chain. We have investigated the variation of the correlation functions with polymer chain length N, solvent density rho, and system size. The measured initial decay rates gamma(p) of the correlation functions were compared with the predictions from a theory of polymer dynamics which uses the Oseen tensor to describe hydrodynamic interactions between monomers. Over the range of chain lengths considered (N = 30-60 monomers), the predicted scaling of gamma(p) proportional to N(-3nu) was observed at high rho, where nu is the polymer scaling exponent. The predicted gamma(p) are generally higher than the measured values. This discrepancy increases with decreasing rho, as a result in the breakdown in the conditions required for the Oseen approximation. The agreement between theory and simulation at high rho improves considerably if the theoretical expression for gamma(p) is modified to avoid sum-to-integral approximations, and if the values of (R(p)2), which are used in the theory, are taken directly from the simulation rather than being calculated using approximate scaling relations. The observed finite-size scaling of gamma(p) is not quantitatively consistent with the theoretical predictions.     

Simulation study on the translocation of a partially charged polymer through a nanopore

The translocation of a partially charged polymer through a neutral nanopore under external electrical field is studied by using dynamic Monte Carlo method on a simple cubic lattice. One monomer in the polymer is charged and it suffers a driving force when it locates inside the pore. Two time scales, mean first passage time τ(FP) with the first monomer restricted to never draw back into cis side and translocation time τ for polymer continuously threading through nanopore, are calculated. The first passage time τ(FP) decreases with the increase in the driving force f, and the dependence of τ(FP) on the position of charged monomer M is in agreement with the theoretical results using Fokker-Planck equation [A. Mohan, A. B. Kolomeisky, and M. Pasquali, J. Chem. Phys. 128, 125104 (2008)]. But the dependence of τ on M shows a different behavior: It increases with f for M < N/2 with N the polymer length. The novel behavior of τ is explained qualitatively from dynamics of polymer during the translocation process and from the free energy landscape.     

Using an incremental mean first passage approach to explore the viscosity dependent dynamics of the unbiased translocation of a polymer through a nanopore

Noting the limitations of the standard characterization of translocation dynamics, an incremental mean first passage process methodology is used to more completely map the unbiased translocation of a polymer through a nanopore. In this approach, the average time t(0) required to reach successively increasing displacements for the first time is recorded - a measure shown to be more commensurate with the mean first passage nature of translocation. Applying this methodology to the results of Langevin dynamics simulations performed in three dimensions across a range of viscosities, a rich set of dynamics spanning regular diffusion at low viscosities to sub-diffusion at higher viscosities is revealed. Further, while the scaling of the net translocation time τ with polymer length N is shown to be viscosity-dependent, common regimes are found across all viscosities: super-diffusive behaviour at short times, an N-independent backbone consistent with τ ～ N(2.0) at low viscosities and τ ～ N(2.2) at higher viscosities for intermediate times, and N-dependent deviations from the backbone near the completion of translocation.