Fundamental Vibrational Analyses of the HCN Monomer, Dimer and Associated Isotopologues

In this work we provide high level ab initio treatments of the structures, vibrational frequencies, and electronic energies of the HCN monomer and dimer systems along with several isotopologues. The plethora of information related to this system within the literature is summarized and serves as a basis for comparison with the results of this paper. The geometry of the dimer and monomer are reported at the all electroncoupled-cluster singles, doubles, and perturbative triples level of theory [AE-CCSD(T)] with the correlation consistent quadruple-zeta quality basis sets with extra core functions (cc-pCVQZ) from Dunning. The theoretical geometries and electronic structures are further analyzed through the use of the Natural Bond Orbital (NBO) method and Natural Resonance Theory (NRT). At the AE-CCSD(T)/cc-pCVQZ level of theory, the full cubic with semi-diagonal quartic force field for nine dimer and four monomer isotopologues (the parent isotopologue along with N-15, C-13, and D derivatives) were obtained to treat the anharmonicity of the vibrations via second order vibrational perturbation theory (VPT2). Lastly, the enthalpy change associated with the formation of the dimer from two monomer units was determined using the focal point analysis. Computations including coupled-cluster through perturbative quadruples as well as basis sets up to six-zeta quality, including core functions (cc-pCVXZ, X=D,T,Q,5,6) were used to extrapolate to the AE-CCSDT(Q)/CBS energy associated with this hydrogen-bond forming process. After appending anharmonic zero-point vibrational, relativistic, and diagonal Born-Oppenheimer corrections, we report a value of -3.93 kcal mol(-1) for the enthalpy of formation. To our knowledge, each set of results (geometries, vibrational frequencies, and energetics) reported in this study represents the highest-level and most reliable theoretical predictions reported for this system.