Cavity correlations and the onset of charge ordering at charged interfaces: a modified Poisson-Fermi approach

TL;DR

This paper introduces a modified Poisson-Fermi approach incorporating charge cavities to accurately model charge layering and ordering at charged interfaces, especially in strong coupling regimes like RTILs.

Contribution

The modified PF approach extends traditional theory by including charge cavities, enabling the prediction of charge layering and ordering in strong coupling conditions.

Findings

Reproduces overscreening and charge ordering phenomena.

Identifies transition to undumped charge layering at high ionic strengths.

Provides a criterion for the stability against charge ordering.

Abstract

Charge layering in the close vicinity of charged interfaces is a well known effect, extensively reported in both experiments and simulations of Room Temperature Ionic Liquids (RTIL) and concentrated aqueous-like electrolytes. The traditional Poisson-Fermi (PF) theory is able to successfully describe overcrowding effects, but fails to reproduce charge ordering even in strong coupling regimes. Simple models, yet capable to investigate the interplay between these important interfacial phenomena, are still lacking. In order to bridge this gap, we herein present a modified PF (mPF) approach which is able to capture layering effects in strong coupling regimes typical of RTIL. The modification is based on the introduction of charge cavities around test-particles, which simply extend the exclusion volume effects to also incorporate the accompanying depletion of charges due to particle insertion. Addition of this simple ingredient is shown to reproduce overscreening and charge ordering, thereby extending the predictive power of the PF approach to strong coupling regimes. Using a linear response theory, we are able to study the emergence of charge ordering based on two characteristic lengths: a wavelength responsible by charge layering, along with a dumping length which screens charge oscillations. At large ionic strengths and strong couplings, the system undergoes a transition to undumped charge layering. The transition takes place when the poles of the Fourier components of the linear potential become real-valued. This criteria allows one to identify the transition line across the parameter space, thus delimiting the region of stability against unscreneed charge ordering.