Abstract:The distribution of ions near a charged surface is an important quantity in many biological and material processes, and has been therefore investigated intensively. However, few theoretical and simulation approaches have included the influence of concentration-induced variations in the local dielectric permittivity of an underlying electrolyte solution. Such local variations have long been observed and known to affect the properties of ionic solution in the bulk and around the charged surface. We propose a hybrid computational model that combines Monte Carlo simulations with continuum electrostatic modeling to investigate such properties. A key component in our hybrid model is a semianalytical formula for the ion-ion interaction energy in a dielectrically inhomogeneous environment. This formula is obtained by solving for the Green's function Poisson's equation with ionic-concentration-dependent dielectric permittivity using a harmonic interpolation method and spherical harmonic series. We also construct a self-consistent continuum model of electrostatics to describe the effect of ionic-concentration-dependent dielectric permittivity and the resulting self-energy contribution. With extensive numerical simulations, we verify the convergence of our hybrid simulation scheme, show the qualitatively different structures of ionic distribution due to the concentration-induced dielectric variations, and compare our simulation results with the self-consistent continuum model. In particular, we study the differences between weakly and strongly charged surfaces and multivalencies of counterions. Our hybrid simulations conform particularly the depletion of ionic concentrations near a charged surface and also capture the charge inversion. We discuss several issues and possible further improvement of our approach for simulations of large charged systems.

Continuum electrostatics for ionic solutions with non-uniform ionic sizes

Continuum theory of electrostatic correlations at charged surfaces

The Calculus of Boundary Variations and the Dielectric Boundary Force in the Poisson-Boltzmann Theory for Molecular Solvation

Classical density functional analysis of the ionic size effects on the properties of charge regulating electric double layers

From solvent free to dilute electrolytes: Essential components for a continuum theory

Diffusive dynamics of charge regulated macro-ion solutions

Hybrid Monte Carlo and Continuum Modeling of Electrolytes with Concentration-Induced Dielectric Variations.

Influence of solvent polarization and non-uniform ion size on electrostatic properties between charged surfaces in an electrolyte solution

Theoretical and practical investigation of ion–ion association in electrolyte solutions

Self-consistent field model for strong electrostatic correlations and inhomogeneous dielectric media

On the equilibrium solutions of electro-energy-reaction-diffusion systems

The Role of Energy Scales for the Structure of Ionic Liquids at Electrified Interfaces: A Theory-Based Approach

Simple Theory of Ionic Activity in Concentrated Electrolytes

Interior Electroneutrality in Nernst–Planck–Navier–Stokes Systems

Thermodynamics of Electrolyte Solutions Near Charged Surfaces: Constant Surface Charge vs. Constant Surface Potential

Effective charges and virial pressure of concentrated macroion solutions

Free and Bound States of Ions in Ionic Liquids, Conductivity, and Underscreening Paradox

Ionic Interactions in Biological and Physical Systems: a Variational Treatment

Finite-Difference Approximations and Local Algorithm for the Poisson and Poisson-Boltzmann Electrostatics

Density-functional Theory of Spherical Electric Double Layers and Zeta Potentials of Colloidal Particles in Restricted-Primitive-model Electrolyte Solutions.

Ion Distributions, Exclusion Coefficients, and Separation Factors of Electrolytes in a Charged Cylindrical Nanopore: A Partially Perturbative Density Functional Theory Study