Driven polymer transport through a periodically patterned channel

Parallel computing for mobilities in periodic geometries

We examine methods for calculating the effective mobilities of molecules driven through periodic geometries in the context of particle-based simulation. The standard formulation of the mobility, based on the long-time limit of the mean drift velocity, is compared to a formulation based on the mean first-passage time of molecules crossing a single period of the system geometry. The equivalence of the two definitions is derived under weaker assumptions than similar conclusions obtained previously, requiring only that the state of the system at subsequent period crossings satisfy the Markov property. Approximate theoretical analyses of the computational costs of estimating these two mobility formulations via particle simulations suggest that the definition based on first-passage times may be substantially better suited to exploiting parallel computation hardware. This claim is investigated numerically on an example system modeling the passage of nanoparticles through the slit-well device. In this case, the traditional mobility formulation is found to perform best when the Péclet number is small, whereas the mean first-passage time formulation is found to converge much more quickly when the Péclet number is moderate or large. The results suggest that, given relatively modest access to modern GPU hardware, this alternative mobility formulation may be an order of magnitude faster than the standard technique for computing effective mobilities of biomolecules through periodic geometries.

Pore shapes effects on polymer translocation

We translocated polymers through pores of different shapes and interaction patterns in three dimensions by Langevin molecular dynamics. There were four simple cylindrical pores of the same length but with different diameters. The results showed that even though decreasing the pore diameter would always decrease the translocation velocity, it was strongly dependent on the shape of the increased pore diameter. Although increasing the pore diameter made the translocation faster in simple cylindrical pores, it was complicated in different pore shapes, e.g. increasing the diameter in the middle decreased the translocation velocity. Investigating polymer shapes through the translocation process and comparing the shapes by the cumulative waiting time for different pore structures reveals the non-equilibrium properties of translocation. Moreover, polymer shape parameters such as gyration radius, polymer center of mass, and average aspect ratio help us to distinguish different pore shapes and/or different polymers.

Simulation study on the migration of diblock copolymers in periodically patterned slits

The forced migration of diblock copolymers (A_NAB_NB) in periodically patterned slits was investigated by using Langevin dynamics simulation. The lower surface of the slit consists of stripe α and stripe β distributed in alternating sequence, while the upper one is formed only by stripe β. The interaction between block A and stripe α is strongly attractive, while all other interactions are purely repulsive. Simulation results show that the migration of the diblock copolymer is remarkably dependent on the driving force and there is a transition region at moderate driving force. The transition driving force f_t, where the transition region occurs, decreases monotonously with increasing length of block B (N_B) but is independent of the polymer length and the periodic length of the slit, which is interpreted from the free energy landscape of diblock copolymer migration. The results also show that periodic slits could be used to separate diblock polymers with different N_B by tuning the external driving force.

Simulation on the translocation of homopolymers through sandwich-like compound channels

The forced translocation of homopolymers through αβα sandwich-like compound channels was investigated by Monte Carlo simulation. The interaction between polymer and part α is strongly attractive, whereas that between polymer and part β is purely repulsive. Simulation results show that the translocation is influenced obviously by the length of part β (Lβ) and the starting position of part β (Lα1). For small Lβ, the translocation is mainly governed by the escaping process, and polymer is trapped near the exit of the channel. However, the translocation time can be tuned by varying Lα1 and the fastest translocation can be achieved at relatively large Lα1. Whereas for large Lβ and small Lα1, the translocation is mainly controlled by the filling process. It is difficult for polymer to enter the channel, and polymer is trapped at the first αβ interface. Finally, the dynamics for the filling process and the escaping process are discussed from the view of free-energy landscape, respectively.

Stretching of DNA confined in nanochannels with charged walls

There is currently a growing interest in control of stretching of DNA inside nanoconfined regions due to the possibility to analyze and manipulate single biomolecules for applications such as DNA mapping and barcoding, which are based on stretching the DNA in a linear fashion. In the present work, we couple Finite Element Methods and Monte Carlo simulations in order to study the conformation of DNA molecules confined in nanofluidic channels with neutral and charged walls. We find that the electrostatic forces become more and more important when lowering the ionic strength of the solution. The influence of the nanochannel cross section geometry is also studied by evaluating the DNA elongation in square, rectangular, and triangular channels. We demonstrate that coupling electrostatically interacting walls with a triangular geometry is an efficient way to stretch DNA molecules at the scale of hundreds of nanometers. The paper reports experimental observations of λ-DNA molecules in poly(dimethylsiloxane) nanochannels filled with solutions of different ionic strength. The results are in good agreement with the theoretical predictions, confirming the crucial role of the electrostatic repulsion of the constraining walls on the molecule stretching.