Fully explicit dissipative particle dynamics simulation of electroosmotic flow in nanochannels

A fully explicit mesoscale simulation of electroosmotic flow (EOF) in nanochannels is presented by an extended dissipative particle dynamics (DPD) method. To avoid formation of ionic pairs through interacting soft-core charges, a Slater-type smearing distribution borrowed from quantum mechanics is utilized to surround each soft DPD ion with a charge cloud. To account for reduced periodicity normal to the walls direction, a corrected version of 3D Ewald sum is implemented in which a dipole moment term is deducted from energy and force terms of non-frozen charges. Simulation box is then elongated normal to walls to dampen spurious interslab interactions by adding vacuum gaps between periodic images. These measures together with the established unit conversions guarantee perfect match to molecular dynamics results. The transition of EOF velocity profile from parabolic (equivalent to overlap of electric double layers) to plug-like shapes is studied across the changing electric field between 0.06 and 0.41 [V/nm], and varying salt concentration from 0.26 to 2.0 [M]. It is found that 1.25 [V] increase in the driving voltage can potentially enhance the electroosmotic flow rate by 8-11 times in the range of ionic concentrations studied. The range of surface zeta potential calculated as 27 < zeta < 52 [mV] in the linear response regime, as identified to occur for 0.24 <= E [V/nm], agrees reasonably with numerical and experimental studies.