Qin-Qin
Hu
a,
Qi-Fa
Chen
a,
Hong-Tao
Zhang
*a,
Jia-Yi
Chen
*b,
Rong-Zhen
Liao
b and
Ming-Tian
Zhang
*a
aCenter of Basic Molecular Science (CBMS), Department of Chemistry, Tsinghua University, Beijing 100084, China. E-mail: mtzhang@mail.tsinghua.edu.cn; zhanght18@tsinghua.org.cn
bKey Laboratory of Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica, Hubei Key Laboratory of Materials Chemistry and Service Failure, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, China. E-mail: jiayi@hust.edu.cn
First published on 5th December 2024
The one-step oxidation of benzene to phenol represents a significant and promising advancement in modern industries focused on the production of high-value-added chemical products. Nevertheless, challenges persist in achieving sufficient catalytic selectivity and preventing over-oxidation. Inspired by copper enzymes, we present a nonsymmetric dicopper complex ([CuII2(TPMAN)(μ-OH)(H2O)]3+, 1) for the selective oxidation of benzene to phenol. Utilizing H2O2 as the oxidant, complex 1 demonstrates remarkable catalytic activity (a TON of 14000 within 29 hours) and selectivity exceeding 97%, comparable to the finest homogeneous catalyst derived from first-row transition metals. It is noteworthy that the significant substituent effect, alongside a negligible kinetic isotope effect (KIE = 1.05), radical trapping experiments, and an inconsistent standard selectivity test of the ˙OH radicals, all contradict the conventional Fenton mechanism and rebound pathway. Theoretical investigations indicate that the active CuII(μ-O˙)CuII–OH species generated through the cleavage of the O–O bond in the CuII(μ-1,1-OOH)CuI intermediate facilitates the hydroxylation of benzene via an electrophilic attack mechanism. The nonsymmetric coordination geometry is crucial in activating H2O2 and in the process of O–O bond cleavage.
Significant advancements have been made in the development of heterogeneous materials11–19 and homogeneous molecular catalysts20–39 that facilitate this one-step phenol production. Nonetheless, these approaches remain distant from achieving large-scale industrial viability due to inadequate phenol yields and selectivity. To enhance the design of more efficient catalysts, extensive research has focused on elucidating the mechanism underlying aromatic hydroxylation, leading to the proposal of several reaction pathways. Notably, one mechanism suggests direct activation of the C–H bond by hydroxyl (˙OH) or hydroperoxyl (˙OOH) radicals produced via the Fenton-type reaction.40–44 Alternative pathways include an electrophilic aromatic substitution pathway and a rebound pathway mediated by highly energetic metal–oxygen species, such as FeIVO species found in cytochrome-P450.45–50 Systematic investigations into molecular catalysts have revealed that the ligand denticity plays a crucial role in modulating the reaction activity and mechanism.51,52 For instance, mononuclear iron complexes, featuring tridentate N3 ligands, preferentially follow a hydroxyl radical pathway due to the instability of the iron–oxygen intermediate, thereby resulting in lower selectivity.53 In contrast, tetradentate N4 coordinated iron catalysts facilitate heterolytic cleavage of the O–O bond in the FeIII–OOH species, aided by a free coordination site, ultimately yielding more reactive FeV
O intermediates that enhance both yield and selectivity.54 However, pentadentate N5 ligands in FeIII–OOH intermediates tend to favor homolytic cleavage of the O–O bond, leading to the generation of FeIV
O and ˙OH, which activates benzene for hydroxylation.55
Inspired by copper-containing metalloenzymes, such as tyrosinase and catechol oxidase,56–61 which perform aromatic ring oxidation, copper-based catalysts have also attracted considerable attention due to their superior catalytic activity compared to other first-row transition metals, particularly in the cases of dicopper catalysts.21–25,31,32 Notably, the Cu2(6-hpa)(μ-OH) catalyst, with 6-hpa representing 1,2-bis[2-[bis(2-pyridylmethyl)aminomethyl]-6-pyridyl]ethane, exhibits a remarkable turnover number (TON) exceeding 12000 for phenol production,23 attributed to bimetallic cooperation catalysis via a μ–η1:η1-O2 type dicopper peroxo intermediate. Research in Cu–oxygen chemistry has shown that ligand denticity can significantly influence the coordination geometries of key Cu–O2 adducts and subsequently alter their reactivity, particularly regarding O–O bond cleavage.61–66 In recent years, a diverse range of dicopper–oxygen adducts, including Cu(μ–η1:η1-O2)Cu, Cu(μ–η2:η2-O2)Cu, and Cu(μ-oxo)2Cu intermediates, have been explored for the hydroxylation reaction of aromatic rings (Fig. 1a and b).57,67–70 Meanwhile, the Cu(μ–η1:η2-O2)Cu type peroxo intermediate remains largely uncharacterized in the context of selective benzene hydroxylation.
Recently, we developed a dicopper water oxidation catalyst [CuII2(TPMAN)(μ-OH)(H2O)]3+(1, Fig. 1c), featuring a nonsymmetric ligand TPMAN (= 1-(7-((bis(pyridin-2-ylmethyl)amino)methyl)-1,8-naphthyridin-2-yl)-N-methyl-N-(pyridin-2-ylmethyl)methanamine), in which two CuII ions respectively possess different free coordination sites. This structure leads to the formation of a μ–η1:η2-O2 type dicopperperoxo intermediate, allowing for enhanced catalytic performance.65 In this study, we employ this structurally well-defined dicopper complex for the direct hydroxylation of benzene with H2O2 as the oxidant, demonstrating exceptional catalytic activity (a TON of 14000 within 29 hours) and selectivity (up to 97%) for phenol production. Mechanistic investigations reveal that the Cu(μ–η1:η2-O2)Cu intermediate can be effectively reduced by H2O2, directly generating the CuII(μ-1,1-OOH)CuI intermediate, which facilitates O–O bond cleavage. The resulting CuII(μ-O˙)CuII–OH species serves as the active entity for the electrophilic hydroxylation of benzene. Compared to previously reported dicopper catalysts,23,57,67–70 our uniquely coordinated dicopper site alters the activation pathway of H2O2, yielding a distinct copper-oxyl intermediate that enhances the selectivity of benzene hydroxylation.
Entry | T (°C) | Conv.a (%) | Yieldb (%) | Selectivity for phenolc (%) |
---|---|---|---|---|
Conditions: 11 mmol of benzene, 44 mmol of 30% H2O2, 0.44 μmol of catalyst 1, and 5 μmol of Et3N in 6.3 mL of CH3CN.a Conversation% = n(consumed benzene)/n(starting benzene).b Yield% = n(phenol)/n(starting benzene).c Selectivity is calculated based on the products detected by GC. Selectivity% = n(phenol)/(n(phenol) + n(benzoquinone)). The reaction time is 20 h.d 24 h.e 12 h. | ||||
1 | 45 | 25 | 3.4 | 78 |
2 | 65 | 61 | 6.7 | 92 |
3 | 70 | 52 | 7.3d | 92 |
4 | 80 | 61 | 7.0e | 92 |
Subsequent investigations focused on the long-term catalytic performance of benzene hydroxylation using catalyst 1, with the time courses of phenol production illustrated in Fig. 2 and S5.† At a temperature of 65 °C, catalyst 1 exhibits a maximum turnover frequency (TOFmax) of approximately 859 h−1 and a maximum turnover number (TON) of 14000 after 29 hours, resulting in a phenol yield of 6.1% and a remarkable selectivity of 97%. This impressive catalytic performance is comparable to the highest TON for benzene hydroxylation using H2O2 in homogeneous catalysis (Table S4†).23 At a reduced reaction temperature of 45 °C, catalyst 1 still displayed a noteworthy TON of 8500 for phenol (Fig. S5†). It is important to note that p-benzoquinone was also generated during the early stages of the reaction alongside phenol, indicating that p-benzoquinone is a product of phenol over-oxidation.33 However, the yield of p-benzoquinone remained consistently low throughout the reaction period (Fig. S3 and S4†). These findings suggest that this nonsymmetric dicopper catalyst is capable of selectively promoting the hydroxylation of benzene to phenol, achieving the highest TON reported to date.
![]() | ||
Fig. 2 Time profile formation of phenol in the hydroxylation of benzene using catalyst 1 (0.05 μmol) with 11 mmol of benzene, 44 mmol of 30% H2O2, and 5 μmol of Et3N in 6.3 mL of CH3CN at 65 °C. |
To further explore the reactivity of the active species, we examined various aromatic substrates, including phenol, toluene, and nitrobenzene, under identical reaction conditions (Table 2 and Fig. S6–S8†). For phenol, p-benzoquinone was the only product detected initially, while catechol was observed after 12 h, resulting in a total TON of 3200 after 23 hours. Using toluene as the substrate yielded o-cresol, p-cresol, and benzaldehyde as the major products, with a total TON of 1700 after 23 hours. The ratio of aromatic to aliphatic oxidations was approximately 2.7:
1.0, and the distribution of cresols was about 1.6
:
1.0
:
0 (o
:
p
:
m), which is in contrast with the ratios typically associated with the Fenton mechanism.71 In the case of nitrobenzene, neither o-nitrophenol nor p-nitrophenol was detected. Consequently, the relative reactivity order based on the TON under the same reaction conditions was determined as follows: phenol > toluene > benzene ≫ nitrobenzene. This trend highlights the propensity for oxidation to occur preferentially on aromatic rings with electron-donating substituents, and the observed sensitivity to the electronic density of the aromatic ring suggests that the oxidation mechanism is likely metal-mediated through an electrophilic intermediate rather than involving radical pathways.
Entrya | Substrate | Products [%] (o/p/m) | TOFc (h−1) | TONd |
---|---|---|---|---|
a 0.44 μmol catalyst, 5 μmol base (Et3N), 44 mmol H2O2 (30%), 6.3 mL CH3CN, 45 °C, and 23 h. The substrate: 1 and 2: 11 mmol benzene; 3: 8.5 mmol phenol; 4: 8.2 mmol toluene; and 5: 8.2 mmol nitrobenzene. b 0.05 μmol catalyst, 5 μmol base (Et3N), 44 mmol H2O2 (30%), 6.3 mL CH3CN, 65 °C, and 23 h. c TOF of all products after 23 h. d TON of all products after 23 h. | ||||
1 | Benzene | Phenol [83] | 196 | 1200 |
p-Benzoquinone [17] | ||||
2b | Benzene | Phenol [97] | 859 | 14![]() |
p-Benzoquinone [3] | ||||
3 | Phenol | Catechol [40] | 332 | 3200 |
p-Benzoquinone [60] | ||||
4 | Toluene | Cresol [73] | 245 | 1700 |
(60![]() ![]() ![]() ![]() |
||||
Benzaldehyde [27] | ||||
5 | Nitrobenzene | No product | — | — |
To investigate the reaction mechanism in depth, 5,5′-dimethyl-1-pyrroline N-oxide (DMPO) was employed as a radical scavenger to assess the formation of hydroxyl (˙OH) or hydroperoxyl (˙OOH) radicals, which are potential active species in the Fenton mechanism. As demonstrated in Fig. 3a, the rates of phenol production in the presence of either 0.5 mM or 5.0 mM DMPO were consistent with the rate observed in the absence of DMPO. Additionally, the introduction of carbon tetrachloride (CCl4) did not yield detectable chlorobenzene or biphenyl species at any point during the reaction. These results indicate that no obvious radical species, including hydroxyl (˙OH) or hydroperoxyl (˙OOH) radicals, exist in the reaction solution to catalyze hydroxylation of benzene through the radical-chain mechanism. Furthermore, the kinetic isotope effect (KIE) was investigated using hexadeuterio-benzene (d6-benzene) as a substrate, resulting in a KIE value of 1.05 determined from the initial rates of phenol production (Fig. 3b). This value suggests that the C–H bond cleavage of benzene is not part of the rate-determining step. This relatively low KIE value differs significantly from the KIE (approximately 1.7–1.8) reported for Fenton-type hydroxylation72 or those existing values greater than 4 for rebound mechanisms.42,73,74 Additionally, a standard selectivity test for the reaction involving the ˙OH radical was conducted using methylcyclohexane under standard conditions.75 The gas chromatography (GC) results reveal that the product ratio of tertiary alcohol to secondary alcohol is approximately 1:
15 (Fig. S9†), which does not align with the expected selectivity ratio of 3
:
10 (ref. 75) for the ˙OH radical mechanism. Furthermore, the lack of hydroxylation reactivity observed with nitrobenzene is inconsistent with the Fenton-type mechanism.72 These findings seem to contradict the Fenton-type mechanism and the rebound process, while they are consistent with an electrophilic aromatic substitution mechanism.51,68,76
Moreover, the reaction kinetics of benzene oxidation catalyzed by 1 was analyzed through the initial rates of phenol production. As depicted in Fig. S10–S12,† the concentrations of benzene, catalyst, and H2O2 each exhibit a first-order linear relationship with the rate of phenol production. Based on these dependencies, the kinetic equation governing the benzene oxidation reaction can be articulated as:
d[Phenol]/dt = kcat[catalyst 1][benzene][H2O2], |
Based on the experimental results outlined above, computational studies were performed to elucidate the details of the catalytic mechanism. The initial [CuII2(TPMAN)(μ-OH)(H2O)] species (labeled as complex 1, Fig. 4a, for detailed structure, see Table S5†), featuring two penta-coordinated copper sites, can react with H2O2 to form the CuII(μ–η1:η1-OOH)CuII intermediate Int1 (Table S5†), accompanied by an energy release of 3.8 kcal mol−1. Following this, the deprotonation of Int1 by Et3N occurs, and the generated dicopperperoxo Int2 ([CuII(μ–η1:η2-O2)CuII], Table S6) is then reduced by 1/2 equivalent of H2O2, leading to the formation of Int3 ([CuII(μ-1,1-OOH)CuI], Table S6†), which has been confirmed to be a vital intermediate in our previous study,64,77 and 1/2 equivalent of O2. The deprotonation process is calculated to be exergonic by 10.0 kcal mol−1, while the reduction of Int2 to Int3 is endergonic by +1.8 kcal mol−1; thus, Int2 has been assigned as the starting point of the energy profile (Fig. 4b). Additionally, the possibility of a direct electrophilic attack of benzene on the peroxide moiety of Int2 can be ruled out, as the calculated energy barrier is +40.7 kcal mol−1 (Fig. S13, and Table S7†).
Subsequently, Int3 undergoes a homolytic O–OH bond cleavage process viaTS1, generating the diradical intermediate Int4′ (Table S8†) with a barrier of +7.6 kcal mol−1. The optimized structure of the transition state TS1 is depicted in Fig. 5. TS1 is a doublet and the Mulliken spin populations on Cu1, O1, Cu2, and O2 are 0.62, 0.42, 0.53, and −0.89, respectively. Frequency analysis of 2TS1 gives only one imaginary frequency of 295.5i cm−1, related to the cleavage of the O1–O2 bond. The length of the breaking O1–O2 bond and the forming hydrogen bond of O1–H1 measures 2.38 Å and 1.91 Å, respectively. The Ab Initio Molecular Dynamics (AIMD) calculations for Int4′ indicate that Int4′ can be regarded as a metastable intermediate since it can rapidly (within one ps, Fig. S14 and S15†) isomerize to the CuII(μ-O˙)CuII–OH intermediate Int4 without any applied potential. Int4 lies at −14.9 kcal mol−1 relative to Int3, with a ground state of the doublet. In 2Int4 (Fig. 5), the distance of the formed Cu2–O2 bond is 1.80 Å, consisting of a β-electron donated by the O2 moiety and an α-electron from Cu2. It should be pointed out that the decrease in Mulliken spin population on O1 (a Mulliken spin population of 0.30) and Cu2 (a Mulliken spin population of −0.11) atoms is due to the interaction of the remaining β-electron on Cu2 with the O1 radical, which possesses an α-electron. Additionally, the formation of the CuII(μ-OH)CuII–O˙ intermediate Int4′′ (Table S9) has also been considered, with a barrier of +23.5 kcal mol−1 (TS1′, Fig. 4b, Table S9†) relative to Int3 and +15.9 kcal mol−1 higher than that of TS1. Compared to Int4′′, the formation of Int4 is kinetically more favorable; therefore, the pathway initiating from Int4 is described in detail here, and other pathways will be briefly discussed when necessary.
The O1 radical in Int4 then electrophilically attacks the benzene substrate directly, forming the essential O1–C1 bond through a quartet transition state TS2, with a significant Mulliken spin population distributed at the O1 (a Mulliken spin population of 0.69) and benzene (a Mulliken spin population of 0.53). The barrier of TS2 is calculated to be +14.0 kcal mol−1 relative to Int4 plus benzene, and the length of the O1–C1 bond in TS2 is 1.92 Å. The generated 4Int5 (Table S10†) is +5.2 kcal mol−1 higher in energy than 2Int4, which contains a benzene radical (a Mulliken spin population of 0.92) and two CuII ion centers. Notably, when nitrobenzene acts as a substrate, the O1–C1 bond formation via the electrophilic attack of Int4 needs to overcome higher barriers of about +15.4–16.7 kcal mol−1 (Tables S11 and S12†), consistent with the above-mentioned experimental findings. If CH3CN as the substrate is activated by Int4 abstracting the hydrogen atom, the energy barrier is also obviously higher than that of the benzene activation (Tables S11 and S12†). Thus, the CH3CN solvent would not interfere with the benzene hydroxylation reaction, and no CH3CN-related products are observed in the gas chromatograph (Fig. S8†). These results further confirm that Int4 activates the benzene via electrophilic attack rather than the abstraction of hydrogen atoms.
The ground state of transition state TS3 (Fig. 5) is a doublet, showing that the electron on the benzene radical is flipped (a Mulliken spin population of −0.43) and coupled with Cu1II (a Mulliken spin population of 0.32) antiferromagnetically. TS3 involves proton transfer from C1 to the OH moiety and electron transfer from the benzene radical to the Cu1II ion. The C1–H2, H2–O2, and Cu2–O1 bond distances in 2TS3 are 1.16, 1.79, and 1.98 Å, respectively; the barrier of TS3 is +6.3 kcal mol−1, and the imaginary frequency is 198.8i cm−1. Downhill from TS3, Int6 (Table S10†) is formed, and it is associated with a large energy release of 49.7 kcal mol−1 relative to Int5. Based on the Mulliken spin population analysis of Int6, Cu1 accepts the β-electron transferred from the benzene radical, resulting in its oxidation state becoming +1; Cu2 remains at +2 and coordinated with the generated H2O. Finally, the oxidant H2O2 reacts with Int6, yielding H2O molecules and the target product phenol and regenerating species Int3 to catalyze the next cycle.
According to the Gibbs free energy diagram derived from the DFT calculations, the electrophilic attack of the catalytically active species Int4 on the benzene (TS2) is the rate-determining step, with a calculated barrier of +14.0 kcal mol−1. This finding aligns with the experimental evidence indicating that the C–H bond cleavage of benzene was not involved in the rate-determining step. The additive base (Et3N) is solely responsible for the deprotonation of Int1 to form Int2, initiating the catalytic process via the generated active intermediate Int3. Moreover, the generation of catalytically active species Int4 (CuII(μ-O˙)CuII–OH) benefits from the nonsymmetric coordination geometry of Int3, which provides vacant coordination sites to receive the OH group.
![]() | (1) |
![]() | (2) |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4dt02872d |
This journal is © The Royal Society of Chemistry 2025 |