Free energy of ion hydration: Interface susceptibility and scaling with the ion size

Redox potentials in ionic liquids: Anomalous behavior?

Redox potentials depend on the nature of the solvent/electrolyte through the solvation energies of the ionic solute species. For concentrated electrolytes, ion solvation may deviate significantly from the Born model predictions due to ion pairing and correlation effects. Recently, Ghorai and Matyushov [J. Phys. Chem. B 124, 3754-3769 (2020)] predicted, on the basis of linear response theory, an anomalous trend in the solvation energies of room temperature ionic liquids, with deviations of hundreds of kJ/mol from the Born model for certain size solutes/ions. In this work, we computationally evaluate ionic solvation energies in the prototypical ionic liquid, 1-butyl-3-methylimidazolium tetrafluoroborate (BMIM/BF4), to further explore this behavior and benchmark several of the approximations utilized in the solvation energy predictions. For comparison, we additionally compute solvation energies within acetonitrile and molten NaCl salt to illustrate the limiting behavior of purely dipolar and ionic solvents. We find that the overscreening effect, which results from the inherent charge oscillations of the ionic liquid, is substantially reduced in magnitude due to screening from the dipoles of the molecular ions. Therefore, for the molten NaCl salt, for which the ions do not have permanent dipoles, modulation of ionic solvation energies from the overscreening effect is most significant. The conclusion is that ionic liquids do indeed exhibit unique solvation behavior due to peak(s) in the electrical susceptibility caused by the ion shell structure; redox potential shifts for BMIM/BF4 are of more modest order ∼0.1 V, but may be larger for other ionic liquids that approach molten salt behavior.

Generalized Helmholtz model describes capacitance profiles of ionic liquids and concentrated aqueous electrolytes

In this work, we propose and validate a generalization of the Helmholtz model that can account for both "bell-shaped" and "camel-shaped" differential capacitance profiles of concentrated electrolytes, the latter being characteristic of ionic liquids. The generalization is based on introducing voltage dependence of both the dielectric constant "ϵr(V)" and thickness "L(V)" of the inner Helmholtz layer, as validated by molecular dynamics (MD) simulations. We utilize MD simulations to study the capacitance profiles of three different electrochemical interfaces: (1) graphite/[BMIm+][BF4-] ionic liquid interface; (2) Au(100)/[BMIm+][BF4-] ionic liquid interface; (3) Au(100)/1M [Na+][Cl-] aqueous interface. We compute the voltage dependence of ϵr(V) and L(V) and demonstrate that the generalized Helmholtz model qualitatively describes both camel-shaped and bell-shaped differential capacitance profiles of ionic liquids and concentrated aqueous electrolytes (in lieu of specific ion adsorption). In particular, the camel-shaped capacitance profile that is characteristic of ionic liquid electrolytes arises simply from combination of the voltage-dependent trends of ϵr(V) and L(V). Furthermore, explicit analysis of the inner layer charge density for both concentrated aqueous and ionic liquid double layers reveal similarities, with these charge distributions typically exhibiting a dipolar region closest to the electrode followed by a monopolar peak at larger distances. It is appealing that a generalized Helmholtz model can provide a unified description of the inner layer structure and capacitance profile for seemingly disparate aqueous and ionic liquid electrolytes.

Solvation of Nanoions in Aqueous Solutions

In recent years it has been increasingly recognized that different classes of large ions with multiple valency have effects conceptually similar to weakly solvated ions in the Hofmeister series, also labeled by the term chaotropic. The term "superchaotropic effect" has been coined because these effects are much more strongly pronounced for nanometer-sized ions, whose adsorption properties often resemble typical surfactants. Despite this growing interest in these nanometer-sized ions, a simple conceptual extension of the Hofmeister series toward nanoions has not been achieved because an extrapolation of the one-dimensional surface charge density scale does not lead to the superchaotropic regime. In this work, we discuss a generic model that is broadly applicable to ions of nearly spherical shape and thus includes polyoxometalates and boron clusters. We present a qualitative classification scheme in which the ion size appears as a second dimension. Ions of different sizes but the same charge density differ in their bulk solvation free energy. As the ions grow bigger at constant surface charge density, they become more stable in solution, but the adsorption behavior is still governed by the surface charge density. A detailed molecular dynamics simulation study of large ions that is based on a shifted Lennard-Jones potential is presented that supports the presented classification scheme.

Electrical noise in electrolytes: a theoretical perspective

Seemingly unrelated experiments such as electrolyte transport through nanotubes, nano-scale electrochemistry, NMR relaxometry and surface force balance measurements, all probe electrical fluctuations: of the electric current, the charge and polarization, the field gradient (for quadrupolar nuclei) and the coupled mass/charge densities. The fluctuations of such various observables arise from the same underlying microscopic dynamics of the ions and solvent molecules. In principle, the relevant length and time scales of these dynamics are encoded in the dynamic structure factors. However, modelling the latter for frequencies and wavevectors spanning many orders of magnitude remains a great challenge to interpret the experiments in terms of physical processes such as solvation dynamics, diffusion, electrostatic and hydrodynamic interactions between ions, interactions with solid surfaces, etc. Here, we highlight the central role of the charge-charge dynamic structure factor in the fluctuations of electrical observables in electrolytes and offer a unifying perspective over a variety of complementary experiments. We further analyze this quantity in the special case of an aqueous NaCl electrolyte, using simulations with explicit ions and an explicit or implicit solvent. We discuss the ability of the standard Poisson-Nernst-Planck theory to capture the simulation results, and how the predictions can be improved. We finally discuss the contributions of ions and water to the total charge fluctuations. This work illustrates an ongoing effort towards a comprehensive understanding of electrical fluctuations in bulk and confined electrolytes, in order to enable experimentalists to decipher the microscopic properties encoded in the measured electrical noise.

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. We propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

Dielectric Behavior near a Spherical Ion

The dielectric response of a polar solvent to an ion is analyzed in terms of the bound charge, the net charge that accumulates near the ion as a consequence of the inhomogeneous polarization of the surrounding solvent. We demonstrate that the total bound charge arising in a full molecular treatment is identical to the total bound charge from standard continuum theory. In continuum theory, the bound charge resides in an infinitely thin layer, while in a molecular description the bound charge is spread over a region of finite width. Near simple atomic ions, the width of the bound charge distribution is roughly 1.3 nm. By simulating a sequence of ion charges from 0.1 to 2 e, where e is the magnitude of the electron charge, we analyze the applicability of linear response theory, which has been used by several authors. With increasing charge, the nonlinear response extends to an increasing distance from the ion. However, outside the region containing bound charge, the response is linear and in accord with continuum theory. Previous attempts to assign a dielectric constant to a solvent in the interfacial region are analyzed.

Electrostatic solvation and mobility in uniform and non-uniform electric fields: From simple ions to proteins

A number of observations related to interfacial electrostatics of polar liquids question the traditional assumption of dielectric theories that bulk dielectric properties can be continuously extended to the dividing surface separating the solute from the solvent. The deficiency of this approximation can be remedied by introducing local interface susceptibilities and the interface dielectric constant. Asymmetries of ionic hydration thermodynamics and of the mobility between cations and anions can be related to different propensities of the water molecules to orient their dipole toward and outward from solutes of opposite charges. This electrostatic asymmetry is reflected in different interface dielectric constants for cations and anions. The interface of water with neutral solutes is spontaneously polarized due to preferential water orientations in the interface. This phenomenon is responsible for a nonzero cavity potential directly related to a nonzero surface charge. This connection predicts that particles allowing a nonzero cavity potential must show mobility in an external electric field even if the net charge of the particle is zero. The theory predicts that a positive cavity potential and a positive surface charge translate to an effectively negative solute charge reported by mobility measurements. Passing of the cavity potential through a minimum found in simulations might be the origin of the maximum of mobility vs the ionic size observed experimentally. Finally, mobility of proteins in the field gradient (dielectrophoresis) is many orders of magnitude greater than predicted by the traditionally used Clausius-Mossotti equation. Two reasons contribute to this disagreement: (i) a failure of Maxwell's electrostatics to describe the cavity-field susceptibility and (ii) the neglect of the protein permanent dipole by the Clausius-Mossotti equation. An analytical relation between the dielectrophoretic susceptibility and dielectric spectroscopy of solutions provides direct access to this parameter, confirming the failure of the Clausius-Mossotti equation in application to protein dielectrophresis.

Interfacial structural crossover and hydration thermodynamics of charged C60 in water

Classical molecular dynamics simulations of the hydration thermodynamics, structure, and dynamics of water in hydration shells of charged buckminsterfullerenes are presented in this study. Charging of fullerenes leads to a structural transition in the hydration shell, accompanied by creation of a significant population of dangling O-H bonds pointing toward the solute. In contrast to the well accepted structure-function paradigm, this interfacial structural transition causes nearly no effect on either the dynamics of hydration water or on the solvation thermodynamics. Linear response to the solute charge is maintained despite significant structural changes in the hydration shell, and solvation thermodynamic potentials are nearly insensitive to the altering structure. Only solvation heat capacities, which are higher thermodynamic derivatives of the solvation free energy, indicate some sensitivity to the local hydration structure. We have separated the solvation thermodynamic potentials into direct solute-solvent interactions and restructuring of the hydration shell and analyzed the relative contributions of electrostatic and nonpolar interactions to the solvation thermodynamics.

Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited

The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between -0.9 and -0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.