![](https://d3ilqtpdwi981i.cloudfront.net/C_RDsobz3-dDWYtvV-NTKqPtjps=/425x550/smart/https://bepress-attached-resources.s3.amazonaws.com/uploads/00/b9/fa/00b9fa86-8d12-4556-8c64-d7b6d7e18b6a/thumbnail_69e2278a-3b1e-4724-950e-30c883148c65.jpg)
In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion"Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.
Available at: http://works.bepress.com/alejandro_garcia1/113/
This article originally appeared in Physical Review Fluids, volume 6, issue 4, 2021, published by the American Physical Society. ©2021 American Physical Society. The article can also be found online at this link.