“…This situation changed with the pioneering papers of Parkin and co-workers [26,27] written by Laikov [28]. We used the PBE functional [29], the TZV2P three-exponential basis set of Gauss-type functions for valence electrons [28] (basis set specification: {3,1}/{5,1} for H, {3,3,2}/{5,5,2} for C, N, and P and {5,5,4}/{9,9,8} for Mo), and the TZV2p (SBK-JC) relativistic pseudopotential [30] for core electrons of all the atoms of the molecular systems for which calculations were performed.…”