Systematic analysis of electronic barrier heights and widths for concerted proton transfer in cyclic hydrogen bonded clusters: (HF)n, (HCl)n and (H2O)n where n=3, 4, 5

The MP2 and CCSD(T) methods are paired with correlation consistent basis sets as large as aug-cc-pVQZ to optimize the structures of the cyclic minima for (HF)(n), (HCl)(n) and (H2O)(n) where n = 3-5, as well as the corresponding transition states (TSs) for concerted proton transfer (CPT). MP2 and CCSD(T) harmonic vibrational frequencies confirm the nature of each minimum and TS. Both conventional and explicitly correlated CCSD(T) computations are employed to assess the electronic dissociation energies and barrier heights for CPT near the complete basis (CBS) limit for all 9 clusters. Results for (HF)(n) are consistent with prior studies identifying C-nh and D-nh point group symmetry for the minima and TSs, respectively. Our computations also confirm that CPT proceeds through C-s TS structures for the C-1 minima of (H2O)(3) and (H2O)(5), whereas the process goes through a TS with D-2d symmetry for the S-4 global minimum of (H2O)(4). This work corroborates earlier findings that the minima for (HCl)(3), (HCl)(4) and (HCl)(5) have C-3h, S-4 and C-1 point group symmetry, respectively, and that the C-nh structures are not minima for n = 4 and 5. Moreover, our computations show the TSs for CPT in (HCl)(3), (HCl)(4) and (HCl)(5) have D-3h, D-2d, and C-2 point group symmetry, respectively. At the CCSD(T) CBS limit, (HF)(4) and (HF)(5) have the smallest electronic barrier heights for CPT (approximate to 15 kcal mol(-1) for both), followed by the HF trimer (approximate to 21 kcal mol(-1)). The barriers are appreciably higher for the other clusters (around 27 kcal mol(-1) for (H2O)(4) and (HCl)(3); roughly 30 kcal mol(-1) for (H2O)(3), (H2O)(5) and (HCl)(4); up to 38 kcal mol(-1) for (HCl)(5)). At the CBS limit, MP2 significantly underestimates the CCSD(T) barrier heights (e.g., by ca. 2, 4 and 7 kcal mol(-1) for the pentamers of HF, H2O and HCl, respectively), whereas CCSD overestimates these barriers by roughly the same magnitude. Scaling the barrier heights and dissociation energies by the number of fragments in the cluster reveals strong linear relationships between the two quantities and with the magnitudes of the imaginary vibrational frequency for the TSs.