Wen-zhi Luoa,
Guang-hui Chen*a,
Song-tao Xiao*b,
Qiang Wangc,
Ze-kun Huangd and
Ling-yu Wangb
aDepartment of Chemistry, Shantou University, Guangdong 515063, China. E-mail: ghchen@stu.edu.cn
bInstitute of Radiochemistry, China Institute of Atomic Energy (CIAE), Beijing, 102413, People's Republic of China. E-mail: xiao200112@163.com
cDepartment of Applied Chemistry, College of Science, Nanjing Tech University, Nanjing 211816, People's Republic of China. E-mail: wangqiang@njtech.edu.cn
dThe Wolfson Department of Chemical Engineering, Israel Institute of Technology, Israel. E-mail: zekun-h@campus.technion.ac.il
First published on 30th July 2019
Enzymatic heme and non-heme Fe(IV)–O species usually play an important role in hydrogen abstraction of biocatalytic reactions, yet duplicating the reactivity in biomimicry remains a great challenge. Based on Xiao et al.'s experimental work [Nat. Chem., 2014, 6(7), 590], we theoretically found that in the presence of the oxidant N2O, the enzyme-like metal organic framework, i.e., magnesium-diluted Fe-MOF-74 [Fe/(Mg)-MOF-74] can activate the C–H bonds of 1,4-cyclohexadiene (CHD) into benzene with a two-step hydrogen abstraction mechanism based on the density functional theory (DFT) level. It is shown that the first transition state about the cleavage of the N–O bond of N2O to form the Fe(IV)–O species is the rate-determining step with activation enthalpy of 19.4 kcal mol−1 and the complete reaction is exothermic by 62.8 kcal mol−1 on quintet rather than on triplet PES. In addition, we proposed a rebound mechanism of cyclic cyclohexane (CHA) hydroxylation to cyclohexanol which has not been studied experimentally. Note that the activation enthalpies on the first hydrogen abstraction for both cyclic CHD and cyclohexane are just 8.1 and 3.5 kcal mol−1, respectively, which are less than that of 13.9 kcal mol−1 for chained ethane. Most importantly, for the hydrogen abstraction of methane catalyzed by M/(Mg)-MOF-74 (M = Cu, Ni, Fe, and Co), we found that the activation enthalpies versus the C–H bond length of methane of TSs, NPA charge of the reacting oxyl atom have linear relationships with different slopes, i.e., shorter C–H bond and less absolute value of NPA charge of oxyl atom are associated with lower activation enthalpy; while for the activation of methane, ethane, propane and CHD catalyzed by Fe/(Mg)-MOF-74, there also exists positive correlations between activation enthalpies, bond dissociation energies (BDEs) and C–H bond lengths in TSs, respectively. We hope the present theoretical study may provide the guideline to predict the performance of MOFs in C–H bond activation reactions.
Metal organic frameworks (MOFs) are important porous crystalline materials that are composed of inorganic metal sites and organic ligands, which has been another vital structure that can be used for heterogeneous catalysis with Fe or other metals.14–19 The inner pores of these MOFs can determine reaction selectivity, and significantly, the structure of the pores can be changed by adjusting the organic ligands and metal sites. Like heme and nonheme enzymes of Fe(IV)–O intermediate, another attractive feature of MOFs is that redox active metals can be utilized as activity center, which can perform behavior similar to that in metalloenzymes, possibly exhibiting powerful enzyme-like catalytic activity. The MOF with isolated terminal Fe(IV)–O is a highly promising area of research. In addition, although the Fe(IV)–O intermediates generally make use of nitrogen-based chelating ligands, the metal center in MOF is also often linked by weak-field ligands, such as aryl oxides and carboxylates, which are surrounded in the coordination environment of framework structure. Thus, these materials may also have special electronic properties and reactivity imparted by their unique coordination environment. Indeed, C–H bond activation of chained hydrocarbons by MOFs have received extensive attentions. For example, Dale R. Pahls et al.20 reported Cu-functionalized NU-1000 MOF material can catalyze methane to methanol; Sarawoot Impeng et al.21 found that the C–H bond activation of ethane on the Fe(IV)–O species in a Zn-based cluster of MOFs; Xiao et al.22 have experimentally reported that a high-spin Fe(IV)–O species can be formed in magnesium-diluted Fe-MOF-74 [Fe/(Mg)-MOF-74] and thus affects the hydroxylation of ethane [Nat. Chem., 2014, 6(7), 590]. The weak-field ligand in magnesium-diluted Fe-MOF-74 obviously provide a suitable environment for a high-spin Fe(IV)–O to induce hydrogen abstraction reaction. Note that Xiao et al.23 have also found that expanded Fe-MOF-74 analogues can experimentally catalyze the reaction of cyclohexane to cyclohexanol but there exists also byproduct cyclohexanone, however, the Mg diluted Fe/(Mg)-MOF-74 in which 5% of the Mg atoms are substituted by Fe can hydroxylate ethane to ethanol to avoid overoxidation.22 And using nitrous oxide as an oxidant, Fe/(Mg)-MOF-74 can catalyze 1,4-cyclohexadiene (CHD) oxidation to produce benzene.22 However, no mechanism has been reported so far about the present reaction of catalysis oxidation of such a type of cyclic hydrocarbon, i.e., CHD. Surely, without carrying out accurate quantum chemistry calculations, we cannot understand the reaction mechanism in detail.
Therefore, we decide to study above reaction mechanism of CHD conversion to benzene catalyzed by Fe/(Mg)-MOF-74 with N2O as oxidant using Kohn–Sham density functional theory (DFT).24 Here we firstly propose the mechanism for oxidation of cyclic CHD, which corresponds to two-step hydrogen abstraction reactions with the primary products being benzene. In addition, we also proposed a rebound mechanism of cyclohexane (CHA) hydroxylation to cyclohexanol using Fe(IV)–O species of Fe/(Mg)-MOF-74 which has not been studied experimentally. Most importantly, we have built up the linear relationship between the activation enthalpies and descriptors, such as, the NPA charge of oxyl atoms, C–H lengths in TSs, and bond dissociation energy (BDE) of C–H bonds for different hydrocarbons with a series of M/(Mg)-MOF-74 (M = Cu, Ni, Fe, and Co), respectively, using the density functional theory (DFT).
Fig. 1 The structure of (a) magnesium-diluted Fe-MOF-74; (b) cluster model; and (c) six atoms model (color code: purple = Fe, red = O, light green = Mg, gray = C, and white = H). |
DFT approaches including dispersion-corrected (DFT-D)26–28 is implemented to treat structures with long range dispersion interactions and give good results for geometry and energy of MOF.29–32 All Kohn–Sham calculations were performed with the M06-L33 exchange–correlation functional including D3 dispersion-corrected,27 which performs well for long range electron correlation effects and transition metal chemistry.22,34–36 Since the B3LYP-D3 has also been used to study the Fe–O chemistry recently,37 we did test calculations to see how it performs compared with M06-L-D3 exchange–correlation functional as shown in Table S1.† From comparison with the bond lengths of the experimental structure and calculated ones, M06-L-D3 exchange–correlation functional performed better, i.e., the calculated relative error remains at about 0.25% compared with the experimental results. To reduce computational cost, we used the following strategy: geometry optimizations were adopted with the LANL2TZ(f)38–40 basis set for the valence electrons of Fe atom and the 6-31G(d, p)41 basis set for the other atoms, while effective core potential (LANL2) was employed for Fe. To get the reliable relative energy, the single-point energies were refined at the optimized geometries with same exchange–correlation functional and a def2-TZVP42 basis set for all atoms. Note that the Fe/(Mg)-MOF-74 cluster initially carved from the crystal structures were fully optimized. All transition states (TSs) were identified with one imaginary frequency, and internal reaction coordinate (IRC) calculations were carried out to prove that the TSs connecting the right reactants and products. The DFT calculations were performed with Gaussian 16 program package.43 The natural population analysis (NPA) were carried out based on the natural bond orbital theory with NBO 3.1 module44,45 as incorporated in Gaussian 16 package, while electron localization function (ELF)46 calculations to determine bonding nature were performed with Multiwfn 3.6 package.47
Scheme 1 The proposed mechanism for the oxidation of CHD to benzene with the catalysis of Fe/(Mg)-MOF-74. |
Fig. 2 The plotted enthalpy profiles of the catalytic oxidation of CHD by Fe/(Mg)-MOF-74 on triplet and quintet surfaces. |
It is found that the relative enthalpies of all species that react on the quintet state surface are generally lower than those of on triplet one. Note that although the barrier of 3TS1 on the triplet surface is just 3.4 kcal mol−1, which is 16.0 kcal mol−1 lower than that of 19.4 kcal mol−1 of 5TS1, which indicates that the triplet state is sometimes lower in energy compared with the quintet to bring about two-state reactivity of a hydroxylation reaction. This is consistent with recent computational work on the hydroxylation of methane using Fe(IV)oxo in MOF-74.25 However, in this case, the adsorption enthalpy from A (the bare cluster) to B (the N2O bound to Fe of A) is less than that of the quintet pathway (2.2 kcal mol−1 vs. 8.5 kcal mol−1), and the relative energy of the reactant is generally higher than that of the quintet surface at least 37.5 kcal mol−1. Note that the activation barrier 5TS2 and 5TS3 on quintet surface are just 3.5 and 4.7 kcal mol−1 which are also less than those of 5.2 and 7.9 kcal mol−1 on triplet, respectively. So, the reaction will dominantly proceed on quintet PES and it is important to analyze the relevant structures of the hydrogen abstraction process between Fe/(Mg)-MOF-74 and CHD on quintet PES.
The important structures for the catalytic cycle reaction on quintet PES is shown in Fig. 3. From Scheme 1 and Fig. 3, it is found that the Fe–O–N angle is 117° with Fe–O bond at 2.48 Å for the structure B which is in agreement with powder neutron diffraction data.22 The three steps including three transition-state structures in the catalytic cycle are discussed in sequence as follows: (1) formation of a Fe(IV)–O intermediate via 5TS1: the first transition state about the cleavage of N–O bond of N2O to form the Fe(IV)–O (A–O) species is the rate-determining step with activation enthalpy of 19.4 kcal mol−1. The Fe–O bond distance is shortened from 2.48 Å in B to 1.87 Å in 5TS1, while the C–H bond distance is elongated from 1.20 to 1.49 Å. Then the C species is formed corresponding to an octahedron with Fe(IV)–O and N2. Note that the process from B to C is highly exothermic (ΔH = −15.6 kcal mol−1); (2) the first hydrogen abstraction of CHD by Fe(IV)–O intermediate via 5TS2: when N2 gets away from C the CHD is bound to Fe(IV)–O to form D species. And the first hydrogen abstraction is accomplished via 5TS2 from D to form the Fe(III)–OH with a nearby CHD radical (E). The activation enthalpy of 5TS2 is 3.5 kcal mol−1, which is less than 13.9 kcal mol−1 reported by Xiao et al.48 for oxidation of ethane. The reason may be that the bond dissociation energy (BDE) of CHD is less than that of ethane (75.1 kcal mol−1 vs. 99.8 kcal mol−1 as listed in Table 3). Note that the Fe–O bond distance is elongated from 1.62 to 1.71 Å from D to 5TS2, while the C–H bond is elongated from 1.10 to 1.23 Å; (3) the second hydrogen abstraction of CHD radical by Fe(III)–OH intermediate to form benzene via 5TS3: the E species undergoes conformational change to form F species that is related to the second hydrogen bound to Fe(III)–OH and the Fe(III)–OH abstracts the second hydrogen of CHD radical to benzene (G) via 5TS3. Note that the Fe–O bond distance is elongated from 1.85 to 1.93 Å from F to 5TS3, while the C–H bond distance of CHD radical is elongated from 1.11 to 1.19 Å, indicating that the second hydrogen of CHD radical is close to O of Fe(III)–OH while the O is away from Fe center. This makes the H2O easily separated from Fe(II) and the catalyst recyclable. At last, the reaction ends up with the products benzene, N2 and H2O, which are in agreement with experiment.22
Fig. 3 The selected important information of the most favorable PES (quintet surface) for oxidation of CHD. |
Project | Species | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
A | B | 5TS1 | C | D | 5TS2 | E | F | 5TS3 | G | |
a Average charge on the five MOF oxygen atoms in the first coordination shell of Fe.b O6 is the oxygen atom transferred from N2O to the Fe center.c Average Fe–O distance for the five oxygen atoms (O1–O5) in the first coordination shell of Fe. | ||||||||||
Oxidation state | II | II | II | IV | IV | III | III | III | III | II |
Fe spin population | 3.65 | 3.65 | 3.78 | 2.90 | 3.12 | 3.57 | 3.92 | 3.90 | 3.65 | 3.66 |
O spin population | — | 0.01 | −0.01 | 0.66 | 0.65 | 0.45 | 0.4 | 0.37 | 0.23 | 0.03 |
Fe NPA charge | 1.09 | 1.03 | 1.15 | 0.95 | 0.95 | 1.10 | 1.21 | 1.19 | 1.17 | 1.02 |
O NPA charge | — | −0.34 | −0.54 | −0.38 | −0.38 | −0.64 | −0.91 | −0.93 | −0.95 | −0.94 |
Charge on O1–O5a | −0.82 | −0.83 | −0.80 | −0.78 | −0.78 | −0.80 | −0.81 | −0.81 | −0.81 | −0.83 |
Fe–O6 distanceb | — | 2.48 | 1.87 | 1.62 | 1.62 | 1.71 | 1.84 | 1.85 | 1.93 | 2.28 |
Fe–O distancec | 2.08 | 2.09 | 2.09 | 2.09 | 2.09 | 2.11 | 2.10 | 2.11 | 2.10 | 2.10 |
Wiberg bond index of Fe–O6 | — | 0.11 | 0.51 | 1.42 | 1.43 | 1.00 | 0.63 | 0.61 | 0.45 | 0.17 |
To gain more electronic structure information about hydrogen abstraction reaction, we also calculated the electron localized function (ELF)46,51 and orbital interaction. We mainly focus on the first hydrogen abstraction step, which possesses the features of significant interest.52 In Fig. 4, the color-scale plot of the electron localization function (ELF) of structure 5TS2, B and C is shown to describe the first hydrogen abstraction of CHD by Fe/(Mg)-MOF-74. ELF is a three-dimensional real space function with a numerical range between 0.0 and 1.0. In short, the areas surrounded by the ELF isosurface is of large ELF values, where electrons are more localized and are thus not easy to get out of it. On the contrary, the weak electron localization area is with small ELF value, indicating that the electrons placed there can be easily delved into other areas. Note that the deep blue area with the small ELF value between O and H atoms indicates the electron localization is weak in B. As the reaction goes on, the H atom gets closer to the O atom, the ELF value slowly increases between O and H in 5TS2. Finally, when H atom transfers to O of Fe(IV)–O nearby, the relatively smooth yellow area of the O–H is formed, where the electrons are more localized between H and O atom in C species.
Fig. 4 Color-scale plot of the electron localization function (ELF) of species 5TS2, B and C for the first hydrogen abstraction of CHD catalyzed by Fe/(Mg)-MOF-74. |
Since we cannot clarify the steric difference of the two transition states for oxidation of CHD, the reason for lower enthalpy of 5TS2 compared with 3TS2 requires an explanation. In Scheme 2, it is shown that on both triplet and quintet PESs, the d-area of Fe(IV)–O gets an electron during the hydrogen abstraction step and the σCH bond orbital of the CHD is converted to a singly occupied orbital of the CHD ϕC. In Fig. S2(a)–(d),† from the plotted spin natural orbitals, natural orbitals and their occupations of 5A–O, 3A–O, 5TS2 and 3TS2, respectively, it is found that in 5TS2, there is a spin-down electron in an ϕC orbital which is located on the CHD and five unpaired spin-up electrons in the d-area. Thus, as shown in Scheme 2(a), the d-area of 5TS2 gets an electron in the unoccupied orbital, while the C–H bond orbital is transferred to a singly occupied spin-down orbital of CHD, which can be called as σ-pathway. On the other hand, in 3TS2, there is a singly occupied spin-up electron in the d-area and another one in the CHD, as shown in Scheme 2(b). The gained d electron is now in dyz orbital which is the main component of the orbital, which can be called as π-pathway. Note that the d-area of 5TS2 gets a spin-up electron in the orbital and 3TS2 gets a spin-down in the orbital, and we have known the orbital-selection principle for the spin dependent orientation selectivity,53–56 therefore, the TSs may adjust the direction and tend to maximize the overlap of the two orbitals participating in the electronic transition. In Scheme 2(a), it is shown that a spin-up electron is transferred to the orbital on quintet, and the favorable geometry of 5TS2 is the overlap. As can be seen from the overlap sketch map in Scheme 2(a), a vertical orientation of 5TS2 will maximize this overlap and suggest the preferred pathway of the transition state. By comparison with quintet, a spin-down electron of C–H on triplet PES is transferred to one of the orbitals of Fe–O, and thus, the overlap sketch map shows that the preferred structure of 3TS2 should involve a sideway attack due to the large steric effect and the unavoidably increased Pauli repulsion, resulting in the attack angle to about 120° and less overlap than the quintet one (Scheme 2(b)). It makes a rationalization that the activation enthalpy of 5TS2 is 3.5 kcal mol−1 which is less than that of 3TS2 with 4.7 kcal mol−1 as shown in Fig. 2.
Scheme 2 Electron shifts diagrams during the first step of hydrogen abstraction on (a) quintet; and (b) triplet. |
In Fig. 5, the favorable quintet PES and important structures for the CHA hydroxylation reaction using enthalpy change are plotted. The calculated ΔH and ΔG for various species are collected in Table S2.† Starting from the H species corresponding to the bound of CHA with Fe(IV)–O, crossing the activation enthalpy 8.1 kcal mol−1 of TS4, which is less than that of 13.9 kcal mol−1 reported by Xiao et al. for oxidation of ethane,48 the species I will be formed corresponding to the broken of C–H bond in CHA with exothermicity by 5.2 kcal mol−1. Note that the Fe–O bond distance is elongated from 1.62 to 1.72 Å from H to 5TS4, while the C–H bond is elongated from 1.10 to 1.32 Å at the same time. The next step corresponds to the reaction of the Fe(III)–OH species with the CHA radical to cyclohexanol via 5TS5. Note that the 5TS5 is formed after rotating a hydrogen atom of the hydroxo group along the Fe(III)–O axis in intermediate I and the CHA radical moves closer to the oxyl on Fe(III) center. Note that from I to 5TS5, the Fe–O bond is elongated from 1.83 to 1.86 Å, with the activation enthalpy just 0.9 kcal mol−1. Finally, the CHA radical abstracted the hydroxyl group to form cyclohexanol (J).
The above reaction mechanism can be called as rebound mechanism since a radical produced by hydrogen abstraction rebound to oxygen atom of hydroxyl group to form alcohol. Note that the activation enthalpy of 8.1 kcal mol−1 for oxidation of CHA is larger than that of 3.5 kcal mol−1 of for hydrogen abstraction step in CHD. And different from the oxidation of CHD to benzene via two step hydrogen abstraction, the oxidation of CHA to cyclohexanol corresponds to a rebound mechanism.
Take the hydrogen abstraction of methane as an example, it is shown that the activation enthalpies and NPA charges in reactions catalyzed by M/(Mg)-MOF-74 (M = Cu, Ni, Fe, and Co) are 8.8, 11.7, 15.3, and 19.4 kcal mol−1 as well as −0.682, −0.687, −0.712, and −0.714e, respectively, as listed in Table 2. In Fig. 6(a), it is shown that there exists the linear relationship between the activation enthalpies and the NPA charges of the reacting oxyl atom, i.e., the less of absolute value of NPA charge of oxyl atom, the lower of the activation enthalpy (R2 = 0.810).
Project | Metal center | |||
---|---|---|---|---|
Cu | Ni | Fe | Co | |
Activation enthalpy (kcal mol−1) | 8.8 | 11.7 | 15.3 | 19.4 |
C–H bond length of TS (Å) | 1.299 | 1.307 | 1.351 | 1.442 |
NPA charge of oxyl (e) | −0.682 | −0.687 | −0.712 | −0.714 |
Fig. 6 Plots of activation enthalpy (ΔH) of methane versus (a) NPA charge of oxyl atom; (b) C–H bond length in TSs for M/(Mg)-MOF-74 (M = Cu, Ni, Fe, and Co). Note that regression line is in red. |
In the following, we turn to analyze the C–H bond lengths of transition state structures. It is shown that the C–H bond lengths of methane in TSs on quintet PES catalyzed by M/Mg-MOF-74 (M = Cu, Ni, Fe, and Co) are 1.299, 1.307, 1.351, and 1.442 Å, respectively, as collected in Table 2. Note that the activation enthalpies in reactions catalyzed by M/Mg-MOF-74 (M = Cu, Ni, Fe, and Co) are 8.8, 11.7, 15.3, and 19.4 kcal mol−1, respectively. And there is a linear correlation between C–H bond lengths of TSs and activation enthalpies of the reaction of different M–O clusters as plotted in Fig. 6(b) (R2 = 0.870). This means shorter C–H bonds in TSs is associated with lower activation enthalpy. Furthermore, the activation enthalpies and C–H bond length of TSs for methane, ethane, propane and CHD catalyzed by Fe/(Mg)-MOF-74 are 15.3, 11.4, 10.3, and 3.5 kcal mol−1, as well as 1.351, 1.310, 1.265, and 1.226 Å, respectively, as listed in Table 3. Obviously, both the activation enthalpies and C–H bond lengths decrease as the number of carbon atom increases. It also has a linear correlation between C–H bond lengths of TSs and activation enthalpies of the reaction of different hydrocarbon as plotted in Fig. 7(a) (R2 = 0.863).
Project | Hydrocarbon | |||
---|---|---|---|---|
Methane | Ethane | Propane | CHD | |
Activation enthalpy (kcal mol−1) | 15.3 | 11.4 | 10.3 | 3.5 |
C–H bond length of TS (Å) | 1.351 | 1.310 | 1.265 | 1.226 |
C–H bond (BDE) (kcal mol−1) | 103.9 | 99.8 | 96.9 | 75.1 |
Fig. 7 Plots of activation enthalpy (ΔH) versus (a) C–H bond length in TSs; and (b) C–H bond dissociation energy (BDE) for C1, C2, C3, C6 species. Note that regression line is in red. |
These results are clearly in agreement with the Hammond hypothesis,59 i.e., a more reactive oxyl (less activation enthalpies) should have a transition state (TS) structure closer to reactants. Note that it is clear that for the hydrogen abstraction reaction of methane with catalysis of Cu/(Mg)-MOF-74, there exist the lowest activation enthalpy as plotted in Fig. 6. The reason may result from that Cu/(Mg)-MOF-74 has the largest NPA charge (−0.682e) and shortest C–H bond length (1.299 Å) in TS compared with those of M/(Mg)-MOF-74 (M = Ni, Fe, and Co). We hope the present theoretical result may stimulate the experimental synthesis of Cu/(Mg)-MOF-74 further.
In addition, inspired by the report on dissociation energy of molecular fragments [–OH, –OCH3, and –O(O)CH] on Pt (1,1,1) surfaces by Eric M. Karp et al.,60 we expect different hydrocarbon catalyzed by Fe/(Mg)-MOF-74 to be related to the reactant's C–H bond dissociation energy (BDE). From Table 3, we found that the activation enthalpies and BDE of C–H bond for methane, ethane, propane and CHD (C1, C2, C3, and C6) are 15.3, 11.4, 10.3, and 3.5 kcal mol−1 as well as 103.9, 99.8, 96.9, and 75.1 kcal mol−1, respectively. It is clear that the BDE of C6 is less than that of C1, C2 and C3, so CHD is the easiest one to be activated. From Fig. 7(b), it is shown that there exists linear correlation between the activation enthalpies and bond dissociation energies (BDEs) of C–H bonds in methane, ethane, propane and CHD (C1, C2, C3, and C6) (R2 = 0.917). According to this, we can estimate activation enthalpy without needing to optimize the TS structures.
In addition, for the catalysis by M/(Mg)-MOF-74 (M = Cu, Ni, Fe, and Co) of hydrogen abstraction of methane, it is found that there exists linear relationships between the activation enthalpies, C–H bond lengths of hydrogen-abstraction transition states, and NPA charge of the reacting oxyl atoms, i.e., shorter C–H bond and less absolute value of NPA charge of oxyl atom are associated with lower activation enthalpy. It should be noted that there is the lowest activation enthalpy for the hydrogen abstraction reaction of methane with the catalysis of Cu/(Mg)-MOF-74. With the catalysis of Fe/(Mg)-MOF-74, it is shown that the C–H bond dissociation energies (BDEs) and the C–H bond lengths in TSs of methane, ethane, propane, and CHD have a positive correlation versus their activation enthalpies, i.e., the lower BDE of hydrocarbon and shorter C–H bond length of TS correspond to lower activation enthalpies. The above descriptors should be useful to predict the catalysis activity of enzyme-like MOF in activating C–H bond. We hope the present theoretical modeling may give a guide toward the synthesis of metalloenzyme-like MOFs for the activation of C–H bond.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra04495g |
This journal is © The Royal Society of Chemistry 2019 |