Theoretical Modeling of Polymer Translocation: From the Electrohydrodynamics of Short Polymers to the Fluctuating Long Polymers

The theoretical formulation of driven polymer translocation through nanopores is complicated by the combination of the pore electrohydrodynamics and the nonequilibrium polymer dynamics originating from the conformational polymer fluctuations. In this review, we discuss the modeling of polymer translocation in the distinct regimes of short and long polymers where these two effects decouple. For the case of short polymers where polymer fluctuations are negligible, we present a stiff polymer model including the details of the electrohydrodynamic forces on the translocating molecule. We first show that the electrohydrodynamic theory can accurately characterize the hydrostatic pressure dependence of the polymer translocation velocity and time in pressure-voltage-driven polymer trapping experiments. Then, we discuss the electrostatic correlation mechanisms responsible for the experimentally observed DNA mobility inversion by added multivalent cations in solid-state pores, and the rapid growth of polymer capture rates by added monovalent salt in α-Hemolysin pores. In the opposite regime of long polymers where polymer fluctuations prevail, we review the iso-flux tension propagation (IFTP) theory, which can characterize the translocation dynamics at the level of single segments. The IFTP theory is valid for a variety of polymer translocation and pulling scenarios. We discuss the predictions of the theory for fully flexible and rodlike pore-driven and end-pulled translocation scenarios, where exact analytic results can be derived for the scaling of the translocation time with chain length and driving force.


Introduction
DNA is the key transmitter of the biological information carrying our genetic heritage. Fast and inexpensive access to this information is essential for various purposes ranging from the treatment of genetic diseases in medicine to the identification of harmful organisms in metagenomic sciences or DNA profiling in forensic sciences [1,2]. During the past three decades, this need has stimulated intensive research work on the development of efficient and low-cost biosequencing techniques such as the field-driven translocation of polymers through nanoscale pores [3]. This biosensing approach consists of mapping the sequence of the polymer portion translocating through the pore from the current perturbations caused by the biopolymer [4][5][6][7][8][9][10]. As the accuracy of this mapping depends sensitively on the duration of the current signal triggered by the presence of the translocating polymer, efficient use of this method requires a high degree of control on the dynamics of the molecule. At this point, one needs theoretical models able to predict the dependence of the polymer translocation dynamics on the experimentally controllable system parameters such as salt concentration, polymer charge and length, pore charge and size, and the external forces driving the translocation process.
The theoretical formulation of polymer translocation is a highly ambitious task. This complex transport process is indeed governed by a combination of effects such as pore electrohydrodynamics resulting from the electrophoretic (EP) and electroosmotic (EO) forces acting on the polymer, direct electrostatic polymer-membrane coupling, and entropic effects originating from conformational polymer fluctuations and steric polymer-membrane interactions. Two rather complementary approaches of distinct nature have so far been adopted for investigating polymer translocation: Approaches based on coarse-grained conformational models and electrohydrodynamic formalisms.
In the case of polymers longer than the translocated pores whose characteristic size L m ∼10-100 nm is comparable with the DNA persistence length, polymer fluctuations are substantial but the electrohydrodynamic forces can be assumed to act locally on DNA, i.e., exclusively on the polymer portion confined to the pore. This scale separation allows to bypass the details of the pore electrohydrodynamics that can be absorbed into the effective force f driving the polymer and the effective pore friction η p on it, enabling coarse-grained modeling of effects associated with non-equilibrium polymer conformations. Such coarse-grained models are easily amenable to molecular dynamics (MD) and Monte Carlo simulations [11][12][13][14][15], but even then it is a challenge to explicitly include electrostatic polymer-membrane interactions and they are usually assumed to be negligible. On the theoretical side, a comprehensive theory for driven polymer translocation dynamics has been developed based on the idea of non-equilibrium tension propagation [16][17][18][19][20][21]. The basic idea in this theory is to focus on the dynamics of a single degree of freedom, the translocation coordinate s(t), and include all the many-body effects arising from the (non-equilibrium) chain conformations on the cis side of the membrane into a time-dependent friction η cis (t). This leads to a Brownian dynamics type of equation for s(t) which makes the problem both analytically and numerically tractable, and allows exact analytic results for the scaling of the translocation time as a function of the chain length. As explained in Section 3, this iso-flux tension propagation (IFTP) theory has been benchmarked for a variety of driven polymer translocation scenarios with excellent agreement with coarse-grained MD simulations and relevant experiments.
In the opposite regime of polymers whose length is comparable to the thickness of the translocated membrane, polymer fluctuations in the pore can be assumed to be negligible but the electrostatic polymer-pore interactions and the electrohydrodynamic pore effects have to be accurately taken into account. A consistent electrohydrodynamic modeling of polymer translocation was initiated by Ghosal in Ref. [22]. Via the coupled solution of the electrostatic Poisson-Boltzmann (PB) and hydrodynamic Navier-Stokes equations, Ghosal derived DNA translocation velocity as the superposition of the EP and EO velocity components. The role played by polymer-pore interactions on the unzipping of a DNA hairpin was investigated in Ref. [23] without pore hydrodynamics. The effect of the EO flow on diffusion-limited polymer capture was studied in Ref. [24] and the predictions of different electrostatic models [25,26] were compared with translocation experiments [27]. Within a Smoluchowski formalism, we incorporated in Ref. [28] mean-field (MF) level electrostatic polymer-membrane interactions into the electrohydrodynamic model of Ref. [22]. This unified polymer translocation theory was extended in Ref. [29] to include electrostatic correlations. The extended theory was applied to the experiments of Ref. [30] to explain the electrohydrodynamic mechanism behind the polyvalent-cation-induced DNA mobility reversal. In the same work, a new mechanism of facilitated polymer capture by charge-inverted EO flow was also identified. Finally, in Ref. [31], we revealed an electrostatic trapping mechanism enabling the extension of the polymer translocation time, which would allow to enhance the duration of the current readout in translocation experiments. Very recently, we have also taken a step towards a unified theory of polymer translocation by incorporating the electrostatic coupling of the membrane with the cis and trans portions of the polymer outside the nanopore into the stiff polymer limit of the IFTP theory [32].
In this article, we present a comparative review of the electrohydrodynamic and coarse-grained approaches described above. In the first part of the manuscript, we discuss in Section 2 the electrohydrodynamic translocation model and its application to various experimental setups. First in Section 2.1, we explain the theoretical framework of the approach. Section 2.2 is devoted to the application of the theory to pressure-voltage-driven translocation experiments in monovalent salt where the system is governed by MF electrohydrodynamics. Section 2.3 is devoted to the translocation experiments with polyvalent salt where the high ion valency results in a departure from MF-level electrohydrodynamics. In Section 2.4, we focus on polymer translocation through α-Hemolysin (αHL) pores of subnanometer radius where the strong confinement results in polarization forces driving the system away from the MF transport regime. In the remaining part of the article, we focus on the regime of long coarse-grained polymers and review the iso-flux tension propagation theory able to account for the conformational polymer fluctuations during translocation. Our main results and prospects are discussed in the Summary and Conclusions section.

Electrohydrodynamic Approach to the Translocation of Short Polymers
In this section, we focus on the translocation of short polymers whose size is comparable to the size of the translocated nanopore. The comparable spatial scale of the polymer and the pore requires a detailed consideration of the pore electrohydrodynamics driving the translocation process. This point is the main motivation behind the electrohydrodynamic translocation model presented in this part. The configuration of the model is illustrated in Figure 1. A cylindrical pore of radius d and total length L m extends along the z axis. The ends of the pore are in contact with an ion reservoir composed of p ionic species, with the species i of valency q i and bulk concentration ρ bi . The pore surface at r = d carries a fixed negative charge distribution of density σ m (r) = −σ m δ(r − d). The translocating polymer on the z axis is modeled as a cylinder of radius a, total length L p , and surface charge density σ p (r) = −σ p δ(r − a). Moreover, the polymer portion located in the pore has length l p . The reaction coordinate of the translocation is the location z p of the lower end of the molecule. The translocation of the polymer from the cis side at z = 0 to the trans side at z = L m is induced by an externally applied hydrostatic pressure gradient ∆P and electric voltage ∆V. In addition to these external driving forces, the polyelectrolyte is also subjected to direct polymer-membrane interactions characterized by the electrostatic potential V p (z p ). Section 2.1 reviews the electrohydrodynamic formalism of polymer translocation introduced in Refs. [28,31]. In Section 2.2, we present the application of this theory to solid-state pores and its comparison with pressure-voltage-driven translocation experiments [33]. In Section 2.3, we discuss the effect of charge correlations on polymer translocation in polyvalent electrolytes and the resulting DNA mobility reversal [31] observed in voltage-driven translocation experiments [30]. Finally, in Section 2.4, we investigate surface polarization effects on polymer translocation through αHL pores of subnanometer confinement.

Electrohydrodynamic Formalism of Polymer Translocation
The translocation process is characterized by the polymer diffusion equation where the polymer current is Schematic of a polymer translocating through a cylindrical pore of radius d, length L m , and negative surface charge density −σ m . The anionic polymer translocating on the z axis is a cylinder of length L p , radius a, and surface charge density −σ p . l p is the length of the polymer portion located in the pore. Translocation is driven by the externally applied voltage ∆V = V t − V c , resulting in the electric field E = −∆V/L mûz , and the hydrostatic pressure ∆P = P c − P t .
In Equations (1) and (2), the functions c(z p , t) and J(z p , t) stand respectively for the polymer density and flux. Moreover, the transverse diffusion coefficient of the cylindrical polymer is D = ln(L p /2a)/(3πηL p β) [34], with the inverse thermal energy β = 1/(k B T) and water viscosity η = 8.91 × 10 −4 Pa s. In Equation (2), the first and second terms correspond respectively to the diffusive flux component, and the convective flux component associated with the polymer velocity v p (z p ).
We restrict ourselves to the steady-state regime where the polymer flux becomes constant and uniform, i.e., J(z p , t) = J st . Introducing the effective polymer potential U p (z p ) defined by Equation (2) can be recast as We integrate now Equation (4) with the absorbing boundary condition (BC) c(z ex ) = 0 at the pore exit and impose the polymer density on the cis side of the reservoir, i.e., c(0) = c cis . This yields The polymer population in the pore is given by the integral of Equation (7), The polymer translocation frequency corresponding to the inverse translocation time is defined as the polymer flux per polymer population in the pore, i.e., τ −1 p = J st /N. Moreover, the polymer capture rate is given by the polymer flux per reservoir concentration, i.e., R c = J st /c cis . Using these definitions together with Equations (6)- (8), the polymer capture rate and translocation time follow as The rate R c corresponds to the average speed at which a successful polymer capture takes place. In the drift regime characterized by weak polymer-pore interactions, the limit V p (z p ) → 0 of Equation (9) yields where the second equality is valid for high voltages and a positive drift velocity. We finally note that for comparison with pressure-voltage trapping experiments, the average translocation velocity will also be needed. The average polymer velocity is defined as

Derivation of the Polymer Velocity v p (z p )
We first note that the evaluation of the polymer capture rate, translocation time, and average velocity defined in Equations (9)-(12) requires the effective polymer potential U p (z p ). The calculation of this potential necessitates in turn the knowledge of the polymer velocity v p (z p ) in Equation (3). To derive the latter, we first express the PB and Stokes equations for the net electrostatic potential φ(r) and liquid velocity u c (r) in the cylindrical nanopore, where we introduced the radial distance r from the pore axis, the Bjerrum length B = βe 2 /(4πε w ) with the solvent (water) permittivity ε w = 80 and the unit charge e, and the density of mobile ions ρ c (r) and fix charges σ(r) = −σ m δ(r − d) − σ p δ(r − a). Next, we eliminate from Equations (13) and (14) the ion density ρ c (r), integrate the resulting equation, and impose the no-slip BC at the pore wall u c (d) = 0 and at the DNA surface u c (a) = v p (z p ). We finally account for Gauss' law φ (a) = 4π B σ p and also the force balance relation on the polymer F el + F dr + F b = 0, with the electric force F el = 2πaL p eE, the hydrodynamic drag force F dr = 2πaL p ηu c (a), and the force F b = −V p (z p ) associated with electrostatic polymer-membrane interactions. This yields the liquid and polymer velocities in the form with the effective polymer diffusion coefficient in the pore medium the coefficient of electrophoretic (EP) polymer mobility µ e = ε w k B T/(eη), and the drift velocity component induced by the external voltage and pressure, where we introduced the geometric factor The first term on the r.h.s. of the drift velocity Equation (18) includes the effect of the voltage-induced EP force on DNA (the first term in the bracket) and the opposing force from the electroosmotic (EO) flow drag (the second term in the bracket). The second term of Equation (18) corresponds in turn to the contribution from the pressure-induced streaming flow to the DNA velocity. Then, the second term of Equation (16) brings the effect of electrostatic polymer-membrane interactions on the polymer velocity. As a result of the no-slip relation v p (z p ) = u c (a), the terms on the r.h.s. of Equation (15) clearly indicate the contribution from the same effects to the convective liquid velocity u c (r). Integrating now Equation (3) with Equation (16), one finally obtains the effective polymer potential in Equations (9)-(12) as 2.1.3. Derivation of the Interaction Potential V p (z p ) We explain next the derivation of the electrostatic polymer-membrane interaction potential V p (z p ) in the MF regime of weak surface charges and physiological monovalent salt concentrations. The extension of this calculation beyond MF electrostatics is rather involved and this generalization can be found in Refs. [29,31]. In the MF linear response regime, the polymer-membrane interaction potential induced by the electrostatic coupling between the membrane potential and the polymer charges Q pol = 2πal p σ p located in the pore reads In Equation (21), the potential φ m (r) induced solely by the membrane charges is obtained from the solution of the PB Equation (13) without the polymer charge, i.e., The calculation of the potentials φ(r) and φ m (r) can be found in Ref. [28]. Moreover, the position-dependent length of the polymer portion in the pore reads where we introduced the auxiliary lengths The terms on the r.h.s. of Equation (23) are associated with the regimes of polymer capture (z p < L − ), translocation at constant length (L − < z p < L + ), and polymer escape (z p > L + ), respectively. Finally, defining the characteristic inverse lengths associated with the drift velocity in Equation (18) and the electrostatic interaction potential in Equation (21), the polymer translocation velocity in Equation (16) and the interaction potential in Equation (20) take the simpler forms The inverse lengths λ d,b in Equation (25) allow to characterize the polymer capture and translocation dynamics within the drift-driven regime λ d λ b of weak polymer-membrane interactions and the barrier-driven regime λ b λ d where the electrostatic polymer-membrane coupling takes over the voltage and/or pressure-induced drift force [31]. The explicit forms of Equations (9)-(12) obtained with the velocity in Equation (26) and the potential profile in Equation (27) are given in Appendix A in terms of the inverse lengths defined in Equation (25). Evaluating now the average polymer velocity of Equation (12) with Equations (26) and (27), one gets where the coefficients J i are also given in Appendix A. Finally, the translocation time in Equation (10) follows as where the explicit forms of the characteristic time for polymer capture τ 1 , translocation τ 2 , and escape τ 3 are given in Appendix B.

Polymer Conductivity of Solid-State Pores: MF Electrohydrodynamics with Monovalent Salt
We consider polymer translocation events in solid-state pores and monovalent salt solutions where the electrohydrodynamic interactions are characterized by MF electrostatics. The DNA surface charge density is fixed to the value σ p = 0.4 e/nm 2 previously obtained in Ref. [35] by fitting experimental current blockade data. In Ref. [31], it was shown that both in the barrier-driven regime λ b λ d and drift-dominated regime λ d λ b , the translocation velocity in Equation (28) can be well approximated Passing to the Debye-Hückel (DH) limit of strong salt, this approximation yields [31] v where we used the DH screening parameter κ 2 = 8π B ρ b and introduced the geometric coefficients g and f m,p given in Appendix A. In Equation (30), the first term on the r.h.s. takes into account the electrophoretic (EP) drift force by the electric field on the polymer charges (positive term) and the electroosmotic (EO) flow drag induced by the counterions attracted by the charged membrane (negative term). In the present case of anionic polymers translocating through like-charged membranes, the EO flow opposing the EP drift reduces the polymer velocity, i.e., σ m ↑ R c ↓ v p ↓. Then, the second term in Equation (30) corresponds to the force exerted on the polymer by the pressure-induced streaming flow through the pore. Finally, the third term including the product σ p σ m > 0 accounts for the electrostatic polymer-membrane interactions. Due to the resulting like-charge repulsion, the negative interaction term acts as an electrostatic barrier hindering the polymer capture by the nanopore.

Comparison with Pressure-Voltage Trapping Experiments
Equation (30) shows that in the drift-driven regime, the average velocity rises linearly with the external voltage ∆V. This linear dependence has been observed in experiments and simulations [36,37]. To understand the pressure dependence of the translocation velocity and time, we now focus on the pressure-driven translocation experiments. Figure 2a compares the average polymer velocity in Equation (28) (solid curve) with the experimental data of Ref. [33] (squares). The numerical values of the model parameters taken from Ref. [33] are the negative external voltage ∆V = −100 mV opposing the drag of the streaming flow, the electrolyte concentration ρ b = 1.6 M, the number of monomers in the DNA sequence N = 615 bps corresponding to the DNA length L p = 180 nm, and the pore radius d = 5 nm. The values of the membrane thickness L m = 200 nm and charge σ m = 0.13 e/nm 2 were adjusted to obtain the optimal agreement with the magnitude of the translocation velocity data.   Figure 2a shows that within the experimental scattering, the theoretical result agrees well with the polymer velocity data. For an analytical insight into the pressure dependence of the experimental data, we recast Equation (30) with the critical pressure where the translocation velocity vanishes and the polymer gets trapped Equation (31) reported in Figure 2a by circles indicates that the average polymer velocity grows linearly with the pressure gradient. Figure 2b compares the theoretical translocation time τ p in Equation (29) (solid curve) with the experimental escape times of Ref. [33]. The theoretical result obtained with the same parameters as in Figure 2a can accurately reproduce the general trend of the experimental data. To identify the scaling of the experimental time data with the pressure gradient, we note that the occurrence of a successful polymer translocation necessitates the polymer of average velocity v p to travel the total distance z ex = L p + L m over the time τ p . This allows to approximate the translocation time as τ p ≈ z ex / v p . Using Equation (31), this yields The approximative formula (33) reported in Figure 2b by open circles indicates that the pronounced rise of the translocation time at low pressures is characterized by the inverse power law scaling

Salt and Polymer Length Dependence of Pressure-Voltage-Driven Translocation Events
We discuss here the influence of the ion density and polymer length on pressure-voltage-driven polymer translocation events. Figure 3a,b displays the ion density dependence of the polymer translocation velocity obtained from Equation (28) (solid curves) and the drift formula (30) (dots at ∆P = ±2 atm). We first focus on Figure 3a where polymer translocation is driven by the streaming flow (∆P > 0) and limited by voltage (∆V < 0). In this case, one notes that added salt increases the translocation velocity (ρ b ↑ v p ↑) and turns the velocity from negative to positive. For an analytical insight into the enhancement of polymer capture by added salt, we Taylor-expand Equation (30) in the high salt density regime κa 1 and κd The screening parameter κ in the first term of Equation (34) indicates that the salt-induced growth of the polymer velocity in Figure 3a results from the enhanced screening of the EP drift force opposing the polymer capture.  In the opposite case of voltage-pressure driven (∆V > 0) and pressure-limited translocation (∆P < 0) displayed in Figure 3b, the salt dependence of the translocation velocity is non-monotonic. More precisely, in the strong salt regime ρ b > 0.1 M, the salt-screening of the EP polymer mobility in Equation (34) is seen to reduce the translocation velocity (ρ b ↑ v p ↓) and switch its sign from positive to negative. The characteristic salt density for polymer trapping follows from Equation (34) as In accordance with Figure 3a,b, Equation (35) predicts a drop of the trapping density with an increasing magnitude of the pressure gradient, |∆P| ↑ ρ * > ↓. We focus now on the dilute salt regime ρ b < 0.1 M of Figure 3b where one notes the enhancement of the translocation velocity with added salt ρ b ↑ v p ↑ and the presence of a second characteristic salt density where the velocity vanishes and the polymer gets trapped. For an insight into these features, we Taylor-expand Equation (30) in the dilute salt regime κa 1 and κd where the expansion coefficients are introduced as Equation (36) shows that in the dilute salt regime of Figure 3b, the negative translocation velocity corresponding to the polymer rejection regime is induced by repulsive polymer-membrane interactions (the third term on the r.h.s.). The screening of these interactions by added dilute salt results in the rise of the translocation velocity (ρ b ↑ v p ↑) and the reversal of its sign from negative to positive. The characteristic dilute salt density for polymer trapping follows from Equation (36) as In agreement with Figure 3b, Equation (38) predicts the increase of the lower characteristic salt density with increasing magnitude of the negative pressure, i.e., |∆P| ↑ ρ * > ↑. Finally, we consider the effect of the polymer length on the translocation dynamics. According to Equation (30), the reduction of the polymer length enhances the repulsive barrier term and reduces the polymer velocity, i.e., L p ↓ v p ↓. The slowing down of the translocation by finite polymer length results from the competition between the externally applied drift force and repulsive polymer-pore interactions; the pressure-voltage-induced drift acts on the entire polymer of length L p , while the electrostatic barrier originates solely from the polymer portion enclosed by the pore. As a result, the net drag force on the translocating polymer drops with the length of the molecule. Due to this balance, the polymer velocity in Equation (30) decreases inversely proportional to the polymer length with the critical molecular length for polymer trapping The characteristic length L * p is plotted in Figure 3c. First, one notes that the competition between the streaming current and the repulsive barrier leads to the decay of the characteristic length in Equation (40) with pressure, i.e., ∆P ↑ L * p ↓. As illustrated in the same figure, the dilute salt expansion of Equation (40) indicates that due to the same balance between the drift force and the electrostatic barrier, the critical length also drops with added salt, i.e., ρ b ↑ L * p ↓. In the next section, we investigate the deviation from the MF polymer transport behavior studied herein by added polyvalent cations.

Correlation-Induced DNA Mobility Inversion by Polyvalent Counterions in Solid-State Pores
In this section, we reconsider the polyvalent cation-induced DNA mobility inversion observed by the experiments of Ref. [30] and theoretically investigated in Ref. [29]. The multivalency of counterions requires the inclusion of charge correlations to the MF potential obtained from Equation (13). The details of this electrostatic one-loop (1l) level correlation-corrected scheme can be found in Ref. [29]. Figure 4 illustrates the polymer mobility µ p = v dr /E against the concentration of quadrivalent spermine (Spm 4+ ) molecules in the NaCl+SpmCl 4 solution of two different NaCl density. The plots compare the theoretical result = µ p = µ e [φ(d) − φ(a)] obtained from Equation (18) (solid curves) with the experimental dynamic light scattering (DLS) and single molecule electrophoresis (SME) experiments of Ref. [30] (squares). The nanopore and polymer charge densities given in the caption are the free parameters of the model that were adjusted to give the best agreement with the experimental data. The pore radius was in turn fixed to the value d = 10 nm corresponding to the characteristic radial size of solid-state nanopores.  Dynamic light scattering (DLS) and single molecule electrophoresis (SME) data of Ref. [30]. ds-DNA molecule has radius a = 1 nm and effective surface charge density σ p = −0.12 e/nm 2 . The nanopore has radius d = 10 nm and fixed surface charge density σ m = −0.006 e/nm 2 . The results are from Ref. [29].
As expected from MF electrophoresis, the dilute Spm 4+ regime of Figure 4 is characterized by a positive DNA mobility µ p > 0 corresponding to the drift of the negatively charged polymer oppositely to the external electric field E, i.e., from the cis to the trans side of the membrane (see Figure 1). However, the increment of the Spm 4+ concentration reduces the DNA mobility and switches its sign from positive to negative, indicating the reversal of the DNA translocation velocity from the cis-trans to the trans-cis direction. This corresponds to a non-MF charge transport picture where the anionic molecule moves parallel with the applied field E. Moreover, the comparison of the top and bottom figures indicates that added monovalent salt weakens charge correlations and rises the DNA mobility (ρ b+ ↑ µ p ↑) and the characteristic spermine density ρ * b4+ for mobility inversion (ρ b+ ↑ ρ * b4+ ↑). Within the experimental scattering, the theory can account for these features with reasonable quantitative accuracy.
The electrohydrodynamic mechanism driving the DNA velocity reversal is illustrated in Figure 5 where we plot the cumulative charge density (top plots) defined as and the convective liquid velocity profile obtained from Equation (15) (bottom plots). Equation (42) corresponds to the net charge of the DNA and its counterions. ρ c (r) is the local mobile charge density of the PB Equation (13). In order to emphasize first the role played by electrophoresis only, in Figure 5a,b, the EO flow was switched off by considering a neutral membrane (σ m = 0). In the dilute spermine regime with ρ b4+ = 0.1 mM (blue curve in Figure 5a), as one approaches the pore wall from the DNA surface, the gradual screening of the DNA charges by the counterions leads to a decreasingly negative total charge density Q cum (r) ≤ 0. This net negative charge coupled to the external field E results in the motion of the DNA and its counterions along the positive z axis (see Figure 1), i.e., u c (r) ≥ 0 and v dr = u c (a) > 0 (blue curve in Figure 5b).  Figure 4. The results are from Ref. [29].
In the larger spermine concentration regime ρ b4+ = 0.6 mM (black curves) and 1.0 mM (red curves), beyond the characteristic distance of ∼1 nm from the DNA surface, electrostatic correlations enhanced by the multivalency of Spm 4+ molecules switch the cumulative charge from negative to positive, indicating the occurrence of DNA charge inversion (CI). Consequently, in the vicinity of the DNA molecule, the charged liquid changes its direction and flows parallel with the external field E, i.e., u c (r) < 0. One notes, however, that at those Spm 4+ densities where CI is not strong enough to invert the electrophoretic force on DNA, the molecule continues to translocate opposite to the external electric field E, i.e., v dr = u c (a) > 0. Upon further increase of the Spm 4+ concentration to the critical value ρ * b4+ = 2.0 mM (purple curves), stronger charge correlations amplify the inverted charge density. As a result, the hydrodynamic drag by the charge inverted liquid on the DNA surface takes over the electric force on the DNA charges, resulting in the reversal of the DNA velocity from positive to negative (v dr < 0) and the direction of the molecule from the trans to the cis side.
These results show that the DNA mobility reversal is driven by a strong enough DNA charge inversion. The additional effect of the EO flow drag on this peculiarity is displayed in Figure 5c,d where we included the finite membrane charge density in Figure 4. The comparison of the left and right plots shows that the attraction of counterions by the anionic membrane charge amplifies the positive liquid density Q cum (r). As the corresponding EO flow positively adds to the hydrodynamic drag force exerted by the charge inverted liquid, the characteristic spermine density for DNA velocity reversal drops with the membrane charge, i.e., |σ m | ↑ ρ * b4+ ↓. Next, we study the transport properties of biological nanopores where the strong pore confinement results in pronounced correlation effects even in monovalent electrolytes.

Polymer Conductivity of Biological Pores: Image-Charge Barrier against Drift Force
In this part, we investigate the voltage-driven polymer transport properties of αHL channels where the strong pore confinement and the low membrane permittivity ε m ≈ 2 ε w gives rise to an image-charge barrier opposing the drift force on the polymer. The nanopore of radius d = 8.5 Å and length L m = 5 nm contains a monovalent KCl salt with bulk concentration ρ ±b = ρ b . The pore confines as well a ss-DNA molecule with radius a ≈ 5 Å and linear charge density τ ≡ −2πaσ p = 0.29 e/Å. We first focus on the experimentally observed rapid rise of polymer translocation rates with added salt. Within a phenomenological approach, this effect was explained in Ref. [7] by image-charge interactions. Within our translocation model, we intend to bring an analytical explanation to this peculiarity. αHL pores are characterized by a non-uniform surface charge distribution with alternating sign [38,39]. Thus, we assume a vanishing average charge density and take σ m = 0. Figure 6a illustrates the polymer translocation rates R c versus the salt concentration in the reservoir for various polymer charge density values. At low ion densities, R c is vanishingly small. Above a critical density ρ * b , R c grows sharply and converges towards the drift velocity v dr . One also sees that the capture of polymers with stronger charge occurs at higher salt concentrations, i.e., |τ| ↑ ρ * b ↑. In order to understand the physical mechanism behind these features, we plot in Figure 6b,c the polymer potential and velocity profiles. At the salt concentration ρ b = 0.04 M where the pore rejects the polymer (purple symbol in Figure 6a), the electrostatic barrier experienced by the molecule reaches the considerably high value of V p (z p )/L p ≈ 9 k B T/nm. To shed light on the origin of this barrier, we note that in a neutral pore where φ m (r) = 0, the electrostatic interaction energy reduces to the polymer self-energy that can be expressed as the following Fourier integral [29], where we introduced the dielectric jump function with the screening parameter p b = κ 2 + q 2 , the dielectric contrast factor γ = ε m /ε w , and the modified Bessel functions K n (x) [40]. The grand potential of Equation (43) giving rise to the electrostatic barrier V p (z p ) corresponds to the interaction energy of the polymer with its electrostatic image. At the dilute salt concentration ρ b = 0.04 M where the highly repulsive image-charge potential V p (z p ) dominates the drift term of Equation (20), the polymer potential U p (z p ) exhibits a minimum followed by an uphill trend and a barrier at z p = L m (dashed purple curve in Figure 6b). Due to this barrier, during the polymer capture regime z p < L m , the polymer velocity v p (z p ) = −βDU p (z p ) drops and switches from positive to negative (purple curve in Figure 6c). The change of the velocity sign indicates polymer trapping by the image-charge barrier at the pore entrance. The system is located in the barrier-driven regime.
In strong salt conditions κl p 1 and κd 1, the polymer grand potential of Equation (43) can be approximated by In Figure 6b, one notes that due to the exponential screening of the image-charge barrier Formula (45), the increment of the salt concentration from ρ b = 0.04 M to 0.1 M reduces the electrostatic potential from V p (z p )/L p ≈ 9 k B T/nm to ≈ 3 k B T/nm. As a result, the drift force on the polymer takes over the electrostatic barrier and the polymer potential U p (z p ) turns to downhill. The polymer is now in the drift-driven regime characterized by Equation (11). Figure 6c shows that this leads to a purely positive velocity, indicating the successful polymer capture and translocation (see also the orange symbol in Figure 6a).
These results indicate that the sharp rise of the polymer capture rates by salt addition originates from the competition between the image-charge barrier and the drift force. The same competition can indeed allow to understand the turnover in the voltage dependence of experimental polymer translocation rates in αHL channels [5,6,27]. Figure 7 illustrates this peculiarity at various salt concentration values. In agreement with the experimental curves of Ref. [41], the translocation rates rise exponentially at low voltages but grow with a weaker slope beyond a crossover voltage ∆V * . One also notes that added salt reduces this critical voltage, i.e., ρ b ↑ ∆V * ↓. The transition in the voltage dependence of the capture rates can be explained in terms of the potential profile U p (z p ) displayed in the inset of Figure 7. At the voltage ∆V = 150 mV (purple curve) located in the exponentially rising regime of the R c − ∆V curves (purple symbols) the image-charge barrier results in a potential trap to be escaped by thermal fluctuations. At the higher voltage ∆V = 250 mV located in the regime where R c increases linearly and gets close to the drift velocity (orange dots), the enhanced drift force takes over the repulsive image-charge barrier and the potential U p (z p ) turns to downhill (orange curve in Figure 7). This indicates that the non-uniform voltage dependence of the translocation rates is a consequence of the transition from the barrier to drift-dominated polymer transport regime.

Limitation of the Stiff Polymer Approximation
The main approximation of the electrohydrodynamic translocation theory presented herein is the modeling of the polymer as a rigid rod. We have shown that despite this approximation, accurate and detailed modeling of the pore electrohydrodynamics enables the theory to quantitatively explain several results obtained in many translocation experiments. This said, because the entropic cost of polymer translocation originating from conformational polymer fluctuations becomes important beyond the DNA persistence length L p 50 nm, the rigid rod approximation limits the quantitative predictive power of the electrohydrodynamic theory to polymer sequences whose size is comparable with the length of solid-state pores. Thus, if one wishes to consider the translocation of long polymer sequences, the inclusion of conformational polymer fluctuations becomes unavoidable. We are currently working in this direction.
The theoretical or even numerical consideration of electrohydrodynamic forces on a fluctuating polymer presents itself as an almost untractable task. However, the rigid rod approximation can be relaxed in the opposite regime of polymers much longer than the nanopore. In this configuration where the force induced by the pore electrohydrodynamics on the polymer can be considered to be local, one can absorb the electrohydrodynamic forces on the DNA into an effective external force f exerted solely on the polymer portion located in the pore, and the effective pore friction η p . This simplification allows to bypass the detailed description of the pore electrohydrodynamics, thereby enabling the accurate consideration of the polymer conformations bringing a major contribution to the translocation of long polymers. The next section of our article is devoted to this type of configurational translocation approach called the iso-flux tension propagation theory.

Iso-Flux Tension Propagation (IFTP) Theory for the Translocation of Long Polymers
This section is devoted to the tension propagation theory of polymer translocation through a nanopore. First, the theoretical model is introduced. Then to show the validity of the tension propagation theory the dynamics of the polymer translocation process is examined at the monomer level by looking at the waiting time distribution, that is the time each bead spends at the pore during the course of translocation. In the next subsection the scaling form of the translocation time τ 2 (cf. Equation (29)), which is the time that the chain needs to completely pass through the nanopore, is obtained for both pore-driven and end-pulled cases. We note that the theory presented here does not include any specific capture or escape processes, but assumes that the translocation starts with the pore being already filled and and stops when the cis side has no monomers left. This corresponds to the assumption that the pore thickness L m L p . Finally, in the last subsection we discuss the application of the theory to semi-flexible and rodlike polymers.

Coarse-Grained Polymer Model
Following our previous works, in this section we denote the polymer contour length by N 0 , and the translocation time τ = τ 2 , since in the theory here τ 1 = τ 3 = 0 corresponding to the thin pore approximation without any specific capture or trapping processes. For brevity, dimensionless units denoted by tilde are used asZ ≡ Z/Z u , with the units of length s u ≡ a, time t u ≡ ηa 2 /(k B T), Here a is the segment length, T is the temperature of the system, k B is the Boltzmann constant, and η is the friction of the solvent per monomer. Variables without tilde are expressed in Lennard-Jones units (for details see Refs. [20,21]).
During the process of polymer translocation the driving force may either act on the monomer(s) inside the pore (pore-driven case) or on the head monomer of the polymer (end-pulled case). For both pore-driven and end-pulled cases when the driving force is switched on a tension front starts to propagate along the backbone of the chain. Consequently, the cis part of the chain can be divided into two parts, mobile and immobile ones (see Figure 8a,c) [16]. Indeed, the part of the chain which experiences tension is mobile and has non-zero net velocity and the rest of the chain is in an immobile equilibrium state with zero average velocity. In the pore-driven case the velocity of the mobile part is towards the pore (see Figure 8a) while for the end-pulled case it is in the direction of the driving force (see Figure 8c). The boundary between the mobile and immobile parts is called the tension front that is located on the cis side for both pore-driven and end-pulled cases (see Figure 8). Both processes comprise two stages, the tension propagation (TP) and post propagation (PP) ones. During the TP stage the tension has not reached the chain end (see Figure 8a,c), while in the PP stage the whole chain has already been influenced by the tension (see Figure 8b,d). i.e.,t <τ TP , for the SS regime. The external driving forcef acts only on the monomer(s) at the pore towards the trans side. N 0 is the contour length of polymer, ands is the number of segments that have already been translocated into the trans side.l +s is the number of beads influenced by the tension force which is less than N 0 during the TP stage.R denotes the location of the tension front. (b) Translocation process for the SS regime in the post propagation (PP) stage when the tension has reached the chain end and after that, i.e.,l +s = N 0 . (c) The same as (a) but for the end-pulled polymer translocation process, where the external driving force acts only on the head monomer in the direction perpendicular to the membrane from cis to trans side. (d) The same as (b) but for the end-pulled case.
The shape of the mobile subchain depends on the strength of the driving force. For the pore-driven case in the limit of weak (N −ν 0 f 1) and moderate forces (1 f N ν 0 ) the mobile subchain is reminiscent of the trumpet (TR) and stem-flower (SF) configurations, respectively, while in the limit of very strong force (N ν 0 f ), in the strong stretching (SS) regime, the mobile subchain is fully straightened [13,19,42,43]. Here, N 0 is the contour length of the polymer, ν is the Flory exponent which is 3/4 and 0.588 for excluded volume chains in a good solvent is 2D and 3D, respectively, andf is the external driving force that acts on the monomers inside the pore for the pore-driven case or on the head monomer for the end-pulled case as depicted in Figure 8. For the end-pulled case as the mobile part extends to the cis and the trans sides, the dynamics of the chain is more complicated than that of the pore-driven case. For the cis side mobile subchain the same scenario as for the pore-driven case is valid here, but instead of the driving force,f , one needs to measure the value of the mediated tension force at the pore,f p . Thus, the shape of the mobile subchain in the cis side fits into the TR and SF regimes if (N 0 −s) −ν f p 1 and 1 f p (N 0 −s) ν , respectively, wheres is the translocation coordinate that is the length of the subchain in the trans side. The cis side mobile subchain is fully straightened if (N 0 −s) ν f p . In the other hand, for the end-pulled case the trans side mobile subchain shape can be either TR, SF, or fully straightened according to the strength of the driving force. For example, the mobile subchain in the trans side is fully straightened if N 0 f [20].
To study polymer translocation through a nanopore, similar to Refs. [13,19,43] the basic framework of Brownian dynamics (BD) in the overdamped limit is employed. According to BD the equation of motion for the translocation coordinates is written as whereΓ(t) is the total effective friction, γ is the surface exponent which is γ ≈ 0.95 and ≈ 0.69 for self-avoiding chains in 2D and 3D, respectively, and γ = 0.5 for ideal chains,ζ(t) is Gaussian white noise that satisfies ζ(t) = 0 and ζ(t)ζ(t ) = 2Γ(t)k B Tδ(t − t ), andf tot is the total force. The effective friction can be written as a sum of the friction due to the mobile part of the chain and the pore frictions. For the pore-driven case the effective friction isΓ(t) =η cis (t) +η p , whereη cis (t) is the friction due to the movement of mobile subchain in the solvent. For the end-pulled case the effective friction is written asΓ(t) =η cis (t) +η TS (t) +η p , whereη cis (t) andη TS (t) are frictions due to the movement of the mobile parts of the chain in the cis and in the trans sides inside the solvent, respectively. The index TS is an abbreviation for the trans side. It should be mentioned that for the pore-driven case of a flexible chain the dynamical trans side friction can be adsorbed into the pore friction as it just contributes a constant factor to it [13,43,44]. We also note that the term proportional to 1 − γ arises from the equilibrium entropy of the chain (for a fixeds) and is small enough to be neglected in the SS regime. It should be noted that it is not fully consistent with the propagation of the tension front on the cis side even for the pore-driven case. Equation (46) gives the time evolution of the translocation coordinate,s, provided that the effective friction,Γ(t), is known. Indeed, the physics of tension propagation theory is embedded inΓ(t). To find the effective friction that is the combination of the mobile subchain and pore frictions, one needs to find the time evolution of the tension front, which gives the dynamics of the friction due to the mobile subchain. To this end, similar to Ref. [42], we assume that the flux of monomers in the mobile domain, φ = ds/dt, is constant in space but evolves in time (the iso-flux assumption). The tension front is located at the distancex = −R from the pore on the cis side. The tension force at the distancex from the pore is obtained by integrating the force-balance equation for a differential element dx located betweenx andx + dx, asf (x,t) =f 0 −φ(t)x. For the pore-driven case the integration is performed from pore tox andf 0 ≡f tot −η pφ (t), while for the end-pulled case the integration is from head monomer to pore and then from pore tox andf 0 ≡f tot −η pφ (t) −η TSφ (t). Here, in the SS regimẽ η TS =s for the end-pulled case. Indeed, by integration of the force balance equation over the whole mobile domain together with the definition of the tension front, where the tension force vanishes, the equation of motion for the monomer flux is written as Then the effective friction is expressed by using Equations (46) and (47) and the definition of the monomer flux,φ ≡ ds/dt, as Γ(t) =R(t) +η p , pore-driven; Γ(t) =R(t) +η p +η TS , end-pulled.
Time evolution ofs is given by Equations (46)-(48), but to have a full solution one still needs the equation of motion for the location of the tension front,R, which is obtained in the TP and PP stages separately. In the TP stage where the tension front has not reached the chain end, for flexible chain one can writeR = A ν N ν , where N =l +s is the number of segments that already influenced by the tension (see Figure 8a) andl is the number of segments in the mobile domain in the cis side. Here we only present the time evolution of the tension front for the strong stretching (SS) regime, where the force is very strong, andl =R. To study the other TR and SF regimes a similar procedure is employed [19]. InsertingÑ inside the equation above forR and performing time derivation, the equation of motion for the tension front for TP stage is obtained aṡR In the PP stage the tension force has already reached the chain and therefore N =l +s = N 0 . By substitutingl =R (in the SS regime) in the above relation and taking the time derivative, the equation for the time evolution of the tension front is written aṡR The time evolution of the tension front for the end-pulled case in the SS regime for TP and PP stages is the same as of pore-driven case, i.e., Equations (49) and (50).
To find a full solution of the iso-flux tension propagation (IFTP) model for the TP stage Equations (46)-(49) must be consistently solved, while for the PP stage one needs to solve Equations (46)-(48) and (50). It should be mentioned that to improve the quantitative accuracy of the IFTP theory the distribution of the initial configurations of the chain can be incorporated into the model through A ν inR = A ν N ν . The details are in Ref. [19]. Moreover, modified versions of IFTP theory have been employed to study the translocation of a semi-flexible or stiff polymer through a nanopore [21] (to be discussed later) as well as polymer translocation through a flickering nanopore under an alternating driving force [45].

Waiting Time Distribution
The waiting time (WT), which is the time that each bead spends at the pore during the course of translocation, is an important quantity that can reveal the dynamics of the process at the monomer level. Figure 9a shows the WT as a function ofs, the translocation coordinate, for the pore-driven polymer translocation. The chain length is N 0 = 128, the external driving force at the pore is f = 5, and the pore friction in the IFTP theory is η p = 3.5. The black curve presents the deterministic case. The red up triangles show the WT when force is chosen randomly but A ν = 1.15 is deterministic. Green left triangles are devoted to the WT when both the force and A ν are stochastic, and finally MD simulation data are shown in blue squares. The stochastic sampling of the initial configurations of the chain smoothens the transition from the TP to the PP stage. In MD simulations (blue squares) the same feature is seen where the initial configurations are sampled by thermalizing the chain before each simulation trajectory. It is clear there is a very good quantitative agreement between the stochastically augmented IFTP theory and MD simulations. Panel (b) is the same as panel (a) but for the end-pulled case. Black curve is for deterministic case while the blue squares show the MD simulation data. Here, the chain length is N 0 = 100, the external driving force is f = 100, and the pore friction in the IFTP theory is η p = 3.
For both the pore-driven and end-pulled cases, it can be seen in Figure 9 that the translocation process is a far-from-equilibrium process in the sense that the conformations of the chain do not correspond to linear response or quasi-equilibrium ones. For the pore-driven case, Figure 9a, in the TP stage, where the tension is still propagating along the backbone of the chain, the number of mobile monomers in the cis side is increasing. Therefore the friction is growing and consequently WT increases monotonically until it gets its maximum when the tension reaches the end of the chain. Then in the second PP stage as the time passes the number of mobile monomers in the cis side decreases, which means the friction due to the mobile part of the chain decreases too, and WT decreases. For the end-pulled case, Figure 9b, in the TP stage similar to the pore-driven case WT increases. However, in the PP stage WT is almost constant. This is because in the SS regime both subchains in the cis and in the trans sides are fully straightened (rodlike) therefore the friction due to the movement of the whole chain in the solvent remains approximately constant. The same as (a) but for the end-pulled case. Black curve is for deterministic case while the blue squares represent the MD simulation data. Here, the chain length is N 0 = 100, the external driving force is f = 100, and the pore friction in the IFTP theory is η p = 3.

Scaling of the Translocation Time for a Flexible Polymer
A fundamental quantity characterizing the polymer translocation process is the translocation time, which is the time that a chain needs to pass through the nanopore. Here, we consider only the case where τ = τ 2 and show how the scaling form of the translocation time can be extracted analytically from the IFTP theory in the SS regime. The same approach can be applied for the TR and SF regimes [19][20][21]45,46].
To obtain an analytical form of the translocation time in the SS regime we use an approximation and we only take into account the contribution of the external driving force to the equation of motion for the translocation coordinate, i.e.,f tot ≈f . Then for the pore-driven case, combiningφ =f / η p +R with the definition of the monomer flux,φ = ds/dt, together with the mass conservation in the TP stage, N =l +s, by integration of N from 0 to N 0 the TP time reads asτ TP = ´N 0 0R (N)dN + η p N 0 /f − ∆τ, where ∆τ = η pR (N 0 ) +R 2 (N 0 )/2 /f . The PP time is obtained by integratingR fromR(N 0 ) to 0 asτ PP = ∆τ. At the end, the whole translocation time,τ =τ PP +τ TP , is written as τ = ´N 0 0R (N)dN +η p N 0 /f with the scaling form To obtain the scaling of the translocation time for the end-pulled case the same procedure is applied, but now the monomer flux isφ =f / η p +R +η TS , whereη TS is the trans side friction and for the SS regime isη TS =s. The whole translocation time for the end-pulled case is given bỹ τ = 1f ´N 0 0R (N)dN +η p N 0 +τ TS , whereτ TS = N 2 0 /(2f ) is the contribution of friction due to the fully straightened trans side subchain. Thus In Figure 10a we plot the translocation exponent α, which is defined asτ ∝ N α 0 , for the pore-driven polymer translocation based on the deterministic IFTP theory as a function of N 0 for various values of the pore friction η p = 1, 5 and 10. Here the external driving force is f = 5, and A ν = 1.15. As can be seen in the short chain limit the values of the effective translocation exponents for different pore friction coefficients are different due to a competition between the two terms in the right hand side of Equation (51). In the long chain limit where the first term due to the mobile subchain friction dominates, the values of the translocation exponents for different pore frictions are the same, as depicted in the figure by 1 + ν. To show the effect of the pore friction on the translocation exponent, the rescaled exponent curves are presented for different values of η p . As can be seen all of them collapse on a single curve as denoted by rescaled data in the figure, i.e., 1 + ν = 1.588. For the pore-driven case the rescaled exponent is defined asτ −η p N 0 /f ∝ N α † 0 and shows only the cis side friction contribution to the dynamics of the translocation process. Figure 10b is similar to Figure 10a but for the end-pulled case, and for various values of the pore friction η p = 1.5, 10 and 20. Here the external driving force is f = 100 acts on the head monomer of the polymer. The green solid line presents the rescaled translocation exponent α † , which is defined asτ −η p N 0 /f ∝ N α † 0 . The black dashed curve shows the rescaled exponent α ‡ = 2 with the definition ofτ − ´N

Scaling of the Translocation Time for a Stiff Polymer
In this subsection the scaling form of the translocation time for a stiff chain is briefly discussed for both pore-driven and end-pulled polymer translocation processes. The complete theory for semi-flexible chains can be found in Ref. [21]. The end-to-end distance of a very stiff chain (rod-like limit) is given byR(N) = N. For both pore-driven and end-pulled translocation processes of the stiff polymer the number of mobile monomers on the cis side isl =R, while on the trans side it is given bỹ s. As the chain is very stiff the total translocation time is much larger than the TP time, i.e.,τ τ tp , and therefore it is a very good approximation to ignore the TP stage. During the PP stage as the tension has already reached the chain end N =s +l = N 0 , and one sets the condition dN/dt = 0. The PP time, which is equal to the total translocation time, is obtained by integratingR fromR(N 0 ) to zero. Then, the total translocation time is written asτ =τ pp =´R N 0 0 dR R +η p +η TS (t)]/f . Knowingη TS (t) =s = N 0 −l together with the fact thatl =R, one can obtain the final scaling form of the total translocation time asτ Here the pore friction term,η p N 0 /f , has a significant correction to asymptotic scaling similar to the flexible case and the effective exponents intermediate values of N 0 will be between unity and two. We note that the scaling form of Equation (52) has recently been derived from an nonequilibrium transport theory in the limit of drift-driven polymer translocation where the electrostatic interactions are weak [32].

Summary and Conclusions
In this article, we have presented a comparative review of the electrohydrodynamic and configurational approaches to polymer translocation. The distinction between these two approaches is based on the relative spatial scales of the polymer and the nanopore. The first part of the article is devoted to the translocation of short polymers whose sizes are comparable to the length of the translocated nanopore. At this scale, an accurate modeling of the translocation process requires the explicit consideration of the electrostatic and hydrodynamic details of the pore medium, such as the electrophoretic and hydrodynamic drag forces on the translocating molecule, and the electrostatic interactions of the molecule with the membrane and the surrounding electrolyte.
First, in Section 2, we have presented the electrohydrodynamic model of polymer translocation where these details are explicitly and consistently included via the coupled solution of the Stokes, Poisson-Boltzmann, and polymer diffusion equations. We then discussed the application of the model to various experimental configurations. In Section 2.2, we presented direct comparisons with pressure-voltage-driven polymer trapping experiments carried out in monovalent salt solutions where the translocation process is governed by MF electrohydrodynamics. For this experimental setup, we showed that the electrohydrodynamic theory can quantitatively reproduce and explain the pressure dependence of the experimental polymer translocation velocity and time data.
Then, in Section 2.3, we focused on polymer translocation experiments conducted in polyvalent electrolytes where the high ion valency drives the system out of the MF electrohydrodynamic regime. We showed that via the inclusion of charge correlations, the electrohydrodynamic theory can predict with quantitative accuracy the experimentally observed inversion of the electrohydrodynamic DNA mobility by added multivalent cations. We elucidated the electrohydrodynamic mechanism underlying the mobility reversal in terms of the DNA charge inversion. Finally, in Section 2.4, we considered translocation experiments conducted with αHL nanopores where the high pore confinement results in strong correlation effects even with monovalent electrolytes. Within the framework of the model, we showed that the experimentally observed salt-induced sharp rise and non-uniform voltage dependence of polymer capture rates can be explained by the competition between the electric drift force and surface polarization effects acting as a repulsive barrier for polymer capture.
Section 3 was devoted to the opposite case of polymers much longer than the nanopore thickness. In this regime, as the pore electrohydrodynamics acts only on a small portion of the polymer, one can make the approximation to absorb the electrohydrodynamic details of the nanopore into an effective force f acting locally on the polymer portion confined to the pore, and the effective (constant) pore friction η p . This simplification allows to bypass the electrohydrodynamic details of the translocation process, enabling accurate modeling of the configurational effects originating from conformational polymer fluctuations. Within the corresponding IFTP theory, we discussed the detailed characterization of the translocation dynamics of polymers with arbitrary length. We presented the predictions of the IFTP theory for the scaling of the polymer translocation time with the polymer length and the variation of the former with the pore friction and the external force driving the polymer. The theory is applicable to a variety of translocation and polymer pulling scenarios, including pore-driven and end-pulling setups discussed here, and is in excellent agreement with MD simulations of the corresponding coarse-grained polymer models.
The grand challenge in polymer translocation consists of amalgamating the regimes of short and long polymers so far separately studied. This requires the consideration of electrohydrodynamic effects and conformational polymer fluctuations on the same footing. In our recent work of Ref. [32], a first attempt in this direction was made by incorporating the electrostatic coupling of the membrane with the cis and trans portions of the polymer outside the nanopore into the stiff polymer limit of the IFTP theory. We are currently working on the relaxation of the stiff polymer constraint. This extension will hopefully allow to better understand the effect of electrostatic interactions on the scaling of the translocation time with the polymer length.