Structure and solvation dynamics of the hydroxide ion in ice-like water clusters: a CCSD(T) and car–parrinello molecular dynamics study

Abstract

<p>Using MP2, CCSD(T) electronic structure theory and <em>ab initio</em> molecular dynamics simulations, we explore the structure, solvation dynamics and vibrational spectra of OH<small>⁻</small>(H<small>₂</small>O)<small>_<em>n</em></small> clusters. Our study reports new cubic and fused cubic global minima structures of OH<small>⁻</small>(H<small>₂</small>O)<small>_<em>n</em></small> for <em>n</em> = 8–26 with surface and interior solvation arrangements. In the case of OH<small>⁻</small>(H<small>₂</small>O)<small>₂₆</small>, we show that MP2 and CCSD(T) calculations predict global minima structures with the hydroxide ion occupying the interior region of a densely packed cubic cluster that is secured by ionic hydrogen bonds. More importantly, results from <em>ab initio</em> molecular dynamics simulations of OH<small>⁻</small>(H<small>₂</small>O)<small>₂₆</small> demonstrate that the hydroxide ion remains in the cluster interior and hexa-coordinated, irrespective of the temperature, up to around 175 K, then incrementally transitions from a surface-exposed penta- (170–200 K), to a tetra- (225 K) to a tri-coordinated OH<small>⁻</small>(H<small>₂</small>O)<small>₃</small> structure at 300 K. Building on our temperature-dependent vibrational power spectra, we are also able to disentangle structure and temperature effects on individual spectral contributions arising from water molecules located in the inner and outer shell of OH<small>⁻</small>(H<small>₂</small>O)<small>₂₆</small>. Some of these theoretical results provide valuable guidance for the interpretation of IRMPD spectra of small hydroxide-water clusters, but there are also several intriguing implications of these results, in particular, for the solvation of the OH<small>⁻</small> ion at the surface of water nanodroplets and aqueous interfaces.</p>