Subnanometric Cu clusters on atomically Fe-doped MoO2 for furfural upgrading to aviation biofuels

Single cluster catalysts (SCCs) are considered as versatile boosters in heterogeneous catalysis due to their modifiable single cluster sites and supports. In this work, we report subnanometric Cu clusters dispersed on Fe-doped MoO2 support for biomass-derived furfural upgrading. Systematical characterizations suggest uniform Cu clusters (composing four Cu atoms in average) are homogeneously immobilized on the atomically Fe-doped ultrafine MoO2 nanocrystals (Cu4/Fe0.3Mo0.7O2@C). The atomic doping of Fe into MoO2 leads to significantly modified electronic structure and consequently charge redistribution inside the supported Cu clusters. The as-prepared Cu4/Fe0.3Mo0.7O2@C shows superior catalytic performance in the oxidative coupling of furfural with C3~C10 primary/secondary alcohols to produce C8~C15 aldehydes/ketones (aviation biofuel intermediates), outperforming the conventionally prepared counterparts. DFT calculations and control experiments are further carried out to interpret the structural and compositional merits of Cu4/Fe0.3Mo0.7O2@C in the oxidative coupling reaction, and elucidate the reaction pathway and related intermediates.

D ownsizing metal active sites into several atoms scale is demonstrated effective in boosting their catalytic performances [1][2][3] . Compared with conventional metal nanoparticles/bulks, atomic site catalysts combine the advantages of both homogeneous and heterogeneous ones, such as extremely high atom utilization, well-defined active sites, and reliable durability, thus showing remarkably enhanced catalytic activities [4][5][6] . Consequently, single cluster catalysts (SCCs) are emerging, which compose of a certain number of metal atoms as active sites (M x , where M and x represent metal and atom number, respectively) 7,8 . Apart from the above-mentioned common merits of atomic site catalysts, the tunable composition and geometry of M x may endow SCCs with other distinctive properties. Thus, a rational decoration of M x could correspondingly modify its charge distributions and consequently affect the catalytic processes 9,10 .
Recently, numerous efforts have been devoted to the synthesis of novel SCCs with unique electronic structures. For instance, tailoring the metal type and atom number of M x could effectively modify its geometrical topology. The distinctively spatial distribution of the metal atoms in M x leads to unexpected size effect and polarized charge distribution, achieving remarkable catalytic performances in various reactions [11][12][13][14] . This protocol is highly dependent on the physicochemical natures of the metals employed (e.g., coordination number, chemical state, and geometry topology), which could only qualitatively change the electronic structure of SCCs. Alternatively, decorating the support of SCCs may modify the electronic structure of M x in a more precise and selective way, i.e., altering the metal-support interactions [15][16][17][18][19][20] . In theory, the even introduction of different dopants into SCC supports could generally and uniformly modulate the electronic property of M x , which, however, is still scarcely investigated in spite of the significance in the pursuing of highly efficient single cluster catalysts for advanced catalytic applications.
In this work, subnanometric Cu clusters immobilized on Fedoped MoO 2 nanocrystals are fabricated via a cation exchangediffusion strategy and subsequent pyrolysis procedure. The obtained nanocomposite features a multi-shelled hollow hierarchical porous octahedron (SHHPO) morphology. Spherical aberration-corrected high-angle annular dark-field scanning transmission electron microscopy (AC HAADF-STEM) and X-ray absorption spectroscopy (XAS) results uncover the homogeneous immobilization of Cu clusters composing four Cu atoms in average (Cu 4 ) upon the atomically Fe-doped ultrafine MoO 2 nanocrystals in the obtained Cu 4 / Fe 0.3 Mo 0.7 O 2 @C. The atomic doping of Fe into MoO 2 leads to significantly modified electronic structure and consequently charge redistribution inside the Cu clusters. Moreover, the asprepared Cu 4 /Fe 0.3 Mo 0.7 O 2 @C shows superior catalytic performance in the oxidative coupling of furfural (FFA) with C 3~C10 primary/secondary alcohols (e.g., isopropanol) to produce the corresponding C 8~C15 aldehydes/ketones (e.g., 4-(2furyl)-buten-2-one, FBO) as aviation biofuel intermediates. The structural and compositional relationship of Cu 4 / Fe 0.3 Mo 0.7 O 2 @C in the oxidative coupling reaction and the possible reaction pathway are also investigated.

Results
Synthesis of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO. The synthesis of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO is illustrated in Fig. 1. First, NENU-5 is assembled by using a reported method 21 . Powder X-ray diffraction (XRD) patterns ( Supplementary Fig. 1a) and scanning electron microscopy (SEM) images ( Supplementary  Fig. 1c-e) reveal the successful synthesis of monodispersed NENU-5 octahedrons in size of 1 µm with smooth external surfaces. N 2 adsorption-desorption isotherms indicate the microporous nature of NENU-5 ( Supplementary Fig. 2b, c). Subsequently, the obtained NENU-5 is immerged into a Fe 3+contained aqueous solution to substitute partial Cu ions with Fe species to yield Fe-NENU-5. Meanwhile, the Fe 3+ -contained aqueous solution could also facilitate the structural evolution of the solid NENU-5 octahedrons through cation exchanging to fabricate multi-shelled nanostructures. The XRD patterns of Fe-NENU-5 match well with those of the parent NENU-5 (Supplementary Fig. 2a), suggesting the well-retained crystalline phase of NENU-5 upon Fe 3+ substitution. The as-prepared Fe-NENU-5 features typical type IV N 2 adsorption-desorption isotherms with an obvious N 2 -uptake increment at high relative pressures (0.9-1.0) ( Supplementary Fig. 2b), revealing the generation of meso-and macropores. Correspondingly, the formation of hierarchical pores in Fe-NENU-5 is also confirmed by the significantly reduced specific surface area (from 615.2 to 399.3 m 2 g −1 , Supplementary Table 1) and enlarged pore-size distribution range ( Supplementary Fig. 2c) in comparison with NENU-5.
Fe-NENU-5 retains octahedral morphology of NENU-5 and uniform particle size (Fig. 2a, b) whereas possesses rougher external surfaces. Transmission electron microscopy (TEM) image ( Fig. 2c) clearly disclose the triple-shelled hollow nanostructure of Fe-NENU-5 with even shell thickness (60-70 nm), intershell spacings (70 nm) and hollow cavities (250 nm). Further, HAADF-STEM and the corresponding elemental mapping images (Fig. 2d, e) show the homogeneous distribution of Cu, Mo, Fe, C, O, and P elements throughout the Fe-NENU-5. Notably, no obvious Cu or Fe aggregations are observed, excluding their residue or adsorption within the nanostructure. These results indicate that the cation exchangediffusion process leads to a partial substitution of Fe to Cu, and simultaneously an evolution of the solid NENU-5 octahedron into multi-shelled hollow hierarchical structure with the generation of meso-and macropores inside.
The morphology evolution of NENU-5 nanocrystals is tracked and characterized by SEM and TEM at different periods of the cation exchange-diffusion process ( Supplementary Fig. 3a-e). At the very beginning, the Fe 3+ cations adsorb on the external surfaces of NENU-5 and in situ substitute some Cu nodes, leading to partial decomposition of NENU-5 nanostructure and the resultant rough external surfaces. With time goes by, the Fe 3+ cations are immerged into the bulk NENU-5 nanocrystals for continuous substitution, resulting in the structural evolution of the solid NENU-5 octahedrons into multi-shelled hollow hierarchical ones within 30 min ( Supplementary Fig. 3a-c). With a further increase in etching time (i.e., 50 and 60 min), doubleand even single-shelled hollow nanostructures ( Supplementary  Fig. 3d, e) are formed.
The content of Fe in the discarded solution decreases while that of Cu correspondingly increases in the cation exchange-diffusion process ( Supplementary Fig. 4), confirming the partial replacement of Cu to Fe in the material. Besides, the weakened intensity of the NENU-5 diffraction peaks in the XRD patterns of the asprepared multi-shelled nanostructures ( Supplementary Fig. 5a) suggest the partially preserved NENU-5 topologies. In addition, the effects of Fe 3+ concentration on the morphology of the obtained Fe-NENU-5 are also investigated. A series of cationexchanging solutions containing 0.05-0.30 g FeCl 3 • 6H 2 O are prepared for the substitution. Characteristic diffraction peaks of NENU-5 are clearly observed in the XRD patterns of the obtained counterparts ( Supplementary Fig. 5b), indicating the partially retained NENU-5 crystalline structure. When the FeCl 3 • 6H 2 O dosage increases from 0 to 0.30 g, the external surfaces of the asprepared Fe-NENU-5 gradually turn from smooth to rough (Supplementary Fig. 3f-j). Specifically, a low FeCl 3 • 6H 2 O loading (<0.05 g) leads to smooth external surfaces and solid yolk-shell nanostructures of Fe-NENU-5 ( Supplementary Fig. 3f, j). Rough external surfaces and triple-shelled hollow nanostructures with irregular intershell gaps are observed at the FeCl 3 • 6H 2 O usage of 0.10 g (Supplementary Fig. 3h). When the FeCl 3 • 6H 2 O loading is higher than 0.20 g, the triple-shelled nanostructure collapses with apparent generation of defects ( Supplementary Fig. 3i, j). These results suggest that a 40 min period of cation exchanging with 0.15 g FeCl 3 • 6H 2 O loading is appropriate for the formation of hollow triple-shelled nanostructure of Fe-NENU-5 ( Fig. 2a-d). An annealing treatment is then carried out under Ar atmosphere to convert Fe-NENU-5 into Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO. The pyrolysis temperature is set at 600°C, higher than the decomposition temperatures of NENU-5 and Fe-NENU-5 determined by thermogravimetric analysis (Supplementary Fig. 6).
Structural characterizations of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO. The as-prepared Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO inherits the octahedral morphology from Fe-NENU-5 but exhibits a shrunken size of ca. 550 nm (Fig. 2f-h). Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO also features type-IV N 2 -sorption isotherms while with a reduced N 2 uptake and an enlarged hysteresis loop ( Supplementary Fig. 1e) as compared with Fe-NENU-5. Besides, the pore-size distribution curves reveal a widened pore-size distribution in 10-100 nm ( Supplementary Fig. 2f), suggesting a large amount of micropores are coalesced into meso-and even macropores during the pyrolysis. Correspondingly, the Brunauer-Emmett-Teller (BET) specific surface area of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C SHHPO is reduced to 169.7 m 2 g −1 with a pore volume of 0.41 cm 3 g −1 (Supplementary Table 1). In comparison, Cu x /MoO 2 @C prepared via the direct pyrolysis of NENU-5 has a slightly larger BET specific surface area (182.9 m 2 g −1 ) with smaller pore volume (0.24 cm 3 g −1 ) due to the more micropores and fewer mesopores.
The multi-shelled hollow hierarchical nanostructure of Cu 4 / Fe 0.3 Mo 0.7 O 2 @C SHHPO is further verified by TEM (Fig. 2h), showing easily identified triple-shelled nanostructures ( Fig. 2i and m) without obvious agglomeration ( Fig. 2j1-j6). AC HAADF-STEM images of the ultrafine Mo oxide nanocrystals (Fig. 2k, l) clearly demonstrate the atomic distribution of Fe and Cu species (isolated bright dots) into the MoO 2 (−1 1 1) lattice planes due to their Z-contrast differences. The support is denoted as Fe 0.3 Mo 0.7 O 2 after Fe incorporation in which the atom ratio of Fe to Mo is determined by ICP-OES (Supplementary Table 1 Fig. 2d). Besides, no characteristic XRD peaks of Cu composites are observed, which is possibly related to their high dispersion and relatively small size. In contrast, the Cu x /MoO 2 @C counterpart shows octahedral morphology with obvious metal agglomerations ( Supplementary Fig. 7), indicating the proposed cation exchange-diffusion strategy is essential for the achievement of both multi-shelled hollow hierarchical nanostructure and high metal dispersion.
X-ray absorption fine structure (XAFS) measurements were performed to investigate the coordination environments of Fe and Cu at atomic level. The Fe K-edge X-ray absorption nearedge structure (XANES) spectra of Cu 4 Table 2). Moreover, no apparent peak at ca.   Figs. 9a, 10a), and the corresponding O 1s peak is also shifted to higher binding energies ( Supplementary Fig. 10b). To exclude the impact of crystal terminations on the XPS spectra, DFT calculations were additionally performed to investigate the electrostatic potential at the surface of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C and Cu x /MoO 2 @C models ( Supplementary Fig. 11). In comparison with Cu x /MoO 2 @C, Cu 4 / Fe 0.3 Mo 0.7 O 2 @C shows significant peak shifts in electrostatic potential and increased work function (5.32 eV vs 4.96 eV). Combining the XPS and DFT results, the observed rigid shift in binding energies is possibly attributed to the electron charge redistributions within the nanocrystals and discrepant surfaces due to the introduction of Fe [24][25][26][27][28] Table 3). The Cu K-edge XANES spectrum of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C is located close to that of Cu foil (Fig. 3i), revealing the metallic state of Cu species 31,32 . A prominent peak at ca. 2.2 Å in the FT k 3 -weighted EXAFS spectrum of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C contributes similarly to Cu foil (Fig. 3j), implying the formation of Cu-Cu bonds, in good agreement with the fitting EXAFS results of R-space spectrum of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C (Fig. 3k) [33][34][35][36] . Besides, the intensity maximum at 7 Å −1 in the WT contour plots of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C (Fig. 3p) is similar to the Cu foil reference, also suggesting the existence of Cu-Cu bonds. Whereas, the dramatic attenuation of post-edge oscillations (Fig. 3i) suggests the Cu species in Cu 4 / Fe 0.3 Mo 0.7 O 2 @C is lacking of long-range orders 33 . In addition, the sharply reduced first-shell intensity in Cu 4 /Fe 0.3 Mo 0.7 O 2 @C (Fig. 3j) indicates a much lower coordination number of Cu in Cu 4 /Fe 0.3 Mo 0.7 O 2 @C than in Cu foil (6.3 to 12), in consistent with the EXAFS fitting data (Supplementary Table 3). These results suggest the Cu species in Cu 4 /Fe 0.3 Mo 0.7 O 2 @C are present in the form of subnanometric clusters with six Cu-Cu bonds, i.e., Cu 4 , of which the topological nanostructure is illustrated in the inset of Fig. 3k. In comparison, the Cu K-edge XANES spectrum of Cu x /MoO 2 @C is located close to that of Cu foil with remarkable post-edge oscillations (Fig. 3i) Supplementary Fig. 14), the prominent band locating at 932.5 eV indicates the existence of Cu 0 species in Cu x /MoO 2 @C (the co-existent Cu 2+ may be related to the surface oxidation). In comparison, Cu 0 species in Cu 4 /Fe 0.3 Mo 0.7 O 2 @C manifest a slight shift towards higher binding energy by ca. 0.5 eV, possibly caused by the electron transfer from Cu to Fe atoms.
To sum up, the cation exchange-diffusion strategy demonstrated here is effective and precise for the fabrication of Cu 4 / Fe 0. 3 [37][38][39][40] . One of the most promising strategies for FFA valorization is the aldol condensation of FFA with aldehydes/ketones (e.g., acetone) for the production of α,β-unsaturated aldehydes/ketones, which are well-known as critical intermediates for the production of aviation biofuels (C 8~C15 alkanes) via a hydrodeoxygenation process [41][42][43][44] . From a green, economic, and sustainable viewpoint, the direct reaction of saturated monobasic alcohols with FFA may be the most preferred and promising, due to the availability and low cost of alcohols in comparison with their oxidation products such as aldehydes and ketones.
In this work, we developed a new reaction route for the production of α,β-unsaturated ketones from FFA, i.e., the oxidative coupling of FFA with secondary alcohols. The onepot cascade oxidative coupling of FFA with isopropanol (IPA) into 4-(2-furyl)-buten-2-one (FBO, Supplementary Fig. 15) is selected as a model reaction. Currently, FBO is produced from the normal aldol condensation of FFA with acetone that is known as a toxic and dangerous solvent. Typically, the reaction is carried out at 120°C under O 2 using IPA as both the reactant and solvent, and the results are summarized in Fig. 5a and Supplementary Table 4. The FFA conversion is below the detecting limit in the absence of catalyst (Fig. 5a

Discussion
Reaction mechanism investigation. On the basis of the intermediates detected in the reaction, the plausible reaction routes are illustrated in Fig. 5c. Initially, FFA is transfer-hydrogenated by IPA into furfuryl alcohol (FFAL) with the generation of equivalent amount of acetone (route I). The obtained FFAL subsequently reacts with another IPA (greatly excessive) molecule through an oxidative-coupling process to produce FBOL (route II), which is finally oxidized into FBO. Another possible reaction route is the conventional aldol condensation of FFA with acetone (route III) 45 , which is also investigated as a potential route in the following reaction mechanism. Some control experiments are carried out to verify the proposed reaction routes (Fig. 5b and Supplementary Table 5). When acetone is employed as the reactant, <15% FBO is obtained (Fig. 5b, columns 1 and 2). Similarly, when the reaction intermediate FFAL is used as the reactant, only 31.9% FBO yield is obtained (Fig. 5b, column 3). These results indicate the FFA transformation over Cu 4 /Fe 0.3 Mo 0.7 O 2 @C is started with the transferring hydrogenation of FFA with IPA, and FFAL could hardly react with acetone under the investigated conditions. In addition, when the reaction atmosphere is changed from O 2 to N 2 , only 5% FBO yield is detected (Fig. 5b, column 4), revealing that the oxidative coupling step (route II) is the dominant route to produce FBO over the Cu 4 /Fe 0.3 Mo 0.7 O 2 @C using FFA and IPA as the substrates.
DFT calculations were performed to deeply elucidate and bring theoretical insight to the reaction mechanism (Figs. 4, 6, and Supplementary Figs. 20, 21). The differential charge densities of Cu 4 /MoO 2 (Fig. 4c, g) and Cu 4 /Fe 0.3 Mo 0.7 O 2 (Fig. 4d, h)  A plausible reaction mechanism is proposed for the oxidative coupling of FFA with IPA to FBO over Cu 4 /Fe 0.3 Mo 0.7 O 2 @C (Fig. 6, for details, see Supplementary Figs. 20, 21). The reaction pathway-1 consists of ten elementary steps. At the very beginning, one FFA molecule bonds with the Cu 0 atom at the top of Cu 4 to form (c1) due to a relatively low binding energy (−0.65 eV). Subsequently, the adsorbed FFA molecule reacts with one adsorbed IPA molecule (c2), overcoming the energy barrier of 1.47 eV (TS-1) and forming FFAL radical and one acetone molecule (illustrated as c3). Afterwards, the generated FFAL radical (c4) further reacts with another IPA molecule to form a (c5) complex, which undergoes a dehydration process with an energy barrier of 0.92 eV (TS-2(p1)) to produce *ph-CHCH-CH(OH)-CH 3 (c6-c7, the ph represents the furan ring). Afterwards, the *ph-CHCH-CH(OH)-CH 3 undergoes a dehydrogenation reaction which tends to take place on the carbon side (to form c9) than at oxygen side (to form c8) due to the former's relatively low energy barrier (0.82 eV in comparison with 1.43 eV). Then, another dehydrogenation reaction occurs on the oxygen side overcoming an energy barrier of 0.61 eV (TS-4) to produce *FBO (c10), which at last desorbs endothermically (0.46 eV) from the catalyst to yield FBO. In addition, the enthalpies and barriers for the first few steps of the reaction mechanism over the Cu 4 /MoO 2 @C were also calculated to interpret the impact of Fe substitution. DFT results (Supplementary Fig. 22) reveal a much higher energy barrier (TS-1) over Cu 4 / MoO 2 than Cu 4 /Fe 0.3 Mo 0.7 O 2 (1.65 eV vs 1.47 eV) within the first few steps, suggesting the atomic Fe species can indeed facilitate the reaction process through regulating the electronic properties.
After the first three elementary steps of pathway-1, the generated acetone may react with the residue FFA. Therefore, the reaction mechanism of FFA and the generated acetone are additionaly investigated in pathway-2. In the fourth elementary step, the energy barrier of the reaction between FFAL radical and acetone to produce (c11) is much higher in comparison with that of pathway-1 (0.91 to 0.78 eV). In addition, the subsequent dehydration reaction (from c12 to c13, TS-2(p2)) requires an energy barrier as high as 1.27 eV, also higher than that of the reaction pathway-1 (0.82-0.92 eV). Therefore, the reaction pathway-2 (the generated acetone react with FFA to give FBO) is not preferential in both thermodynamic and kinetic ways as compared with the reaction pathway-1.
To investigate the general applicability of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C, C 3~C10 secondary and primary alcohols were also employed to synthesize the corresponding C 8~C15 ketones and aldehydes as aviation biofuel intermediates under the optimized reaction conditions (Fig. 5d, Supplementary Tables 6, 7). Good to excellent yields (81.7-100%) of the C 8~C15 ketones and aldehydes are achieved through slightly modifying the reaction conditions (e.g., extending the reaction time or increasing the reaction temperature). These reaction results demonstrate the high practical potentials of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C in the valorization of the biomass-derived FFA to produce biofuels.
After reaction, the catalyst was easily isolated from the reaction mixture and directly reused after washing and drying. No obvious activity loss is observed in the continuous recycling tests ( Supplementary Fig. 23a), implying the good stability and recyclability of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C. Besides, the Cu, Fe and Mo contents in the reaction solution collected by hot filtration (Supplementary Fig. 23b) are below the detecting limits of ICP-OES, indicating the loss of metals is negligible, which could account for the maintenance of catalytic activity. Furthermore, XRD ( Supplementary Fig. 24f) and XPS ( Supplementary  Fig. 24i-k) results confirm the well-preserved crystallization and composition of the Cu 4 /Fe 0.3 Mo 0.7 O 2 @C. TEM ( Supplementary  Fig. 24a), HAADF-STEM ( Supplementary Fig. 24b), AC HAADF-STEM ( Supplementary Fig. 24c), the corresponding EDX mapping (Supplementary Fig. 24d) images and elemental line scan profiles (Supplementary Fig. 24e) and N 2 adsorption-desorption isotherms ( Supplementary Fig. 24g, h) demonstrate the well-retained multi-shelled hollow hierarchical nanostructure consisting of Cu 4 clusters on the Fe-doped MoO 2 (Cu 4 /Fe 0.3 Mo 0.7 O 2 @C).
In summary, we have demonstrated that the electronic property of M x in SCC can be efficiently modulated through atomically decorating the support. As an example, subnanometric Cu clusters are immobilized onto Fe-doped MoO 2 support to fabricate the Cu 4 /Fe 0.3 Mo 0.7 O 2 @C nanocomposite. Systematical characterizations reveal the uniform subnanometric Cu clusters composing of four Cu atoms in average, atomically Fe-doped MoO 2 support and distinctive metal-support interactions between them, which eventually lead to charge redistribution inside the subnanometric Cu clusters. The obtained Cu 4 / Fe 0.3 Mo 0.7 O 2 @C shows excellent catalytic performance in the newly developed reaction route of FFA to FBO, i.e., oxidative coupling of FFA with IPA, outperforming the conventionally prepared counterparts. DFT and control experiments uncover the reaction pathways and suggest the superior catalytic performance of Cu 4 /Fe 0.3 Mo 0.7 O 2 @C is mostly originated from the charge redistribution inside the Cu clusters. Moreover, this catalyst is also highly efficient for the preparation of C 8~C15 aviation biofuel intermediates using the corresponding secondary and even primary alcohols as feedstock. The proposed strategy of modulating the electronic structures of metal single clusters through atom-decoration of support may present a new dimension for the development of advanced catalysts with precisely adjusted electronic properties for various frontier applications.

Methods
Synthesis of NENU-5 truncated octahedron. All chemicals were purchased from commercial sources and directly used without further purification. In a typical synthesis, Cu(OAc) 2   Synthesis of Fe-NENU-5. Typically, the as-synthesized NENU-5 (0.05 g) was dispersed into a 60 mL aqueous solution containing FeCl 3 • 6H 2 O (0.15 g). After stirring for 40 min, the resulting yellow-green precipitates were collected by centrifugation, washed twice with deionized water and ethanol, and dried under vacuum at 50°C overnight to obtain Fe-NENU-5. Synthesis of Cu x /MoO 2 @C. Typically, the as-prepared NENU-5 was placed in a tubular furnace (BTF-1200C, Anhui BEQ Equipment Technology Co., Ltd.) and heated to 600°C at a ramp rate of 1°C min −1 and kept for 120 min under Ar atmosphere to yield Cu x /MoO 2 @C.

Synthesis of
Synthesis of Cu-BTC and Cu@C. In a typical synthesis, Cu(OAc) 2 • H 2 O (1 mmol, 0.2 g) was dissolved in 40 mL deionized water. After stirring for 1 h under ambient condition, a 40 mL ethanol solution containing H 3 BTC (0.67 mmol, 0.14 g) was added to the above solution under vigorous stirring. After stirring for 14 h, the resulting blue precipitates were collected by centrifugation, washed twice with deionized water and ethanol, respectively, and finally dried under vacuum at 50°C overnight to obtain Cu-BTC. The as-prepared Cu-BTC was placed in a tubular furnace (BTF-1200C, Anhui BEQ Equipment Technology Co., Ltd.), which was heated to 600°C at a ramp rate of 1°C min −1 and kept for 120 min under Ar atmosphere to yield Cu@C.
Synthesis of Fe-Cu-BTC and Fe-Cu@C. The as-synthesized Cu-BTC (0.05 g) was dispersed into a 60 mL aqueous solution containing FeCl 3 • 6H 2 O (0.15 g). After stirring for 40 min, the resulting precipitates were collected by centrifugation, washed twice with deionized water and ethanol, respectively, and finally dried under vacuum at 50°C overnight to yield Fe-Cu-BTC. The as-prepared Fe-Cu-BTC was placed in a tubular furnace (BTF-1200C, Anhui BEQ Equipment Technology Co., Ltd.), which was heated to 600°C at a ramp rate of 1°C min −1 and kept for 120 min under Ar atmosphere to yield Fe-Cu@C.
Synthesis of Fe/MoO 2 . In a typical synthesis, the commercial MoO 2 (3.75 g) was added to 60 mL aqueous solution containing FeCl 3 • 6H 2 O (0.15 g) under vigorous stirring. After stirring for 14 h under ambient condition, the resulting precipitates were collected and then heated to 600°C at a ramp rate of 1°C min −1 and kept for 120 min under Ar atmosphere to yield Fe/MoO 2 .
Synthesis of Cu-Fe-MoO 2 . In a typical synthesis, Cu(OAc) 2 • H 2 O (0.1 g), FeCl 3 • 6H 2 O (0.08 g), and 2.0 g commercial MoO 2 were added to 50 mL deionized water. The mixture was dried and subsequently subjected to pyrolysis at 600°C with a ramp rate of 1°C·min −1 and kept for 120 min under Ar atmosphere to yield Cu-Fe-MoO 2 .
Materials characterization. The size and morphology of materials were studied by scanning electron microscopy (SEM) and transmission electron microscopy (TEM). SEM was carried out on a JEOL-6700 instrument. High-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and spherical aberration correction HAADF-STEM (AC HAADF-STEM) were recorded on a FEI Titan Cubed Themis G2 300S/TEM with a probe corrector and a monochromator at 200 kV. Powder X-ray diffraction (PXRD) patterns of the samples were obtained on a Rigaku diffractometer (D/max-IIIA, 3 kW) with Cu K α radiation (λ = 1.5406 Å) at a voltage of 40 kV and a current of 10 mA at room temperature. Brunauer-Emmett-Teller (BET) surface area and pore-size measurements were performed on a Micromeritics ASAP 2020 M instrument at 77 K. Before the analysis, the samples were degassed at 50°C for 12 h. X-ray photoelectron spectroscopy (XPS) was collected on a Thermo Scientific K-Alpha system with the C 1s peak (284.6 eV) as reference. The metal contents of the samples were determined by ICP-OES on an Optima 8300 instrument. The C and O elemental contents of the samples were measured on a Euro Vector EA3000 instrument. Thermogravimetric analysis (TGA) was performed on a NETZSCH STA449C instrument loaded with 5 mg sample using a heating rate of 5°C/min under argon atmosphere. The X-ray absorption experiments were carried out at the XAS station (BL14W1) of the Shanghai Synchrotron Radiation Facility (SSRF). The electron storage ring was operated at 3.5 GeV. Si (311) double-crystal was used as the monochromator, and the data was collected using solid-state detector under ambient conditions. The beam size was limited by the horizontal and vertical slits with the area of 1 × 4 mm 2 during XAS measurements.
Calculation details. All the calculations are performed in the framework of the density functional theory with the projector augmented plane-wave method, as implemented in the Vienna ab initio simulation package. The generalized gradient approximation proposed by Perdew, Burke, and Ernzerhof is selected for the exchange-correlation potential. The cut-off energy for plane wave is set to 400 eV. The energy criterion is set to 10 −5 eV in iterative solution of the Kohn-Sham equation. All the structures are relaxed until the residual forces on the atoms have declined to <0.05 eV/Å. The weak interaction was described by DFT + D3 method using empirical correction in Grimme's scheme. The vacuum space was set to be more than 12 Å, which was enough to avoid the interaction between periodical images. For MoO 2 , a four layer of (−1 1 1) surface was used, and 33% Mo atoms was replaced by Fe for the structure of Fe-doped MoO 2 . The Brillouin zones of structures were sampled with 3 × 3 × 1 k points. The minimum energy pathway for transition state searching process was determined by using a climbing image nudged elastic band method (CINEB).
One-pot cascade oxidative coupling of FFA with IPA. The one-pot cascade oxidative coupling of FFA with IPA was carried out in a 25 mL high-pressure reactor (NSG25-P5-T3-SS1-SV, Anhui Chem-n Instrument Co., Ltd.) equipped with a magnetic stirrer. In a typical run, FFA (0.5 mmol), catalyst (Cu 3.6 mol% relative to FFA), K 2 CO 3 (0.1 mmol), and IPA (5 mL) were sealed in the reactor. The reactor was evacuated, refilled with 2 bar O 2 , and heated to 120°C for 16 h under magnetic stirring. After reaction, the reactor was cooled to room temperature. The catalyst was isolated from the solution by centrifugation and directly reused after washing and drying. The product was quantified by a GC-MS spectrometer (Agilent, 7890B GC/5977A MS) equipped with a DB-35 MS UI capillary column (0.25 mm × 30 m).
The conversions and selectivities were calculated using the following equations: Data availability