Andreas Härtel

The electric double layer (EDL) has a pivotal role in screening charges on surfaces as in supercapacitor electrodes or colloidal and polymer solutions. Its structure is determined by correlations between the finite-sized ionic charge carriers of the underlying electrolyte, and, this way, these correlations affect the properties of the EDL and of applications utilizing EDLs. We study the structure of EDLs within classical density functional theory (DFT) in order to uncover whether a structural transition in the first layer of the EDL that is driven by changes in the surface potential depends on specific particle interactions or has a general footing. This transition has been found in full-atom simulations. Thus far, investigating the in-plane structure of the EDL for the primitive model (PM) using DFT has proved a challenge. We show here that the use of an appropriate functional predicts the in-plane structure of EDLs in excellent agreement with molecular dynamics simulations. This provides the playground to investigate how the structure factor within a layer parallel to a charged surface changes as a function of both the applied surface potential and its separation from the surface. We discuss pitfalls in properly defining an in-plane structure factor and fully map out the structure of the EDL within the PM for a wide range of electrostatic electrode potentials. However, we do not find any signature of a structural crossover and conclude that the previously reported effect is not fundamental but rather occurs due to the specific force field of ions used in the simulations.

The electric double layer (EDL) has a pivotal role in screening charges on surfaces as in supercapacitor electrodes or colloidal and polymer solutions. Its structure is determined by correlations between the finite-sized ionic charge carriers of the underlying electrolyte and, this way, these correlations affect the properties of the EDL and of applications utilizing EDLs. We study the structure of EDLs within classical density functional theory (DFT) in order to uncover whether a structural transition in the first layer of the EDL that is driven by changes in the surface potential depends on specific particle interactions or has a general footing. This transition has been found in full-atom simulations. Thus far, investigating the in-plane structure of the EDL for the primitive model (PM) using DFT proved a challenge. We show here that the use of an appropriate functional predicts the in-plane structure of EDLs in excellent agreement with molecular dynamics (MD) simulations. This provides the playground to investigate how the structure factor within a layer parallel to a charged surface changes as function of both the applied surface potential and its separation from the surface. We discuss pitfalls in properly defining an in-plane structure factor and fully map out the structure of the EDL within the PM for a wide range of electrostatic electrode potentials. However, we do not find any signature of a structural crossover and conclude that the previously reported effect is not fundamental but rather occurs due to the specific force field of ions used in the simulations.

Underscreening is a collective term for charge correlations in electrolytes decaying slower than the Debye length. Anomalous underscreening refers to phenomenology that cannot be attributed alone to steric interactions. Experiments with concentrated electrolytes and ionic fluids report anomalous underscreening, which so far has not been observed in simulation. We present Molecular Dynamics simulation results exhibiting anomalous underscreening that can be connected to cluster formation. A theory that accounts for ion pairing confirms the trend. Our results challenge the classic understanding of dense electrolytes impacting the design of technologies for energy storage and conversion.

The primitive model describes ions by point charges with an additional hard-core interaction. In classical density-functional theory (DFT) the mean-field electrostatic contribution can be obtained from the first order of a functional perturbation of the pair potential for an uncharged reference system of hard spheres. This mean-field electrostatic term particularly contributes at particle separations that are forbidden due to hard-core overlap. In this work we modify the mean-field contribution such that the pair potential is constant for distances smaller than the contact distance of the ions. We motivate our modification by the underlying splitting of the potential, which is similar to the splitting of the Weeks–Chandler–Andersen potential and leads to higher-order terms in the respective expansion of the functional around the reference system. The resulting formalism involves weighted densities similar to the ones found in fundamental measure theory. To test our modifications, we analyze and compare density profiles, direct and total correlation functions, and the thermodynamic consistency of the functional via a widely established sum rule and the virial pressure formula for our modified functional, for established functionals, and for data from computer simulations. We found that our modifications clearly show improvements compared to the standard mean-field functional, especially when predicting layering effects and direct correlation functions in high concentration scenarios; for the latter we also find improved consistency when calculated via different thermodynamic routes. In conclusion, we demonstrate how modifications toward higher order corrections beyond mean-field functionals can be made and how they perform, by this providing a basis for systematic future improvements in classical DFT for the description of electrostatic interactions.

Recently, the reversible heat production during the electric double layer (EDL) buildup in a sodium chloride solution was measured experimentally [Janssen et al., Phys. Rev. Lett. 119, 166002 (2017)] and matched with theoretical predictions from density functional theory and molecular dynamics simulations [Glatzel et al., J. Chem. Phys. 154, 064901 (2021)]. In the latter, it was found that steric interactions of ions with the electrode’s walls, which result in the so-called Stern layer, are sufficient to explain the experimental results. As only symmetric ion sizes in a restricted primitive model were examined, it is instructive to investigate systems of unequal ion sizes that lead to modified Stern layers. In this work, we explore the impact of ion asymmetry on the reversible heat production for each electrode separately. In this context, we further study an extension of the primitive model where hydration shells of ions can evade in the vicinity of electrode’s walls. We find a strong dependence on system parameters such as particle sizes and the total volume taken by particles. Here, we even found situations where one electrode was heated and the other electrode was cooled at the same time during charging, while, in sum, both electrodes together behaved very similarly to the already mentioned experimental results. Thus, heat production should also be measured in experiments for each electrode separately. By this, the importance of certain ingredients that we proposed to model electrolytes could be confirmed or ruled out experimentally, finally leading to a deeper understanding of the physics of EDLs.

Electrolytes play an important role in a plethora of applications ranging from energy storage to biomaterials. Notwithstanding this, the structure of concentrated electrolytes remains enigmatic. Many theoretical approaches attempt to model the concentrated electrolyte by introducing the idea of ion pairs, with ions either being tightly “paired” with a counter-ion or “free” to screen charge. In this study, we reframe the problem into the language of computational statistics and test the null hypothesis that all ions share the same local environment. Applying the framework to molecular dynamics simulations, we find that this null hypothesis is not supported by data. Our statistical technique suggests the presence of two distinct local ionic environments at intermediate concentrations, whose differences surprisingly originate in like charge correlations rather than unlike charge attraction. Through considering the effect of these “aggregated” and “non-aggregated” states on bulk properties including effective ion concentration and dielectric constant, we identify a scaling relation between the effective screening length and theoretical Debye length, which applies across different dielectric constants and ion concentrations.

Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT) and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one is based on the simplest mean-field treatment of Coulomb interactions (MFC), while the other two employ the Mean Spherical Approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function, whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that are in excellent agreement with MD simulations in the near field, in the far field, we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT, we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.