Heats of formation on the way from B₂H₆ to B₂₀H₁₆: thermochemical consequences of multicenter bonding in <i>ab initio</i> 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 approach. 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 B₅H₉, by applying the intrinsic bond orbital (IBO) method.