Heats of Formation on the Way from B2H6 to B20H16: Thermochemical Consequences of Multicenter Bonding in Ab initio and DFT Methods

The objective of this study is to evaluate the effectiveness of various computational methods in reproducing the experimental heats of formation of boron hydrides using the atomization energy ap- proach. The results have demonstrated that the empirical dispersion combined with the BJ damping function provided too large intramolecular dispersion energies, thereby compromising the accuracy of the outcomes produced by the DFT-D3 methods. Additionally, the CCSD(T) method has reproduced the experimental values only when combined with a basis set optimized for an accurate description of the core-valence correlation effect. Furthermore, the number of multicenter bonds present in the molecules under examination has also reflected their stability, as indicated by the heats of formation. Finally, a five-center two-electron (5c-2e) bond has emerged in B5H9, by applying the Intrinsic Bond Orbital (IBO) method.