Molecular Force Field Development for Aqueous Electrolytes: 2. Polarizable Models Incorporating Crystalline Chemical Potential and Their Accurate Simulations of Halite, Hydrohalite, Aqueous Solutions of NaCl, and Solubility

Molecular Modeling of Water-in-Salt Electrolytes: A Comprehensive Analysis of Polarization Effects and Force Field Parameters in Molecular Dynamics Simulations

Accurate modeling of highly concentrated aqueous solutions, such as water-in-salt (WiS) electrolytes in battery applications, requires proper consideration of polarization contributions to atomic interactions. Within the force field molecular dynamics (MD) simulations, the atomic polarization can be accounted for at various levels. Nonpolarizable force fields implicitly account for polarization effects by incorporating them into their van der Waals interaction parameters. They can additionally mimic electron polarization within a mean-field approximation through ionic charge scaling. Alternatively, explicit polarization description methods, such as the Drude oscillator model, can be selectively applied to either a subset of polarizable atoms or all polarizable atoms to enhance simulation accuracy. The trade-off between simulation accuracy and computational efficiency highlights the importance of determining an optimal level of accounting for atomic polarization. In this study, we analyze different approaches to include polarization effects in MD simulations of WiS electrolytes, with an example of a Na-OTF solution. These approaches range from a nonpolarizable to a fully polarizable force field. After careful examination of computational costs, simulation stability, and feasibility of controlling the electrolyte properties, we identify an efficient combination of force fields: the Drude polarizable force field for salt ions and non-polarizable models for water. This cost-effective combination is sufficiently flexible to reproduce a broad range of electrolyte properties, while ensuring simulation stability over a relatively wide range of force field parameters. Furthermore, we conduct a thorough evaluation of the influence of various force field parameters on both the simulation results and technical requirements, with the aim of establishing a general framework for force field optimization and facilitating parametrization of similar systems.

Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides

To assess the salt solubilities of six alkali halides in aqueous systems, we proposed a thermodynamic cycle and an efficient molecular modeling methodology. The Gibbs free energy changes for vaporization, dissociation, and dissolution were calculated using the experimental data of ionic thermodynamic properties obtained from the NBS tables. Additionally, the Marcus' and Tissandier's solvation free energy data for Li+, Na+, K+, Cl-, and Br- ions were compared with the conventional solvation free energies by substituting into our self-consistent thermodynamic cycle. Furthermore, Tissandier's absolute solvation free energy data were used as the training set to refit the Lennard-Jones parameters of OPLS-AA force field for ions. To predict salt solubilities, an assumption of a pseudo-solvent was proposed to characterize the coupling work of a solute with its environment from infinitely diluted to saturated solutions, indicating that the Gibbs energy change of solvation process is a function of ionic strength. Following the self-consistency of the cycle, the newly derived formulas were used to determine the salt solubilities by interpolating the intersection of Gibbs free energy of dissolution and the zero free energy line. The refined ion parameters can also predict the structure and thermodynamic properties of aqueous electrolyte solutions, such as densities, pair correlation functions, hydration numbers, mean activity coefficients, vapor pressures, and the radial dependences of the net charge at 298.15 K and 1 bar. Our method can be used to characterize the solid-liquid equilibria of ions or charged particles in aqueous systems. Furthermore, for highly concentrated strong electrolyte systems, it is essential to introduce accurate water models and polarizable force fields.

Molecular dynamics of preferential adsorption in mixed alkali–halide electrolytes at graphene electrodes

Understanding the microscopic behavior of aqueous electrolyte solutions in contact with graphene and related carbon surfaces is important in electrochemical technologies, such as capacitive deionization or supercapacitors. In this work, we focus on preferential adsorption of ions in mixed alkali–halide electrolytes containing different fractions of Li+/Na+ or Li+/K+ and/or Na+/K+ cations with Cl− anions dissolved in water. We performed molecular dynamics simulations of the solutions in contact with both neutral and positively and negatively charged graphene surfaces under ambient conditions, using the effectively polarizable force field. The simulations show that large ions are often intuitively attracted to oppositely charged electrodes. In contrast, the adsorption behavior of small ions tends to be counterintuitive. In mixed-cation solutions, one of the cations always supports the adsorption of the other cation, while the other cation weakens the adsorption of the first cation. In mixed-cation solutions containing large and small cations simultaneously, adsorption of the larger cations varies dramatically with the electrode charge in an intuitive way, while adsorption of the smaller cations changes oppositely, i.e., in a counterintuitive way. For (Li/K)Cl mixed-cation solutions, these effects allow the control of Li+ adsorption by varying the electrode charge, whereas, for LiCl single-salt solutions, Li+ adsorption is nearly independent of the electrode charge. We rationalize this cation–cation lever effect as a result of a competition between three driving forces: (i) direct graphene–ion interactions, (ii) the strong tendency of the solutions to saturate the network of non-covalent intermolecular bonds, and (iii) the tendency to suppress local charge accumulation in any region larger than typical interparticle distances. We analyze the driving forces in detail using a general method for intermolecular bonding based on spatial distribution functions and different contributions to the total charge density profiles. The analysis helps to predict whether an ion is more affected by each of the three driving forces, depending on the strength of the ion solvation shells and the compatibility between the contributions of the charge density profiles due to the ion and water molecules. This approach is general and can also be applied to other solutions under different thermodynamic conditions.

Atomistic simulation framework for molten salt vapor–liquid equilibrium prediction and its application to NaCl

Knowledge of the vapor–liquid equilibrium (VLE) properties of molten salts is important in the design of thermal energy storage systems for solar power and nuclear energy production applications. The high temperatures involved make their experimental determination problematic, and the development of both macroscopic thermodynamic correlations and predictive molecular-based methodologies are complicated by the requirement to appropriately incorporate the chemically reacting vapor-phase species. We derive a general thermodynamic-based atomistic simulation framework for molten salt VLE prediction and show its application to NaCl. Its input quantities are temperature-dependent ideal-gas free energy data for the vapor phase reactions and density and residual chemical potential data for the liquid. If these are not available experimentally, the former may be predicted using standard electronic structure software, and the latter may be predicted by means of classical atomistic simulation methodology. The framework predicts the temperature dependence of vapor pressure, coexisting phase densities, vapor phase composition, and vaporization enthalpy. It also predicts the concentrations of vapor phase species present in minor amounts (such as the free ions), quantities that are extremely difficult to measure experimentally. We furthermore use the results to obtain an approximation to the complete VLE binodal dome and the critical properties. We verify the framework for molten NaCl, for which experimentally based density and chemical potential data are available in the literature. We then apply it to the analysis of NaCl simulation data for two commonly used atomistic force fields. The framework can be readily extended to molten salt mixtures and to ionic liquids.

Single-Particle Dynamics at the Intrinsic Surface of Aqueous Alkali Halide Solutions

The distribution of ions in the proximity of the liquid-vapor interface of their aqueous solution has been the subject of an intense debate during the last decade. The effects of ionic polarizability have been one of its salient aspects. Much less has been said about the corresponding dynamical properties, which are substantially unexplored. Here, we investigate the single-particle dynamics at the liquid-vapor interface of several alkali halide solutions, using molecular dynamics simulations with polarizable and nonpolarizable force fields and intrinsic surface analysis. We analyze the diffusion coefficient, residence time, and velocity autocorrelation function of water and ions and investigate how these properties depend on the molecular layer where they reside. While anions are found in the first molecular layer for relatively long times, cations are only making quick excursions into it, thanks to thermal fluctuations. The in-layer residence time of ions and their molar fraction in the layer turned out to be linearly dependent on each other. We interpret this unexpected result using a simple two-state model. In addition, we found that, unlike water and other neat molecular liquids that show a different diffusion mechanism at the surface than in the bulk of their liquid phase, ions do not enjoy enhanced mobility in the surface layer of their aqueous solution. This result indicates that ions in the surface layer are shielded by their nearest water neighbors from being exposed to the vapor phase as much as possible. Such positions are available for the ions at the negatively curved troughs of the molecularly rugged liquid surface.