Surface of aqueous solutions of alkali halides: layer by layer analysis

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.

Surface Affinity of Alkali and Halide Ions in Their Aqueous Solution: Insight from Intrinsic Density Analysis

The surface tension of all aqueous alkali halide solutions is higher than that of pure water. According to the Gibbs adsorption equation, this indicates a net depletion of these ions in the interfacial region. However, simulations and experiments show that large, soft ions, such as I^-, can accumulate at the liquid/vapor interface. The presence of a loose hydration shell is usually considered to be the reason for this behavior. In this work, we perform computer simulations to characterize the liquid-vapor interface of aqueous alkali chloride and sodium halide solutions systematically, considering all ions from Li⁺ to Cs⁺ and from F^- to I^-. Using computational methods for the removal of surface fluctuations, we analyze the structure of the interface at a dramatically enhanced resolution, showing that the positive excess originates in the very first molecular layer and that the next 3-4 layers account for the net negative excess. With the help of a fictitious system with charge-inverted ion pairs, we also show that it is not possible to rationalize the surface affinity of ions in solutions in terms of the properties of anions and cations separately. Moreover, the surface excess is generally dominated by the smaller of the two ions.

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

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.

A force field of Li+, Na+, K+, Mg2+, Ca2+, Cl-, and SO4 2- in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions

In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ and the anions Cl^- and SO₄ ^2-. These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.