Optimizing Open Iron Sites in Metal–Organic ... - ACS Publications

Mar 31, 2017 - organic frameworks (MOFs) with open Fe sites are promising candidates for .... linkers that could be used to build MOFs and analyzed ho...
0 downloads 0 Views 2MB Size
Forum Article www.acsami.org

Optimizing Open Iron Sites in Metal−Organic Frameworks for Ethane Oxidation: A First-Principles Study Peilin Liao,†,‡ Rachel B. Getman,§ and Randall Q. Snurr*,† †

Department of Chemical & Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States School of Materials Engineering, Purdue University, West Lafayette, Indiana 47907, United States § Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States ‡

S Supporting Information *

ABSTRACT: Activation of the C−H bonds in ethane to form ethanol is a highly desirable, yet challenging, reaction. Metal− organic frameworks (MOFs) with open Fe sites are promising candidates for catalyzing this reaction. One advantage of MOFs is their modular construction from inorganic nodes and organic linkers, allowing for flexible design and detailed control of properties. In this work, we studied a series of single-metal atom Fe model systems with ligands that are commonly used as MOF linkers and tried to understand how one can design an optimal Fe catalyst. We found linear relationships between the binding enthalpy of oxygen to the Fe sites and common descriptors for catalytic reactions, such as the Fe 3d energy levels in different reaction intermediates. We further analyzed the three highest-barrier steps in the ethane oxidation cycle (including desorption of the product) with the Fe 3d energy levels. Volcano relationships are revealed with peaks toward higher Fe 3d energy and stronger electron-donating group functionalization of linkers. Furthermore, we found that the Fe 3d energy levels positively correlate with the electron-donating strength of functional groups on the linkers. Finally, we validated our hypotheses on larger models of MOF-74 iron sites. Compared with MOF-74, functionalizing the MOF-74 linkers with NH2 groups lowers the enthalpic barrier for the most endothermic step in the reaction cycle. Our findings provide insight for catalyst optimization and point out directions for future experimental efforts. KEYWORDS: metal−organic frameworks, DFT, catalyst screening, ethane, ethanol, nitrous oxide complex,12−14 which can convert methane to methyl bisulfate (CH3OSO3H) with good yield. The methyl bisulfate can then be reacted with water to form methanol. Challenges for this chemistry include the need to recover the SO2 byproduct and corrosion. Inspired by enzymes containing iron or copper in the active sites,15,16 another direction of active research is to introduce active sites containing transition metal elements to zeolites, such as Fe-ZSM-5,17,18 Cu-ZSM-5,19,20 and Co-ZSM5,21,22 as well as to investigate various iron-oxo structures23,24 and tricopper clusters.25,26 Some remaining challenges for these

1. INTRODUCTION Methane is the major component in natural gas (>84 vol %).1 It is highly desirable to develop efficient chemical processes that convert methane to industrially useful products, by functionalizing methane or forming longer chain organic molecules.2−7 However, activating the C−H bond in methane is very difficult because the C−H bond has a high bond dissociation energy of 431 kJ/mol.8 Compared to methane, the C−H bond in ethane is weaker (410 kJ/mol),8 and ethane is the second most abundant component (>2 vol %) found in natural gas.1 Thus, many researchers have focused on the slightly less difficult problem of C−H activation in ethane as a stepping stone toward eventual activation of methane. One way to activate the C−H bond in methane or ethane is through oxidation.9−11 One promising class of catalysts with high turnover frequencies for this reaction is the Pt bidiazine © XXXX American Chemical Society

Special Issue: Hupp 60th Birthday Forum Received: February 14, 2017 Accepted: March 31, 2017

A

DOI: 10.1021/acsami.7b02195 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces

Figure 1. (a) Ligand structures with complete and abbreviated names. Dashed lines in structures mark where atoms coordinate to Fe. Functional groups are indexed by numbers, where no prime symbol, a subsequent prime symbol, and a subsequent double prime symbol indicate para, meta, and ortho substitutions, respectively. Only single substitution on the linker is tested; therefore, when no ligand is specified for Rn (n = 1, 2, 3), Rn is the H atom. (b) Top view (looking from the missing ligand position toward the Fe site) for sample five-coordinated pyramidal A structures (see Figure 3) with labels underneath. Without and with a prime symbol after the ligand names correspond to anion ligands opposite to each other within the pyramidal plane (“PP” arrangement) and one of the anion ligands located at the pyramidal top (“PT” arrangement), respectively. Gray dashed circles indicate functional group substitutions. Color scheme: H, white; C, gray; N, blue; O, red; and Fe, purple. Figure S3 shows additional sample structures.

catalysts are difficulty in extracting methanol and the need to reactivate the catalysts, for example, with steam. Recently, Xiao et al.27 reported using two MOFs, Fe2(dobdc) (also known as Fe-MOF-74, dobdc4− = 2,5-dioxido-1,4benzenedicarboxylate) and Fe0.1Mg1.9(dobdc), with open Fe sites to catalyze the conversion of ethane to ethanol with N2O as the oxygen source (C2H6 + N2O → C2H5OH + N2). There are several advantages to using MOFs for catalysis. First, they provide a solid framework for supporting well dispersed catalytic Fe sites. Second, the combination of inorganic nodes and organic linkers enables the modular design of catalysts. Last but not least, specifically for the open Fe sites, Xiao et al. proposed that the MOF-74 framework structure stabilizes the high spin state of Fe, which is the spin state found in Fe enzyme catalysts for methane activation.28,29 Verma et al.30 published a follow-up computational study, where they applied Kohn− Sham density functional theory (KS-DFT) to study reactive intermediates and analyze the full reaction pathway. On the basis of these previous findings, it is promising to look further into open Fe sites in MOFs. Vogiatzis et al.31 searched through a database of MOFs32 to find those with coordinatively unsaturated Fe sites. In their stepwise screening procedure, many of the resulting 46 structures were manually discarded based on criteria such as likely instability. For the remaining 7 structures, DFT calculations were performed to assess the formation energy and stability of the high-spin iron(IV)-oxo intermediate, which is known to play a key role in the reaction mechanism. On the basis of their predictions, one of the materials was synthesized and tested experimentally. In this work, we adopted another, more “bottom-up” approach in which we created structures on the computer using common linkers that could be used to build MOFs and analyzed how the

different linkers affect the reaction energetics as a means to gain physical insights for improving the catalysts.

2. COMPUTATIONAL DETAILS The Gaussian 09 Rev. D.01 program33 was used to perform DFT calculations. Multireference methods may be necessary to provide accurate predictions for Fe with an unfilled 3d shell, ranging from 4 to 6 3d electrons for Fe(IV) to Fe(II). Prior work by Verma et al.30 indicates that KS-DFT with the M06-L exchange-correlation functional34 gives reasonable agreement with the correlated wave function method (multistate complete active space second-order perturbation theory (MS-CASPT2)) for Fe-containing MOFs. So we adopted the computational scheme used by Verma et al. for our DFT calculations, which also enables a direct comparison with their results. We used the M06-L exchange-correlation functional34 and the def2-TZVP35 allelectron basis sets with density fitting for all atoms. An ultrafine grid was employed to evaluate integrals. Previously, Verma et al.30 predicted an overall quintet S = 2 spin state as most stable in the Fe-MOF-74 cluster. We also found the S = 2 spin states to be the most stable for a few representative clusters (see Table S1) and set S = 2 for all structures reported in this forum article. We adopted N2O as the oxidant as in the previous experimental and computational work of Xiao et al.27 and Verma et al.30 Verma et al.30 used two structural models: one with 88 atoms and the other with 26 atoms. The larger 88-atom model, obtained by cutting a segment from MOF-74 with an 8 Å radius around the central Fe ion, contains three metal ions, while the smaller 26-atom model comprises one Fe and the first coordination sphere around it, i.e., methanol, formate, and 3-hydroxyprop-2-enoic acid. The baseline model that we use in this work comprises one Fe cation coordinated with five ligands, which consist of methanol and formate ions. We validated our findings using this model by comparing reaction energetics with those calculated using the same 88-atom model as Verma et al. (see section 3.5). B

DOI: 10.1021/acsami.7b02195 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces To explore the effect of different organic linkers, we constructed a series of Fe catalyst structures by coordinating different common MOF linkers32 to the Fe site. The construction started with building what we call “A” structures (see Figures 1 and 3), which are the initial states in the catalytic cycles. The Fe cation in the A structure is in the +2 oxidation state. Therefore, to ensure charge neutrality, building the five-coordinated A structure requires two −1 charged ligands and three neutral ligands. Inspired by the 26-atom model employed by Verma et al.,30 our two simplest structures, labeled Fm and Fm′ (see Figure 1), used methanol as the neutral ligand and formate anion as the ligand with −1 charge. We tested two possible arrangements for the two formate groups: one arrangement with the two groups directly opposite to each other within the plane of the pyramid (“PP” arrangement, see Figure 2) and the other with one of the anion ligands

similar trends when we tested the effects of allowing Fe and its NN atoms to relax (see section 1 and Figure S2 in the Supporting Information). In section 3.5, we calculated the reaction pathway for the functionalized 88-atom MOF-74 model in ref 30. In that case, as in ref 30, the reactive species, the Fe atom, and its NN O atoms were allowed to relax while keeping all other atoms fixed.

3. RESULTS AND DISCUSSION 3.1. Reaction Pathway of Ethane Oxidation. The reaction pathway studied in this work is shown in Figure 3.

Figure 2. Schematic illustration for the PP (“pyramid plane”) and PT (“pyramid top”) arrangements of ligands around an Fe center. located at the top of the pyramid (“PT” arrangement). To differentiate the different arrangements, we append a prime “ ′ ” after the names with the “PT” arrangement (see, for example, Figure 1, where Fm′ denotes the Fm structure with a PT arrangement for the formate groups). Then, we replaced either two of the methanol groups or the two formate groups in Fm and Fm′ with different common linkers. We built the Py series by replacing methanol with pyridine, and the BA and Im series by replacing formate with benzoate or imidazolate anions, respectively (see Figure 1 and Figure S3). We also introduced functional groups with various electron donating/withdrawing strengths to these linkers, to examine the effect of electron donating/withdrawing groups on the electronic structure of Fe and, hence, the reaction energetics. The functional groups tested are listed next to the illustration of linkers in Figure 1. The electron donating groups activate the ortho and para positions of aromatic rings, while the electron withdrawing groups activate the meta position. Accordingly, we placed electron donating groups at the para position and also tested a couple of structures at the ortho positions on the BA linker. One exception is for the COOH group. Even though it is slightly electron withdrawing, it is placed at the para position since it acts as a linking group in many MOFs. The other electron withdrawing group tested is −CN, which is placed at the meta position. For the Im group, we only tested one substituted location for functional groups. The index of the functional groups is ordered such that smaller index numbers label stronger electron donating groups. In MOFs, the linkers are somewhat constrained by the crystal structure. To mimic this in a simple way, we fixed the Fe-nearest neighbor (NN) bond lengths to be 2.06 Å, which is the average Fe−O bond length of the 26-atom MOF-74 model as determined by Verma et al.30 Each Fe-NN bond is set to be perpendicular to all other Fe-NN bonds. This is one simple approximate model. We tested the effects of Fe-NN bond length relaxation. We found that while individual Fe-NN bond distances vary, the average bond distance has a maximal difference of 3%. The relative spin stability is unchanged for the case of Im3 tested. Details can be found in the Supporting Information. After this initial construction of A structures, geometry relaxations were carried out on all atoms except Fe and its NN atoms. In all other structures in the reaction cycle, the optimized A structure is held fixed, allowing only the reactive molecular species on top of the open Fe sites to relax. This optimization scheme was used to provide fast screening and allow comparison across different ligands. We obtained

Figure 3. Catalytic reaction cycle for ethane oxidation with N2O on an open Fe site. The bottom dashed arrow indicates where steps that were included by Verma et al.30 were neglected in this work. Oxygen atoms are shown as nearest neighbor atoms to Fe, although for the Py/ Py′ and Im/Im′ series of structures, some of these O atoms would be changed to N atoms. Structure labels are inside boxes next to each structure.

N2O adsorbs on the Fe site, oxidizes Fe to iron oxo, and desorbs as N2, whereas ethane is oxidized by iron oxo to form ethanol. Figure 4 presents sample enthalpy changes along the reaction pathway (full energetics for all catalysts can be found in Table S2). Initial adsorption of N2O is exothermic (from “A” to “B” in Figure 3). After that, the Fe(II) cation is oxidized to form the iron(IV) oxo (FeO) (from “B” to “C”, via TS1). Desorption of N2 with adsorption of ethane causes very little change in enthalpy (from “C” to “D”, via “A_O” which corresponds to the isolated cluster without adsorption of N2 or ethane). The iron oxo then oxidizes ethane to ethanol (from “D” to “E”, via TS2). The last step is ethanol desorption to regenerate the “A” structure. Compared to the pathway reported by Verma et al.,30 the pathway in Figure 3 skips several steps between “D” and “E” that correspond to the formation of Fe−OH and a ·C2H5 radical, rotation of the OH, and rebound of ·C2H5 to form ethanol. According to Verma et al.,30 the enthalpies of activation and reaction for each of these steps are BA2 and BA3″ > BA3) when we moved the functional group from para to ortho positions (see Figure 8c). Overall, this provides additional support for our conclusion and suggests that we can derive design principles from simple models and extend them to more realistic systems.48



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.7b02195. Discussion of the effects of geometry optimization, energies of different spin states, reaction enthalpies, and additional structures (PDF) G

DOI: 10.1021/acsami.7b02195 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces



(16) Sazinsky, M. H.; Lippard, S. J. Methane Monooxygenase: Functionalizing Methane at Iron and Copper. In Sustaining Life on Planet Earth: Metalloenzymes Mastering Dioxygen and Other Chewy Gases; Kroneck, P. M. H., Sosa Torres, M. E., Eds.; Springer International Publishing: Switzerland, 2015; pp 205−256. (17) Hammond, C.; Dimitratos, N.; Jenkins, R. L.; Lopez-Sanchez, J. A.; Kondrat, S. A.; Hasbi ab Rahim, M.; Forde, M. M.; Thetford, A.; Taylor, S. H.; Hagen, H.; Stangland, E. E.; Kang, J. H.; Moulijn, J. M.; Willock, D. J.; Hutchings, G. J. Elucidation and Evolution of the Active Component within Cu/Fe/ZSM-5 for Catalytic Methane Oxidation: From Synthesis to Catalysis. ACS Catal. 2013, 3, 689−699. (18) Hammond, C.; Forde, M. M.; Ab Rahim, M. H.; Thetford, A.; He, Q.; Jenkins, R. L.; Dimitratos, N.; Lopez-Sanchez, J. A.; Dummer, N. F.; Murphy, D. M.; Carley, A. F.; Taylor, S. H.; Willock, D. J.; Stangland, E. E.; Kang, J.; Hagen, H.; Kiely, C. J.; Hutchings, G. J. Direct Catalytic Conversion of Methane to Methanol in an Aqueous Medium by Using Copper-Promoted Fe-ZSM-5. Angew. Chem., Int. Ed. 2012, 51, 5129−5133. (19) Vanelderen, P.; Hadt, R. G.; Smeets, P. J.; Solomon, E. I.; Schoonheydt, R. A.; Sels, B. F. Cu-ZSM-5: A Biomimetic Inorganic Model for Methane Oxidation. J. Catal. 2011, 284, 157−164. (20) Woertink, J. S.; Smeets, P. J.; Groothaert, M. H.; Vance, M. A.; Sels, B. F.; Schoonheydt, R. A.; Solomon, E. I. A [Cu2O]2+ Core in CuZSM-5, the Active Site in the Oxidation of Methane to Methanol. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 18908−18913. (21) Beznis, N. V.; Weckhuysen, B. M.; Bitter, J. H. Partial Oxidation of Methane Over Co-ZSM-5: Tuning the Oxygenate Selectivity by Altering the Preparation Route. Catal. Lett. 2010, 136, 52−56. (22) Kung, M. C.; Lin, S. S.-Y.; Kung, H. H. In Situ Infrared Spectroscopic Study of CH4 Oxidation Over Co−ZSM-5. Top. Catal. 2012, 55, 108−115. (23) Hohenberger, J.; Ray, K.; Meyer, K. The Biology and Chemistry of High-Valent Iron-Oxo and Iron-Nitrido Complexes. Nat. Commun. 2012, 3, 720. (24) Ray, K.; Pfaff, F. F.; Wang, B.; Nam, W. Status of Reactive NonHeme Metal−Oxygen Intermediates in Chemical and Enzymatic Reactions. J. Am. Chem. Soc. 2014, 136, 13942−13958. (25) Grundner, S.; Markovits, M. A. C.; Li, G.; Tromp, M.; Pidko, E. A.; Hensen, E. J. M.; Jentys, A.; Sanchez-Sanchez, M.; Lercher, J. A. Single-Site Trinuclear Copper Oxygen Clusters in Mordenite for Selective Conversion of Methane to Methanol. Nat. Commun. 2015, 6, 7546. (26) Liu, C.-C.; Mou, C.-Y.; Yu, S. S.-F.; Chan, S. I. Heterogeneous Formulation of the Tricopper Complex for Efficient Catalytic Conversion of Methane into Methanol at Ambient Temperature and Pressure. Energy Environ. Sci. 2016, 9, 1361−1374. (27) Xiao, D. J.; Bloch, E. D.; Mason, J. A.; Queen, W. L.; Hudson, M. R.; Planas, N.; Borycz, J.; Dzubak, A. L.; Verma, P.; Lee, K.; Bonino, F.; Crocellà, V.; Yano, J.; Bordiga, S.; Truhlar, D. G.; Gagliardi, L.; Brown, C. M.; Long, J. R. Oxidation of Ethane to Ethanol by N2O in a Metal−organic Framework with Coordinatively Unsaturated iron(II) Sites. Nat. Chem. 2014, 6, 590−595. (28) Price, J. C.; Barr, E. W.; Glass, T. E.; Krebs, C.; Bollinger, J. M. Evidence for Hydrogen Abstraction from C1 of Taurine by the HighSpin Fe(IV) Intermediate Detected during Oxygen Activation by Taurine:α-Ketoglutarate Dioxygenase (TauD). J. Am. Chem. Soc. 2003, 125, 13008−13009. (29) Sinnecker, S.; Svensen, N.; Barr, E. W.; Ye, S.; Bollinger, J. M.; Neese, F.; Krebs, C. Spectroscopic and Computational Evaluation of the Structure of the High-Spin Fe(IV)-Oxo Intermediates in Taurine: α-Ketoglutarate Dioxygenase from Escherichia Coli and Its His99Ala Ligand Variant. J. Am. Chem. Soc. 2007, 129, 6168−6179. (30) Verma, P.; Vogiatzis, K. D.; Planas, N.; Borycz, J.; Xiao, D. J.; Long, J. R.; Gagliardi, L.; Truhlar, D. G. Mechanism of Oxidation of Ethane to Ethanol at Iron(IV)−Oxo Sites in Magnesium-Diluted Fe2(dobdc). J. Am. Chem. Soc. 2015, 137, 5770−5781. (31) Vogiatzis, K. D.; Haldoupis, E.; Xiao, D. J.; Long, J. R.; Siepmann, J. I.; Gagliardi, L. Accelerated Computational Analysis of

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Randall Q. Snurr: 0000-0003-2925-9246 Notes

The authors declare the following competing financial interest(s): R.Q.S. has a financial interest in the start-up company NuMat Technologies, which is seeking to commercialize metal-organic frameworks.



ACKNOWLEDGMENTS This work was supported as part of the Inorganometallic Catalyst Design Center, an Energy Frontier Research Center funded by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under Award DESC0012702. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, and the Quest high performance computing facility at Northwestern University.



REFERENCES

(1) Faramawy, S.; Zaki, T.; Sakr, A. A.-E. Natural Gas Origin, Composition, and Processing: A Review. J. Nat. Gas Sci. Eng. 2016, 34, 34−54. (2) Lunsford, J. H. Catalytic Conversion of Methane to More Useful Chemicals and Fuels: A Challenge for the 21st Century. Catal. Today 2000, 63, 165−174. (3) Havran, V.; Duduković, M. P.; Lo, C. S. Conversion of Methane and Carbon Dioxide to Higher Value Products. Ind. Eng. Chem. Res. 2011, 50, 7089−7100. (4) Karakaya, C.; Kee, R. J. Progress in the Direct Catalytic Conversion of Methane to Fuels and Chemicals. Prog. Energy Combust. Sci. 2016, 55, 60−97. (5) Zakaria, Z.; Kamarudin, S. K. Direct Conversion Technologies of Methane to Methanol: An Overview. Renewable Sustainable Energy Rev. 2016, 65, 250−261. (6) Ismagilov, Z. R.; Matus, E. V.; Tsikoza, L. T. Direct Conversion of Methane on Mo/ZSM-5 Catalysts to Produce Benzene and Hydrogen: Achievements and Perspectives. Energy Environ. Sci. 2008, 1, 526−541. (7) Labinger, J. A.; Bercaw, J. E. Understanding and Exploiting C-H Bond Activation. Nature 2002, 417, 507−514. (8) Dean, J. A. Lange’s Handbook of Chemistry, 15th ed.; McGrawHill, Inc.: New York, 1998. (9) Guo, Z.; Liu, B.; Zhang, Q.; Deng, W.; Wang, Y.; Yang, Y. Recent Advances in Heterogeneous Selective Oxidation Catalysis for Sustainable Chemistry. Chem. Soc. Rev. 2014, 43, 3480−3524. (10) Tang, P.; Zhu, Q.; Wu, Z.; Ma, D. Methane Activation: The Past and Future. Energy Environ. Sci. 2014, 7, 2580−2591. (11) Schwarz, H. Chemistry with Methane: Concepts Rather than Recipes. Angew. Chem., Int. Ed. 2011, 50, 10096−10115. (12) Zimmermann, T.; Soorholtz, M.; Bilke, M.; Schüth, F. Selective Methane Oxidation Catalyzed by Platinum Salts in Oleum at Turnover Frequencies of Large-Scale Industrial Processes. J. Am. Chem. Soc. 2016, 138, 12395−12400. (13) Palkovits, R.; Antonietti, M.; Kuhn, P.; Thomas, A.; Schüth, F. Solid Catalysts for the Selective Low-Temperature Oxidation of Methane to Methanol. Angew. Chem., Int. Ed. 2009, 48, 6909−6912. (14) Periana, R. A.; Mironov, O.; Taube, D.; Bhalla, G.; Jones, C. J. Catalytic, Oxidative Condensation of CH4 to CH3COOH in One Step via CH Activation. Science 2003, 301, 814−818. (15) Tinberg, C. E.; Lippard, S. J. Dioxygen Activation in Soluble Methane Monooxygenase. Acc. Chem. Res. 2011, 44, 280−288. H

DOI: 10.1021/acsami.7b02195 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Forum Article

ACS Applied Materials & Interfaces Metal−Organic Frameworks for Oxidation Catalysis. J. Phys. Chem. C 2016, 120, 18707−18712. (32) Chung, Y. G.; Camp, J.; Haranczyk, M.; Sikora, B. J.; Bury, W.; Krungleviciute, V.; Yildirim, T.; Farha, O. K.; Sholl, D. S.; Snurr, R. Q. Computation-Ready, Experimental Metal-Organic Frameworks: A Tool To Enable High-Throughput Screening of Nanoporous Crystals. Chem. Mater. 2014, 26, 6185−6192. (33) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, N. J.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö .; Foresman, J. B.; Ortiz, J. V; Cioslowski, J.; Fox, D. J. Gaussian 09, revision D.01; Gaussian, Inc.: Wallingford, CT, 2009. (34) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Function. Theor. Chem. Acc. 2008, 120, 215−241. (35) Weigend, F. Accurate Coulomb-Fitting Basis Sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057−1065. (36) Norskov, J. K.; Bligaard, T.; Rossmeisl, J.; Christensen, C. H. Towards the Computational Design of Solid Catalysts. Nat. Chem. 2009, 1, 37−46. (37) Man, I. C.; Su, H.-Y.; Calle-Vallejo, F.; Hansen, H. A.; Martínez, J. I.; Inoglu, N. G.; Kitchin, J.; Jaramillo, T. F.; Nørskov, J. K.; Rossmeisl, J. Universality in Oxygen Evolution Electrocatalysis on Oxide Surfaces. ChemCatChem 2011, 3, 1159−1165. (38) Peterson, A. A.; Nørskov, J. K. Activity Descriptors for CO2 Electroreduction to Methane on Transition-Metal Catalysts. J. Phys. Chem. Lett. 2012, 3, 251−258. (39) Calle-Vallejo, F.; Loffreda, D.; M, K. T.; Sautet, P. Introducing Structural Sensitivity into Adsorption−energy Scaling Relations by Means of Coordination Numbers. Nat. Chem. 2015, 7, 403−410. (40) Trasatti, S. Work Function, Electronegativity, and Electrochemical Behaviour of Metals: III. Electrolytic Hydrogen Evolution in Acid Solutions. J. Electroanal. Chem. Interfacial Electrochem. 1972, 39, 163−184. (41) Greeley, J. Computational Screening for Improved Heterogeneous Catalysts and Electrocatalysts. In Heterogeneous Catalysis at Nanoscale for Energy Applications; Tao, F., Schneider, W. F., Kamat, P. V., Eds.; John Wiley & Sons, Inc.: New York, 2015; pp 139−159. (42) Boisen, A.; Dahl, S.; Nørskov, J. K.; Christensen, C. H. Why the Optimal Ammonia Synthesis Catalyst Is Not the Optimal Ammonia Decomposition Catalyst. J. Catal. 2005, 230, 309−312. (43) Suntivich, J.; May, K. J.; Gasteiger, H. A.; Goodenough, J. B.; Shao-Horn, Y. A Perovskite Oxide Optimized for Oxygen Evolution Catalysis from Molecular Orbital Principles. Science 2011, 334, 1383− 1385. (44) Ye, S.; Geng, C.-Y.; Shaik, S.; Neese, F. Electronic Structure Analysis of Multistate Reactivity in Transition Metal Catalyzed Reactions: The Case of C-H Bond Activation by Non-Heme Iron(iv)-Oxo Cores. Phys. Chem. Chem. Phys. 2013, 15, 8017−8030. (45) Andrikopoulos, P. C.; Michel, C.; Chouzier, S.; Sautet, P. In Silico Screening of Iron-Oxo Catalysts for CH Bond Cleavage. ACS Catal. 2015, 5, 2490−2499. (46) Bligaard, T.; Nørskov, J. K.; Dahl, S.; Matthiesen, J.; Christensen, C. H.; Sehested, J. The Brønsted−Evans−Polanyi

Relation and the Volcano Curve in Heterogeneous Catalysis. J. Catal. 2004, 224, 206−217. (47) Kazaryan, A.; Baerends, E. J. Ligand Field Effects and the High Spin−High Reactivity Correlation in the H Abstraction by Non-Heme Iron(IV)−Oxo Complexes: A DFT Frontier Orbital Perspective. ACS Catal. 2015, 5, 1475−1488. (48) Pellizzeri, S.; Jones, I. A.; Doan, H. A.; Snurr, R. Q.; Getman, R. B. Using Gas-Phase Clusters to Screen Porphyrin-Supported Nanocluster Catalysts for Ethane Oxidation to Ethanol. Catal. Lett. 2016, 146, 2566−2573.

I

DOI: 10.1021/acsami.7b02195 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX