Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power-law dependence on the length, as is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, while the rate of translocated beads at each time step is linearly correlated with the number of resident beads in the pore. Finally, analysis of the statistics over the translocation work unravels the importance of the hydrodynamic interactions in the process. © 2009 IOP Publishing Ltd.

We investigate the process of biopolymer translocation through a narrow pore using a multiscale approach which explicitly accounts for the hydrodynamic interactions of the molecule with the surrounding solvent. The simulations confirm that the coupling of the correlated molecular motion to hydrodynamics results in significant acceleration of the translocation process. Based on these results, we construct a phenomenological model which incorporates the statistical and dynamical features of the translocation process and predicts a power-law dependence of the translocation time on the polymer length with an exponent alpha approximate to 1.2. The actual value of the exponent from the simulations is alpha=1.28 +/- 0.01, which is in excellent agreement with experimental measurements of DNA translocation through a nanopore, and is not sensitive to the choice of parameters in the simulation. The mechanism behind the emergence of such a robust exponent is related to the interplay between the longitudinal and transversal dynamics of both translocated and untranslocated segments. The connection to the macroscopic picture involves separating the contributions from the blob shrinking and shifting processes, which are both essential to the translocation dynamics.

We present a detailed description of biopolymer translocation through a nanopore in the presence of a solvent using an innovative multiscale methodology that treats the biopolymer at the microscopic scale as combined with a self-consistent mesoscopic description for the solvent fluid dynamics. We report evidence for quantized current blockade depending on the folding configuration and offer detailed information on the role of hydrodynamic correlations in speeding up the translocation process.

We present a multiscale approach to the modeling of polymer dynamics in the presence of a fluid solvent. The approach combines Langevin molecular dynamics ( MD) techniques with a mesoscopic lattice Boltzmann (LB) method for the solvent dynamics. A unique feature of the present approach is that hydrodynamic interactions between the solute macromolecule and the aqueous solvent are handled explicitly, and yet in a computationally tractable way due to the dual particle-field nature of the LB solver. The suitability of the present LB-MD multiscale approach is demonstrated for the problem of polymer fast translocation through a nanopore. We also provide an interpretation of our results in the context of DNA translocation through a nanopore, a problem that has attracted much theoretical and experimental attention recently.

We propose a discrete lattice version of the Fokker-Planck kinetic equation in close analogy with the lattice-Boltzmann scheme. Our work extends an earlier one-dimensional formulation to arbitrary spatial dimension D. A generalized Hermite-Gauss procedure is used to construct a discretized kinetic equation and a Chapman-Enskog expansion is applied to adapt the scheme so as to correctly reproduce the macroscopic continuum equations. The linear stability of the algorithm with respect to the finite time step Delta t is characterized by the eigenvalues of the collision matrix. A heuristic second-order algorithm in Delta t is applied to investigate the time evolution of the distribution function of simple model systems, and compared to known analytical solutions. Preliminary investigations of sedimenting Brownian particles subjected to an orthogonal centrifugal force illustrate the numerical efficiency of the Lattice-Fokker-Planck algorithm to simulate nontrivial situations. Interactions between Brownian particles may be accounted for by adding a standard Bhatnagar-Gross-Krook collision operator to the discretized Fokker-Planck kernel.

A lattice version of the Fokker-Planck equation is introduced. The resulting numerical method is illustrated through the calculation of the electric conductivity of a one-dimensional charged fluid at zero and finite-temperature.

A lattice version of the Fokker-Planck equation, accounting for dissipative interactions, not resolved on the molecular scale, is applied to the study of electrorheological transport of a one-dimensional charged fluid, and is found to yield quantitative agreement with a recent analytical solution.