We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by the 3d periodic images of the walls is constant inside the simulation cell, with the field produced by the transverse images of the charged plates canceling out. The non-neutral confined electrolyte in an external potential can be simulated using 3d Ewald summation with a suitable renormalization of the electrostatic energy, to remove a divergence, and a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting a particle-particle particle-mesh approach.

BrazilWe study, using Density Functional theory and Monte Carlo simulations, aqueous electrolyte solutions between charged infinite planar surfaces, in a contact with a bulk salt reservoir. In agreement with recent experimental observations [Z. Luo et al., Nat. Comm. 6, 6358 (2015)], we find that the confined electrolyte lacks local charge neutrality. We show that a Density Functional Theory (DFT) based on a bulk-HNC expansion properly accounts for strong electrostatic correlations and allows us to accurately calculate the ionic density profiles between the charged surfaces, even for electrolytes containing trivalent counterions. The DFT allows us to explore the degree of local charge neutrality violation, as a function of plate separation and bulk electrolyte concentration, and to accurately calculate the interaction force between the charged surfaces.

We explore, using recently developed efficient Monte Carlo simulation method, interaction of anionic polyelectrolyte solution with a like-charged membrane. In addition to polyions, solution also contains salt with either monovalent, divalent, or trivalent counterions. In agreement with recent experimental observations, we find that multivalent counterions can lead to strong adsorption of polyions to a like charged surface.On the other hand, addition of 1:1 electrolyte diminishes the adsorption induced by the multivalent counterions. Dielectric discontinuity at the interface is found to play only a marginal role for the polyion adsorption.

We present an efficient method for simulating Coulomb systems confined by metal electrodes. The approach relies on Green function techniques to obtain the electrostatic potential for an infinite periodically replicated system. This avoids the use of image charges or an explicit calculation of the induced surface charge, both of which dramatically slows down the simulations. To demonstrate the utility of the new method, we use it to obtain the ionic density profiles and the differential capacitances, which are of great practical and theoretical interest, for a lattice model of an ionic liquid.

We study, using Monte Carlo simulations, the interaction between charged colloidal particles confined to the air-water interface. The dependence of force on ionic strength and counterion valence is explored. For 1:1 electrolyte, we find that the electrostatic interaction at the interface is very close to the one observed in the bulk. On the other hand, for salts with multivalent counterions, an interface produces an enhanced attraction between like charged colloids. Finally, we explore the effect of induced surface charge at the air-water interface on the interaction between colloidal particles.

We present an efficient approach for simulating Coulomb systems confined by planar polarizable surfaces. The method is based on the solution of the Poisson equation using periodic Green functions. It is shown that the electrostatic energy arising from the surface polarization can be decoupled from the energy due to the direct Coulomb interaction between the ions. This allows us to combine an efficient Ewald summation method, or any other fast method for summing over the replicas, with the polarization contribution calculated using Green function techniques. We apply the method to calculate density profiles of ions confined between the charged dielectric and metal surfaces.

Focus Article: Oscillatory and long-range monotonic exponential decays of electrostatic interactions in ionic liquids and other electrolytes: The significance of dielectric permittivity and renormalized charges

We study ionic liquids interacting with electrified interfaces. The ionic fluid is modeled as a Coulomb lattice gas. We compare the ionic density profiles calculated using a popular modified Poisson-Boltzmann equation with the explicit Monte Carlo simulations. The modified Poisson-Boltzmann theory fails to capture the structural features of the double layer and is also unable to correctly predict the ionic density at the electrified interface. The lattice Monte Carlo simulations qualitatively capture the coarse grained structure of the double layer in the continuum. We propose a convolution relation that semiquantitatively relates the ionic density profiles of a continuum ionic liquid and its lattice counterpart near an electrified interface.2