The intimate relationship between the dielectric response and the decay of intermolecular correlations and surface forces in electrolytes

Structure and Screening in Confined Electrolytes: The Role of Ion Association

The effect of ionic association on the structure and property of confined electrolytes is investigated using the classical density functional theory. We find that ionic association strongly affects the ion distribution, surface force, and screening behavior of confined electrolytes. The decay length ξ, which can describe the screening effect of high-concentration electrolytes, satisfies a scaling relationship ξ/λD ∼ (σ/λD)n, with λD being the Debye length and σ representing the ion diameter. We find that n = 1.5 in the nonassociation model, which is contributed by the charge correlation, but n = 3 in the association model, which is contributed by the density correlation. The ion association changes the concentration-dependent characteristics of the screening length by promoting the shift of the decay behavior from the charge-dominated regime to the density-dominated regime. Our result reveals the importance of ion association for electrolyte structure and screening behaviors.

Nonlocal structural effects of water on DNA homology recognition

The mechanism behind mutual recognition of homologous DNA sequences prior to genetic recombination is one of the remaining puzzles in molecular biology. Leading models of homology recognition, based on classical electrostatics, neglect the short-range nonlocal screening effects arising from structured water around DNA, and hence may only provide insight for relatively large separations between interacting DNAs. We elucidate the role of the effects of the nonlocal dielectric response of water on DNA-DNA interaction and show that these can dramatically enhance the driving force for recognition.

Measuring the Capacitance of Carbon in Ionic Liquids: From Graphite to Graphene

The physical electrochemistry of the carbon/ionic liquids interface underpins the processes occurring in a vast range of applications spanning electrochemical energy storage, iontronic devices, and lubrication. Elucidating the charge storage mechanisms at the carbon/electrolyte interface will lead to a better understanding of the operational principles of such systems. Herein, we probe the charge stored at the electrochemical double layer formed between model carbon systems, ranging from single-layer graphene to graphite and the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI). The effect of the number of graphene layers on the overall capacitance of the interface is investigated. We demonstrate that in pure EMIM-TFSI and at moderate potential biases, the electronic properties of graphene and graphite govern the overall capacitance of the interface, while the electrolyte contribution to the latter is less significant. In mixtures of EMIM-TFSI with solvents of varying relative permittivity, the complex interplay between electrolyte ions and solvent molecules is shown to influence the charge stored at the interface, which under certain conditions overcomes the effects of relative permittivity. This work provides additional experimental insights into the continuously advancing topic of electrochemical double-layer structure at the interface between room temperature ionic liquids and carbon materials.

Exploring how cation entropy influences electric double layer formation and electrochemical reactivity

Electric double layers are crucial to energy storage and electrocatalytic device performance. While double layer formation originates in electrostatic interactions, electric double layer properties are governed by a balance of both electrostatic and entropic driving forces. Favorable ion-surface electrostatic interactions attract counterions to charged surfaces to compensate, or "screen," potentials, but the confinement of these same ions from a bulk reservoir to the interface incurs an entropic penalty. Here, we use a dicationic imidazolium ionic liquid and its monovalent analogue to explore how cation valence and entropy influence double layer formation and electrochemical reactivity using CO2 electroreduction as a model reaction. We find that divalent and monovalent cations display similar CO2 reduction kinetics but differ vastly in steady-state reactivity due to rapid electrochemically induced precipitation of insulating dicationic (bi)carbonate films. Using in situ surface-enhanced Raman scattering spectroscopy, we find that potential-dependent cation reorientation occurs at similar potentials between the two ionic liquids, but the introduction of a covalent link in the divalent cation imparts a more ordered double layer structure that favors (bi)carbonate precipitation. In mixed monovalent-divalent electrolytes, we find that the divalent cations dominate interfacial properties by preferentially accumulating at surfaces even at very low relative concentrations. Our findings confirm that ion entropy plays a key role in modulating local electrochemical environments. Furthermore, we highlight how double layer properties are sensitive to the properties of counterions that pay the lowest entropic penalty to accumulate at interfaces. Overall, we illustrate that ion entropy provides a new knob to tune reaction microenvironments and unveil how entropy plays a major role in modulating electrochemical reactivity in mixed ion electrolytes.

Colloidal Systems in Concentrated Electrolyte Solutions Exhibit Re-entrant Long-Range Electrostatic Interactions due to Underscreening

Surface force measurements have revealed that at very high electrolyte concentrations as well as in neat and diluted ionic liquids and deep eutectic solvents, the range of electrostatic interactions is far greater than the Debye length. Here, we explore the consequences of this underscreening for soft-matter and colloidal systems by investigating the stability of nanoparticle dispersions, the self-assembly of ionic surfactants, and the thickness of soap films. In each case, we find clear evidence of re-entrant properties due to underscreening at high salt concentrations. Our results show that underscreening in concentrated electrolytes is a general phenomenon and is not dependent on confinement by macroscopic surfaces. The stability of systems at very high salinity due to underscreening may be beneficially applied to processes that currently use low-salinity water.

Mesoscopic Inhomogeneities in Concentrated Electrolytes

A mesoscopic theory for water-in-salt electrolytes combining density functional and field-theoretic methods is developed in order to explain the unexpectedly large period of the oscillatory decay of the disjoining pressure observed in recent experiments for the lithium bis(trifluoromethylsulfonyl)-imide (LiTFSI) salt [T. S. Groves et al., J. Phys. Chem. Lett. 2021, 12, 1702]. We assumed spherical ions with different diameters and implicit solvent, inducing strong, short-range attraction between ions of the same sign. For this highly simplified model, we calculated correlation functions. Our results indicate that mesoscopic inhomogeneities can occur when the sum of the Coulomb and the water-mediated interactions between like ions is attractive at short and repulsive at large distances. We adjusted the attractive part of the potential to the water-in-LiTFSI electrolyte and obtained both the period and the decay rate of the correlations, in semiquantitative agreement with the experiment. In particular, the decay length of the correlations increases nearly linearly with the volume fraction of ions.

Ionic screening in bulk and under confinement

Recent experiments have shown that the repulsive force between atomically flat, like-charged surfaces confining room-temperature ionic liquids or concentrated electrolytes exhibits an anomalously large decay length. In our previous publication [J. Zeman, S. Kondrat, and C. Holm, Chem. Commun. 56, 15635 (2020)], we showed by means of extremely large-scale molecular dynamics simulations that this so-called underscreening effect might not be a feature of bulk electrolytes. Herein, we corroborate these findings by providing additional results with more detailed analyses and expand our investigations to ionic liquids under confinement. Unlike in bulk systems, where screening lengths are computed from the decay of interionic potentials of mean force, we extract such data in confined systems from cumulative charge distributions. At high concentrations, our simulations show increasing screening lengths with increasing electrolyte concentration, consistent with classical liquid state theories. However, our analyses demonstrate that-also for confined systems-there is no anomalously large screening length. As expected, the screening lengths determined for ionic liquids under confinement are in good quantitative agreement with the screening lengths of the same ionic systems in bulk. In addition, we show that some theoretical models used in the literature to relate the measured screening lengths to other observables are inapplicable to highly concentrated electrolytes.

Overscreening and Underscreening in Solid-Electrolyte Grain Boundary Space-Charge Layers

Polycrystalline solids can exhibit material properties that differ significantly from those of equivalent single-crystal samples, in part, because of a spontaneous redistribution of mobile point defects into so-called space-charge regions adjacent to grain boundaries. The general analytical form of these space-charge regions is known only in the dilute limit, where defect-defect correlations can be neglected. Using kinetic Monte Carlo simulations of a three-dimensional Coulomb lattice gas, we show that grain boundary space-charge regions in nondilute solid electrolytes exhibit overscreening-damped oscillatory space-charge profiles-and underscreening-decay lengths that are longer than the corresponding Debye length and that increase with increasing defect-defect interaction strength. Overscreening and underscreening are known phenomena in concentrated liquid electrolytes, and the observation of functionally analogous behavior in solid electrolyte space-charge regions suggests that the same underlying physics drives behavior in both classes of systems. We therefore expect theoretical approaches developed to study nondilute liquid electrolytes to be equally applicable to future studies of solid electrolytes.

Primitive model electrolytes in the near and far field: Decay lengths from DFT and simulations

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.

Correlation Length in Concentrated Electrolytes: Insights from All-Atom Molecular Dynamics Simulations

We study the correlation length of the charge-charge pair correlations in concentrated electrolyte solutions by means of all-atom, explicit-solvent molecular dynamics simulations. We investigate LiCl and NaI in water, which constitute highly soluble, prototypical salts for experiments, as well as two more complex, molecular electrolyte systems of lithium bis(trifluoromethane)sulfonimide (LiTFSI), a salt commonly employed in electrochemical storage systems, in water, and in an organic solvent mixture of dimethoxyethane and dioxolane. Our simulations support the recent experimental observations as well as theoretical predictions of a nonmonotonic behavior of the correlation length with increasing salt concentration. We observe a Debye-Hückel like regime at low concentration, followed by a minimum reached when d/λ_D ≃ 1, where λ_D is the Debye correlation length and d is the effective ionic diameter, and an increasing correlation length with salt concentration in very concentrated electrolytes. As in the experiments, we find that the screening length in the concentrated regime follows a universal scaling law as a function d/λ_D for all studied salts. However, the scaling exponent is significantly lower than the experimentally measured one and lies in the range of the theoretical predictions based on much simpler electrolyte models.