Molecular simulation of liquid-vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Electrokinetic properties of NaCl solution via molecular dynamics simulations with scaled-charge electrolytes

Scaling ionic charges has become an alternative to polarizable force fields for representing indirect charge transfer effects in molecular simulations. In our work, we apply molecular dynamics simulations to investigate the properties of NaCl aqueous solutions in homogeneous and confined media. We compare classical integer- and scaled-charge force fields for the ions. In the bulk, we validate the force fields by computing equilibrium and transport properties and comparing them with experimental data. Integer-charge ions overestimate dielectric saturation and ionic association. Both force fields present an excess in ion-ion correlation, which leads to a deviation in the ionic conductivity at higher ionic strengths. Negatively charged quartz is used to simulate the confinement effect. Electrostatic interactions dominate counter-ion adsorption. Full-charge ions have stronger and more defined adsorption planes. We obtain the electroosmotic mobility of the solution by combining the shear plane location from non-equilibrium simulations with the ionic distribution from equilibrium simulations. From the Helmholtz-Smoluchowski equation, the zeta potential and the streaming potential coupling coefficient are computed. From an atomic-scale perspective, our molecular dynamics simulations corroborate the hypothesis of maximum packing of the Stern layer, which results in a stable and non-zero zeta potential at high salinity. The scaled-charge model representation of both properties is in excellent qualitative and quantitative agreement with experimental data. With our work, we demonstrate how useful and precise simple scaled-charge models for electrolytes can be to represent complex systems, such as the electrical double layer.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na⁺ and Cl^- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.

Entropy scaling of viscosity for molecular models of molten salts

Entropy scaling relates dynamic and thermodynamic properties by reducing the viscosity to a function of only the residual entropy. Molecular simulations are used to investigate the entropy scaling of the viscosity of three models of sodium chloride and five monovalent salts. Even though the correlation between the potential energy and the virial is weak, entropy scaling applies at liquid densities for all models and salts investigated. At lower densities, entropy scaling breaks down due to the formation of ion pairs and chains. Entropy scaling can be used to develop more extendable correlations for the dynamic properties of molten salts.

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.

Freezing point depression of salt aqueous solutions using the Madrid-2019 model

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl₂, and Li₂SO₄) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered.

Anomalous Phase Behaviors of Monolayer NaCl Aqueous Solutions Induced by Effective Coulombic Interactions within Angstrom-Scale Slits

Interests in subnanofluidic devices have called for molecular dynamics (MD) simulation studies of the thermodynamic behavior of monolayer salt solution within angstrom-scale slits. However, it still remains a grand challenge to accurately describe the Coulombic interactions by incorporating the effects of charge transfer and electronic dielectric screening. Herein, by using the electronic continuum model, where the effective ion charges are fine-tuned with a scaling factor of λ, we present simulation evidence that the effective Coulombic interactions among Na⁺/Cl^- ions can strongly affect the behavior of monolayer ionic aqueous solution. Our microsecond-scale MD simulations show that only the counterions with moderate effective charges (0.3 ≤ λ ≤ 0.8) can dissolve in monolayer water, whereas the high effective charges (λ ≥ 0.85) induce ions to assemble into monolayer nanocrystals, and ions with the low effective charges (λ ≤ 0.2) exhibit gas-like nanobubble. These findings could provide deeper insights into the physical chemistry behind subnanofluidic iontronic devices.

Vapor pressures and vapor phase compositions of choline chloride urea and choline chloride ethylene glycol deep eutectic solvents from molecular simulation

Despite the widespread acknowledgment that deep eutectic solvents (DESs) have negligible vapor pressures, very few studies in which the vapor pressures of these solvents are measured or computed are available. Similarly, the vapor phase composition is known for only a few DESs. In this study, for the first time, the vapor pressures and vapor phase compositions of choline chloride urea (ChClU) and choline chloride ethylene glycol (ChClEg) DESs are computed using Monte Carlo simulations. The partial pressures of the DES components were obtained from liquid and vapor phase excess Gibbs energies, computed using thermodynamic integration. The enthalpies of vaporization were computed from the obtained vapor pressures, and the results were in reasonable agreement with the few available experimental data in the literature. It was found that the vapor phases of both DESs were dominated by the most volatile component (hydrogen bond donor, HBD, i.e., urea or ethylene glycol), i.e., 100% HBD in ChClEg and 88%-93% HBD in ChClU. Higher vapor pressures were observed for ChClEg compared to ChClU due to the higher volatility of ethylene glycol compared to urea. The influence of the liquid composition of the DESs on the computed properties was studied by considering different mole fractions (i.e., 0.6, 0.67, and 0.75) of the HBD. Except for the partial pressure of ethylene glycol in ChClEg, all the computed partial pressures and enthalpies of vaporization showed insensitivity toward the liquid composition. The activity coefficient of ethylene glycol in ChClEg was computed at different liquid phase mole fractions, showing negative deviations from Raoult's law.

New Features of the Open Source Monte Carlo Software Brick-CFCMC: Thermodynamic Integration and Hybrid Trial Moves

We present several new major features added to the Monte Carlo (MC) simulation code Brick-CFCMC for phase- and reaction equilibria calculations (https://gitlab.com/ETh_TU_Delft/Brick-CFCMC). The first one is thermodynamic integration for the computation of excess chemical potentials (μ^ex). For this purpose, we implemented the computation of the ensemble average of the derivative of the potential energy with respect to the scaling factor for intermolecular interactions (). Efficient bookkeeping is implemented so that the quantity is updated after every MC trial move with negligible computational cost. We demonstrate the accuracy and reliability of the calculation of μ^ex for sodium chloride in water. Second, we implemented hybrid MC/MD translation and rotation trial moves to increase the efficiency of sampling of the configuration space. In these trial moves, short Molecular Dynamics (MD) trajectories are performed to collectively displace or rotate all molecules in the system. These trajectories are accepted or rejected based on the total energy drift. The efficiency of these trial moves can be tuned by changing the time step and the trajectory length. The new trial moves are demonstrated using MC simulations of a viscous fluid (deep eutectic solvent).