Physical Chemistry

A thermalized electrokinetics model including stochastic reactions suitable for multiscale simulations of reaction-advection-diffusion systems



We introduce a scheme to simulate the spatial and temporal evolution of the densities of charged species, taking into account diffusion, thermal fluctuations, coupling to a carrier fluid, and chemical reactions. To this end, the diffusive fluxes in the electrokinetic model by Capuani et al. [1] are supplemented with thermal fluctuations. Chemical reactions are included via an additional source term in the mass balance equation. The diffusion-reaction model is then coupled to a solver for fluctuating hydrodynamics based on the lattice Boltzmann method. This combination is particularly useful for soft matter simulations, due to the ability to couple particles to the lattice-Boltzmann fluid. These could, e.g., be charged colloids or polymers, which then interact with an ion distribution. We describe one implementations based on the automatic code generation tools pystencils and lbmpy, and another one that is contained in the molecular dynamics package ESPResSo and that allows for an easy coupling of particles to the density fields. We validate our implementations by comparing to several known analytic results. Our method can be applied to coarse-grained catalysis problems as well as to many other multi-scale problems that require the coupling of explicit-particle simulations to flow fields, diffusion, and reaction problems in arbitrary geometries.

Version notes

Better description of the used LB. Explained the solving of the Poisson equation via Fourier transformation. Added numerical works showing the grid independence. Added numerical works towards the reaction-diffusion-advection channel.


Thumbnail image of chemRXive_V3.pdf

Supplementary material

Thumbnail image of fluctuating_reactive_ek.ipynb
Jupyter notebook to demonstrate the themalized reaction-advection-diffusion algorithm which is being implemented using pystencils. To fully run this script the other 3 files ('', 'ideal_gas_md_hist.dat' and 'coulomb_gas_md_hist.dat') are required.
Thumbnail image of
A script to solve the poisson equation using pystencils. This script is being imported in the 'fluctuating_reactive_ek.ipynb' notebook.
Thumbnail image of ideal_gas_md_hist.dat
Density probability histogram data of an ideal gas MD simulation. This data is being compared to the data obtained by the 'fluctuating_reactive_ek.ipynb' notebook.
Thumbnail image of coulomb_gas_md_hist.dat
Density probability histogram data of an coulomb gas MD simulation. This data is being compared to the data obtained by the 'fluctuating_reactive_ek.ipynb' notebook.