Performance of density functional theory for describing hetero‐metallic active‐site motifs for methane‐to‐methanol conversion in metal‐exchanged zeolites

Methane‐to‐methanol conversion (MMC) can be facilitated with high methanol selectivities by copper‐exchanged zeolites. There are however two open questions regarding the use of these zeolites to facilitate the MMC process. The first concerns the possibility of operating the three cycles in the stepwise MMC process by these zeolites in an isothermal fashion. The second concerns the possibility of improving the methanol yields by systematic substitution of some copper centers in these active sites with other earth‐abundant transition metals. Quantum‐mechanical computations can be used to compare methane activation by copper oxide species and analogous mixed‐metal systems. To carry out such screening, it is important that we use theoretical methods that are accurate and computationally affordable for describing the properties of the hetero‐metallic catalytic species. We have examined the performance of 47 exchange‐correlation density functionals for predicting the relative spin‐state energies and chemical reactivities of six hetero‐metallic [M‐O‐Cu]2+ and [M‐O2‐Cu]2+, (where MCo, Fe, and Ni), species by comparison with coupled cluster theory including iterative single, double excitations as well as perturbative treatment of triple excitations, CCSD(T). We also performed multireference calculations on some of these systems. We considered two types of reactions (hydrogen addition and oxygen addition) that are relevant to MMC. We recommend the use of τ‐HCTH and OLYP to determine the spin‐state energy splittings in the hetero‐metallic motifs. ωB97, ωB97X, ωB97X‐D3, and MN15 performed best for predicting the energies of the hydrogen and oxygen addition reactions. In contrast, local, and semilocal functionals do poorly for chemical reactivity. Using [Fe‐O‐Cu]2+ as a test, we see that the nonlocal functionals perform well for the methane CH activation barrier. In contrast, the semilocal functionals perform rather poorly. © 2018 Wiley Periodicals, Inc. We have examined the performance of exchange‐correlation density functionals to predict the relative spin‐state energies, methane C‐H activation energies and chemical reactivities of hetero‐metallic analogues of copper oxide cores implicated in methane activation in zeolites. CCSD(T) and multireference calculations are used for benchmark data. τ‐HCTH and OLYP performwell for predicting the spinstate energy splittings. ωB97, ωB97X, ωB97XD3 and MN15 performed best for predicting the methane C‐H activation energy as wel