2019/7/29 スイス連邦工科大学チューリッヒ校(ETH) （チューリッヒ工科大学）
・ ETH と仏・Total社が、CO2 と水素をメタノールに効率的に直接変換する、酸化インジウムとパラジウムを使用した触媒を新たに開発。
・ ETH は数年前に酸化インジウム触媒によるメタノール生成（副生物はほぼ水のみ）の実証に成功し、高度な安定性を確認しているが、不十分な触媒活性度を補うための、酸化インジウムの多量利用が商用化の可能性を阻んでいた。
・ 現在、産業規模のメタノール生成は化石燃料由来であり、これに伴ってカーボンフットプリントも高くなっている。CO2 からメタノールを生成する今回の技術では、メタノールから燃料を合成・燃焼後に CO2 を再利用するため、クローズドなカーボンサイクルの実現が可能と考える。
・ ETH と Total 社は、本技術について特許を出願。Total社は同技術のスケールアップと実証ユニットでの同技術の実証を数年内に実施予定。
（関連情報) Nature Communications 掲載論文(フルテキスト）
Atomic-scale engineering of indium oxide promotion by palladium for methanol production via CO2 hydrogenation
–scale engineering of indium oxide promotion by palladium for methanol production via CO2 hydrogenation
Metal promotion is broadly applied to enhance the performance of heterogeneous catalysts to fulfill industrial requirements. Still, generating and quantifying the effect of the promoter speciation that exclusively introduces desired properties and ensures proximity to or accommodation within the active site and durability upon reaction is very challenging. Recently, In2O3 was discovered as a highly selective and stable catalyst for green methanol production from CO2. Activity boosting by promotion with palladium, an efficient H2-splitter, was partially successful since palladium nanoparticles mediate the parasitic reverse water–gas shift reaction, reducing selectivity, and sinter or alloy with indium, limiting metal utilization and robustness. Here, we show that the precise palladium atoms architecture reached by controlled co-precipitation eliminates these limitations. Palladium atoms replacing indium atoms in the active In3O5 ensemble attract additional palladium atoms deposited onto the surface forming low-nuclearity clusters, which foster H2 activation and remain unaltered, enabling record productivities for 500 h.
Metal promotion is a key strategy in heterogeneous catalysis to enhance activity, selectivity, and stability of catalytic centers by tuning of their electronic and geometric properties. Numerous industrial catalysts rely on promoters to comply with requirements related to reactor scale, operating conditions, and atom economy, in order to minimize costs, waste production, and energy input, ensuring higher profits and environmental friendliness1. To maximize benefits, the added element should carry the required functionality without introducing detrimental properties, be located in high vicinity to or within the active sites, and remain unperturbed upon use in the targeted reaction. Fulfilling all these goals is a demanding task that necessitates precision synthesis methods defined based on information provided by in-depth characterization, mechanistic elucidation, and testing about the most suitable promoter speciation2,3,4. The latter is a broad-term embracing identity, oxidation state, and molecular structure of the chemical species.
In the frame of research efforts devoted to mitigating global warming and overcoming our reliance on fossil feedstocks5,6,7,8,9,10,11,12,13,14,15,16, promotion by palladium has been attempted on a breakthrough catalyst, In2O3, recently discovered for CO2 hydrogenation to methanol17,18. Indeed, although this reducible oxide is very selective and durable, especially when supported on ZrO2, owing to the presence of frustrated Lewis pairs19,20 and single-ensembles21,22, its H2-splitting ability is limited23. Palladium was shown to generate bulk intermetallic compounds with indium24,25, which were superior to In2O3 in liquid-phase CO2 hydrogenation26, but were exclusively active for the competitive reverse water–gas shift (RWGS) reaction when applied in the gas phase as nanoparticles supported on silica27. A synthetic method devised to deposit palladium nanoparticles on In2O3 minimizing alloy formation28 led to a system displaying enhanced methanol formation, which was ascribed to the higher availability of activated hydrogen, fostering the desired hydrogenation reaction and vacancy formation on the oxide surface. These contrasting findings highlight that the true relevant speciation and action mechanism of palladium are poorly understood. More importantly, to maximize the promoting effect and thus process-level advantages, the H2-splitting ability of palladium should be exploited while minimizing its high efficacy for the RWGS reaction as a stand-alone metallic phase29,30,31. As single palladium atoms supported on carbon nitride are active in alkyne semi-hydrogenation and photocatalytic CO2 reduction32,33,34 and possess distinct electronic properties compared to agglomerated palladium, calculations indicate that the (R)WGS reaction on Pd(111) requires at least three palladium atoms35, and the selectivity of platinum on MoS2 in CO2 hydrogenation was dictated by the cluster mono-/binuclearity36, we conceived atomic dispersion as a strategy to curtail detrimental effects while preserving beneficial attributes. This approach would simultaneously grant improved metal utilization and robustness, if palladium atoms are well anchored to the In2O3 structure. Aiming at deriving synthesis–structure–performance relations, we emphasized experimental methods as well as theoretical means, which comprise a prominent tool to uncover the impact of structure on activity and selectivity37,38,39 but have rarely been applied to link preparation to structure.
Here, we introduce a co-precipitation method effective in incorporating isolated palladium atoms in the In2O3 lattice forming low-nuclearity palladium clusters anchored to the active site, as elucidated by state-of-the-art characterization and density functional theory (DFT) investigations. Kinetic and long-term tests show how this speciation overcomes selectivity and stability limitations associated with palladium nanoparticles, rationalizing the unparalleled performance of this catalyst for CO2-based methanol synthesis.
Impact of promoter nature and loading, and synthesis method
We prepared palladium-promoted In2O3 by co-precipitation and a sol–gel method targeting an atomic palladium distribution in the oxide, and applied spray deposition and dry and wet impregnation as well as simultaneous chemical reduction of palladium and indium salts to attain benchmark materials comprising palladium nanoparticles and a PdIn intermetallic compound. All catalysts possessed sufficiently high surface area and pore volume (SBET > 100 m2 g−1, Vpore > 0.3 cm3 g−1) and the nominal palladium loading of 0.75 wt.% was closely matched (Supplementary Table 1). Their catalytic performance was assessed in CO2 hydrogenation to methanol at 553 K and 5 MPa and compared with that of pure In2O3. The presence of palladium increased the methanol space–time yield (STY) in all cases, except for the intermetallic solid. The co-precipitated catalyst (CP) was the best material (STY = 0.61 gMeOH h−1 gcat−1), closely followed by that produced by dry impregnation (DI). Reference solids comprising palladium introduced into the inert TiO2 support by CP and DI mainly produced CO, as expected from the reported RWGS activity of palladium nanoparticles31. After confirming the superiority of palladium with respect to other typical hydrogenation metals (Supplementary Fig. 1), we unraveled the impact of its loading by testing CP and DI materials containing 0.25–10 wt.% Pd (Fig. 1a). The initial loading of 0.75 wt.% turned out to be optimal, as the STY was inferior at lower contents and remained equal or dropped for higher amounts. Interestingly, the CP catalyst with this palladium content retained a STY of 0.60 gMeOH h−1 gcat−1 for 93 h on stream after a minimal initial deactivation, whereas the DI counterpart lost half as of its productivity after 74 h (Fig. 1b). Comparing the methanol selectivity (SMeOH) over CP and DI catalysts and In2O3 at a CO2 conversion of 3% (97, 78, and 89%, respectively, Fig. 1c) revealed that palladium has a beneficial effect on the product distribution when added through CP but not via DI.
Catalytic performance of palladium-promoted In2O3 catalysts. a Impact of Pd loading on methanol space–time yield (STY) after 1 h on stream for palladium-promoted In2O3 catalysts obtained by co-precipitation (CP) and dry impregnation (DI), b temporal evolution of the methanol STY for CP and DI catalysts, with pure In2O3 as a reference, and c methanol selectivity over the same materials after 48 h on stream, when a CO2 conversion of 3% is attained by operation at the indicated WHSV. In b, data of an additional CP catalyst is presented that was aged for 200 instead of 1 h. Reaction conditions: T = 553 K, P = 5 MPa, H2:CO2 = 4, and WHSV = 24,000 cm3STP h−1 gcat−1 unless otherwise stated
Experimental description of palladium sites
CP and DI samples were studied through a battery of characterization techniques to elucidate the palladium speciation. X-ray diffraction (XRD) confirmed the bixbyite (i.e., a defective fluorite-type) structure of In2O3 in both systems (Supplementary Fig. 2), the lowest-energy surface of which presents a protruding In3O5(O) substructure, which generates the active site in acetylene and CO2 hydrogenation after removal of the (O) atom21, and an adjacent sunken region. The average particle size of In2O3 was calculated at 9 nm for the fresh materials and increased to 15 and 25 nm for used CP and DI catalysts, respectively. Moderate sintering was also observed earlier for unpromoted In2O317,22. A weak reflection at 39.5° 2θ specific to Pd0 was detected only for the used DI material. High angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images and indium and palladium maps acquired by energy-dispersive X-ray spectroscopy (EDX) indicate that palladium is extremely well dispersed in the fresh samples (Fig. 2). The CP material used in the reaction for 93 h appeared virtually unaltered. In contrast, In2O3 particles increased their size to ca. 30 nm and palladium nanoparticles formed using the DI catalyst for 74 h (Supplementary Fig. 2), the latter having an average size of 2.8 nm based on high resolution transmission electron microscopy (HRTEM, Supplementary Figs. 3 and 4). A determination of the palladium dispersion in CP and DI catalysts was attempted by H2 and CO chemisorption. The reduction behavior of the fresh catalysts was studied (Supplementary Fig. 5) through temperature-programmed reduction with H2 and CO (H2-TPR and CO-TPR) and temperature-programmed desorption of CO (CO-TPD) to identify suitable conditions for the reductive pre-treatment and the volumetric analysis. Barely measurable H2 and CO uptakes were registered for used DI samples (Supplementary Fig. 5). We put forward strong metal support interactions (SMSI) and/or Pd–In alloying as plausible reasons for the dramatic change in electronic properties of palladium particles in promoted In2O3 samples. The similar palladium dispersion determined for a Pd/TiO2 material with the same promoter loading and average particle size as DI 74 h (Supplementary Fig. 4), a catalyst for which the SMSI phenomenon has been widely documented40, supports our hypothesis. No and a very weak CO chemisorption was detected by diffuse-reflectance infrared Fourier transform spectroscopy (CO-DRIFTS, Supplementary Fig. 6) for fresh and used CP samples and DI 74 h, respectively. The first evidence agrees with studies reporting that isolated charged palladium species do not adsorb CO32 and isolated metallic Pd atoms weakly interact with CO at close-to-liquid-nitrogen temperatures41, and the second with the volumetric analysis. Based on electron microscopy, atomic dispersion of surface palladium species is plausible for fresh CP and DI catalysts and CP 1 h, while a value of 42% was estimated for DI 74 h considering the TEM-determined average particle size and a spherical geometry42.
Microscopy analysis of palladium-promoted In2O3 catalysts. HAADF-STEM images and corresponding EDX maps of indium (blue) and palladium (green) for the catalysts produced by a CP and b DI with 0.75 wt.% Pd in fresh form and after the tests depicted in Fig. 1b. The scale bar in all panels is 50 nm
According to X-ray photoelectron spectroscopy (XPS), the fresh CP and DI catalysts exclusively contained Pd2+ species (Fig. 3a), which chiefly reduced to Pd0 upon use in the reaction for 1 h. Residual Pd2+ species are likely owing to the short exposure of the samples to air for the ex situ analysis, as no uptake possibly related to the reduction of oxidized palladium was evidenced by H2-TPR for CP used for 1 h at H2:CO2 ratios of 4 and 2 (Supplementary Fig. 7). Despite the same bulk metal loading, palladium is almost 3-times more abundant on the surface of the DI catalyst (1.1 at.%, Supplementary Table 2) than on the CP sample (0.4 at.%), which supports Pd incorporation throughout the oxide crystals in the latter. After removing ca. 4 nm of the surface by sputtering with Ar+ ions, metallic palladium was still detected for DI 1 h. As this material only contains palladium deposited onto the In2O3 surface, the signal of the promoter should vanish. Still, information on the bulk is obviously mixed with that of the surface, as we are dealing with a polycrystalline material with an average particle size of ca. 8 nm, and XPS has a penetration depth of ca. 6 nm. For CP 1 h, palladium in oxidized form was evidenced for the sputtered sample, suggesting that the promoter in the interior of the oxide is not reduced upon reaction. The Pd3d signal of CP 16 h exhibited unaltered shape and intensity, as expected based on the stability of this material. In contrast, the surface Pd content for the DI material almost halved after 16 h on stream, in agreement with its sintering.
Experimental and theoretical elucidation of palladium-promoted In2O3 catalysts. a XPS Pd3d core-level spectra, and b Pd–K-edge EXAFS spectra of the CP and DI catalysts in fresh form and after use in CO2 hydrogenation for 1 h. EXAFS spectra simulated for the latter are added in red. The dashed lines mark expected positions of signals of the species indicated52,53,54. c Abundance of Pd-containing fragments emitted from the surface of used CP and DI materials in TOF-SIMS analysis. d H2-TPR profiles at 5 MPa of the fresh CP and DI catalysts with In2O3 as a reference. e XPS O1s core-level spectra of CP and DI catalysts after use for 1 and 16 h. The contribution of lattice oxygen, oxygen close to a defect, and oxygen in a hydroxyl group retrieved from deconvolution are marked in red, green, and blue, respectively, and labeled with their relative abundance in %. f In115 NMR spectra collected for pure In2O3 and the CP and DI catalysts in fresh form and after use for 16 h. g (left) Structure of a regularly terminated In2O3(111) surface featuring, owing to the strong anisotropy, a depression (pocket) where Pd atoms preferentially accommodate in the fresh DI catalyst. The coloring of In atoms at the surface represents the energy associated with their replacement by palladium (Eseg, in eV), which takes place upon CP. 1–3 mark the lattice positions most preferred by the promoter. (right) Initial palladium-promoted In2O3(111) surface models for the CP and DI catalysts and their evolution to account for structural modifications occurring upon their use for in the reaction. The energy difference between states for each transition (ΔE, in eV) and the Bader charges of the Pd atoms (qB in |e‒|) are indicated in blue. Atoms are shaded with a progressively lighter color upon increasing distance from the surface towards the bulk
To further elucidate the palladium speciation in the two systems, X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS, Fig. 3b) spectroscopy were measured at the Pd K-edge. The fresh CP material contains isolated Pd2+ species, which possess four O neighbors (Pd–O scattering path, Supplementary Fig. 8 and Supplementary Table 3) in square geometry at a distance of ca. 2.01 Å in the first coordination shell and are embedded in In2O3, as evidenced by the second coordination shell of Pd following the Pd–O-In scattering path. Accordingly, they are strongly bound to the oxide surface or can be in bulk lattice positions if oxygen vacancies are present in their first coordination shell. For the fresh DI catalyst, the O shell around Pd is similar to that in the CP material and in PdO, but the interaction with the In2O3 lattice is much weaker, as hinted by the low coordination number of the Pd–O-In scattering path. Consequently, palladium is likely present as an oxide highly dispersed on In2O3. After 1 h on stream, palladium becomes metallic in both catalysts forming nanostructures smaller than in the intermetallic reference, which contains nanoparticles of ca. 2 nm. In the CP material, these actually comprise clusters of 4–5 atoms, which do not grow significantly upon longer use in the reaction. Due to the similar scattering factors of palladium and indium, their composition can be described as PdxIn4‒x tetramers or Pdx+1In4‒x pentamers. The resemblance of the XANES spectrum of the used CP sample to that of the intermetallic reference suggests that these clusters include both elements. As the Pd K-edge was probed, one atom must be Pd. As Pd–O-In scattering is detected and the palladium amount is low, it is unlikely that the remaining 3–4 atoms are all Pd atoms. Considering that palladium modifies the intrinsic activity of the catalytic site in indium oxide, it is reasonable that palladium replaces one of the In atoms in the In3O5 ensemble. Hence, most likely configurations are Pd2In2 and Pd3In2. These species can only be generated by CP and with low palladium loading, as the Pd–O-In scattering is absent in the EXAFS spectra of the CP catalyst with 3.5 wt.% Pd and the DI solid. Confirmation of the low clustering degree of palladium was attained by time of flight-secondary ions mass spectrometry (TOF-SIMS). Ionic species emitted from the surface of CP 1 h predominantly contained a single Pd atom, and no ions containing > 2 Pd atoms. Relevantly, an equivalent result was attained for CP 16 h, in line with the high resilience of palladium to aggregate in the CP catalyst (Fig. 3c).
To investigate the promoter impact on vacancy formation, H2-TPR at the process pressure (5 MPa, Fig. 3d and Supplementary Fig. 9) was conducted. Surface reduction of In2O3 occurred at much lower temperatures for all promoted catalysts than for the unpromoted oxide (292–313 versus 495 K), highlighting the assistance of palladium in H2 splitting. The signal of the DI sample is more intense, hinting at a higher number of vacancies. For the same solid, an additional signal at 249 K related to palladium reduction was detected, which supports the segregation of the promoter from In2O3. The findings are further corroborated by deconvolution of the XPS O1s signal (Fig. 3e): CP and DI used for 1 h contain a comparable amount of vacancies and after 16 h on stream their concentration remains practically unaltered for the CP catalyst, whereas it increases by 29% for the DI system. The latter sample also features a greater amount of surface hydroxides (25% more than CP), likely caused by dissociated hydrogen spilled from palladium particles to In2O3. A greater abundance of vacancies is also supported by the higher quantity of CO2 desorbed from this catalyst than from CP 1 h (Supplementary Fig. 10). To further unravel changes of the oxide by palladium addition, 115In solid-state nuclear magnetic resonance (NMR) spectra were acquired to selectively probe indium. Notably, the large quadrupolar interaction of the NMR active spin 9/2 In nucleus makes this a challenging analysis, which has never been reported for a bulk indium-based material. Two magnetically inequivalent types of In atoms were detected in the spectrum of cubic In2O3, one causing the central transition with maximum at 280 ppm and the other originating two peaks at 3750 and ‒6250 ppm (Fig. 3f). The asymmetry of these functions is owing to quadrupolar coupling and some disorder, owing to the small size of the oxide particles. The spectra of fresh CP and DI catalysts resemble that of In2O3, but are noisier, mostly because of the diamagnetic character of the PdO present. Additional signals between 5000 and 8000 and −7000 and 11,000 ppm were observed for CP 1 h and DI 1 h owing to the formation of oxygen vacancies, which modify the electronic environment of adjacent In nuclei. Although these are equivalent in all used samples, the particles size of In2O3 in DI 16 h is larger, implying a higher density of vacancies per surface area in this material. The analysis further indicates that the In2O3 structure must be preserved in models of promoted catalysts selected for theoretical modeling.
Theoretical description of palladium sites and mechanism
Different palladium-promoted In2O3(111) surfaces were analyzed to unravel the most likely environment of the noble metal in CP and DI catalysts (Fig. 3g). As palladium and indium cations are atomically mixed upon CP, palladium can replace indium in the In2O3 lattice. However, DFT simulations indicate that palladium is substantially more stable at a surface rather than a bulk position, i.e., segregation is favored towards position 1 in the In3O5(O) ensemble (− 0.87 eV). The replacement of a second or third In atom in the ensemble is unfavored, with islanding energies taking positive values (up to 0.75 eV, Supplementary Table 4). The Pd coordination is square planar, with a Pd–O distance of 2.05 Å, in good agreement with the experimental EXAFS data for the fresh sample. The substitution of Pd atoms in the In3O5(O) protrusion promotes the formation of two oxygen vacancies upon reaction (Supplementary Fig. 11), creating a PdIn2O4 site. Hence, Pd species become more reduced, as evident from the core-level Pd3d spectra and the Bader charges (Fig. 3g). The PdxIn4‒x and Pdx+1In4‒x clusters identified by EXAFS in the used CP sample can be rationalized by the spontaneous migration of Pd atoms from the small depression (pocket) to the adjacent PdIn2O4 sites. Still, an ensemble with 2 Pd atoms likely is more representative of the actual catalyst than one with 3 Pd atoms, as the fresh catalyst contains only few extra-lattice Pd atoms (Supplementary Table 2), the addition of a Pd atom to a palladium dimer is thermoneutral, and configurational entropy favors dispersion43. Notably, the electronic structure of Pd2In4O4 ensembles strongly differs from that of In2O3-supported Pd nanoparticles28 and PdIn intermetallics27,44 based on local densities of states calculations (Supplementary Fig. 12). However, CP models with 0–4 Pd atoms were considered for the reaction mechanism study detailed below.
In the case of DI, Pd atoms are deposited onto an already formed oxide surface preferentially allocating in its depressions. Upon exposure to the reaction mixture, the oxidic palladium species lose the ligands becoming metallic and tend to generate nanoparticles since (i) clustering energies are favorable (−0.4 to −1.0 eV), (ii) Pd atoms on the surface are more abundant than in the CP material (Supplementary Table 2), and (iii) the barriers for the diffusion of Pd atoms between pockets are mild (Supplementary Fig. 13). This phase separation leads to a DI catalyst exhibiting In3O5(O) ensembles and Pd nanoparticles as active sites. The latter easily activate H2 and, by spillover, promote the formation of additional vacancies on the oxide surface (Fig. 3g). Simulated EXAFS spectra of the active structures in CP and DI solids are in reasonable agreement with experimental findings (Fig. 3b). Figure 4 presents the most significant CP and DI catalysts’ models and their features, highlighting the remarkable matching of information derived from experiments and from theory.
Experimental and theoretical description of palladium sites. Most relevant surface models representing the CP and DI catalyst in fresh, activated, and used forms and their main features derived from experiments and theory. The color code for the atoms in the models is the same as described in caption to Fig. 3
Gibbs energy profiles for methanol and CO formation from CO2 were calculated based on the elementary steps found optimal for In2O3. The activation of the CP catalyst starts with H2 dissociation onto two O atoms located nearby the lattice palladium species, followed by transfer of one H atom to form water, which desorbs creating a vacancy. A second vacancy is generated similarly (Supplementary Fig. 11). This can be followed by the creation of active sites with higher palladium nuclearities, i.e., PdxIn2O4 (x > 1), through palladium migration from pockets to PdIn2O4 ensembles. For the DI catalyst, vacancy formation is equivalent, but H2 is split on Pd atoms in the pockets. Notably, H2 dissociation on In2O3 is heterolytic, whereas it is homolytic in the presence of intra- or extra-lattice Pd atoms, which reduce the H2 splitting barrier. CO2 hydrogenation to methanol on CP and DI surface models with 1–4 Pd atoms comprises CO2 adsorption (Supplementary Figs. 14–17), followed by H2 activation and sequential transfer of the two adsorbed H atoms to the C atom. Further H2 activation enables the hydrogenation of one O atom leading to H2COOH. Cleavage of the OH group generates a CH2O species, which is hydrogenated at the C and the O centers, in this order, producing methanol. The latter then desorbs. In the RWGS reaction, CO2 adsorbs on the surface once a third vacancy is created, and is hydrogenated at one O atom. The formed OH group is cleaved and moves to a vacancy, whereas CO desorbs. Based on the span model, CP surfaces with 1–3 Pd atoms and the DI surface with 1 Pd atom are selective towards methanol and to a greater extent than bulk In2O3, whereas, for the CP surface with 4 Pd atoms and the DI surface with 3 Pd atoms, the intrinsic activity of palladium for the RWGS reaction dominates and CO is the favored product. This agrees with the prerequisite of 3 extra-lattice Pd atoms to mediate the water–gas shift reaction on TiO2-supported palladium35. The fact that the boost in methanol selectivity is limited for the CP model with a sole palladium species and remarkable for those with 2–3 Pd atoms supports the higher relevance of Pd2In2O4 and Pd3In2O4 as active sites. Regarding the equilibrated DI catalyst, CO2 hydrogenation can be described as a superposition of RWGS on palladium, i.e., the Pd(111) surface35,45, and methanol and CO formation on neighboring undoped In2O322, the latter being facilitated by hydrogen spillover from adjacent palladium nanoparticles. To account for coverage effects in this complex surface reaction, a microkinetic model was built (vide infra).
Experimental selectivity-promoter nuclearity correlation
Intrigued by the selectivity-nuclearity correlation uncovered by the theoretical studies for CP catalysts, we conducted a complementary experiment to link palladium agglomeration and its impact on methanol selectivity. When the DI catalyst was tested keeping the CO2 conversion at 2–3% (Supplementary Fig. 18), the methanol selectivity remained at ca. 83% for 6 h, followed by a linear decline owing to a progressively greater impact of the RWGS reaction till a value of 58% after 9 h on stream. Characterization of the specimen retrieved at the end of the test confirmed severe sintering of palladium and In2O3 (particle sizes of 12 and 28 nm, respectively). A sample generated by a repeated experiment ending before the onset of the selectivity change still showed an extremely high-palladium dispersion based on STEM-EDX, as expected. This supports that agglomerates of only few Pd atoms are required to start favoring the RWGS reaction and highlights the importance of controlling the promoter nanostructure upon synthesis and preserving it upon reaction.
To tackle the former aspect, we explored the impact of the aging time in CP, as continuous dissolution and re-precipitation occur at this stage, which could lead to palladium segregation from In2O3. For aging times up to 48 h, the methanol selectivity of CP catalysts was unaltered at ca. 80%, hinting that palladium clusters formed still were very small (Supplementary Fig. 19). A sensible difference was observed for the material aged for 200 h (64%), which was found to comprise large palladium particles (up to 10 nm) as well as finely distributed palladium (Supplementary Fig. 20), with its performance originating from the superposition of these two distinct promoter speciations. The In2O3 particle size was 12 nm, corroborating dynamic restructuring upon aging. The longer term behavior of this sample evidenced less-pronounced deactivation than for the DI catalyst and a sustained methanol STY intermediate between that of the CP catalyst attained with an aging of 1 h and the DI solid.
Evaluation of reaction kinetics and catalyst durability
To gain insights into the origin of palladium promotion, the kinetics of CO2 hydrogenation was investigated in detail. An experimental analysis was performed over the CP catalyst, the DI material used for 48 h, and the Pd/TiO2 catalyst already serving as a reference for characterization purposes, under the assumption of first order kinetics, as determined for In2O322. Moreover, a microkinetic model was developed and applied to all CP and DI models35,45. Based on experimental tests at variable temperature (Fig. 5a and Supplementary Fig. 21), the apparent activation energy (Ea) for CO2 hydrogenation to methanol for the CP catalyst is ca. 20% lower than that previously estimated for pure In2O322, whereas that of the RWGS reaction is equivalent. The microkinetic model (Supplementary Table 6) indicates that the apparent activation energy for methanol production is lowered by 56 or 4 kJ mol–1 when 1 or 2 exposed Pd atoms are present, respectively, and the RWGS reaction is hindered by a barrier higher by at least + 80 kJ mol–1 compared with In2O3. This means that competing transformation mainly occurs on unpromoted In3O522 and is effectively blocked by the existence of 1 or 2 exposed Pd atoms at promoted sites. When 3 or more exposed Pd atoms are present, the active site starts behaving like metallic palladium, as the barrier of the RWGS reaction is decreased by –7 kJ mol–1 and methanol production is hindered. For the DI sample, methanol synthesis was facilitated to a comparable extent, but the RWGS reaction was also favored (activation energy decreased by ca. 10%, Fig. 5a), supporting the dual beneficial-detrimental action of palladium nanoparticles. In line with these results, a low activation energy for the parasitic reaction was also determined for Pd/TiO2. Similarly, In2O3 surfaces containing three or more exposed Pd atoms form CO rather than methanol and the Pd(111) surface is associated with the lowest activation barrier for the RWGS reaction of all models (Supplementary Table 6). The experimental measurements additionally indicate that a fourfold higher methanol STY (1.01 gMeOH h−1 gcat−1, CO2 conversion = 9.7%, methanol selectivity = 78%) can be attained over the CP material at a 20-K lower temperature than for In2O3, implying higher productivity with lower energy input. The CP and equilibrated DI catalysts were further tested at variable partial pressures of reactants and products to derive their reaction orders (Table 1, Supplementary Figs. 22 and 23). The values obtained for CO2 and H2 for both catalysts are close to 0, thus being equal and substantially lower than for In2O322, respectively. The reaction order for methanol was barely altered by palladium addition through either synthesis method, whereas the lower reaction order for H2 indicates superior H2 splitting ability over both systems. The negative rate dependence on the water pressure was greatly alleviated for the CP catalyst, whereas no significant change compared to In2O3 was determined for the DI sample. Facilitated water desorption represents a strong contributor to catalyst stability, as water-induced sintering is a major deactivation mechanism for the pure oxide22, whereby we speculate that surface OH groups condensation likely drives the coalescence of In2O3 crystals. This was corroborated by an extended test over the CP material aged for 200 h (Fig. 1b). Owing to the portion of palladium in the lattice, the material did not exhibit pronounced performance deterioration over 80 h (Fig. 1b). As expected, the methanol STY was overall inferior (0.54 versus 0.61 gMeOH h−1 gcat−1) to the CP catalyst aged for 1 h since the methanol selectivity was reduced by the presence of the palladium particles (64% versus 80%, Supplementary Fig. 19). As the methanol production rate was almost constant upon altering the H2:CO2 ratio, the CP catalyst could be practically operated with a sub-stoichiometric H2 feed, thus greatly lowering the size of the recycle stream needed in the process, which will substantially reduce investment and operating costs for reactor, heat exchangers, and compressors46. Applying the optimal conditions identified in the kinetic analysis (H2:CO2 = 4, P = 5 MPa, T = 553 K, and WHSV = 48,000 cm3STP h−1 gcat−1), catalyst durability was assessed in a 500-h test (Fig. 6). The CO2 conversion slightly dropped (from 9.7 to 9.2%) and the methanol selectivity increased (from 75 to 78%) in the first 200 h, whereas remaining unaltered until the end of the run. The average particle size of In2O3 increased by only 2 nm compared with CP 96 h (17 and 15 nm, respectively). Owing to the addition of TiO2 as a diluent in this test, which could not be separated from the used catalyst, changes in porous properties and vacancy density could not be investigated. The methanol STY was 0.96 gMeOH h−1 gcat−1, which is, to our knowledge, the highest sustained productivity reported for CO2-based methanol synthesis over a heterogeneous catalyst (Table 2). Specifically, our catalyst reaches a 15%-higher methanol STY than the best performer reported and at a ca. 60% shorter residence time. The latter enables a reduction in reactor size by 60%, which is a strong gain for a prospective industrial process. Moreover, CP is a synthesis method that has been successfully implemented for catalyst manufacture at the large scale for numerous applications, whereas the polymer-assisted route used for the benchmark appears of limited industrial amenability.
Activation energies and palladium role for palladium-promoted In2O3 catalysts. a Arrhenius plots and calculated apparent activation energies (Ea) for CO2 hydrogenation to methanol (green) and the RWGS reaction (red) over In2O3, CP, DI, and Pd-TiO2 catalysts. Data regression was performed in the region outlined by the dashed lines. Solid lines indicate the thermodynamic equilibria of the two independent reactions. b Schematic representation of the distinct role of palladium in equilibrated CP and DI systems with In2O3 and Pd-TiO2 as references. The thickness of the arrows qualitatively indicates the products formation rates. Reaction conditions: P = 5 MPa, H2:CO2 = 4, and WHSV = 48,000 cm3STP h−1 gcat −1
Table 1 Reaction orders of reactants and products
Extended stability test of the best palladium-promoted In2O3 catalyst. Evolution of the catalytic performance of the CP catalyst during a 500-h test conducted under optimized reaction conditions: T = 553 K, P = 5 MPa, H2:CO2 = 4, and WHSV = 48,000 cm3STP h−1 gcat−1
Table 2 Comparison of heterogeneous catalysts for CO2 hydrogenation to methanol
We successfully produced a palladium-promoted In2O3 catalyst featuring palladium atomically dispersed in the oxide matrix by CP and a reference system comprising palladium nanoparticles by DI, the first showing a stable activity enhancement compared with bulk In2O3 exclusively directed to methanol production and the second displaying a decaying improvement of both methanol and CO formation. Extensive characterization by a battery of techniques, including solid-state 115In NMR, here introduced as a novel technique to elucidate vacancy formation and short-range crystallinity in (palladium-promoted) In2O3 catalysts, and theoretical modeling provided sound understanding of the structure and functioning of the active sites in the two systems. In the CP catalyst, embedding one Pd atom in the In2O3 matrix enables a controlled growth of extra-lattice atoms, leading to the stabilization of low-nuclearity palladium clusters. The number of vacancies does not increase but their electronic properties are strongly modified, permitting improved H2 dissociation. The small size of the palladium cluster is crucial to curtail the RWGS reaction on palladium sites, which is relevant already for clusters of 3 extra-lattice promoter atoms. In the DI catalyst, palladium steadily agglomerates with time on stream forming particles. These aid methanol production by supplying activated hydrogen that fosters the formation of surface vacancies on In2O3 and the hydrogenation of the CO2 therein adsorbed. Still, they consume a significant part of the H2 split by converting CO2 into CO on their own surface (Fig. 5b). These contrasting palladium speciations also explain the profoundly diverse time-on-stream behaviors of the two catalysts. Palladium embedded into indium oxide attracts the scarce palladium in pockets, avoiding its clustering, and facilitates water desorption minimizing In2O3 sintering, whereas palladium deposited onto In2O3 easily agglomerates and leads to catalyst sintering by spilling too much activated H2 to the oxide, causing over-reduction, and by excessively retaining water, provoking crystals coalescence. Hence, anchoring of palladium clusters to the In2O3 lattice is an unprecedented approach to effectively enhance activity and selectivity, and more importantly, to maintain the promotional effect in the long term, an aspect fully disregarded in previous studies on palladium-promoted In2O3 systems. In addition, low-nuclearity palladium clusters grant the possibility to operate the catalyst at reduced temperature, with hydrogen-lean feeds, and in the presence of greater amounts of water, altogether implying strong economic and ecologic benefits at a process level. Overall, our strategy to engineer promotion at the atomic scale marks a step ahead toward green methanol production and holds great general potential for tailoring new or existing promoted systems in current and emerging applications of heterogeneous catalysis.
To obtain palladium-promoted (0.25–10 wt.% Pd) In2O3 catalysts by co-precipitation, In(NO3)3·xH2O (3.51, 3.50, 3.48, 3.47, 3.46, 3.45, 3.40, or 3.20 g, Sigma-Aldrich, 99.99%, x = 6.9) and Pd(NO3)2·xH2O (4.63, 11.6, 23.15, 34.8, 46.30 69.6, 162.5, or 463 mg, Sigma-Aldrich, > 99.99% metals basis, x = 5.5) were dissolved in deionized water (50 cm3) in a round-bottomed flask (250 cm3). Na2CO3 (10.0 g) was dissolved in deionized water (100 cm3) and added dropwise (ca. 25 cm3, 3 cm3 min−1) to the metals solution under stirring at room temperature to reach a pH value of 9.2. After aging the resulting slurry for 1 h, deionized water was added (50 cm3). Then, the precipitate was recovered by high-pressure filtration, washed with deionized water (three times, 500 cm3 each time), and dried in a vacuum oven (1.5 kPa, 323 K, 1.5 h). To produce palladium-containing (0.25–3.5 wt.% Pd) In2O3 catalysts by dry impregnation, a round-bottomed flask (25 cm3) was loaded with Pd(NO3)2·xH2O (7.8, 23.5, or 109.4 mg) and deionized water (0.39 g). After adding In2O3 (1.00 g) and five stainless-steel spheres (radius = 3 mm), the flask was rotated (ca. 45 r.p.m.) using a Büchi R-114 rotary evaporator at room temperature and pressure. After 12 h, the pressure was lowered to 2 kPa and the temperature raised to 333 K for 1 h to allow the evaporation of the solvent. Catalysts were calcined for 3 h at 573 K (2 K min−1) in static air prior to characterization or catalytic testing. The syntheses of all other materials reported in this study are detailed in the Supplementary Methods.
The metal content of the catalysts was determined by inductively coupled plasma optical emission spectrometry. When Na2CO3 was used as the precipitating agent, the absence of sodium in calcined materials was confirmed by X-ray fluorescence spectroscopy as well as by STEM-EDX and XPS for selected materials. The water amount in In and Pd precursors was derived from thermogravimetric analysis. Porous and structural properties of the catalysts were assessed by N2 sorption and powder XRD. Palladium dispersion, speciation and coordination geometry, and reducibility were accessed by volumetric chemisorption of H2 and CO, CO-DRIFTS, STEM-EDX, and HRTEM, X-ray absorption spectroscopy (XAS), TOF-SIMS, and by TPR. Surface concentration and oxidation states of C, O, Pd, and In were determined by XPS. Oxygen vacancies and CO2 adsorption capacity were characterized by XPS, 115In NMR, and TPR and TPD methods. Details to the characterization techniques are provided in the Supplementary Methods and in Supplementary Fig. 24.
CO2 hydrogenation to methanol was conducted in a high-pressure continuous-flow fixed-bed reactor setup comprising mass flow controllers to feed gases (Bronkhorst, El-Flow F-201CV), pneumatic valves, a reactor with an inner diameter of 2.2 mm housed in an electrically-heated aluminum brass furnace, a pressure transducer, a burst plate calibrated at 6 MPa, a syringe pump to feed liquids (Chemyx Nexus 6000), an online gas chromatograph (Agilent 7890 A equipped with Agilent DB-1 and GS-GasPro columns), and a computer control by a custom protocol within the LabView software. In a typical test, the reactor was loaded with 100 mg of catalyst with a particle size of 100–125 μm, which was held in place by a bed of quartz wool and heated from ambient temperature to 553 K (5 K min−1) at 5 MPa under a He flow (20 cm3 min−1). After 3 h, the gas flow was switched to the reactant mixture (40 cm3 min−1) comprising H2 and CO2 (Messer, 99.997% and 99.999%, respectively) in a molar ratio of 4:1. For kinetic tests, temperature (473–653 K) and inlet partial pressures of reagents (pH2 = 3.5–4.5 MPa, pCO2 = 0.5–1.5 MPa) and products (pMeOH and pH2O = 0.05–0.25 MPa) were varied at a doubled WHSV attained by loading only 50 mg of catalyst diluted in 50 mg of TiO2 (100–125 μm, Sigma-Aldrich, > 99.9%). Water (ABCR-Chemicals, HPLC grade) or methanol (Sigma-Aldrich, > 99.9%, anhydrous) were fed to the reactor inlet by a high-pressure syringe pump (Chemyx, Nexus 6000). The effluent stream was sampled every 12 min and analyzed by online gas chromatography. Data evaluation procedures are presented in the Supplementary Methods and details to all catalytic tests are compiled in Supplementary Table 7.
DFT simulations were performed using the Vienna Ab initio Simulation Package (VASP)47,48 using the Perdew–Burke–Ernzerhof density functional49. Valence electrons were described with a plane-wave energy cutoff of 500 eV, whereas core electrons were represented by projector augmented-wave pseudopotentials50. The In2O3(111) surface was modeled as a p(1 × 1) slab containing five O−In−O trilayers and optimized using a Γ-centered (3 × 3 × 1) k-point mesh. In each surface, the three outermost O−In−O trilayers were allowed to relax, whereas the two bottommost layers were fixed in their bulk positions. The initial structure of the CP catalyst was modeled by substituting an In atom in the In3O5(O) ensemble by a Pd atom. Three further structures include the addition of 1–3 Pd atoms bonded to the Pd atom in the oxide lattice. For the DI catalyst, 1–3 Pd atoms were added in the concave zone of the In2O3(111) surface. For both types of catalysts, O vacancy formation and palladium migration were explored. Energy profiles for CO2 hydrogenation to methanol and CO were calculated using the most feasible paths identified in our previous study22 on all CP and DI models, including the water-assisted step in the RWGS. To reduce the number of steps, the adsorption of the third H2 molecule was moved after H2COOH formation. The reaction network occurring on the Pd(111) surface, representing the palladium nanoparticles in the equilibrated DI sample, was also taken from previous studies35,45. The Gibbs free energies of the adsorbed species considered only vibrational contributions. For gas phase, translational and rotational terms were also included and obtained from Gaussian 0951. The reaction network was interrogated through a microkinetic model based on a differential reactor under steady-state conditions. The kinetic parameters were modeled from transition state theory, whereas adsorption processes followed the Hertz–Knudsen equation46. The resulting system of equations was resolved in Maple using a floating-point precision of 96 digits. The temperatures and pressures of operation were chosen to match the experimental ones. For relevant systems, the XPS core-level Pd3d shifts were computed in VASP, using metallic Pd as reference52, and were compared with the Bader charges (Supplementary Table 9). EXAFS spectra were obtained from XAS simulations attained using feff853 after processing with Athena54 in the same way as for the experimental data and subsequent normalization.
The authors declare that the data supporting the findings of this study are available within the article and its Supplementary Information file. All DFT data and the corresponding structures are accessible at the ioChem-BD database55 following the link: https://doi.org/10.19061/iochem-bd-1-106. All other relevant source data are available from the corresponding author upon reasonable request.
Bartholomew, C. H. & Farrauto, R. J. Fundamentals of Industrial Catalytic Processes (Wiley, Hoboken, 2006).
Hutchings, G. J. Promotion in heterogeneous catalysis: a topic requiring a new approach? Catal. Lett. 75, 1–12 (2001).
van Santen, R. A. Chemical basis of metal catalyst promotion. Surf. Sci. 251-252, 6–11 (1991).
Sheng, H. & Lobo, R. F. Iron-promotion of silica-supported copper catalysts for furfural hydrodeoxygenation. ChemCatChem 8, 3402–3408 (2016).
Kondratenko, E. V., Mul, G., Baltrusaitis, J., Larrazábal, G. O. & Pérez-Ramírez, J. Status and perspectives of CO2 conversion into fuels and chemicals by catalytic, photocatalytic and electrocatalytic processes. Energ. Environ. Sci. 6, 3112–3135 (2013).
Olah, G. A. Beyond oil and gas: the methanol economy. Angew. Chem. Int. Ed. 44, 2636–2639 (2005).
Ali, K. A., Abdullah, A. Z. & Mohamed, A. R. Recent development in catalytic technologies for methanol synthesis from renewable sources: a critical review. Renew. Sust. Energ. Rev. 44, 508–518 (2015).
Behr, A. Methanol: the basic chemical and energy feedstock of the future. Angew. Chem. Int. Ed. 53, 12674–12674 (2014).
Álvarez, A. et al. Challenges in the greener production of formates/formic acid, methanol, and DME by heterogeneously catalyzed CO2 hydrogenation processes. Chem. Rev. 117, 9804–9838 (2017).
Wang, J. et al. A highly selective and stable ZnO-ZrO2 solid solution catalyst for CO2 hydrogenation to methanol. Sci. Adv. 3, e1701290 (2017).
Vogt, C. et al. Unravelling structure sensitivity in CO2 hydrogenation over nickel. Nat. Catal. 1, 127–134 (2018).
Wang, L. et al. Incorporating nitrogen atoms into cobalt nanosheets as a strategy to boost catalytic activity toward CO2 hydrogenation. Nat. Energy 2, 869–876 (2017).
Behrens, M. Promoting the synthesis of methanol: understanding the requirements for an industrial catalyst for the conversion of CO2. Angew. Chem. Int. Ed. 55, 14906–14908 (2016).
Xu, J. et al. Methanol synthesis from CO2 and H2 over Pd/ZnO/Al2O3: catalyst structure dependence of methanol selectivity. Appl. Catal. A 514, 51–59 (2016).
Yang, H. B. et al. Atomically dispersed Ni(I) as the active site for electrochemical CO2 reduction. Nat. Energy 3, 140–147 (2018).
Loiland, J. A., Wulfers, M. J., Marinkovic, N. S. & Lobo, R. F. Fe/γ-Al2O3 and Fe-K/γ-Al2O3 as reverse water-gas shift catalysts. Catal. Sci. Technol. 6, 5267–5279 (2016).
Martin, O. et al. Indium oxide as a superior catalyst for methanol synthesis by CO2 hydrogenation. Angew. Chem. Int. Ed. 55, 6261–6265 (2016).
Gao, P. et al. Direct conversion of CO2 into liquid fuels with high selectivity over a bifunctional catalyst. Nat. Chem. 9, 1019–1024 (2017).
Zhang, S. et al. Solid frustrated-Lewis-pair catalysts constructed by regulations on surface defects of porous nanorods of CeO2. Nat. Commun. 8, 15266 (2017).
Ghuman, K. K. et al. Photoexcited surface frustrated Lewis pairs for heterogeneous photocatalytic CO2 reduction. J. Am. Chem. Soc. 138, 1206–1214 (2016).
Albani, D. et al. Semihydrogenation of acetylene on indium oxide: proposed single-ensemble catalysis. Angew. Chem. Int. Ed. 56, 10755–10760 (2017).
Frei, M. S. et al. Mechanism and microkinetics of methanol synthesis via CO2 hydrogenation on indium oxide. J. Catal. 361, 313–321 (2018).
Ye, J., Liu, C.-j, Mei, D. & Ge, Q. Methanol synthesis from CO2 hydrogenation over a Pd4/In2O3 model catalyst: a combined DFT and kinetic study. J. Catal. 317, 44–53 (2014).
Lorenz, H. et al. Pd–In2O3 interaction due to reduction in hydrogen: consequences for methanol steam reforming. Appl. Catal. A 374, 180–188 (2010).
Rameshan, C. et al. Impregnated and co-precipitated Pd–Ga2O3, Pd–In2O3 and Pd–Ga2O3–In2O3 catalysts: influence of the microstructure on the CO2 selectivity in methanol steam reforming. Catal. Lett. 148, 3062–3071 (2018).
García-Trenco, A. et al. PdIn intermetallic nanoparticles for the hydrogenation of CO2 to methanol. Appl. Catal. B 220, 9–18 (2018).
Ye, J., Ge, Q. & Liu, C.-J. Effect of PdIn bimetallic particle formation on CO2 reduction over the PdIn/SiO2 catalyst. Chem. Eng. Sci. 135, 193–201 (2015).
Rui, N. et al. CO2 hydrogenation to methanol over Pd/In2O3: effects of Pd and oxygen vacancy. Appl. Catal. B 218, 488–497 (2017).
Arena, F. et al. Solid-state interactions, adsorption sites and functionality of Cu-ZnO/ZrO2 catalysts in the CO2 hydrogenation to CH3OH. Appl. Catal. A 350, 16–23 (2008).
Erdöhelyi, A., Pásztor, M. & Solymosi, F. Catalytic hydrogenation of CO2 over supported palladium. J. Catal. 98, 166–177 (1986).
Pettigrew, D. J., Trimm, D. L. & Cant, N. W. The effects of rare earth oxides on the reverse water-gas shift reaction on palladium/alumina. Catal. Lett. 28, 313–319 (1994).
Vilé, G. et al. A stable single-site palladium catalyst for hydrogenations. Angew. Chem. Int. Ed. 54, 11265–11269 (2015).
Mitchell, S., Vorobyeva, E. & Pérez-Ramírez, J. The multifaceted reactivity of single‐atom heterogeneous catalysts. Angew. Chem. Int. Ed. 130, 15538–15552 (2018).
Gao, G., Jiao, Y., Waclawik, E. R. & Du, A. Single atom (Pd/Pt) supported on graphitic carbon nitride as an efficient photocatalyst for visible-light reduction of carbon dioxide. J. Am. Chem. Soc. 138, 6292–6297 (2016).
Li, Q., García-Muelas, R. & López, N. Microkinetics of alcohol reforming for H2 production from a fair density functional theory database. Nat. Commun. 9, 526 (2018).
Li, H. et al. Synergetic interaction between neighbouring platinum monomers in CO2 hydrogenation. Nat. Nanotechnol. 13, 411–417 (2018).
Nørskov, J. K., Bligaard, T., Rossmeisl, J. & Christensen, C. H. Towards the computational design of solid catalysts. Nat. Chem. 1, 37–46 (2009).
Medford, A. J. et al. From the Sabatier principle to a predictive theory of transition-metal heterogeneous catalysis. J. Catal. 328, 36–42 (2015).
Jimenez-Izal, E. & Alexandrova, A. N. Computational design of clusters for catalysis. Annu. Rev. Phys. Chem. 69, 377–400 (2018).
Tauster, S. Strong metal-support interactions. Acc. Chem. Res. 20, 389–394 (1987).
Abbet, S. et al. Identification of defect sites on MgO(100) thin films by decoration with Pd atoms and studying CO adsorption properties. J. Am. Chem. Soc. 123, 6172–6178 (2001).
Bergeret, G. & Gallezot, P. in Handbook of Heterogeneous Catalysis Vol. 2. (eds Ertl, G., Knözinger, H., Schüth, F. & Weitkamp, J.) 738–765 (Wiley-VCH, Weinheim, 2008).
Capdevila-Cortada, M. & López, N. Entropic contributions enhance polarity compensation for CeO2(100) surfaces. Nat. Mater. 16, 328–334 (2017).
Snider, J. L. et al. Revealing the synergy between oxide and alloy phases on the performance of bimetallic In–Pd catalysts for CO2 hydrogenation to methanol. ACS Catal. 9, 3399–3412 (2019).
Garcı́a-Muelas, R., Li, Q. & López, N. Density functional theory comparison of methanol decomposition and reverse reactions on metal surfaces. ACS Catal. 5, 1027–1036 (2015).
Chorkendorff, I. & Niemantsverdriet, J. W. in Concepts of Modern Catalysis and Kinetics Vol. 1 (eds Chorkendorff, I. & Niemantsverdriet) 79–128 (Wiley-VCH, Weinheim, 2003)
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558–561 (1993).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865–3868 (1996).
Blöchl, P. E. Projector augmented-wave method, Phys. Rev. B 50, 17953–17979 (1994).
Frisch, M. J. et al. Gaussian 09, Revision A.02. (Gaussian, Inc., Wallingford CT, 2016).
Köhler, L. & Kresse, G. Density functional study of CO on Rh(111). Phys. Rev. B 70, 165405 (2004).
Newville, M. EXAFS analysis using FEFF and FEFFIT. J. Synchrotron Radiat. 8, 96–100 (2001).
Ravel, B. & Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat. 12, 537–541 (2005).
Álvarez-Moreno, M. et al. Managing the computational chemistry big data problem: the ioChem-BD platform. J. Chem. Inf. Model. 55, 99–103 (2015).
Liao, F. et al. A promising low pressure methanol synthesis route from CO2 hydrogenation over Pd@Zn core-shell catalysts. Green. Chem. 19, 270–280 (2017).
Li, M. M. J. et al. CO2 hydrogenation to methanol over catalysts derived from single cationic layer CuZnGa LDH precursors. ACS Catal. 8, 4390–4401 (2018).
Duyar, M. S. et al. A highly active molybdenum phosphide catalyst for methanol synthesis from CO and CO2. Angew. Chem. Int. Ed. 57, 15045–15050 (2018).
Total Research & Technology Feluy is thanked for sponsoring this project. R.G.-M. thanks the AGAUR 2017 SGR 90 project for financial support. BSC-RES is acknowledged for generous computational resources. S. Mitchell is thanked for the electron microscopy measurements, L. Piveteau and R. Verel for help with the acquisition and analysis of NMR spectra, the Scientific Center for Optical and Electron Microscopy (ScopeM) at the ETH Zurich and the NMR-Service of the Institute for Chemical and Bioengineering at the ETH Zurich, for the use of their facilities, and Marçal Capdevila-Cortada for useful discussions.
Institute for Chemical and Bioengineering, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 1, 8093, Zurich, Switzerland
Matthias S. Frei, Cecilia Mondelli, Klara S. Kley, Begoña Puértolas & Javier Pérez-Ramírez
Institute of Chemical Research of Catalonia (ICIQ), The Barcelona Institute of Science and Technology, Av. Països Catalans 16, 43007, Tarragona, Spain
Rodrigo García-Muelas & Núria López
Paul Scherrer Institute, Forschungsstrasse 111, 5232, Villigen, Switzerland
Olga V. Safonova
Total Research & Technology Feluy, Zone Industrielle Feluy C, 7181, Seneffe, Belgium
Joseph A. Stewart & Daniel Curulla Ferré
J.P.-R. and C.M. conceived and coordinated all stages of this research. M.S.F., K.S.K. and B.P. prepared and characterized the catalysts and conducted the catalytic tests. O.V.S. coordinated acquisition and evaluation of X-ray absorption spectroscopy data. R.G.-M. and N.L. conducted the computational studies. J.A.S. and D.C.F. contributed in setting industrial targets for the experimental program. All authors contributed to the writing of the manuscript.
Correspondence to Javier Pérez-Ramírez.
The authors declare no competing interests.
Peer review information: Nature Communications thanks Yuhan Sun, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.