Non-convertional hydrogen bonding interaction of BH3NH3 complexes: a comparative theoretical study

A new type of hydrogen bond, called a dihydrogen bond, has recently been introduced. In this bond hydrogen is donated to (hydridic) hydrogen. In this paper, ab initio HF, MP2 and DFT(B3LYP) levels of theory with different basis sets in combination with counterpoise procedure for basis set superposition error correction have been applied to BH3NH3 dimer and BH3NH3 complexes of methane, hydrogen cyanide, ammonia, water, methanol and hydrogen fluoride to understand the features of dihydrogen bond. The optimized geometric parameters and interaction energies for various isomers at different levels are estimated. The structures obtained at different computational levels are in agreement with each other. Dihydrogen bond does not occur in both BH3NH3⋯CH4 and BH3NH3⋯NH3 complexes. Apart from the B-H⋯H-N dihydrogen bond found in the BH3NH3 crystal and dimmer, the B-H⋯H-X (XC, O, F) dihydrogen bonds have been observed in the BH3NH3⋯HCN, BH3NH3⋯H2O, BH3NH3⋯CH3OH and BH3NH3⋯HF complexes, while the classic H bonds also exist in the last three complexes. As for the complexes in which only dihydrogen bonds appear the strength of dihydrogen bonds ranges from 17.9 to 18.9 kJ mol−1 at B3LYP/6-311++g(d,p) level. Binding energies obtained from the MP2 and B3LYP optimized structures are more sensitive to basis sets than those from the HF method. Larger basis functions generally tend to produce slightly longer intermolecular distances, and the B3LYP and MP2 methods generate shorter intermolecular distances though they usually produce longer bond lengths compared with those at the HF level. The infrared spectrum frequencies, IR intensities and the vibrational frequency shifts are reported. Finally the solution phase studies on BH3NH3⋯HF complex are also carried out using the Onsager reaction field model with a range of dielectric constants from 2 to 80 at B3LYP/6-311++g(d,p) level.