“…[19] The DFT calculations have been carried out using the Gaussian 09 program [20] employing the hybrid B3LYP exchange-correlation functional [21] with the unrestricted Kohn-Sham solution, [22] and TZVP basis sets. [23] The unrestricted B3LYP/TZVP level of theory proved reliable in previous studies of cationic, [24] anionic, [25] and neutral [10] AlVO 5 À , [27] V 2 SiO 8 À , [28] V 2 SiO 7 + , [11d] and V 3 PO 10 +[11a,c] and their gas-phase reactions with small alkanes. For the optimization of transition structures (TS), we employed either the Berny algorithm [29] or the synchronous transitguided quasi-Newton (STQN) method.…”