# 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

## Abstract

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}more »

- Authors:

- Publication Date:

- Research Org.:
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States)

- Sponsoring Org.:
- USDOE

- OSTI Identifier:
- 1094928

- Report Number(s):
- PNNL-SA-97024

Journal ID: ISSN 0021-9606; JCPSA6; KC0301020

- DOE Contract Number:
- AC05-76RL01830

- Resource Type:
- Journal Article

- Journal Name:
- Journal of Chemical Physics

- Additional Journal Information:
- Journal Volume: 139; Journal Issue: 11; Journal ID: ISSN 0021-9606

- Publisher:
- American Institute of Physics (AIP)

- Country of Publication:
- United States

- Language:
- English

### Citation Formats

```
Miliordos, Evangelos, Aprà, Edoardo, and Xantheas, Sotiris S.
```*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*. United States: N. p., 2013.
Web. doi:10.1063/1.4820448.

```
Miliordos, Evangelos, Aprà, Edoardo, & Xantheas, Sotiris S.
```*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*. United States. doi:10.1063/1.4820448.

```
Miliordos, Evangelos, Aprà, Edoardo, and Xantheas, Sotiris S. Tue .
"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". United States. doi:10.1063/1.4820448.
```

```
@article{osti_1094928,
```

title = {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},

author = {Miliordos, Evangelos and Aprà, Edoardo and Xantheas, Sotiris S.},

abstractNote = {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.},

doi = {10.1063/1.4820448},

journal = {Journal of Chemical Physics},

issn = {0021-9606},

number = 11,

volume = 139,

place = {United States},

year = {2013},

month = {1}

}