A test of ab initio-generated, radial intermolecular potential energy functions for five axially-symmetric, hydrogen-bonded complexes BMIDLINE HORIZONTAL ELLIPSISHF, where B = N2, CO, PH3, HCN and NH3

The radial potential energy functions for five axially symmetric, hydrogen-bonded complexes BMIDLINE HORIZONTAL ELLIPSISHF (B = N-2, CO, PH3, HCN and NH3) have been calculated ab initio at the explicitly correlated level of theory CCSD(T)(F12c)/cc-pVTZ-F12 as a function of the hydrogen-bond distance r(ZMIDLINE HORIZONTAL ELLIPSISH), where Z is the hydrogen-bond acceptor atom of B. The remaining geometric parameters were optimised at each point in the calculation. The functions so generated were then used to estimate the spectroscopic constants omega(sigma), omega(sigma)x(sigma),alpha(sigma), and the v(sigma) = 1 <- 0 transition wavenumber associated with the intermolecular stretching mode nu(sigma) of each BMIDLINE HORIZONTAL ELLIPSISHF by using two equivalent approaches. Both involved the assumption that the vibrational modes of the B and HF molecules were sufficiently stiff relative to the intermolecular stretching mode that BMIDLINE HORIZONTAL ELLIPSISHF could be treated in a pseudo-diatomic approximation. One approach used derivatives of the potential evaluated at the distance r = r(e) while the other used the potential constants obtained by non-linear regression fits of three analytical functions (Morse, Rydberg and Hulburt-Hirschfelder) to the ab initio calculated points. The two approaches would give exactly the same results if the functions were a perfect fit. The H-H function gave the best fit. The determined spectroscopic constants were found to be in reasonable agreement with the limited number available by experiment. The efficacy of the approach was tested for the diatomic molecule (HCl)-Cl-35 by taking advantage of both an accurate RKR-type potential and an accurate set of spectroscopic constants. It was also established that the relationship D-e = k(sigma)/(2a(2)) between two measures (D-e and k(sigma), both calculated ab initio) of the strength of the hydrogen bond in the BMIDLINE HORIZONTAL ELLIPSISHF complexes (and required if the H-H function were an accurate representation of the BMIDLINE HORIZONTAL ELLIPSISHF potential functions) holds to an excellent level of approximation, and supports the conclusion that this function is appropriate to represent the hydrogen bond in the complexes investigated.

Automated summary

The study calculated the radial potential energy functions for five axially symmetric, hydrogen-bonded complexes, estimating spectroscopic constants and transition wavenumbers in reasonable agreement with experimental data. It was also confirmed that the H-H function effectively represents the hydrogen bond in the investigated complexes.