Diffusion of ions in supercritical water: Dependence on ion size and solvent density and roles of voids and necks

Hydrogen Bond Kinetics, Ionic Dynamics, and Voids in the Binary Mixtures of Protic Ionic Liquids with Alkanolamines

Classical molecular dynamics simulations were used to investigate the structural and dynamical properties of the mixtures of ionic liquids (ILs) with the conjugate forms of the cation in a 1:1 molar ratio. The experimental studies suggested the combination of ethanolamines and ILs as novel absorbents for acidic gases such as CO₂ and H₂S, which provide the advantage of efficient absorption of gases at low pressures. However, the microscopic properties of the ionic mixtures are not studied. From our computational investigations, the densities of mixtures are reported and compared with the experimental results. The structural evolution of mixtures is reported by radial distribution functions, coordination numbers, void analysis, and spatial distribution functions. The mixtures' dynamic properties were studied by analyzing the hydrogen bond, ion-pair, and ion-cage lifetimes of the system. Monoethanolammonium and triethanolammonium ILs show different types of spatial distribution functions. The cations have lesser effect on dynamics compared with anions. The charge on the anion greatly affects the dynamics of mixtures. The dianion mixtures show slower dynamics than the monoanionic mixtures. The hydrogen bonding between cations and anions is stronger than that between cations and neutral molecules due to strong coulombic attractive forces. The cations spend more time around the dianions as compared to monoanions. The distributions of voids show that the void sizes are smaller in triethanolamine-based mixtures. The sulfobenzoate-based mixtures show voids smaller than those of pyridine-3-carboxylate-based mixtures due to more available free space between the entities, which facilitates the overall dynamics.

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

ReaxFF molecular dynamics simulations of electrolyte-water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal-chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li⁺ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na⁺ and K⁺ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

Nanostructure domains, voids, and low-frequency spectra in binary mixtures of N,N-dimethylacetamide and ionic liquids with varying cationic size

Classical MD simulations were carried out on binary mixtures of DMA with hydroxide based ammonium ILs, TEAH, TPAH and TBAH, at three different mole fractions of IL.