A Brownian dynamics-finite element method for simulating DNA electrophoresis in nonhomogeneous electric fields

The objective of this work is to develop a numerical method to simulate DNA electrophoresis in complicated geometries. The proposed numerical scheme is composed of three parts: (1) a bead-spring Brownian dynamics (BD) simulation, (2) an iterative solver-enhanced finite element method (FEM) for the electric field, and (3) the connection algorithm between FEM and BD. A target-induced searching algorithm is developed to quickly address the electric field in the complex geometry which is discretized into unstructured finite element meshes. We also develop a method to use the hard-sphere interaction algorithm proposed by Heyes and Melrose [J. Non-Newtonian Fluid Mech. 46, 1 (1993)] in FEM. To verify the accuracy of our numerical schemes, our method is applied to the problem of lambda-DNA deformation around an isolated cylindrical obstacle for which the analytical solution of the electric field is available and experimental data exist. We compare our schemes with an analytical approach and there is a good agreement between the two. We expect that the present numerical method will be useful for the design of future microfluidic devices to stretch and/or separate DNA.