Finding Order in the Disordered Hydration Shell of Rapidly Exchanging Water Molecules around the Heaviest Alkali Cs+ and Fr+

Tracing mechanistic pathways and reaction kinetics toward equilibrium in reactive molten salts

In the dynamic environment of multi-component reactive molten salts, speciation unfolds as a complex process, involving multiple competing reaction pathways that are likely to face free energy barriers before reaching the reaction equilibria. Herein, we unravel intricate speciation in the AlCl3-KCl melt compositions with rate theory and ab initio molecular dynamics simulations. We find that the compositions with 100 and 50 mol% AlCl3 exclusively comprise neutral Al2Cl6 dimers and charged AlCl4- monomers, respectively. In intermediate AlCl3-KCl compositions, the chemical speciation proves to be a very complex process, requiring over 0.5 nanosecond to reach an equilibrium distribution of multiple species. It is a consequence of the competitive formation and dissociation of additional species, including charged Al dimers, trimers, and tetramers. Here, the species formation occurs through ion exchange events, which we explain by computing free energy landscapes and employing a Marcus-like rate theory. We show that both interspecies and intraspecies ion exchanges are probable and are dictated by the local structural reorganization reflected in the change of local coulombic fields. The species distributions are validated by comparing computed Raman spectra and neutron structure factors with the available experimental data. We find an excellent simulation-experiment agreement in both cases. Nevertheless, Raman spectroscopy turns out to be particularly advantageous for distinguishing between unique species distributions because of the distinct vibrational signatures of different species. The mechanistic insight into reaction dynamics gained in this study will be essential for the advancement of molten salts as reactive media in high-temperature energy applications.

Femtosecond Electron-Transfer Dynamics across the D2O/Cs+/Cu(111) Interface: The Impact of Hydrogen Bonding

Hydrogen bonding is essential in electron-transfer processes at water-electrode interfaces. We study the impact of the H-bonding of water as a solvent molecule on real-time electron-transfer dynamics across a Cs+-Cu(111) ion-metal interface using femtosecond time-resolved two-photon photoelectron spectroscopy. We distinguish in the formed water-alkali aggregates two regimes below and above two water molecules per ion. Upon crossing the boundary of these regimes, the lifetime of the excess electron localized transiently at the Cs+ ion increases from 40 to 60 fs, which indicates a reduced alkali-metal interaction. Furthermore, the energy transferred to a dynamic structural rearrangement due to hydration is reduced from 0.3 to 0.2 eV concomitantly. These effects are a consequence of H-bonding in the water-water interaction and the beginning formation of a nanoscale water network. This finding is supported by real-space imaging of the solvatomers and vibrational frequency shifts of the OH stretching and bending modes calculated for these specific interfaces.

Molecular Dynamics Simulation of the Synergistic Effect of Alkali/Surfactant/Polymer on the Formation and Stabilization of Water-Based Foam Systems

The stable maintenance effect of a chemical oil displacement agent on a foam liquid film usually creates problems with the oilfields surface system. To achieve comprehensive insights into the influence mechanism of these chemical agent components on the foam liquid film, an "SDBS/HPAM/OH^-" water-based foam simulation system and corresponding control systems were constructed by adjusting the categories and quantities of component molecules by molecular dynamics (MD) simulation. The simulated results indicated that the foam stability follows the order of "SDBS/HPAM/OH^-" system > "SDBS/HPAM" system > "SDBS" system. The smaller the inclination angle of the SDBS molecular tail chain, the greater the tendency of the SDBS molecular configuration to be "upright" at the gas-liquid interface, which is not conducive to preventing the aggregation and penetration of gas molecules at the gas-liquid interface. Although the presence of HPAM molecules can significantly enhance the stability of the liquid film by restricting the liquid film's drainage and the diffusion of gas molecules, the addition of HPAM molecules would weaken the formation ability of the foam liquid film. Through decreasing the aggregation of cations around the co-adsorption layer, OH^- not only enhances the interfacial activity of SDBS molecules, but also reduces the electrostatic repulsion between -COO^- groups on the HPAM molecular chain, which makes the foam more stable. With an increase in the pH, SDBS concentration, and HPAM concentration, the stability of foam liquid film was strengthened. These results are helpful in facilitating new insights into the formation and stabilization mechanism of water-based foams. In particular, they provide support for the development and application of new defoaming technologies.

A Holistic Approach for Elucidating Local Structure, Dynamics, and Speciation in Molten Salts with High Structural Disorder

To examine ion solvation, exchange, and speciation for minority components in molten salts (MS) typically found as corrosion products, we propose a multimodal approach combining extended X-ray absorption fine structure (EXAFS) spectroscopy, optical spectroscopy, ab initio molecular dynamics (AIMD) simulations, and rate theory of ion exchange. Going beyond conventional EXAFS analysis, our method can accurately quantify populations of different coordination states of ions with highly disordered coordination environments via linear combination fitting of the EXAFS spectra of these coordination states computed from AIMD to the experimental EXAFS spectrum. In a case study of dilute Ni(II) dissolved in the ZnCl₂+KCl melts, our method reveals heterogeneous distributions of coordination states of Ni(II) that are sensitive to variations in temperature and melt composition. These results are fully explained by the difference in the chloride exchange dynamics at varied temperatures and melt compositions. This insight will enable a better understanding and control of ion solubility and transport in MS.

Dynamic ionic radius of alkali metal ions in aqueous solution: a pulsed-field gradient NMR study

The dynamic behavior of alkali metal ions, Li⁺, Na⁺, K⁺, Rb⁺ and Cs⁺ in aqueous solutions is one of the most important topics in solution chemistry. Since these alkali metals contain nuclear magnetic resonance (NMR) active nuclei, it is possible to directly measure the diffusion constants of the alkali metal ions using the pulsed field gradient (PFG) NMR method. In this paper, the ⁷Li, ²³Na, ⁸⁷Rb, ¹³³Cs and ¹H resonances are observed for diffusion constants in aqueous solution and the solvent H₂O. Until now, the values of the diffusion constant have been lacking when discussing hydration effects around alkali metal ions. It is known that the static ionic radius (R_ion) increases with increasing the atomic number, and the experimental diffusion constants also increase with increasing the atomic number, which is opposite to the Stokes–Einstein (SE) relation. It suggests that alkali metal ions diffuse through a space of 10⁻⁶ m accompanying the hydrated spheres with a time interval of 10⁻³ s. For each alkali metal ion, the dynamic ionic radius is evaluated.

Stokes radius (dynamic ionic radius) of the alkali metal ions versus the ionic radius (R_ion) at 303 K. The dotted line is a guide for the 1 : 1 relation.

X-ray scattering reveals ion clustering of dilute chromium species in molten chloride medium

Enhancing the solar energy storage and power delivery afforded by emerging molten salt-based technologies requires a fundamental understanding of the complex interplay between structure and dynamics of the ions in the high-temperature media. Here we report results from a comprehensive study integrating synchrotron X-ray scattering experiments, ab initio molecular dynamics simulations and rate theory concepts to investigate the behavior of dilute Cr3+ metal ions in a molten KCl-MgCl2 salt. Our analysis of experimental results assisted by a hybrid transition state-Marcus theory model reveals unexpected clustering of chromium species leading to the formation of persistent octahedral Cr-Cr dimers in the high-temperature low Cr3+ concentration melt. Furthermore, our integrated approach shows that dynamical processes in the molten salt system are primarily governed by the charge density of the constituent ions, with Cr3+ exhibiting the slowest short-time dynamics. These findings challenge several assumptions regarding specific ionic interactions and transport in molten salts, where aggregation of dilute species is not statistically expected, particularly at high temperature.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

Hydration structure and water exchange kinetics at xenotime-water interfaces: implications for rare earth minerals separation

Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral-water interfaces. Here, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5-180 picoseconds), dictated by ion-water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.

Structure and dynamics of the molten alkali-chloride salts from an X-ray, simulation, and rate theory perspective

Molten salts are of great interest as alternative solvents, electrolytes, and heat transfer fluids in many emerging technologies.