Forced translocation of a polymer: Dynamical scaling versus molecular dynamics simulation

Forced and spontaneous translocation dynamics of a semiflexible active polymer in two dimensions

Polymer translocation is a fundamental topic in non-equilibrium physics and is crucially important to many biological processes in life. In the present work, we adopt two-dimensional Langevin dynamics simulations to study the forced and spontaneous translocation dynamics of an active filament. The influence of polymer stiffness on the underlying dynamics is explicitly analyzed. For the forced translocation, the results show a robust stiffness-induced inhibition, and the translocation time exhibits a dual-exponent scaling relationship with the bending modulus. Tension propagation (TP) is also examined, where we find prominent modifications in terms of both activity and stiffness. For spontaneous translocation into a pure solvent, the translocation time is almost independent of the polymer stiffness. However, when the polymer is translocated into a porous medium, an intriguing non-monotonic alteration of translocation time with increasing chain stiffness is demonstrated. The semiflexible chain is beneficial for translocation while the rigid chain is not conducive. Stiffness regulation on the diffusion dynamics of the polymer in porous media shows a consistent scenario. The interplay of activity, stiffness, and porous crowding provides a new mechanism for understanding the non-trivial translocation dynamics of an active filament in complex environments.

Polymer translocation driven by longitudinal and transversal time-dependent end-pulling forces

In this article, we simulate the translocation of a semiflexible homopolymer through an extended pore, driven by both a constant and a time-dependent end-pulled force, employing a model introduced in previous studies. The time dependence is simplistically modeled as a cosine function, and we distinguish between two scenarios for the driving--longitudinal force and transversal force-depending on the relative orientation of the force, parallel or perpendicular, respectively, with respect to the pore axis. Besides some key differences between the two drivings, the mean translocation times present a large minimum region as a function of the frequency of the force that is typical of the resonant activation effect. The presence of the minimum is independent on the elastic characteristics of the polymeric chains and reveals a linear relation between the optimum mean translocation time and the corresponding period of the driving. The mean translocation times show different scaling exponents with the polymer length for different flexibilities. Lastly, we derive an analytical expression of the mean translocation time for low driving frequency, which clearly agrees with the simulations.

Translocation of Hydrophobic Polyelectrolytes under Electrical Field: Molecular Dynamics Study

We studied the translocation of polyelectrolyte (PE) chains driven by an electric field through a pore by means of molecular dynamics simulations of a coarse-grained HP model mimicking high salt conditions. Charged monomers were considered as polar (P) and neutral monomers as hydrophobic (H). We considered PE sequences that had equally spaced charges along the hydrophobic backbone. Hydrophobic PEs were in the globular form in which H-type and P-type monomers were partially segregated and they unfolded in order to translocate through the narrow channel under the electric field. We provided a quantitative comprehensive study of the interplay between translocation through a realistic pore and globule unraveling. By means of molecular dynamics simulations, incorporating realistic force fields inside the channel, we investigated the translocation dynamics of PEs at various solvent conditions. Starting from the captured conformations, we obtained distributions of waiting times and drift times at various solvent conditions. The shortest translocation time was observed for the slightly poor solvent. The minimum was rather shallow, and the translocation time was almost constant for medium hydrophobicity. The dynamics were controlled not only by the friction of the channel, but also by the internal friction related to the uncoiling of the heterogeneous globule. The latter can be rationalized by slow monomer relaxation in the dense phase. The results were compared with those from a simplified Fokker-Planck equation for the position of the head monomer.

Polymer Translocation Time

The force- and flow-induced translocation processes of linear and ring polymers are studied using a combination of multiparticle collision dynamics and molecular dynamics, focusing on the behavior of the polymer translocation time. We compare the force- and flow-induced translocations of linear and ring polymers. It is found that when the translocation time (τ*) is characterized by scaling exponents, δ, δ', and α, via the relations τ* ∼ f^δN^α and τ* ∼ J^δ'N^α, the scaling exponents are not constants. For long chains tested, α = 1.0 for both force- and flow-induced translocations. The difference between the force- and flow-induced translocations stems from different monomer crowding effects due to distinct flow patterns outside the channel. Furthermore, general relations for polymer translocation time are derived for these two translocation processes, which are in good agreement with the simulation results. Our results provide clear molecular pictures for the force- and flow-induced translocations, which shed light on the understanding of translocation dynamics and provide guidance for practical applications such as molecular sequencing and ultrafiltration.

Effects of active crowder size and activity-crowding coupling on polymer translocation

Polymer translocation in complex environments is crucially important to many biological processes in life. In the present work, we adopted two-dimensional Langevin dynamics simulations to study the forced and unbiased polymer translocation dynamics in active and crowded media. The translocation time and probability are analyzed in terms of active force F_a, volume fraction φ and also the crowder size. The non-trivial active crowder size effect and activity-crowding coupling effect as well as the novel mechanism of unbiased translocation between two active environments with different active particle sizes are clarified. Firstly, for forced translocation, we reveal an intriguing non-monotonic dependence of the translocation time on the crowder size in the case of large activity. In particular, crowders of intermediate size similar to the polymer segment are proven to be the most favorable for translocation. Moreover, a facilitation-inhibition crossover of the translocation time with increasing volume fraction is observed, indicating a crucial activity-crowding coupling effect. Secondly, for unbiased translocation driven by different active crowder sizes, the translocation probability demonstrates a novel turnover phenomenon, implying the appearance of an opposite directional preference as the active force exceeds a critical value. The translocation time in both directions decreases monotonically with the active force. The asymmetric activity effect together with the entropic driving scenario provides a reasonable picture for the peculiar behavior observed in unbiased translocation.

Pulling a DNA molecule through a nanopore embedded in an anionic membrane: tension propagation coupled to electrostatics

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that in the realistic case of a single-stranded DNA molecule, dilute salt conditions characterized by weak charge screening, and a negatively charged membrane, the translocation dynamics is unexpectedly accelerated despite the presence of large repulsive electrostatic interactions between the polymer coil on thecisside and the charged membrane. This is due to the rapid release of the electrostatic potential energy of the coil during translocation, leading to an effectively attractive force that assists end-driven translocation. The speedup results in non-monotonic polymer length and membrane charge dependence of the exponentαcharacterizing the translocation timeτ∝N0αof the polymer with lengthN0. In the regime of long polymersN0 ≳   500, the translocation exponent exceeds its upper limitα= 2 previously observed for the same system without electrostatic interactions.

Translocation of Charged Polymers through a Nanopore in Monovalent and Divalent Salt Solutions: A Scaling Study Exploring over the Entire Driving Force Regimes

Langevin dynamics simulations are performed to study polyelectrolytes driven through a nanopore in monovalent and divalent salt solutions. The driving electric field E is applied inside the pore, and the strength is varied to cover the four characteristic force regimes depicted by a rederived scaling theory, namely the unbiased (UB) regime, the weakly-driven (WD) regime, the strongly-driven trumpet (SD(T)) regime and the strongly-driven isoflux (SD(I)) regime. By changing the chain length N, the mean translocation time is studied under the scaling form 〈 τ 〉 ∼ N α E - δ . The exponents α and δ are calculated in each force regime for the two studied salt cases. Both of them are found to vary with E and N and, hence, are not universal in the parameter's space. We further investigate the diffusion behavior of translocation. The subdiffusion exponent γ p is extracted. The three essential exponents ν s , q, z p are then obtained from the simulations. Together with γ p , the validness of the scaling theory is verified. Through a comparison with experiments, the location of a usual experimental condition on the scaling plot is pinpointed.

Sequencing of semiflexible polymers of varying bending rigidity using patterned pores

We study the translocation of a semiflexible polymer through extended pores with patterned stickiness, using Langevin dynamics simulations. We find that the consequence of pore patterning on the translocation time dynamics is dramatic and depends strongly on the interplay of polymer stiffness and pore-polymer interactions. For heterogeneous polymers with periodically varying stiffness along their lengths, we find that variation of the block size of the sequences and the orientation results in large variations in the translocation time distributions. We show how this fact may be utilized to develop an effective sequencing strategy. This strategy involving multiple pores with patterned surface energetics can predict heteropolymer sequences having different bending rigidity to a high degree of accuracy.

Polymer translocation under a pulling force: Scaling arguments and threshold forces

DNA translocation through nanopores is one of the most promising strategies for next-generation sequencing technologies. Most experimental and numerical works have focused on polymer translocation biased by electrophoresis, where a pulling force acts on the polymer within the nanopore. An alternative strategy, however, is emerging, which uses optical or magnetic tweezers. In this case, the pulling force is exerted directly at one end of the polymer, which strongly modifies the translocation process. In this paper, we report numerical simulations of both linear and structured (mimicking DNA) polymer models, simple enough to allow for a statistical treatment of the pore structure effects on the translocation time probability distributions. Based on extremely extended computer simulation data, we (i) propose scaling arguments for an extension of the predicted translocation times τ∼N^{2}F^{-1} over the moderate forces range and (ii) analyze the effect of pore size and polymer structuration on translocation times τ.

Driven translocation of a semi-flexible polymer through a nanopore

We study the driven translocation of a semi-flexible polymer through a nanopore by means of a modified version of the iso-flux tension propagation theory, and extensive molecular dynamics (MD) simulations. We show that in contrast to fully flexible chains, for semi-flexible polymers with a finite persistence length [Formula: see text] the trans side friction must be explicitly taken into account to properly describe the translocation process. In addition, the scaling of the end-to-end distance R N as a function of the chain length N must be known. To this end, we first derive a semi-analytic scaling form for R N, which reproduces the limits of a rod, an ideal chain, and an excluded volume chain in the appropriate limits. We then quantitatively characterize the nature of the trans side friction based on MD simulations. Augmented with these two factors, the theory shows that there are three main regimes for the scaling of the average translocation time τ ∝ N α . In the rod [Formula: see text], Gaussian [Formula: see text] and excluded volume chain [Formula: see text] ≫ 10 6 limits, α = 2, 3/2 and 1 + ν, respectively, where ν is the Flory exponent. Our results are in good agreement with available simulations and experimental data.

Translocation Time through a Nanopore with an Internal Cavity Is Minimal for Polymers of Intermediate Length

The translocation of polymers through nanopores with large internal cavities bounded by two narrow pores is studied via Langevin dynamics simulations. The total translocation time is found to be a nonmonotonic function of polymer length, reaching a minimum at intermediate length, with both shorter and longer polymers taking longer to translocate. The location of the minimum is shown to shift with the magnitude of the applied force, indicating that the pore can be dynamically tuned to favor different polymer lengths. A simple model balancing the effects of entropic trapping within the cavity against the driving force is shown to agree well with simulations. Beyond the nonmonotonicity, detailed analysis of translocation uncovers rich dynamics in which factors such as going to a high force regime and the emergence of a tail for long polymers dramatically change the behavior of the system. These results suggest that nanopores with internal cavities can be used for applications such as selective extraction of polymers by length and filtering of polymer solutions, extending the uses of nanopores within emerging nanofluidic technologies.

Dynamics of Polymer Translocation: A Short Review with an Introduction of Weakly-Driven Regime

As emphasized in a recent review (by V.V. Palyulin, T. Ala-Nissila, R. Metzler), theoretical understanding of the unbiased polymer translocation lags behind that of the (strongly) driven translocation. Here, we suggest the introduction of a weakly-driven regime, as described by the linear response theory to the unbiased regime, which is followed by the strongly-driven regime beyond the onset of nonlinear response. This provides a concise crossover scenario, bridging the unbiased to strongly-driven regimes.

Conformation Change, Tension Propagation and Drift-Diffusion Properties of Polyelectrolyte in Nanopore Translocation

Using Langevin dynamics simulations, conformational, mechanical and dynamical properties of charged polymers threading through a nanopore are investigated. The shape descriptors display different variation behaviors for the cis- and trans-side sub-chains, which reflects a strong cis-trans dynamical asymmetry, especially when the driving field is strong. The calculation of bond stretching shows how the bond tension propagates on the chain backbone, and the chain section straightened by the tension force is determined by the ratio of the direct to the contour distances of the monomer to the pore. With the study of the waiting time function, the threading process is divided into the tension-propagation stage and the tail-retraction stage. At the end, the drift velocity, diffusive property and probability density distribution are explored. Owing to the non-equilibrium nature, translocation is not a simple drift-diffusion process, but exhibits several intermediate behaviors, such as ballistic motion, normal diffusion and super diffusion, before ending with the last, negative-diffusion behavior.

Driven polymer translocation in good and bad solvent: Effects of hydrodynamics and tension propagation

We investigate the driven polymer translocation through a nanometer-scale pore in the presence and absence of hydrodynamics both in good and bad solvent. We present our results on tension propagating along the polymer segment on the cis side that is measured for the first time using our method that works also in the presence of hydrodynamics. For simulations we use stochastic rotation dynamics, also called multiparticle collision dynamics. We find that in the good solvent the tension propagates very similarly whether hydrodynamics is included or not. Only the tensed segment is by a constant factor shorter in the presence of hydrodynamics. The shorter tensed segment and the hydrodynamic interactions contribute to a smaller friction for the translocating polymer when hydrodynamics is included, which shows as smaller waiting times and a smaller exponent in the scaling of the translocation time with the polymer length. In the bad solvent hydrodynamics has a minimal effect on polymer translocation, in contrast to the good solvent, where it speeds up translocation. We find that under bad-solvent conditions tension does not spread appreciably along the polymer. Consequently, translocation time does not scale with the polymer length. By measuring the effective friction in a setup where a polymer in free solvent is pulled by a constant force at the end, we find that hydrodynamics does speed up collective polymer motion in the bad solvent even more effectively than in the good solvent. However, hydrodynamics has a negligible effect on the motion of individual monomers within the highly correlated globular conformation on the cis side and hence on the entire driven translocation under bad-solvent conditions.

Polyelectrolyte Threading through a Nanopore

Threading charged polymers through a nanopore, driven by electric fields E, is investigated by means of Langevin dynamics simulations. The mean translocation time 〈 τ 〉 is shown to follow a scaling law Nα, and the exponent α increases monotonically from 1.16 (4) to 1.40 (3) with E. The result is double-checked by the calculation of mean square displacement of translocation coordinate, which asserts a scaling behavior tβ (for t near τ) with β complying with the relation αβ = 2. At a fixed chain length N, 〈τ〉 displayed a reciprocal scaling behavior E-1 in the weak and also in the strong fields, connected by a transition E-1.64(5) in the intermediate fields. The variations of the radius of gyration of chain and the positions of chain end are monitored during a translocation process; far-from-equilibrium behaviors are observed when the driving field is strong. A strong field can strip off the condensed ions on the chain when it passes the pore. The total charges of condensed ions are hence decreased. The studies for the probability and density distributions reveal that the monomers in the trans-region are gathered near the wall and form a pancake-like density profile with a hump cloud over it in the strong fields, due to fast translocation.

Evaluating the applicability of the Fokker-Planck equation in polymer translocation: a Brownian dynamics study

Brownian dynamics (BD) simulations are used to study the translocation dynamics of a coarse-grained polymer through a cylindrical nanopore. We consider the case of short polymers, with a polymer length, N, in the range N = 21-61. The rate of translocation is controlled by a tunable friction coefficient, γ0p, for monomers inside the nanopore. In the case of unforced translocation, the mean translocation time scales with polymer length as <τ1> ∼ (N - Np)(α), where Np is the average number of monomers in the nanopore. The exponent approaches the value α = 2 when the pore friction is sufficiently high, in accord with the prediction for the case of the quasi-static regime where pore friction dominates. In the case of forced translocation, the polymer chain is stretched and compressed on the cis and trans sides, respectively, for low γ0p. However, the chain approaches conformational quasi-equilibrium for sufficiently large γ0p. In this limit the observed scaling of <τ1> with driving force and chain length supports the Fokker-Planck (FP) prediction that <τ> ∝ N/fd for sufficiently strong driving force. Monte Carlo simulations are used to calculate translocation free energy functions for the system. The free energies are used with the FP equation to calculate translocation time distributions. At sufficiently high γ0p, the predicted distributions are in excellent agreement with those calculated from the BD simulations. Thus, the FP equation provides a valid description of translocation dynamics for sufficiently high pore friction for the range of polymer lengths considered here. Increasing N will require a corresponding increase in pore friction to maintain the validity of the FP approach. Outside the regime of low N and high pore friction, the polymer is out of equilibrium, and the FP approach is not valid.

Using a Péclet number for the translocation of a polymer through a nanopore to tune coarse-grained simulations to experimental conditions

Coarse-grained simulations are often employed to study the translocation of DNA through a nanopore. The majority of these studies investigate the translocation process in a relatively generic sense and do not endeavor to match any particular set of experimental conditions. In this manuscript, we use the concept of a Péclet number for translocation, P(t), to compare the drift-diffusion balance in a typical experiment vs a typical simulation. We find that the standard coarse-grained approach overestimates diffusion effects by anywhere from a factor of 5 to 50 compared to experimental conditions using double stranded DNA (dsDNA). By defining a Péclet control parameter, λ, we are able to correct this and tune the simulations to replicate the experimental P(t) (for dsDNA and other scenarios). To show the effect that a particular P(t) can have on the dynamics of translocation, we perform simulations across a wide range of P(t) values for two different types of driving forces: a force applied in the pore and a pulling force applied to the end of the polymer. As P(t) brings the system from a diffusion dominated to a drift dominated regime, a variety of effects are observed including a non-monotonic dependence of the translocation time τ on P(t) and a steep rise in the probability of translocating. Comparing the two force cases illustrates the impact of the crowding effects that occur on the trans side: a non-monotonic dependence of the width of the τ distributions is obtained for the in-pore force but not for the pulling force.

Polymer translocation: the first two decades and the recent diversification

Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous-infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis.

Criteria for minimal model of driven polymer translocation

While the characteristics of the driven translocation for asymptotically long polymers are well understood, this is not the case for finite-sized polymers, which are relevant for real-world experiments and simulation studies. Most notably, the behavior of the exponent α, which describes the scaling of the translocation time with polymer length, when the driving force fp in the pore is changed, is under debate. By Langevin dynamics simulations of regular and modified translocation models using the freely jointed-chain polymer model we find that a previously reported incomplete model, where the trans side and fluctuations were excluded, gives rise to characteristics that are in stark contradiction with those of the complete model, for which α increases with fp. Our results suggest that contribution due to fluctuations is important. We construct a minimal model where dynamics is completely excluded to show that close alignment with a full translocation model can be achieved. Our findings set very stringent requirements for a minimal model that is supposed to describe the driven polymer translocation correctly.

Driven polymer transport through a periodically patterned channel

We study the driven transport of polymers in a periodically patterned channel using Langevin dynamics simulations in two dimensions. The channel walls are patterned with periodically alternating patches of attractive and non-attractive particles that act as trapping sites for the polymer. We find that the system shows rich dynamical behavior, observing giant diffusion, negative differential mobility, and several different transition mechanisms between the attractive patches. We also show that the channel can act as an efficient high-pass filter for polymers longer than a threshold length Nthr, which can be tuned by adjusting the length of the attractive patches and the driving force. Our findings suggest the possibility of fabricating polymer filtration devices based on patterned nanochannels.

Dynamics of polymer ejection from capsid

Polymer ejection from a capsid through a nanoscale pore is an important biological process with relevance to modern biotechnology. Here, we study generic capsid ejection using Langevin dynamics. We show that even when the ejection takes place within the drift-dominated region there is a very high probability for the ejection process not to be completed. Introducing a small aligning force at the pore entrance enhances ejection dramatically. Such a pore asymmetry is a candidate for a mechanism by which viral ejection is completed. By detailed high-resolution simulations we show that such capsid ejection is an out-of-equilibrium process that shares many common features with the much studied driven polymer translocation through a pore in a wall or a membrane. We find that the ejection times scale with polymer length, τ ∼ N(α). We show that for the pore without the asymmetry the previous predictions corroborated by Monte Carlo simulations do not hold. For the pore with the asymmetry the scaling exponent varies with the initial monomer density (monomers per capsid volume) ρ inside the capsid. For very low densities ρ ≤ 0.002 the polymer is only weakly confined by the capsid, and we measure α = 1.33, which is close to α=1.4 obtained for polymer translocation. At intermediate densities the scaling exponents α = 1.25 and 1.21 for ρ = 0.01 and 0.02, respectively. These scalings are in accord with a crude derivation for the lower limit α = 1.2. For the asymmetrical pore precise scaling breaks down, when the density exceeds the value for complete confinement by the capsid, ρ ⪆ 0.25. The high-resolution data show that the capsid ejection for both pores, analogously to polymer translocation, can be characterized as a multiplicative stochastic process that is dominated by small-scale transitions.

Driven translocation of a semi-flexible chain through a nanopore: a Brownian dynamics simulation study in two dimensions

We study translocation dynamics of a semi-flexible polymer chain through a nanoscopic pore in two dimensions using Langevin dynamics simulation in presence of an external bias F inside the pore. For chain length N and stiffness parameter κb considered in this paper, we observe that the mean first passage time <τ> increases as <τ(κb)>~<τ(κb=0)>lp(aN) , where κb and lp are the stiffness parameter and persistence length, respectively, and aN is a constant that has a weak N dependence. We monitor the time dependence of the last monomer xN(t) at the cis compartment and calculate the tension propagation time (TP) ttp directly from simulation data for <xN(t)> ~ t as alluded in recent nonequlibrium TP theory [T. Sakaue, Phys. Rev. E 76, 021803 (2007)] and its modifications to Brownian dynamics tension propagation theory [T. Ikonen, A. Bhattacharya, T. Ala-Nissila, and W. Sung, Phys. Rev. E 85, 051803 (2012); and J. Chem. Phys. 137, 085101 (2012)] originally developed to study translocation of a fully flexible chain. We also measure ttp from peak position of the waiting time distribution W(s) of the translocation coordinate s (i.e., the monomer inside the pore), and explicitly demonstrate the underlying TP picture along the chain backbone of a translocating chain to be valid for semi-flexible chains as well. From the simulation data, we determine the dependence of ttp on chain persistence length lp and show that the ratio ttp∕<τ> is independent of the bias F.

Polymer translocation dynamics in the quasi-static limit

Monte Carlo (MC) simulations are used to study the dynamics of polymer translocation through a nanopore in the limit where the translocation rate is sufficiently slow that the polymer maintains a state of conformational quasi-equilibrium. The system is modeled as a flexible hard-sphere chain that translocates through a cylindrical hole in a hard flat wall. In some calculations, the nanopore is connected at one end to a spherical cavity. Translocation times are measured directly using MC dynamics simulations. For sufficiently narrow pores, translocation is sufficiently slow that the mean translocation time scales with polymer length N according to <τ> ∝ (N - N(p))(2), where N(p) is the average number of monomers in the nanopore; this scaling is an indication of a quasi-static regime in which polymer-nanopore friction dominates. We use a multiple-histogram method to calculate the variation of the free energy with Q, a coordinate used to quantify the degree of translocation. The free energy functions are used with the Fokker-Planck formalism to calculate translocation time distributions in the quasi-static regime. These calculations also require a friction coefficient, characterized by a quantity N(eff), the effective number of monomers whose dynamics are affected by the confinement of the nanopore. This was determined by fixing the mean of the theoretical distribution to that of the distribution obtained from MC dynamics simulations. The theoretical distributions are in excellent quantitative agreement with the distributions obtained directly by the MC dynamics simulations for physically meaningful values of N(eff). The free energy functions for narrow-pore systems exhibit oscillations with an amplitude that is sensitive to the nanopore length. Generally, larger oscillation amplitudes correspond to longer translocation times.

Through the eye of the needle: recent advances in understanding biopolymer translocation

In recent years polymer translocation, i.e., transport of polymeric molecules through nanometer-sized pores and channels embedded in membranes, has witnessed strong advances. It is now possible to observe single-molecule polymer dynamics during the motion through channels with unprecedented spatial and temporal resolution. These striking experimental studies have stimulated many theoretical developments. In this short theory-experiment review, we discuss recent progress in this field with a strong focus on non-equilibrium aspects of polymer dynamics during the translocation process.

Cis-trans dynamical asymmetry in driven polymer translocation

During polymer translocation driven by, e.g., voltage drop across a nanopore, the segments in the cis side are incessantly pulled into the pore, which are then pushed out of it into the trans side. This pulling and pushing of polymer segments are described in the continuum level by nonlinear transport processes known, respectively, as fast and slow diffusions. By matching solutions of both sides through the mass conservation across the pore, we provide a physical basis for the cis and trans dynamical asymmetry, a feature repeatedly reported in recent numerical simulations. We then predict how the total driving force is dynamically allocated between cis (pulling) and trans (pushing) sides, demonstrating that the trans-side event adds a weak finite-chain length effect to the dynamical scaling.

Translocation dynamics of freely jointed Lennard-Jones chains into adsorbing pores

Polymer translocation into adsorbing nanopores is studied by using the Fokker-Planck equation of chain diffusion along the energy landscape calculated with Monte Carlo simulations using the incremental gauge cell method. The free energy profile of a translocating chain was found by combining two independent sub-chains, one free but tethered to a hard wall, and the other tethered inside an adsorbing pore. Translocation dynamics were revealed by application of the Fokker-Planck equation for normal diffusion. Adsorption of polymer chains into nanopores involves a competition of attractive adsorption and repulsive steric hindrance contributions to the free energy. Translocation times fell into two regimes depending on the strength of the adsorbing pore. In addition, we found a non-monotonic dependence of translocation times with increasing adsorption strength, with sharp peak associated with local free energy minima along the translocation coordinate.

Influence of non-universal effects on dynamical scaling in driven polymer translocation

We study the dynamics of driven polymer translocation using both molecular dynamics (MD) simulations and a theoretical model based on the non-equilibrium tension propagation on the cis side subchain. We present theoretical and numerical evidence that the non-universal behavior observed in experiments and simulations are due to finite chain length effects that persist well beyond the relevant experimental and simulation regimes. In particular, we consider the influence of the pore-polymer interactions and show that they give a major contribution to the non-universal effects. In addition, we present comparisons between the theory and MD simulations for several quantities, showing extremely good agreement in the relevant parameter regimes. Finally, we discuss the potential limitations of the present theories.

Influence of polymer-pore interaction on the translocation of a polymer through a nanopore

The translocation of a bond fluctuation polymer through an interacting nanopore is studied using dynamic Monte Carlo simulation. A driving force F is applied only for monomers inside the pore. The influence of polymer-pore interaction on the scaling relation τ~N(α) is studied for both unbiased and biased translocations, with τ the translocation time and N the polymer length. Results show that the exponent α is dependent on the polymer-pore interaction. For a noninteracting pore, we find α=2.48 for unbiased translocation and α=1.35 for strong biased translocation; for strong attraction, we find α=2.35 for unbiased translocation and α=1.22 for strong biased translocation. The unbiased translocation corresponds to the low-NF regime whereas the strong biased translocation corresponds to the high-NF regime.

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.

Process time distribution of driven polymer transport

We discuss the temporal distribution of dynamic processes in driven polymer transport inherent to flexible chains due to stochastic tension propagation. The stochasticity originates from the disordered initial configuration of an equilibrium polymer coil, which results in random paths for tension propagation. We consider the process time for when translocation occurs across a fixed pore and when stretching occurs by pulling the chain end. A scaling argument for the mean and standard deviation of the process time is provided using the two-phase picture for stochastic propagation. The two cases are found to differ remarkably. The process time distribution of the translocation exhibits substantial spreading even in the long-chain limit, unlike that found for the dynamics of polymer stretching. In addition, the process time distribution in the driven translocation is shown to have a characteristic asymmetric shape.

Unifying model of driven polymer translocation

We present a Brownian dynamics model of driven polymer translocation, in which nonequilibrium memory effects arising from tension propagation (TP) along the cis side subchain are incorporated as a time-dependent friction. To solve the effective friction, we develop a finite chain length TP formalism, based on the idea suggested by Sakaue [Phys. Rev. E 76, 021803 (2007)]. We validate the model by numerical comparisons with high-accuracy molecular dynamics simulations, showing excellent agreement in a wide range of parameters. Our results show that the dynamics of driven translocation is dominated by the nonequilibrium TP along the cis side subchain. Furthermore, by solving the model for chain lengths up to 10^{10} monomers, we show that the chain lengths probed by experiments and simulations are typically orders of magnitude below the asymptotic limit. This explains both the considerable scatter in the observed scaling of translocation time with respect to chain length, and some of the shortcomings of present theories. Our study shows that for a quantitative theory of polymer translocation, explicit consideration of finite chain length effects is required.

Event distributions of polymer translocation

We present event distributions for the polymer translocation obtained by extensive Langevin dynamics simulations. Such distributions have not been reported previously and they provide new understanding of the stochastic characteristics of the process. We extract at a high length scale resolution distributions of polymer segments that continuously traverse through a nanoscale pore. The obtained log-normal distributions together with the characteristics of polymer translocation suggest that it is describable as a multiplicative stochastic process. In spite of its clear out-of-equilibrium nature the forced translocation is surprisingly similar to the unforced case. We find forms for the distributions almost unaltered with a common cut-off length. We show that the individual short-segment and short-time movements inside the pore give the scaling relations τ∼Nα and τ∼f-β for the polymer translocation.