Coupled molecular-dynamics and finite-element-method simulations for the kinetics of particles subjected to field-mediated forces

Phys Rev E. 2019 Jun;99(6-1):063307. doi: 10.1103/PhysRevE.99.063307.

Abstract

A computational approach that couples molecular-dynamics (MD) and the-finite-element-method (FEM) technique is here proposed for the theoretical study of the dynamics of particles subjected to electromechanical forces. The system consists of spherical particles (modeled as micrometric rigid bodies with proper densities and dielectric functions) suspended in a colloidal solution, which flows in a microfluidic channel in the presence of a generic nonuniform variable electric field generated by electrodes. The particles are subjected to external forces (e.g., drag or gravity) which satisfy a particlelike formulation that is typical of the MD approach, along with an electromechanical force that, in turn, requires the three-dimensional self-consistent solutions of correct continuum field equations during the integration of the equations of motion. In the MD-FEM method used in this work, the finite element method is applied to solve the continuum field equations while the MD technique is used for the stepwise explicit integration of the equations of motion. Our work shows the potential of coupled MD-FEM simulations for the study of electromechanical particles and opens a double perspective for implementing (a) MD away from the field of atomistic simulations and (b) the continuum-particle approach to cases where the conventional force evaluation used in MD is not applicable.