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 BH 3 NH 3 dimer and BH 3 NH 3 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 BH 3 NH 3 ⋯CH 4 and BH 3 NH 3 ⋯NH 3 complexes. Apart from the B–H⋯H–N dihydrogen bond found in the BH 3 NH 3 crystal and dimmer, the B–H⋯H–X (XC, O, F) dihydrogen bonds have been observed in the BH 3 NH 3 ⋯HCN, BH 3 NH 3 ⋯H 2 O, BH 3 NH 3 ⋯CH 3 OH and BH 3 NH 3 ⋯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.9kJmol −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 BH 3 NH 3 ⋯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.