Simulations of electroosmotic flow in charged nanopores using Dissipative Particle Dynamics with Ewald summation

Spiers Memorial Lecture: Towards understanding of iontronic systems: electroosmotic flow of monovalent and divalent electrolyte through charged cylindrical nanopores

In many practical applications, ions are the primary charge carrier and must move through either semipermeable membranes or through pores, which mimic ion channels in biological systems. In analogy to electronic devices, the "iontronic" ones use electric fields to induce the charge motion. However, unlike the electrons that move through a conductor, motion of ions is usually associated with simultaneous solvent flow. A study of electroosmotic flow through narrow pores is an outstanding challenge that lies at the interface of non-equilibrium statistical mechanics and fluid dynamics. In this paper, we will review recent works that use dissipative particle dynamics simulations to tackle this difficult problem. We will also present a classical density functional theory (DFT) based on the hypernetted-chain approximation (HNC), which allows us to calculate the velocity of electroosmotic flows inside nanopores containing 1 : 1 or 2 : 1 electrolyte solution. The theoretical results will be compared with simulations. In simulations, the electrostatic interactions are treated using the recently introduced pseudo-1D Ewald summation method. The zeta potentials calculated from the location of the shear plane of a pure solvent are found to agree reasonably well with the Smoluchowski equation. However, the quantitative structure of the fluid velocity profiles deviates significantly from the predictions of the Smoluchowski equation in the case of charged pores with 2 : 1 electrolyte. For low to moderate surface charge densities, the DFT allows us to accurately calculate the electrostatic potential profiles and the zeta potentials inside the nanopores. For pores with 1 : 1 electrolyte, the agreement between theory and simulation is particularly good for large ions, for which steric effects dominate over the ionic electrostatic correlations. The electroosmotic flow is found to depend very strongly on the ionic radii. In the case of pores containing 2 : 1 electrolyte, we observe a reentrant transition in which the electroosmotic flow first reverses and then returns to normal as the surface change density of the pore is increased.

Widom insertion method in simulations with Ewald summation

We discuss the application of the Widom insertion method for calculation of the chemical potential of individual ions in computer simulations with Ewald summation. Two approaches are considered. In the first approach, an individual ion is inserted into a periodically replicated overall charge neutral system representing an electrolyte solution. In the second approach, an inserted ion is also periodically replicated, leading to the violation of the overall charge neutrality. This requires the introduction of an additional neutralizing background. We find that the second approach leads to a much better agreement with the results of grand canonical Monte Carlo simulation for the total chemical potential of a neutral ionic cluster.

Reversal of Electroosmotic Flow in Charged Nanopores with Multivalent Electrolyte

We study the reversal of electroosmotic flow in charged cylindrical nanopores containing multivalent electrolyte. Dissipative particle dynamics is used to simulate the hydrodynamics of the electroosmotic flow. The electrostatic interactions are treated using 3D Ewald summation, corrected for a pseudo-one-dimensional geometry of the pore. We observe that, for sufficiently large surface charge density, condensation of multivalent counterions leads to the reversal of the pore's surface charge. This results in the reversal of electroosmotic flow. Our simulations show that the Smoluchowski equation is able to quantitatively account for the electroosmotic flow through the nanopore, if the shear plane is shifted from the position of the Stern contact surface.

Electroosmotic Flow in Polarizable Charged Cylindrical Nanopores

We present a simulation method to study electroosmotic flow in charged nanopores with dielectric contrast between their interior and the surrounding medium. To perform simulations, we separate the electrostatic energy into the direct Coulomb and the polarization contributions. The polarization part is obtained using periodic Green functions and can be expressed as a sum of fast converging modified Bessel functions. On the other hand, the direct Coulomb part of the electrostatic energy is calculated using fast converging three-dimensional (3D) Ewald summation method, corrected for a pseudo one-dimensional (1D) geometry. The effects of polarization are found to be particularly important for systems with multivalent counterions and narrow nanopores. Depending on the surface charge density, polarization can increase the volumetric flow rate by 200%. For systems with 3:1 electrolyte, we observe that there is a saturation of the volumetric flow rate. In this case, for polarizable pores, the flow rate is 100% higher than for nonpolarizable pores.