Communication: Kinetic and pairing contributions in the dielectric spectra of electrolyte solutions

A screening of results on the decay length in concentrated electrolytes

Screening of electrostatic interactions in room-temperature ionic liquids and concentrated electrolytes has recently attracted much attention as surface force balance experiments have suggested the emergence of unanticipated anomalously large screening lengths at high ion concentrations. Termed underscreening, this effect was ascribed to the bulk properties of concentrated ionic systems. However, underscreening under experimentally relevant conditions is not predicted by classical theories and challenges our understanding of electrostatic correlations. Despite the enormous effort in performing large-scale simulations and new theoretical investigations, the origin of the anomalously long-range screening length remains elusive. This contribution briefly summarises the experimental, analytical and simulation results on ionic screening and the scaling behaviour of screening lengths. We then present an atomistic simulation approach that accounts for the solvent and ion exchange with a reservoir. We find that classical density functional theory (DFT) for concentrated electrolytes under confinement reproduces ion adsorption at charged interfaces surprisingly well. With DFT, we study confined electrolytes using implicit and explicit solvent models and the dependence on the solvent's dielectric properties. Our results demonstrate how the absence vs. presence of solvent particles and their discrete nature affect the short and long-range screening in concentrated ionic systems.

Solvent Effects on Structure and Screening in Confined Electrolytes

Using classical density functional theory, we investigate the influence of solvent on the structure and ionic screening of electrolytes under slit confinement and in contact with a reservoir. We consider a symmetric electrolyte with implicit and explicit solvent models and find that spatially resolving solvent molecules is essential for the ion structure at confining walls, excess ion adsorption, and the pressure exerted on the walls. Despite this, we observe only moderate differences in the period of oscillations of the pressure with the slit width and virtually coinciding decay lengths as functions of the scaling variable σ_{ion}/λ_{D}, where σ_{ion} is the ion diameter and λ_{D} the Debye length. Moreover, in the electrostatic-dominated regime, this scaling behavior is practically independent of the relative permittivity and its dependence on the ion concentration. In contrast, the crossover to the hard-core-dominated regime depends sensitively on all three factors.

An Atomic Insight into the Chemical Origin and Variation of the Dielectric Constant in Liquid Electrolytes

The dielectric constant is a crucial physicochemical property of liquids in tuning solute-solvent interactions and solvation microstructures. Herein the dielectric constant variation of liquid electrolytes regarding to temperatures and electrolyte compositions is probed by molecular dynamics simulations. Dielectric constants of solvents reduce as temperatures increase due to accelerated mobility of molecules. For solvent mixtures with different mixing ratios, their dielectric constants either follow a linear superposition rule or satisfy a polynomial function, depending on weak or strong intermolecular interactions. Dielectric constants of electrolytes exhibit a volcano trend with increasing salt concentrations, which can be attributed to dielectric contributions from salts and formation of solvation structures. This work affords an atomic insight into the dielectric constant variation and its chemical origin, which can deepen the fundamental understanding of solution chemistry.

Dielectric constant and density of aqueous alkali halide solutions by molecular dynamics: A force field assessment

The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation.

Dielectric Decrement for Aqueous NaCl Solutions: Effect of Ionic Charge Scaling in Nonpolarizable Water Force Fields

We investigate the dielectric constant and the dielectric decrement of aqueous NaCl solutions by means of molecular dynamic simulations. We thereby compare the performance of four different force fields and focus on disentangling the origin of the dielectric decrement and the influence of scaled ionic charges, as often used in nonpolarizable force fields to account for the missing dynamic polarizability in the shielding of electrostatic ion interactions. Three of the force fields showed excessive contact ion pair formation, which correlates with a reduced dielectric decrement. In spite of the fact that the scaling of charges only weakly influenced the average polarization of water molecules around an ion, the rescaling of ionic charges did influence the dielectric decrement, and a close-to-linear relation of the slope of the dielectric constant as a function of concentration with the ionic charge was found.

Preferential solvation and ion association properties in aqueous dimethyl sulfoxide solutions

We study the solvation and the association properties of ion pairs in aqueous dimethyl sulfoxide (DMSO) solution by atomistic molecular dynamics (MD) simulations. The ion pair is composed of two lithium and a single sulfonated diphenyl sulfone ion whose properties are studied under the influence of different DMSO concentrations. For increasing mole fractions of DMSO, we observe a non-ideal behavior of the solution as indicated by the derivatives of the chemical activity. Our findings are complemented by dielectric spectra, which also verify a complex DMSO-water mixing behavior. In agreement with these results, further simulation outcomes reveal an aqueous homoselective solvation of the ion species which fosters the occurrence of pronounced ion association constants at higher DMSO mole fractions. The consequences of this finding are demonstrated by lower ionic conductivities for increasing concentrations of DMSO.

A coarse-grained polarizable force field for the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate

We present a coarse-grained polarizable molecular dynamics force field for the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIm][PF₆]). For the treatment of electronic polarizability, we employ the Drude model. Our results show that the new explicitly polarizable force field reproduces important static and dynamic properties such as mass density, enthalpy of vaporization, diffusion coefficients, or electrical conductivity in the relevant temperature range. In situations where an explicit treatment of electronic polarizability might be crucial, we expect the force field to be an improvement over non-polarizable models, while still profiting from the reduction of computational cost due to the coarse-grained representation.

Revealing the signature of dipolar interactions in dynamic spectra of polydisperse magnetic nanoparticles

We investigate, via a modified mean field approach, the dynamic magnetic response of a polydisperse dipolar suspension to a weak, linearly polarised, AC field. We introduce an additional term into the Fokker-Planck equation, which takes into account dipole-dipole interaction in the form of the first order perturbation, and allows for particle polydispersity. The analytical expressions, obtained for the real and imaginary dynamic susceptibilities, predict three measurable effects: the increase of the real part low-frequency plateaux; the enhanced growth of the imaginary part in the low-frequency range; and the shift of the imaginary part maximum. Our theoretical predictions find an experimental confirmation and explain the changes in the spectrum.

Temperature-dependent dynamic correlations in suspensions of magnetic nanoparticles in a broad range of concentrations: a combined experimental and theoretical study

We study the effects of temperature and concentration on the dynamic spectra of polydisperse magnetic nanoparticle suspensions.

Dielectric response of the water hydration layer around spherical solutes

We calculate the local dielectric function ɛ(r) inside the hydration layer around a spherical solute (i) from molecular dynamics simulations including explicit solutes and (ii) theoretically using the nonlocal dielectric function of bulk water which includes the radial electric field, but not the explicit solute. The observed agreement between the two concepts shows that while ɛ(r) is strongly different from bulk, this difference is not due to restructuring of the hydrogen bond network but is mostly a consequence of the field geometry. The dielectric response differs for anions and cations, yielding a natural explanation for the well-known charge asymmetry of ionic solvation in agreement with experimental data.

Kinetic dielectric decrement revisited: phenomenology of finite ion concentrations

With the help of a recently developed non-equilibrium approach, we investigate the ionic strength dependence of the Hubbard-Onsager dielectric decrement. We compute the depolarization of water molecules caused by the motion of ions in sodium chloride solutions from the dilute regime (0.035 M) up close to the saturation concentration (4.24 M), and find that the kinetic decrement displays a strong non-monotonic behavior, in contrast to the prediction of available models. We introduce a phenomenological modification of the Hubbard-Onsager continuum theory, which takes into account the screening due to the ionic cloud at the mean-field level and, which is able to describe the kinetic decrement at high concentrations including the presence of a pronounced minimum.

Dissecting ion-specific dielectric spectra of sodium-halide solutions into solvation water and ionic contributions

Using extensive equilibrium molecular dynamics simulations we determine the dielectric spectra of aqueous solutions of NaF, NaCl, NaBr, and NaI. The ion-specific and concentration-dependent shifts of the static dielectric constants and the dielectric relaxation times match experimental results very well, which serves as a validation of the classical and non-polarizable ionic force fields used. The purely ionic contribution to the dielectric response is negligible, but determines the conductivity of the salt solutions. The ion-water cross correlation contribution is negative and reduces the total dielectric response by about 5%-10% for 1 M solutions. The dominating water dielectric response is decomposed into different water solvation shells and ion-pair configurations, by this the spectral blue shift and the dielectric decrement of salt solutions with increasing salt concentration is demonstrated to be primarily caused by first-solvation shell water. With rising salt concentration the simulated spectra show more pronounced deviations from a single-Debye form and can be well described by a Cole-Cole fit, in quantitative agreement with experiments. Our spectral decomposition into ionic and different water solvation shell contributions does not render the individual contributions more Debye-like, this suggests the non-Debye-like character of the dielectric spectra of salt solutions not to be due to the superposition of different elementary relaxation processes with different relaxation times. Rather, the non-Debye-like character is likely to be an inherent spectral signature of solvation water around ions.

A new lattice Monte Carlo simulation for dielectric saturation in ion-containing liquids

We develop a new, rapid method for the lattice Monte Carlo simulation of ion-containing liquids that accounts for the effects of the reorganization of solvent dipoles under external electrostatic fields. Our results are in reasonable agreement with the analytical solutions to the dielectric continuum theory of Booth for single ions, ion pairs, and ionic cross-links. We also illustrate the substantial disparity between the dielectric functions for like and unlike charges on the nanometer scale. Our simulation rationalizes the experimental data for the dependence of the bulk dielectric value of water on ion concentrations in terms of saturated dipoles near ions.