Pd/CsPMA: structural features and H 2 -reduction/re-oxidation properties.
Pd/CsPMA was prepared by precipitating PMA anions ([PMo12O40]3−) with Cs+ and Pd2+ at sub-ambient temperatures. The separated yellow solids, as observed by scanning electronic microscopy (SEM) and transmission electronic microscopy (TEM), are sphere particles with a size approximately ranging from 50 to 300 nm (Figures S1 and S2). X-ray diffraction (XRD) results suggest that the Pd/CsPMA material has a cubic crystal structure distinct from those of MoO3 and PMA (Fig. 1a). From integrated differential phase contrast imaging (iDPC), the crystals are assembled from individual PMA units (Fig. 1b, S3-S5), based on which a structural model of Pd/CsPMA material was constructed (Figure S6). Pd/CsPMA and PMA gave similar Raman spectra, suggesting the preservation of the anion structure during catalyst synthesis (Fig. 1c). The actual Pd loading of Pd/CsPMA was determined to be 0.23 wt.% by inductively coupled plasma optical emission spectroscopy. Pd species distribute evenly on Pd/CsPMA, as shown in the high-angle annular dark-field (HAADF) imaging results under STEM mode (Fig. 1d). The X-ray photoelectron spectroscopy (XPS) analysis (Fig. 1e) and in situ diffuse reflectance infrared spectroscopy (DRIFTS) with NO as the probe molecule (Fig. 1f) confirm that Pd exists as charged Pd cations, likely situated at the 4-fold hollow sites of PMA by replacing Cs+, as indicated in previous studies on single-atom POM systems.38,41
When Pd/CsPMA was exposed to forming gas (5% H2/N2) flow at room temperature (24 ℃), a rapid color change from yellow to dark blue was observed within twenty minutes. In situ H2-DRIFT spectrum of Pd/CsPMA present a notable peak at above 3200 cm− 1 for Pd/CsPMA (Figure S7), indicating the formation of large amount of surface hydroxyl groups, which is not the case for CsPMA. This implies that Pd plays a crucial role in dissociating H2 into H atoms and transferring H species from the noble metal sites to the support.42,43 The H atoms donate electrons to PMA and combine with PMA oxygen to form − OH groups on the surface. In situ NO-DRIFTS analysis confirms that Pd species remain highly dispersed during H2 treatment at room temperature (Figure S8). We then prepared a Pd/CsPMA-H material, by treating the Pd/CsPMA under H2 (4 bar) in water for 15 min. As shown in Fig. 1g, the color of the catalyst also changed from yellow (Pd/CsPMA) to dark blue (Pd/CsPMA-H). UV-Visible (UV-Vis) adsorption spectra of the Pd/CsPMA-H suspension show an obvious decrease of Mo(VI) signal at 320–340 nm and an increase of Mo(V) and Mo(IV) in the range of 600–900 nm, illustrating the reduction of Mo by hydrogen spillover (Fig. 1h).44,45 Determined by redox titration using acidified FeCl3 solution (Figure S9), each POM anion on average accepts 7.5 electrons after 15-minute H2 treatment, corroborating the UV-Vis adsorption data that Mo(VI) was only partially reduced. Additional H2 treatment does not increase the degree of hydrogen spillover.
The reduced Pd/CsPMA-H powder readily re-oxidizes when exposed to air at room temperature (represented as Pd/CsPMA-H-O), as judged from the recovered yellowish color and the UV-Vis absorption results (Fig. 1g and 1h). This highlights that the reduced Pd/CsPMA can activate O2 at ambient temperature. Pd/CsPMA-H-O exhibits similar morphology, crystalline structure, and POM anion characteristic Raman peaks to those of the pristine Pd/CsPMA (Figure S1, S10, 1a and 1c), which prove that the repeated reduction-oxidation steps would not change the structure of Pd/CsPMA, thus hinting at its potential as a catalyst in redox reactions.
Catalytic performance in partial methane oxidation. The methane oxidation reaction was carried out at room temperature in an aqueous solution with CH4 and O2.46 The unactivated Pd/CsPMA catalyst exhibited negligible methane conversion activity (Table 1, Entry 1), and after reaction the Mo species in Pd/CsPMA maintained an oxidation state of + 6 (Figure S11). To validate adding electrons to the catalyst is conducive to aerobic oxidation of methane, Pd/CsPMA was first treated by H2 in the liquid phase at room temperature to form the reduced catalyst, Pd/CsPMA-H. When Pd/CsPMA-H was used for methane oxidation with O2 under identical conditions for 30 mins, we observed methane conversion with methanol as the sole product (5.2 µmolgcat−1h− 1) (Table 1, Entry 2). Shorting the reaction to 5 minutes provided proportionally increased methanol production rate (28.8 µmolgcat−1h− 1), suggesting that Pd/CsPMA-H is only active in the initial stage of reaction (Table 1, Entry 3). By adjusting the partial pressure of CH4 and O2, it was possible to improve methanol productivity, with this metric rising to 67.4 µmolgcat−1h− 1 (Table 1, Entry 4). The color of the spent catalyst changed to greenish yellow, suggesting that they were re-oxidized by O2 during the course of the reaction (Figure S11). These results highlight two important findings. First, H2-reduced Pd/CsPMA selectively transform methane and O2 into methanol at room temperature. Second, Pd/CsPMA-H does not maintain its reduced state under an oxidizing atmosphere, thus quickly losing its initial activity.
To sustain a reduced catalyst state, we then attempted the methane oxidation reaction using a non-explosive gas mixture of H2, O2, CH4 and N2, while H2 pretreatment was no longer applied. Under an optimized H2/O2 ratio, Pd/CsPMA presents a methanol productivity of 28.5 µmolgcat−1h− 1 (a noble metal specific productivity of 12.4 mmolgPd−1h− 1), again with methanol as the only liquid product detected together with negligible gas-phase product in a 30-minute reaction (Table 1, Entry 5). The activity and selectivity are exceptional among most reported noble-metal-containing catalysts for thermal catalytic aerobic oxidation of methane to methanol at room temperature (Table S1). Interestingly, we observed that in a typical 30-mins reaction, the reaction proceeds at higher rates in the 2nd 15 mins compared to the 1st, likely due to an induction period required to form the reduced catalyst that is active for methane oxidation (Figure S12). The ratio of H2 and CH4 almost remained constant during the reaction, suggesting that H2 is not preferentially consumed by O2 as compared to CH4 (Table S2).
Table 1
Oxidation of methane to methanol using Pd/CsPMA catalyst. Reaction condition: 2 mL D2O, 10 mg catalyst, room temperature, 800 rpm. For Entry 1, 10 bar CH4, 1 bar O2 and 9 bar N2 were applied. For Entry 2, 3, 7 and 8, the catalyst was treated with 4 bar of H2 for 15 min at room temperature before running the reaction under the same condition. For Entry 4, 20 bar of CH4 was used. For Entry 5, 10 bar CH4, 4 bar H2 and 0.3 bar O2 balanced with 25.7 bar N2 were applied. For Entry 6, the solvent used was H2O instead of D2O. For Entry 7 and 8, 1 µmol of ascorbic acid and Na2S were added, respectively, after H2 pretreatment but before charging methane and O2.
Entry | Condition (number in parenthesis refer to gas pressure. unit: bar) | CH3OH prod. (µmol gcat−1 h− 1) | CH3OH Sel. | Specific activity (mmol gPd−1 h− 1) | Color after reaction |
1 | Single step: CH4(10) + O2(1), 30 min | 0 | N/A | 0 | Yellow |
2 | Step 1: H2(4); Step 2: CH4(10) + O2(1), 30 min | 5.2 | 100.0% | 2.3 | Yellow |
3 | Step 1: H2(4); Step 2: CH4(10) + O2(1), 5 min | 28.8 | 100.0% | 12.5 | Green |
4 | Step 1: H2(4); Step 2: CH4(20) + O2(0.3) 5 min | 67.4 | 100.0% | 29.3 | Green |
5 | Single step: CH4(10) + H2(4) + O2(0.3), 30 min | 28.5 | 100.0% | 12.4 | Dark blue |
6 | Single step: CH4 + H2O2(400µmol), 30 min | 0 | N/A | 0 | Yellow |
7 | Single step: H2(4); Ascorbic acid, CH4(10) + O2 (1), 5 min | 42.2 | 100.0% | 18.4 | Green |
8 | Single step: H2(4); Na2S, CH4(10) + O2(1), 5 min | 48.0 | 100.0% | 20.9 | Blue |
Xiao and colleagues previously reported29 selective methane oxidation with H2 and O2 using AuPd particles confined in a zeolite support. In that case, H2O2 was produced locally, which then acts as the oxidant for methane oxidation. In contrast, H2O2 is unlikely to be the key oxidant in our system, because (i) no H2O2 was detected when supplying a mixture of H2 and O2 over the catalyst (detection limit: 0.5 ppm, Figure S13) and (ii) when external H2O2 was added, no methanol was detected after 30 mins despite 54% of the H2O2 was decomposed (Table 1, Entry 6 and Table S3). Methanol was produced only when H2O2 was used together with H2 (Table S3), but the best activity (e.g., 14.4 µmolgcat−1h− 1) was still not as high as the case when H2 and O2 are supplied. Molecular O2 works as a better oxidant than H2O2 in our case, possibly because as a weaker oxidant, O2 allows the catalyst to maintain the desired reduced state under H2.
Methane oxidation was performed under varied partial pressure of CH4, H2, and O2. The reaction is first order with respect to CH4 between 0–15 bar (Fig. 2a), indicating methane activation participates in the rate-determining step. On the other hand, the methanol yield as a function of the H2 to O2 ratio follows a volcano-type curve, with the optimal pressure of H2 at 4 bar (Fig. 2b). A threshold of 2 bar H2 partial pressure must be reached to create a state with sufficient degree of Mo reduction for C − H bond activation. When O2 is eliminated from the reaction system, methanol was only detected in trace amount (2 µmolgcat−1h− 1), likely due to residual O2 in water.
Several control catalyst samples were studied to confirm the critical role of reduced Mo species as catalyst active sites. CsPMA, Rh/CsPMA and Pt/CsPMA, showing no reduction of PMA by H2 at room temperature, exhibited no methanol productivity (Fig. 2c). Rh/CsPMA and Pt/CsPMA were then pretreated with H2 at elevated temperature (50 ℃ and 100 ℃, respectively). After the pretreatment, both catalysts turned dark blue signaling Mo reduction, and in the following step both catalysts produced methanol from methane and oxygen at room temperature (Table S4). Replacing PMA by phosphotungstate (PTA), Pd/CsPTA also did not show methane activation under the same conditions (Figure S7). To explore whether our finding, i.e., H2-reduced Mo in Pd/CsPMA readily converts methane and O2 to methanol, is generalizable, 1 wt.% Pd/MoO3, 1 wt.% Pd/WO3 and 10 wt.% Pd/C were evaluated under our methane oxidation conditions (Fig. 2d, Table S5). The H2-pretreated 1 wt.% Pd/MoO3-H indeed exhibited a methanol productivity of 14.4 µmolgcat−1h− 1 in a 5-min reaction, while the other two catalysts were inactive. Interestingly, MoO2, with Mo(IV) as the only Mo species, showed no methane activation properties with O2. The above results indicate that Pd is not indispensable for methane activation, but Mo-species under a partially reduced chemical state is a critical element for methane activation with oxygen.
Methanol formation pathway on reduced PMA: experimental evidence. Since radical mechanism is one of the mainstream mechanisms to activate methane C − H bond in aqueous solution, EPR experiments using 5,5-dimethyl-1-pyrroline N-oxide (DMPO) radical trapping agent were performed to determine the presence of reactive oxygen species (Figure S14).47 The appearance of characteristic quadruple peaks (1:2:2:1) for H2O2 aqueous solution indicates trapped ·OH radicals, while a triple peak was observed for the aqueous suspension of pristine Pd/CsPMA, probably coming from the oxidative dimerization of DMPO.48 However, these are the cases where no methanol productivity was observed. By contrast, for the H2 pretreated Pd/CsPMA-H catalyst, no similar signals were detected. To provide further evidence that ·OH radicals do not operate in our system, methane oxidation reactions were carried out over the Pd/CsPMA-H catalyst in the presence of ·OH radical scavengers (Na2S or ascorbic acid) (Table 1, Entry 6–7). In these cases, the methanol productivity of the Pd/CsPMA-H catalyst did not decrease, but rather increased to 42.2 and 48.0 µmolgcat−1h− 1, respectively, compared to the 28.8 µmolgcat−1h− 1 without the scavengers. This suggests that methane is not activated by free ·OH intermediates. The increased productivity may be ascribed to the elongated lifetime for reduced PMA species under O2, due to the reductivity of scavengers.
To further rule out the role of Pd in C − H activation, selective poisoning tests were conducted by adding benzyl mercaptan (BzM). BzM binds strongly with Pd, thus preventing Pd-reactant interaction (Fig. 3a and S15).40 As expected, after introducing 2 equivalents (relative to Pd) BzM to block Pd in Pd/CsPMA before H2 treatment, PMA was not able to be reduced to enable methane oxidation (BzM-H2-CH4/O2 in Fig. 2e). In contrast, when the same amount of BzM was added to pre-reduced Pd/CsPMA-H catalyst (H2-BzM-CH4/O2 in Fig. 2e), the methanol productivity mimics the performance without BzM. This experiment, together with the earlier presented fact that reduced Rh/CsPMA and Pt/CsPMA are also active in methane oxidation to methanol at room temperature, provide compelling evidence that methane activation proceeds not on Pd sites but on reduced PMA.
The reduced catalyst Pd/CsPMA-H has considerable solubility in H2O, likely due to the strong interaction between spilt H-induced surface OH groups and water molecules.51,52 This provides an opportunity to determine whether methanol formation requires an extended surface, or individual PMA anions. Hence, H2-reduced Pd/CsPMA-H aqueous suspension was separated by ultrafiltration centrifugation (Nominal Molecular Weight Limit: 10 kDa) and evaluated in methane oxidation. The resulting solution and solid residue were both found to offer activity (Figure S16), suggesting that H2-reduced PMA is the active sites for methane conversion, regardless of staying in the solution or solid phase.
To enable catalyst recycling, the spent reaction solution was first re-oxidized by 1 bar O2 at room temperature, followed by the addition of the Cs+ salt as a precipitation agent. Over 99% of the Pd/CsPMA catalyst can be regenerated (Table S6). A reductive atmosphere (2 bar H2, 0.3 bar O2, 10 bar CH4 diluted with N2 to a total pressure of 40 bar) was used for reaction sessions during successive recycling tests to ensure sufficient activation of Mo in each cycle. No loss in methanol formation activity was observed over five successive reactions (Fig. 2f), and furthermore the recovered catalyst exhibited identical XRD patterns compared with the fresh one (Figure S17), indicating that PMA structure remained intact during recycling.
In view of the important role of reduced PMA on methane activation, we investigated the coordination structure around Mo before and after hydrogen reduction by pseudo in situ electron paramagnetic resonance (EPR). As expected, no signal was detected for Pd/CsPMA before H2 treatment since the Mo species predominantly exist as Mo(VI), which is EPR inactive (Figure S18).49 Measured under the H2 atmosphere, Pd/CsPMA-H exhibits a broad, featureless EPR spectrum centred at a g value of approximately 1.93 (Fig. 3a), which can be ascribed to the emergence of Mo(V) species.50–52 The unresolved hyperfine splitting and the asymmetric line shape which almost resemble a broad isotropic line, would point to a Mo(V) centre with six O atoms at close bond lengths where free movement or rotation is relatively unhindered.53 Afterwards, the EPR tube containing Pd/CsPMA-H was vacuumized and sealed (denoted as Pd/CsPMA-H-vac). The treatment resulted in an anisotropic EPR spectrum with much narrower line shape and resolution of the splitting arising from the hyperfine interaction of the unpaired electron in Mo(V) with the 95Mo and 97Mo nuclei, both with I = 5/2 and overall natural abundance of 25.5%. This spectral change would imply a large extent of anisotropic distortion of MoOx polyhedra in vacuo (Fig. 3a). Based on simulation results (Fig. 3b, Figure S19, Table S7), the spectrum consists of the signals of three Mo(V) components, with the dominant one containing an O vacancy on the bridge site of two MoOx polyhedra, signified as F (g1 = 1.955, g2 = 1.954, g3 = 1.857, weight = 62%) (Figure S20).54,55 The remaining signals originate from pentacoordinated Mo(V) without molybdenyl oxygen (structure C) and hexacoordinated Mo(V) center (structure A), respectively. Although the results cannot give information on EPR-silent, deep reduced Mo(IV) species, the insights about Mo(V) already provide knowledge on the largely weakened Mo − O bonds and the tendency of O deficiency around Mo centers under reduced state that may facilitate interaction with oxygen or methane molecules.
When Pd/CsPMA-H was exposed to CH4 and then analyzed by EPR (Fig. 3a), the spectrum exhibited a similar pattern as Pd/CsPMA-H-vac, which illustrates the reduced Mo does not directly coordinate with methane molecules. However, when Pd/CsPMA-H was exposed to air, the EPR signals disappeared immediately. Fully consistent with EPR data, in situ Raman experiments reveal that Mo-Ob-Mo and Mo-Oc-Mo are significantly weakened by H2 reduction, which is recovered after treating with O2 (Figure S21). Based on these, we posit that the MoOx sites with O vacancies first activate O2, creating surface active O to then enable C − H activation by abstracting H from methane.
Theoretical insights into the formation of methanol over Pd/CsPMA-H. DFT calculations were performed to further verify potential mechanisms for CH4 activation and CH3OH formation over the Pd/CsPMA-H catalyst (Fig. 4 and Figure S22). For the creation of active sites, calculations show H2 activation occurs heterolytically over the Pd single atom site with a relatively low barrier of 0.40 eV, resulting in a OH and PdH species (intermediate III). Subsequently, spillover of H to two neighboring oxygen sites is thermodynamically favorable (intermediate VI), followed by H2O generation with a slightly higher barrier of 0.50 eV (intermediate V). Upon formation, the H2O species remains adsorbed bridging two Mo sites, and subsequently desorb to the gas phase which is energetically uphill by 0.81 eV (intermediate VI). This results in a lattice oxygen vacancy in the bridging O site which is unreactive towards methane but can readily adsorb O2 from the gas phase with an adsorption energy of -0.14 eV (intermediate VII). In this adsorbed state, one oxygen from O2 is situated in the lattice O position, while the second oxygen extends away from the surface and does not chemically bond with the surface. The O-O bond length of the adsorbed species is also significantly elongated from the DFT gas-phase length of 1.23 Å to 1.32 Å, suggesting the O2 has become highly activated and could be reactive to methane. In accordance with this expectation, the C − H activation on the adsorbed O2 occurs with a barrier of 0.61 eV to form a O2H species and a CH3 radical (intermediate VIII). Subsequently, the methane radical can be captured by the same O2H species via a rebound adsorption mechanism to form methanol, simultaneously cleaving the O − O bond and regenerating the lattice oxygen site (intermediate IX).
In this mechanism, methane activation is enabled by the formation of activated O2 species resulting from H spillover and vacancy formation. The rate-limiting step is the C − H activation, in agreement with experimentally observed first-order kinetics for methane. We also find this mechanism to agree with the volcano-type curve for H2 to O2 ratio, as a balance must be struck to ensure the survival of the activated O2 species (intermediate VI). An overabundance of gaseous O2 will react with hydroxyl groups and prevent vacancy formation, while an overabundance of H2 will likely react with the adsorbed O2 before it can activate the methane. The mechanism was calculated on oxygen/vacancy sites which are non-adjacent to the Pd, indicating that the presence of Pd was not required for methanol formation beyond the initial H2 activation step, in accordance with the poisoning studies. Alternative mechanisms were also considered which were found to be unfavorable or inconsistent with experimental observations (Figure S23-25, as well as discussions provided there).