Optimal geometries and harmonic vibrational frequencies of the global minima of water clusters (H2O)n, n = 2–6, and several hexamer local minima at the CCSD(T) level of theory
We report the first optimum geometries and harmonic vibrational frequencies for the ring pentamer and several water hexamer (prism, cage, cyclic and two book) at the CCSD(T)/aug-cc-pVDZ level of theory. All five hexamer isomer minima previously reported by MP2 are also minima on the CCSD(T) potential energy surface (PES). In addition, all CCSD(T) minimum energy structures for the n=2-6 cluster isomers are quite close to the ones previously obtained by MP2 on the respective PESs, as confirmed by a modified Procrustes analysis that quantifies the difference between any two cluster geometries. The CCSD(T) results confirm the cooperative effect of the homodromic ring networks (systematic contraction of the nearest-neighbor (nn) intermolecular separations with cluster size) previously reported by MP2, albeit with O-O distances shorter by ~0.02 Å, indicating that MP2 overcorrects this effect. The harmonic frequencies at the minimum geometries were obtained by the double differentiation of the CCSD(T) energy using an efficient scheme based on internal coordinates that reduces the number of required single point energy evaluations by ~15% when compared to the corresponding double differentiation using Cartesian coordinates. Negligible differences between MP2 and CCSD(T) are found for the librational modes, while uniform increases of ~15 and ~25 cm-1 are observed for the bending and “free” OH harmonic frequencies. The largest differences between MP2 and CCSD(T) are observed for the harmonic hydrogen bonded frequencies. The CCSD(T) red shifts from the monomer frequencies (Δω) are smaller than the MP2 ones, due to the fact that the former produces shorter elongations (ΔR) of the respective hydrogen bonded OH lengths from the monomer value with respect to the latter. Both the MP2 and CCSD(T) results for the hydrogen bonded frequencies were found to closely follow the relation - Δω = s · ΔR, with a rate of s = 20.3 cm-1 / 0.001 Å. The CCSD(T) harmonic frequencies, when corrected using the MP2 anharmonicities obtained from second order vibrational perturbation theory (VPT2), produce anharmonicCCSD(T) estimates that are within < 60 cm-1 from the measured infrared (IR) active bands of the n=2-6 clusters and furthermore trace the observed red shifts with respect to the monomer (Δν) quite accurately. The energetic order between the various hexamer isomers on the PES (prism has the lowest energy) previously reported at MP2 was found to be preserved at the CCSD(T) level, whereas the inclusion of anharmonic corrections further stabilizes the cage among the hexamer isomers.
- Research Organization:
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC05-76RL01830
- OSTI ID:
- 1094928
- Report Number(s):
- PNNL-SA-97024; JCPSA6; KC0301020
- Journal Information:
- Journal of Chemical Physics, Vol. 139, Issue 11; ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)
- Country of Publication:
- United States
- Language:
- English
Similar Records
High-level ab-initio calculations for the four low-lying families of minima of (H2O)20: II. Spectroscopic signatures of the dodecahedron, fused cubes, face-sharing pentagonal prisms, and edge-sharing pentagonal prisms hydrogen bonding networks
Development of Transferable Interaction Models for Water: II. Accurate Energetics of the First Few Water Clusters from First Principles