Ab Initio Study of Water Polarization in the Hydration Shell of Aqueous Hydroxide: Comparison between Polarizable and Nonpolarizable Water Models

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger-Weber potentials. The model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. The model enables simulation of early stage Al(OH)₃ nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

Slowing Down of the Molecular Reorientation of Water in Concentrated Alkaline Solutions

It is generally accepted that the hydroxide ion (OH^-) is a strong hydrogen bond acceptor and that its anomalously high diffusion constant in water results from a Grotthuss-like structural diffusion mechanism. However, the spatial extent over which OH^- ions influence the dynamics of the hydrogen-bond network of water remained largely unclear. Here, we measure the ultrafast dynamics of OH groups of HDO molecules interacting with the deuterated hydroxide ion OD^-. For solutions with OD^- concentrations up to 4 M, we find that HDO molecules that are not directly interacting with the ions have a reorientation time constant of ∼2.7 ps, similar to that of pure liquid water. When the concentration of OD^- ions is increased, the reorientation time constant increases, indicating a strong slowing down of the structural dynamics of the solution.

Investigation of the Structure of Concentrated NaOH Aqueous Solutions by Combining Molecular Dynamics and Wide-Angle X-ray Scattering

Classical molecular dynamics has been performed with explicit polarization on NaOH aqueous solutions from 0.5 mol L^-1 up to 9.7 mol L^-1. We adapted a force field of OH^- for polarizable simulation in order to reproduce the NaOH structural and thermodynamics properties in aqueous solutions. A good agreement between theoretical and experimental results has been found. Wide-angle X-ray scattering (WAXS) intensities issued from molecular dynamics are compared to experimental ones measured on Synchrotron facilities. The structure of the first coordination shell of Na⁺ has been studied to determine the variation of the oxygen number and hydroxide oxygen around the cation. In addition, Na⁺-OH^- McMillan-Mayer potential issued from molecular dynamics simulations has been calculated and allows for calculating Na⁺-OH^- pair association constant of 0.1 L mol^-1, which is in good agreement with the experiments.

Hydroxide Hydrogen Bonding: Probing the Solvation Structure through Ultrafast Time Domain Raman Spectroscopy

The mechanism of charge transport in aqueous media is critical in molecular, materials, and life sciences. The structure of the solvated hydroxide ion has been an area of some controversy. Polarization-resolved ultrafast time domain polarizability relaxation is used here to resolve the terahertz frequency Raman spectrum of hydroxide solutions. The measurements reveal the totally symmetric hydrogen-bond stretching (HO(-)···HOH) mode of the solvated hydroxide, permitting an experimental measurement of the bond force constant. The observed polarized Raman spectra are compared with those obtained from DFT calculations performed on HO(-)(H2O)n clusters. Good agreement between the observed frequency and the polarization dependence is found for the n = 3 or 4 clusters, particularly for those in which the solvating water molecules adopt a planar structure. The frequency of the symmetric stretch increases with concentration, consistent with an effect of ionic strength on either the H-bond or the structure of the cluster.

Proton transfer in concentrated aqueous hydroxide visualized using ultrafast infrared spectroscopy

While it is generally recognized that the hydroxide ion can rapidly diffuse through aqueous solution due to its ability to accept a proton from a neighboring water molecule, a description of the OH(-) solvation structure and mechanism of proton transfer to the ion remains controversial. In this report, we present the results of femtosecond infrared spectroscopy measurements of the O-H stretching transition of dilute HOD dissolved in NaOD/D(2)O. Pump-probe, photon echo peak shift, and two-dimensional infrared spectroscopy experiments performed as a function of deuteroxide concentration are used to assign spectral signatures that arise from the OH(-) ion and its solvation shell. A spectral feature that decays on a ∼110 fs time scale is assigned to the relaxation of transiently formed configurations wherein a proton is equally shared between a HOD molecule and an OD(-) ion. Over picosecond waiting times, features appear in 2D IR spectra that are indicative of the exchange of population between OH(-) ions and HOD molecules due to deuteron transfer. The construction of a spectral model that includes spectral relaxation, chemical exchange, and thermalization processes, and self-consistently treats all of our data, allows us to qualitatively explain the results of our experiments and gives a lower bound of 3 ps for the deuteron transfer kinetics.