Transferable Ion Force Fields in Water from a Simultaneous Optimization of Ion Solvation and Ion-Ion Interaction

A Coarse-Grained SPICA Makeover for Solvated and Bare Sodium and Chloride Ions

Aqueous ionic solutions are pivotal in various scientific domains due to their natural prevalence and vital roles in biological and chemical processes. Molecular dynamics has emerged as an effective methodology for studying the dynamic behavior of these systems. While all-atomistic models have made significant strides in accurately representing and simulating these ions, the challenge persists in achieving precise models for coarse-grained (CG) simulations. Our study introduces two optimized models for sodium and chloride ions within the nonpolarizable surface property fitting coarse-grained force field (SPICA-FF) framework. The two models represent solvated ions, such as the original FF model, and unsolvated or bare ions. The nonbonded Lennard-Jones interactions were reparameterized to faithfully reproduce bulk properties, including density and surface tension, in sodium chloride solutions at varying concentrations. Notably, these optimized models replicate experimental surface tensions at high ionic strengths, a property not well-captured by the ions of the original model in the SPICA-FF. The optimized unsolvated model also proved successful in reproducing experimental osmotic pressure. Additionally, the newly reparameterized ion models capture hydrophobic interactions within sodium chloride solutions and show qualitative agreement when modeling structural changes in phospholipid bilayers, aligning with experimental observations. For aqueous solutions, these optimized models promise a more precise representation of the ion behavior.

Force field for halide and alkali ions in water based on single-ion and ion-pair thermodynamic properties for a wide range of concentrations

A classical non-polarizable force field for the common halide (F-, Cl-, Br-, and I-) and alkali (Li+, Na+, K+, and Cs+) ions in SPC/E water is presented. This is an extension of the force field developed by Loche et al. for Na+, K+, Cl-, and Br- (JPCB 125, 8581-8587, 2021): in the present work, we additionally optimize Lennard-Jones parameters for Li+, I-, Cs+, and F- ions. Li+ and F- are particularly challenging ions to model due to their small size. The force field is optimized with respect to experimental solvation free energies and activity coefficients, which are the necessary and sufficient quantities to accurately reproduce the electrolyte thermodynamics. Good agreement with experimental reference data is achieved for a wide range of concentrations (up to 4 mol/l). We find that standard Lorentz-Berthelot combination rules are sufficient for all ions except F-, for which modified combination rules are necessary. With the optimized parameters, we show that, although the force field is only optimized based on thermodynamic properties, structural properties are reproduced quantitatively, while ion diffusion coefficients are in qualitative agreement with experimental values.

Sodium Triflate Water-in-Salt Electrolytes in Advanced Battery Applications: A First-Principles-Based Molecular Dynamics Study

Offering a compelling combination of safety and cost-effectiveness, water-in-salt (WiS) electrolytes have emerged as promising frontiers in energy storage technology. Still, there is a strong demand for research and development efforts to make these electrolytes ripe for commercialization. Here, we present a first-principles-based molecular dynamics (MD) study addressing in detail the properties of a sodium triflate WiS electrolyte for Na-ion batteries. We have developed a workflow based on a machine learning (ML) potential derived from ab initio MD simulations. As ML potentials are typically restricted to the interpolation of the data points of the training set and have hardly any predictive properties, we subsequently optimize a classical force field based on physics principles to ensure broad applicability and high performance. Performing and analyzing detailed MD simulations, we identify several very promising properties of the sodium triflate as a WiS electrolyte but also indicate some potential stability challenges associated with its use as a battery electrolyte.

Force Fields for High Concentration Aqueous KOH Solutions and Zincate Ions

Motivated by increasing interest in electrochemical devices that include highly alkaline electrolytes, we investigated two force fields for potassium hydroxide (KOH) at high concentrations in water. The "FNB" model uses the SPC/E water model, while the "FHM" model uses the TIP4P/2005 water model. We also developed parameters to describe zincate ions in these solutions. The density and viscosity of KOH using the FHM model are in better agreement with experiment than the values from the FNB model. Comparing the properties of the zincate solutions to the available experimental data, we find that both force fields agree reasonably well, although the FHM parameters give a better prediction of the viscosity. The developed force field parameters can be used in future simulations of zincate/KOH solutions in combination with other species of interest.

Multiscale Modeling of Aqueous Electric Double Layers

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Anomalous Concentration Dependence of Viscosity: Hidden Role of Cross-Correlations in Aqueous Electrolyte Solutions

The viscosity of aqueous electrolyte solutions exhibits well-known composition-dependent anomalies that show certain definitive trends and universal features. The viscosity of LiCl and NaCl solutions increases with concentration in a monotonic fashion, while solutions of KCl, RbCl, and CsCl exhibit a more complex behavior. Here, the viscosity first decreases and then increases with increasing concentration, with a rather broad minimum at intermediate concentrations (ca. 1-3 m). To unearth the origin of such puzzling behavior, we carried out detailed molecular-level analyses by interrogating the exact Green-Kubo expression of viscosity in terms of the stress-stress time correlation function (SS-TCF). The total SS-TCF can be decomposed into a collection of three self- and three cross-SS-TCFs arising from the three constituent components (water, cations, and anions). Mode coupling theory (MCT) analysis for the friction on ions and the viscosity of the solution suggests the possible importance of two-particle static and time-dependent cross-correlations between water and the ions. We calculate the viscosity and other dynamical properties for all five electrolyte (LiCl, NaCl, KCl, RbCl, and CsCl) solutions over a range of concentrations, using two models of water (SPC/E and TIP4P/2005). The total viscosity derives non-negligible contributions from all of the terms. The cross-correlations are found to be surprisingly large and seen to play a hidden role in the concentration dependence. However, the importance of cross-correlations is often not discussed. Our study leads to a theoretical understanding of the microscopic origin of the observed anomalies in the composition dependence of viscosity across all five electrolytes.

Molecular Simulation of the Impact of Defects on Electrolyte Intrusion in Zeolites

We have investigated through molecular simulation the intrusion of electrolytes in two representative pure-silica zeolites, silicalite-1 and chabazite, in which point defects were introduced in varying amounts. We distinguish between two types of defects, considering either "weak" or "strong" silanol nest defects, resulting in different hydration behaviors. In the presence of weak defects, the hydration process occurs through a homogeneous nucleation process, while with strong defects, we observe an initial adsorption followed by a filling of the nanoporous volume at a higher pressure. However, we show that electrolytes do not penetrate the zeolites, and these defects appear to have only marginal influence on the thermodynamics of electrolyte intrusion. While replacing pure water by the electrolyte solution shifts the intrusion pressure toward higher values because of the drop of water saturation vapor pressure, an increase in hydrophilicity of the framework due to point defects has the opposite effect, showing that controlling the amount of defects in zeolites is crucial for storage energy applications.

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.

Finding Infinities in Nanoconfined Geothermal Electrolyte Static Dielectric Properties and Implications on Ion Adsorption/Pairing

Infinities should naturally occur in the dielectric responses of ionic solutions relevant to many geochemical, energy storage, and electrochemical applications at a strictly zero frequency. Using molecular dynamics simulations cross-referenced with coarse-grained Monte Carlo models, using nanoslit pore models at hydrothermal conditions, and treating confined mobile charges as polarization, we demonstrate the far reaching consequences. The dielectric permittivity profile perpendicular to the slit (ϵ⊥(z)) increases, not decreases, with ionic concentration, unlike in the more widely studied megahertz-to-gigahertz frequency range. In confined electrolytes, the divergences in ϵ⊥(z) correctly describe crossovers between bulk- and surface-dominated dielectric behavior. Nanoconfinement at low ionic concentrations changes monovalent ion energetics by 1-2 kJ/mol, but no dielectric property studied so far is universally correlated to ion adsorption or ion-ion interactions. We caution that infinities signal violation of the "electrical insulator" dielectric assumption.

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

We develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

Interfacial Tension-Temperature-Pressure-Salinity Relationship for the Hydrogen-Brine System under Reservoir Conditions: Integration of Molecular Dynamics and Machine Learning

Hydrogen (H2) underground storage has attracted considerable attention as a potentially efficient strategy for the large-scale storage of H2. Nevertheless, successful execution and long-term storage and withdrawal of H2 necessitate a thorough understanding of the physical and chemical properties of H2 in contact with the resident fluids. As capillary forces control H2 migration and trapping in a subsurface environment, quantifying the interfacial tension (IFT) between H2 and the resident fluids in the subsurface is important. In this study, molecular dynamics (MD) simulation was employed to develop a data set for the IFT of H2-brine systems under a wide range of thermodynamic conditions (298-373 K temperatures and 1-30 MPa pressures) and NaCl salinities (0-5.02 mol·kg-1). For the first time to our knowledge, a comprehensive assessment was carried out to introduce the most accurate force field combination for H2-brine systems in predicting interfacial properties with an absolute relative deviation (ARD) of less than 3% compared with the experimental data. In addition, the effect of the cation type was investigated for brines containing NaCl, KCl, CaCl2, and MgCl2. Our results show that H2-brine IFT decreases with increasing temperature under any pressure condition, while higher NaCl salinity increases the IFT. A slight decrease in IFT occurs when the pressure increases. Under the impact of cation type, Ca2+ can increase IFT values more than others, i.e., up to 12% with respect to KCl. In the last step, the predicted IFT data set was used to provide a reliable correlation using machine learning (ML). Three white-box ML approaches of the group method of data handling (GMDH), gene expression programming (GEP), and genetic programming (GP) were applied. GP demonstrates the most accurate correlation with a coefficient of determination (R2) and absolute average relative deviation (AARD) of 0.9783 and 0.9767%, respectively.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Competition between β-Sheet and Coacervate Domains Yields Diverse Morphologies in Mixtures of Oppositely Charged Homochiral Polypeptides

Biomolecular assembly processes involving competition between specific intermolecular interactions and thermodynamic phase instability have been implicated in a number of pathological states and technological applications of biomaterials. As a model for such processes, aqueous mixtures of oppositely charged homochiral polypeptides such as poly-l-lysine and poly-l-glutamic acid have been reported to form either β-sheet-rich solid-like precipitates or liquid-like coacervate droplets depending on competing hydrogen bonding interactions. Herein, we report studies of polypeptide mixtures that reveal unexpectedly diverse morphologies ranging from partially coalescing and aggregated droplets to bulk precipitates, as well as a previously unreported re-entrant liquid-liquid phase separation at high polypeptide concentration and ionic strength. Combining our experimental results with all-atom molecular dynamics simulations of folded polypeptide complexes reveals a concentration dependence of β-sheet-rich secondary structure, whose relative composition correlates with the observed macroscale morphologies of the mixtures. These results elucidate a crucial balance of interactions that are important for controlling morphology during coacervation in these and potentially similar biologically relevant systems.

NMR Investigation of Water in Salt Crusts: Insights from Experiments and Molecular Simulations

The evaporation of water from bare soil is often accompanied by the formation of a layer of crystallized salt, a process that must be understood in order to address the issue of soil salinization. Here, we use nuclear magnetic relaxation dispersion measurements to better understand the dynamic properties of water within two types of salt crusts: sodium chloride (NaCl) and sodium sulfate (Na₂SO₄). Our experimental results display a stronger dispersion of the relaxation time T₁ with frequency for the case of sodium sulfate as compared to sodium chloride salt crusts. To gain insight into these results, we perform molecular dynamics simulations of salt solutions confined within slit nanopores made of either NaCl or Na₂SO₄. We find a strong dependence of the value of the relaxation time T₁ on pore size and salt concentration. Our simulations reveal the complex interplay between the adsorption of ions at the solid surface, the structure of water near the interface, and the dispersion of T₁ at low frequency, which we attribute to adsorption-desorption events.

Pair-Reaction Dynamics in Water: Competition of Memory, Potential Shape, and Inertial Effects

When described by a one-dimensional reaction coordinate, pair-reaction rates in a solvent depend, in addition to the potential barrier height and the friction coefficient, on the potential shape, the effective mass, and the friction relaxation spectrum, but a rate theory that accurately accounts for all of these effects does not exist. After a review of classical reaction-rate theories, we show how to extract all parameters of the generalized Langevin equation (GLE) and, in particular, the friction memory function from molecular dynamics (MD) simulations of two prototypical pair reactions in water, the dissociation of NaCl and of two methane molecules. The memory exhibits multiple time scales and, for NaCl, pronounced oscillatory components. Simulations of the GLE by Markovian embedding techniques accurately reproduce the pair-reaction kinetics from MD simulations without any fitting parameters, which confirms the accuracy of the approximative form of the GLE and of the parameter extraction techniques. By modification of the GLE parameters, we investigate the relative importance of memory, mass, and potential shape effects. Neglect of memory slows down NaCl and methane dissociation by roughly a factor of 2; neglect of mass accelerates reactions by a similar factor, and the harmonic approximation of the potential shape gives rise to slight acceleration. This partial error cancellation explains why Kramers' theory, which neglects memory effects and treats the potential shape in harmonic approximation, describes reaction rates better than more sophisticated theories. In essence, all three effects, friction memory, inertia, and the potential shape nonharmonicity, are important to quantitatively describe pair-reaction kinetics in water.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg²⁺ and SO₄ ^2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO₄ ^2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO₄ ^2- ions around 1100 cm^-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO₄ and for Na₂SO₄ in the presence of Mg²⁺ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na₂SO₄ solutions in the absence of Mg²⁺ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO₄ ^2--Mg²⁺ ion pairs.

Hydration Structure of Na+ and K+ Ions in Solution Predicted by Data-Driven Many-Body Potentials

The hydration structure of Na⁺ and K⁺ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na⁺ and K⁺ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

Nature of Cations Critically Affects Water at the Negatively Charged Silica Interface

Understanding the collective behavior of ions at charged surfaces is of paramount importance for geological and electrochemical processes. Ions screen the surface charge, and interfacial fields break the centro-symmetry near the surface, which can be probed using second-order nonlinear spectroscopies. The effect of electrolyte concentration on the nonlinear optical response has been semi-quantitatively explained by mean-field models based on the Poisson-Boltzmann equation. Yet, to explain previously reported ion-specific effects on the spectroscopic response, drastic ion-specific changes in the interfacial properties, including surface acidities and dielectric permittivities, or strong ion adsorption/desorption had to be invoked. Here, we use sum-frequency generation (SFG) spectroscopy to probe the symmetry-breaking of water molecules at a charged silica surface in contact with alkaline metal chloride solutions (LiCl, NaCl, KCl, and CsCl) at various concentrations. We find that the water response varies with the cation: the SFG response is markedly enhanced for LiCl compared to CsCl. We show that within mean-field models, neither specific ion-surface interactions nor a reduced dielectric constant of water near the interface can account for the variation of spectral intensities with cation nature. Molecular dynamics simulations confirm that the decay of the electrochemical potential only weakly depends on the salt type. Instead, the effect of different salts on the optical response is indirect, through the reorganization of the interfacial water: the salt-type-dependent alignment of water directly at the interface can explain the observations.

Is an Inductive Effect Explicit Account Required for Atomic Charges Aimed at Use within the Force Fields?

Polarization and inductive effects are the concepts that have been widely used in qualitative and even quantitative descriptions of experimentally observed properties in chemistry. The polarization effect has proven to be important in cases of biomolecular modeling though still the vast majority of molecular simulations use the classical non-polarizable force fields. In the last few decades, a lot of effort has been put into promoting the polarization effect and incorporating it into modern force fields and charge calculation methods. In contrast, the inductive effect has not attracted such attention and is effectively absent in both classic and modern force fields. Thus, a question is whether this difference corresponds to the difference in the physical significance of the effects and their explicit account, or is an artifact that should be corrected in the next generation of force fields. The significance of the electronic effects is studied in this paper through the prism of performance of specific models for atomic charge calculation that take into explicit account a nested set of effects: the formal charge, the nearest neighbors, the inductive effect, and finally the model, which takes into account all effects, which are possible to account for using atomic charges. The specific choice for the methods is the following: formal charges, MMFF94 bond charge increments, Dynamic Electronegativity Relaxation (DENR), and RESP. We propose a special scheme for the separate estimation of each particular effect contribution. By pairwise comparing the residual molecular electrostatic potential (MEP) errors of those charge models (aimed at best reproducing the quantum chemical reference MEP), we sequentially revealed how the account of each effect contributes to the better-quality MEP reproduction. The following relative importance of effects was estimated; thus, the natural hierarchy of the effects was established. First, the account of formal charges is of primordial importance. Second, the nearest neighbors account is the next in significance. Third, the explicit account of inductive effect in empirical charge calculation schemes was shown to significantly─both qualitatively and quantitatively─improve the quality of MEP reproduction. Fourth, the contribution of polarization is indirectly assessed. Surprisingly, it is of the order of magnitude of the inductive effect even for the molecular systems, for which it is anticipated to be more significant. Finally, the relative importance of anisotropic effects in neutral molecules was additionally reviewed.

Shining (Infrared) Light on the Hofmeister Series: Driving Forces for Changes in the Water Vibrational Spectra in Alkali-Halide Salt Solutions

The Hofmeister series is frequently used to rank ions based on their behavior from chaotropes ("structure breakers"), which weaken the surrounding hydrogen-bond network, to kosmotropes ("structure makers"), which enhance it. Here, we use fluctuation theory to investigate the energetic and entropic driving forces underlying the Hofmeister series for aqueous alkali-halide solutions. Specifically, we exploit the OH stretch infrared (IR) spectrum in isotopically dilute HOD/D₂O solutions as a probe of the effect of the salt on the water properties for different concentrations and choice of halide anion. Fluctuation theory is used to calculate the temperature derivative of these IR spectra, including decomposition of the derivative into different energetic contributions. These contributions are used to determine the thermodynamic driving forces in terms of effective internal energy and entropic contributions. This analysis implicates entropic contributions as the key factor in the Hofmeister series behavior of the OH stretch IR spectra, while the effective internal energy is nearly ion-independent.

Transport of thin water films: from thermally activated random walks to hydrodynamics

Under ambient atmospheric conditions, a thin film of water wets many solid surfaces, including insulators, ice, and salt. The film thickness as well as its transport behavior sensitively depend on the surrounding humidity. Understanding this intricate interplay is of highest relevance for water transport through porous media, particularly in the context of soil salinization induced by evaporation. Here, we use molecular simulations to evaluate the transport properties of thin water films on prototypical salt and soil interfaces, namely NaCl and silica solid surfaces. Our results showtwo distinct regimes for water transport: at low water coverage, the film permeance scales linearly with the adsorbed amount, in agreement with the activated random walk model.For thicker water films, the permeance scales as the adsorbed amount to the power of 3, in line with the Stokes equation. By comparing results obtained for silica and NaCl surfaces, we find that, at low water coverage, water permeance at the silica surface is considerably lower than at the NaCl surface, which we attribute to difference in hydrogen bonding. We also investigate the effect of atomic surface defects on the transport properties. Finally, in the context of water transport through porous material, we determine the humidity-dependent crossover between a vapor dominated and a thin film dominated transport regimes depending on the pore size.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.

Where Lennard-Jones Potentials Fail: Iterative Optimization of Ion-Water Pair Potentials Based on Ab Initio Molecular Dynamics Data

The use of the Lennard-Jones (LJ) potential in computer simulations of aqueous electrolyte solutions is widespread. The standard approach is to parametrize LJ potential parameters against thermodynamic solution properties, but problems in representing the local structural and dynamic properties of ion hydration shells remain. The r^-12-term in the LJ potential is responsible for this as it leads to overly repulsive ion-water interactions at short range. As a result, the LJ potential predicts blue-shifted vibrational peaks of the cations' rattling mode and too large negative ion hydration entropies. We demonstrate that cation-water effective pair potentials derived from ab initio MD data have softer short-range repulsions and represent hydration shell properties significantly better. Our findings indicate that replacing the LJ potential with these effective pair potentials offers a promising route to represent thermodynamic solution properties and local interactions of specific ions with nonpolarizable force field models.

Molecular dynamics simulations of the evaporation of hydrated ions from aqueous solution

Although important for atmospheric processes and gas-phase catalysis, very little is known about the hydration state of ions in the vapor phase. Here we study the evaporation energetics and kinetics of a chloride ion from liquid water by molecular dynamics simulations. As chloride permeates the interface, a water finger forms and breaks at a chloride separation of ≈ 2.8 nm from the Gibbs dividing surface. For larger separations from the interface, about 7 water molecules are estimated to stay bound to chloride in saturated water vapor, as corroborated by continuum dielectrics and statistical mechanics models. This ion hydration significantly reduces the free-energy barrier for evaporation. The effective chloride diffusivity in the transition state is found to be about 6 times higher than in bulk, which reflects the highly mobile hydration dynamics as the water finger breaks. Both effects significantly increase the chloride evaporation flux from the quiescent interface of an electrolyte solution, which is predicted from reaction kinetic theory.

Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E

Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D

Magnesium is essential in many vital processes. To correctly describe Mg2+ in physiological processes by molecular dynamics simulations, accurate force fields are fundamental. Despite the importance, force fields based on the commonly used 12-6 Lennard-Jones potential showed significant shortcomings. Recently progress was made by an optimization procedure that implicitly accounts for polarizability. The resulting microMg and nanoMg force fields (J. Chem. Theory Comput. 2021, 17, 2530-2540) accurately reproduce a broad range of experimental solution properties and the binding affinity to nucleic acids in TIP3P water. Since countless simulation studies rely on available water models and ion force fields, we here extend the optimization and provide Mg2+ parameters in combination with the SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D water models. For each water model, the Mg2+ force fields reproduce the solvation free energy, the distance to oxygens in the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, and the binding affinity and distance to the phosphate oxygens on nucleic acids. We present two parameter sets: MicroMg yields water exchange on the microsecond time scale and matches the experimental exchange rate. Depending on the water model, nanoMg yields accelerated water exchange in the range of 106 to 108 exchanges per second. The nanoMg parameters can be used to enhance the sampling of binding events, to obtain converged distributions of Mg2+, or to predict ion binding sites in biomolecular simulations. The parameter files are freely available at https://github.com/bio-phys/optimizedMgFFs.

Specific Ion Effects in Different Media: Current Status and Future Challenges

We discuss the current state of research as well as the future challenges for a deeper understanding of specific ion effects in protic and aprotic solvents as well as various additional media. Despite recent interest in solute or interfacial effects, we focus exclusively on the specific properties of ions in bulk electrolyte solutions. Corresponding results show that many mechanisms remain unknown for these simple media, although theoretical, computational, and experimental studies have provided some insights into explaining individual observations. In particular, the importance of local interactions and electronic properties is emphasized, which enabled a more consistent interpretation of specific ion effects over the past years. Despite current insufficient knowledge, we also discuss future challenges in relation to dynamic properties as well as the influence of different concentrations, different solvents, and solute contributions to gain a deeper understanding of specific ion effects for technological applications.

Molecular Simulations of Aqueous Electrolytes: Role of Explicit Inclusion of Charge Transfer into Force Fields

We describe here simulations of aqueous salt solutions that are performed using an explicit charge transfer force field. The emphasis of the discussion is on the calculation of a dynamical property of the solutions: self-diffusion of water. While force fields that are based on pairwise additive potentials or on potentials with explicit inclusion of polarization or with scaled charges can provide at best a qualitative agreement with experiments, force fields with explicit inclusion of charge transfer can produce quantitative agreement with experiment for NaCl and KCl solutions. We argue that a force field with explicit charge transfer contains new physics absent in the previously used force fields described in recent reviews of molecular simulations of aqueous electrolytes.

Determining Dielectric Constants for Complex Solvent Mixtures by Microwave Sensing and Model Prediction

The frequency-dependent dielectric constant is a basic fluid property that is currently challenging to determine for complex liquid mixtures. Here, we report the determination of effective dielectric constants for various solvent mixtures under flow conditions using a simple in-line microwave Fabry-Pérot interferometer cable sensor. An ideal solution model-based mixing rule has been established and demonstrated for significantly improved prediction of dielectric constants for single-phase solvent mixtures. However, the current mixing rules exhibit large deviations for immiscible water/oil dispersions apparently because of the effects of strong interfacial polarizations on the overall mixture polarizability that are not accounted for by the models.