“…Careful consideration was given to the possible occurrence of different spin states by checking the stability of the wave function and accounting for spin contamination (refer to Tables S1 and S2 in the Supporting Information for the expectation value of the total spin, ⟨ S 2 ⟩). The metal-zeolite systems were modeled by the M06-L meta-GGA functional , along with the 6-31g(d,p) basis for the Al, C, H, and O atoms and the Lanl2DZ basis set and effective core potential for the Ni and Si atoms (for computational efficiency). , Effective core potential offers significant computational savings and has been reported in the literature for understanding chemical reactivity in zeolite-supported metal clusters. ,− The M06-L functional has been parametrized for transition metal bonding , and employed to model metal clusters in zeolites. , In benchmarking calculations, Truhlar and co-workers have demonstrated that, for binding in zeolite clusters, the M06-L functional gives a mean unsigned error of only 0.87 kcal/mol . Furthermore, the M06-L functional was parametrized to include medium-range dispersion interactions. , All the atoms in the zeolite cluster, including the metal species, were allowed to relax, except for the capping H atoms that were frozen in Cartesian space (see the Supporting Information for details regarding the optimized structures).…”