Nanometric Cu-ZnO Particles Supported on N-Doped Graphitic Carbon as Catalysts for the Selective CO2 Hydrogenation to Methanol

The quest for efficient catalysts based on abundant elements that can promote the selective CO2 hydrogenation to green methanol still continues. Most of the reported catalysts are based on Cu/ZnO supported in inorganic oxides, with not much progress with respect to the benchmark Cu/ZnO/Al2O3 catalyst. The use of carbon supports for Cu/ZnO particles is much less explored in spite of the favorable strong metal support interaction that these doped carbons can establish. This manuscript reports the preparation of a series of Cu-ZnO@(N)C samples consisting of Cu/ZnO particles embedded within a N-doped graphitic carbon with a wide range of Cu/Zn atomic ratio. The preparation procedure relies on the transformation of chitosan, a biomass waste, into N-doped graphitic carbon by pyrolysis, which establishes a strong interaction with Cu nanoparticles (NPs) formed simultaneously by Cu2+ salt reduction during the graphitization. Zn2+ ions are subsequently added to the Cu–graphene material by impregnation. All the Cu/ZnO@(N)C samples promote methanol formation in the CO2 hydrogenation at temperatures from 200 to 300 °C, with the temperature increasing CO2 conversion and decreasing methanol selectivity. The best performing Cu-ZnO@(N)C sample achieves at 300 °C a CO2 conversion of 23% and a methanol selectivity of 21% that is among the highest reported, particularly for a carbon-based support. DFT calculations indicate the role of pyridinic N doping atoms stabilizing the Cu/ZnO NPs and supporting the formate pathway as the most likely reaction mechanism.


Introduction
Methanol is among the most valuable products that can derive hydrogenation from CO 2 [1][2][3].Being in the liquid state at ambient conditions, methanol has other important advantages compared to alternative products formed in CO 2 hydrogenation, including water solubility, non-corrosiveness, high volumetric energy content, and easy transformation into gasoline [4,5] and aromatics [6,7], among other chemicals [8,9].Although methanol is currently produced on a large multi-ton scale, and there was an estimated 100 millions of metric tons produced in 2020, there is still the possibility to considerably increase methanol production [10], particularly if application of methanol as a fuel or hydrogen carrier is finally implemented [11][12][13].In any case, market forecasts indicate that methanol production will at least double by 2030 [12].
Since formic acid has a low H 2 content and is corrosive, another important advantage of methanol vs. formic acid is its much higher mass and volumetric energy content, about four times higher than that of formic acid [14,15].Methanol can be used directly as transportation fuel in combustion engines [15] and it can also be directly used as liquid fuel in proton exchange membrane fuel cells [16].Methanol is also considered a liquid organic hydrogen carrier [17,18], with a H 2 storage capacity of about 19 wt.% [18,19].Although methanol reforming will still emit unwanted CO 2 , the cycle could have a zero CO 2 footprint if methanol is in turn formed from CO 2 [20].Equation (1) corresponds to methanol synthesis from CO 2 hydrogenation.CO As indicated in Equation ( 1), the partial CO 2 hydrogenation to methanol is an exothermic reaction, with the equilibrium towards methanol formation being more favorable at low temperatures and high pressures [21,22].However, the slow reaction kinetics determine that in order to achieve measurable reaction rates, heating of the system and the use of suitable catalysts are required to form methanol.
At high temperatures, CO 2 conversion can be limited by equilibrium composition.In addition, besides methanol, CO appears generally as a competing product (Equation ( 2)).Typical CO 2 hydrogenation mixtures are composed of methanol and CO in various proportions accompanied by lesser amounts of methane, thus decreasing methanol selectivity.Formation of CO and CH 4 prevails in CO 2 hydrogenation at high temperatures.
Due to these constraints, CO 2 hydrogenation to methanol is carried out at temperatures in a range from 250 to 300 • C and high pressures, for which a compromise between thermodynamics and kinetics requirements can be reached.Also, to overcome this thermodynamic limitation, photocatalytic CO 2 reduction to methanol is gaining importance, including the use of g-C 3 N 4 [22][23][24].
After the discovery by BASF of copper chromite (Cu-CrO 3 ) as a catalyst for CO 2 hydrogenation to methanol [25], the most widely used catalyst is copper and zinc oxide supported on alumina (Cu-ZnO/Al 2 O 3 ), which is considered as the current benchmark catalyst for the process [25].The typical atomic Cu/Zn proportion is 2:1 and loading on Al 2 O 3 can be over 30 wt.% [26][27][28][29][30].Although the catalyst is generally denoted as Cu-ZnO, in situ studies suggest that Cu-ZnO is a precursor of the active species formed under the reaction conditions by Cu restructuration, ZnO chemical reduction to Zn metal, and Cu-Zn alloying [31,32].Thus, even though under ambient conditions Zn is present as oxide and Cu and Zn are in different phases, the as-prepared Cu-ZnO/Al 2 O 3 should be considered a precursor of the active sites that have been proposed to be Cu NPs decorated by Zn atoms [32].Besides alumina, zirconia in different crystallographic phases is considered also as a suitable support [33][34][35].
In spite of the fact that Cu-ZnO/Al 2 O 3 was reported many years ago and considering the current intense research on catalytic CO 2 hydrogenation to methanol [23,33,36], progress in the development of alternative, more efficient, catalysts that could promote methanol synthesis from CO 2 at lower temperatures is still unsatisfactory [33] or based on less abundant elements [37][38][39].Among the non-containing-Cu catalysts, indium oxide can promote CO 2 hydrogenation to methanol with a high selectivity [37][38][39].Besides facets of In 2 O 3 , oxygen vacancies are the active sites [40].Pd doping increases the activity of In 2 O 3 without much negative influence on the selectivity [41].However, indium is considered as a scarce element, particularly compared with abundant Cu and Zn, and it is included in the list of critical raw materials to be avoided [42].
In a series of articles, we have been showing that N-doped graphitic carbons are suitable supports to develop highly selective Fe-Co catalysts for various CO 2 hydrogenation reactions, including the Sabatier reaction [43], the reverse water gas shift (RWGS) (Equation ( 2)) [44], and for the formation of C 2+ products [45].In this context, it is also of interest to expand the use of N-doped graphitic carbon matrices as supports for Cu-ZnO NPs [Cu-ZnO@(N)C] and to determine their catalytic activity for methanol synthesis under operation conditions compatible with the thermodynamic limitations of the process.In this way, the materials here prepared based on chitosan derived from biomass wastes represent a clear example of waste valorization, applying circular economy principles.The transformation of chitosan into a N-doped graphitic carbon, supporting metal NPs to be used as catalysts, considerably increases the value of the biomass waste.
In comparison to inorganic oxides, carbonaceous supports for Cu-ZnO have been significantly less studied, with most of these studies being limited to carbon nanotubes for methanol synthesis from CO/H 2 [46] and steam reforming of methanol [47].In one of the few precedents on the use of carbon supports, Cu-CuO and ZnO were formed in a porous carbon FDU-15 obtained by pyrolysis of resol, observing a similar performance for CO 2 hydrogenation to methanol than the benchmark BASF catalyst [48].Therefore, the catalytic activity of Cu-ZnO@(N)C samples still appears to be worth exploring.As it will be shown below, the experimental data support that Cu-ZnO@(N)C is an efficient, selective, and stable catalyst for partial CO 2 hydrogenation to methanol, resulting in a notable methanol productivity of 83 g CH3OH kg catalyst −1 h −1 .

Materials and Methods
2.1.Synthesis of Samples Cu@(N)C and Cu-ZnO@(N)C Cu@(N)C (sample 1) and Cu-ZnO@(N)C (samples 2-6) were obtained by dissolving 1 g chitosan with 625 µL acetic acid in 50 mL Milli-Q water.After chitosan dissolved completely, the solution was introduced dropwise, with a syringe (0.8 mm diameter needle), in an aqueous solution of sodium hydroxide (500 mL, 2 M).The hydrogel microspheres were formed immediately and immersed in NaOH solution for 2 h and then profusely washed with distilled water to pH 7. Afterwards, the resulting hydrogel microspheres were washed by a series of ethanol/water baths with an increasing concentration of ethanol (10,30,50,70,90, 100 vol.%, respectively) for 15 min in each and immersed in 100 mL Cu(OAc) 2 -ethanol solution with different concentrations, as indicated in Table S1, for 2 days with a slow stirring, then washed with anhydrous ethanol, and subsequently dried by supercritical CO 2 .Drying using supercritical CO 2 ensures high porosity and large surface area of the aerogel microspheres in comparison with alternative drying procedures [44].The resulting aerogel microspheres were pyrolyzed under Ar flow (200 mL/min), increasing the temperature at a rate of 2 • C/min up to 200 • C for 2 h and then to 900 • C for 2 h.The as-prepared samples did not exhibit pyrophoric properties.The resulting Cu@(N)C was ground into powder and immersed in 30 mL Zn(OAc) 2 -ethanol solution with different concentrations for 2 days with slow stirring.After removal of ethanol at 60 • C overnight, the Zn 2+ -containing Cu@(N)C was heated at a rate of 2 • C/min up to 200 • C for 2 h to obtain the final Cu-ZnO@(N)C.

Preparation of Cu-ZnO/Al 2 O 3
Copper and zinc oxide supported on alumina (Cu-Zn/Al 2 O 3 ) were prepared by impregnation method in two or one steps.Al 2 O 3 powder (270 mg) was dispersed into an ethanol solution (20 mL) of Cu(OAc) 2 (254.5 mg) and the suspension was stirred at room temperature overnight until the solvent evaporated.Then, the obtained powder was pyrolyzed under Ar flow (200 mL/min), increasing the temperature at a rate of 2 • C/min up to 200 • C for 2 h and then to 900 • C for 2 h.After cooling at room temperature, the resulting Cu/Al 2 O 3 powder was impregnated in ethanol solution (20 mL) with Zn(OAc) 2 (10 mg) and the suspension was stirred at room temperature overnight until the solvent evaporated.The resulting powder was annealed at a rate of 2 • C/min up to 200 • C for 2 h under Ar flow (200 mL/min).Another analogous Cu-Zn/Al 2 O 3 sample in which Cuimpregnated Al 2 O 3 was not submitted to pyrolysis was also prepared (Cu-Zn/Al 2 O 3 -wp, wp meaning without pyrolysis).A third sample was prepared as Cu-Zn/Al 2 O 3-wp except that the impregnation of Cu(OAc) 2 and Zn(OAc) 2 was carried out with the same amounts and times but in a single step (Cu-Zn/Al 2 O 3 -imp).

Sample Characterization
X-ray diffraction (XRD) patterns were obtained in a Philips XPert diffractometer (Košice, Slovakia) (40 kV and 45 mA) equipped with a graphite monochromator employing Ni-filtered Cu Kα radiation (1.541178 Å).Raman spectra were collected with a Horiba Jobin Yvon-Labram HR UV-visible-NIR (Kyoto, Japan) (200-1600 nm).Raman microscope spectrometer model had a 514 nm laser.The carbon and nitrogen content of the samples was determined by combustion chemical analysis by using a CHNS FISONS elemental analyzer (Spain).The chemical analysis was determined by ICP-OES (iCAP 7400, Thermo Scientific, Waltham, MA, USA) from the mother liquor after digesting the Cu-ZnO@(N)C samples in aqua regia at 60 • C for one day.High-resolution field emission scanning electron microscopy (HR-FESEM) images were acquired by using a Zeiss GeminiSEM500 apparatus (Jena, Germany).High-resolution transmission electron microscope (HR-TEM) images were recorded in a JEOL JEM 2100F (Košice, Slovakia) under an accelerating voltage of 200 kV coupled with an X-Max 80 energy-dispersive X-ray detector (Oxford instruments, Abingdon, UK).This HR-TEM is equipped with dark-field and high-angle field image detectors that facilitate the observation of phase contrast with different atomic numbers.Samples for measurement were prepared by dropping a few drops of the suspended material in ethanol or dichloromethane on a carbon-coated nickel grid and drying at room temperature overnight.The average metal particle size was determined by measuring the diameter of a statistically relevant number of metal NPs in dark field TEM images using the program J-image.The results of this measurement are presented in the furthest-right column of Table 1.
Table 1.List of samples under study and their main analytical and physicochemical parameters.

Sample
No.

Computational Models and Methods
Periodic DFT calculations were conducted using the Cambridge Serial Total Energy Package (CASTEP) module with the exchange-correlation functional described by Perdew-Burke-Ernzerhof (revised version for solids) within the generalized gradient approximation (GGA-PBEsol).Tkatchenko and Scheffler (TS) dispersion corrections, scheme, were incorporated along with the exchange and correlation functional to improve the structural and vibrational properties.Slabs were separated by 15 Å along [001], perpendicular to the surface.A self-consistent field method (tolerance 1.0 × 10 −7 eV/atom) was employed in conjunction with plane-wave basis sets with a cutoff energy of 500 eV in reciprocal space.All structures were geometry optimized until energy was converged to 1.0 × 10 −6 eV/atom, maximum force to 0.025 eV/Å, and maximum displacement to 5.0 × 10 −3 Å.
The transition states of every elementary reaction in the CO 2 hydrogenation process were obtained by the complete linear synchronous (LST) and quadratic synchronous transit (QST) methods.The adsorption energy of species over the ZnCu/N-C models was calculated as E b = E total − E A − E sur, where E total represents the total energy of the catalytic surface with the adsorbed molecule and E A and E sur are the energies of isolated adsorbate molecule and the clean surface, respectively.The energy of an isolated molecule (E A ) is computed by placing it in the same lattice box (about 15 × 15 × 18 Å 3 ).The activation energy barriers (E a ) of every step of the elementary reactions are defined as E a = E TS − E IS , where E IS and E TS are the total energy of the initial state (IS) and transition state (TS), respectively.

Catalyst Testing
Samples 1-6 were tested for CO 2 hydrogenation in a PID Eng&Tech setup (Microactivity Spain) equipped with a stainless steel (316 SS) fixed-bed tube microreactor (Autoclave Engineers) featured with an inner K-type thermocouple in contact with the catalyst packed bed.Flow rates of reagent gases (H 2 5.0-Linde and CO 2 4.5-Linde) were controlled with two calibrated mass flow controllers (EL-FLOW Select-Bronkhorst, Nijverheidsstraat, The Netherlands).At least two independent catalytic tests were performed for each sample.
For each test, 40 mg of catalyst powder was introduced in the reactor, air was removed by flushing the system at room temperature for 15 min with ten times the flow rates used during the catalytic tests, followed by 10 min at the flow rates used during the experiments (3.2 mL•min −1 H 2 , 0.8 mL•min −1 CO 2 ).The reactor was afterwards slowly pressurized at 40 bars without changing the flow rates of the gas reagents.Two successive GC injections were performed with the gases passing through the reactor at room temperature to check the stability of the gas-phase composition inside the pressurized reactor.Four reaction temperatures from 150 to 300 • C with steps of 50 • C were investigated.For each temperature, GC injection was performed 90 min after the stabilization of the temperature to achieve a steady-state regime of the reactor setup.
Gas samples from the reactor output were passed through a transfer line kept thermostated at 110 • C to be analyzed using an Agilent 7890A gas chromatography instrument with TCD detection.Product separation was performed using capillary PLOT columns: molecular sieve (RT-Molsieve 5A-Restek, for CO, CO 2 , and light hydrocarbons) and divinylbenzene (SupelQ-Supelco, for methanol) with H 2 carrier gas (linear velocities between 30 and 41 cm/s).Each gas sample was injected through a remote-controlled 6-way valve (A4C6WE-Vici, thermostated at the same temperature as the transfer line), with a 90 µL injection loop.The system is operated under continuous flow and some of the experiments are ran for several tens of hours to ensure that the products formed do not derive from the catalyst.

Sample Preparation and Characterization
Initial attempts to prepare Cu-Zn@(N)C were carried out following similar preparation procedures as those described in related precedents on Fe-Co@(N)C [43,44] consisting of the pyrolysis of chitosan powders embedding simultaneously Co and Fe metal ions.During the pyrolysis, chitosan, a polysaccharide of glucosamine, becomes converted in turbostratic N-doped graphitic carbon that can be completely exfoliated to single-or few-layer defective graphene upon sonication [49].The strong reductive conditions of the hot carbon material during the pyrolysis under inert atmosphere are responsible for the reduction of metal ions to the metallic state [50].In the present case, chitosan solutions were impregnated with a mixture of Cu(OAc) 2 and Zn(OAc) 2 in various molar ratios.However, due to the low boiling point of Zn metal and the flushing Ar flow, these attempts with Cu-Zn@(N)C were met with failure.Complete Zn evaporation occurred under the conditions of pyrolysis, with only Cu on the (N)C support remaining after pyrolysis.
To overcome this limitation, ZnO was incorporated after the pyrolysis of Cu 2+ salts adsorbed on chitosan.While it would be possible also to simultaneously incorporate Cu and Zn after chitosan pyrolysis and formation of the N-doped graphitic carbon, it was anticipated that the strong interaction between the Cu metal NPs and the graphene sheets of the carbon resulting from the pyrolysis of Cu 2+ -chitosan [51] would be lost if the sample was prepared by impregnation of N-doped graphitic carbon by Cu(OAc) 2 .Data in the literature have shown that pyrolysis at 900 • C of Cu 2+ ions adsorbed on chitosan renders a material in which the resulting Cu NPs become strongly grafted on the (N)C support, as deduced from the relatively small particle size, the flat morphology of the Cu NPs, their preferential (111) facet orientation matching the graphene structure, and binding energy shifts in XPS [51][52][53].Therefore, in the present study, we proceeded to incorporate Cu(OAc) 2 to the chitosan before pyrolysis to obtain Cu@(N)C and, subsequently, to impregnate the desired Zn(OAc) 2 amount on the preformed Cu@(N)C.It was reasoned that in this way the interaction between the Cu NPs and defective graphene sheets of the carbon matrix as a substrate would be strong [51] and subsequent deposition of Zn(OAc) 2 could still reconstruct the Cu-Zn alloy under the reaction conditions.Scheme 1 illustrates the steps performed in the preparation of Cu-ZnO@(N)C.Further details can be found in the experimental section.
could still reconstruct the Cu-Zn alloy under the reaction conditions.Scheme 1 illustrates the steps performed in the preparation of Cu-ZnO@(N)C.Further details can be found in the experimental section.
A series of samples with different atomic Cu/Zn ratios were prepared, trying to cover a wide range of Cu/Zn ratios, including a Cu/Zn ratio of around 3 that is close to the composition of the benchmark Cu-ZnO/Al2O3 catalyst.Table S1 in the supporting information indicates the exact weights of Cu(OAc)2 and Zn(OAc)2 used in the preparation of samples 1-6.Unavoidably, the total Cu + Zn loading on the material varied from sample to sample due to poor control of the weight loss during pyrolysis, although for samples 3-6 the Cu + Zn loading was close to 10%.Note that sample 1 containing only Cu was the material to which Zn(OAc)2 was not added and, therefore, it has a lower total metal percentage.Table 1 summarizes the set of samples under study, the main Cu-Zn analytical data, and the average metal particle size.It is worth noting that although NaOH was used in the formation of Cu(OAc)2-impregnated chitosan beads, due to their water solubility and the sublimation of any possible residual Na during the pyrolysis, the Na content of the final Cu-ZnO@(N)G samples was negligible.The percentages of Cu and Zn in the samples were determined by ICP-OES elemental analysis after treating the Cu@(N)C and Cu-ZnO@(N)C samples with aqua regia, quantifying the metal content of the digested liquor.These data summarized in Table 1 show that the Cu/Zn ratio ranges from ∞ for sample 1, which does not contain ZnO, to 0.53 for sample 2, which is the sample with the highest ZnO proportion.This variation in the Cu/Zn ratio allows one to gain information on the influence of this parameter on methanol selectivity.A series of samples with different atomic Cu/Zn ratios were prepared, trying to cover a wide range of Cu/Zn ratios, including a Cu/Zn ratio of around 3 that is close to the composition of the benchmark Cu-ZnO/Al 2 O 3 catalyst.Table S1 in the supporting information indicates the exact weights of Cu(OAc) 2 and Zn(OAc) 2 used in the preparation of samples 1-6.Unavoidably, the total Cu + Zn loading on the material varied from sample to sample due to poor control of the weight loss during pyrolysis, although for samples 3-6 the Cu + Zn loading was close to 10%.Note that sample 1 containing only Cu was the material to which Zn(OAc) 2 was not added and, therefore, it has a lower total metal percentage.Table 1 summarizes the set of samples under study, the main Cu-Zn analytical data, and the average metal particle size.It is worth noting that although NaOH was used in the formation of Cu(OAc) 2 -impregnated chitosan beads, due to their water solubility and the sublimation of any possible residual Na during the pyrolysis, the Na content of the final Cu-ZnO@(N)G samples was negligible.
The percentages of Cu and Zn in the samples were determined by ICP-OES elemental analysis after treating the Cu@(N)C and Cu-ZnO@(N)C samples with aqua regia, quantifying the metal content of the digested liquor.These data summarized in Table 1 show that the Cu/Zn ratio ranges from ∞ for sample 1, which does not contain ZnO, to 0.53 for sample 2, which is the sample with the highest ZnO proportion.This variation in the Cu/Zn ratio allows one to gain information on the influence of this parameter on methanol selectivity.
High-resolution TEM images show that the Cu@(N)C and Cu-ZnO@(N)C samples contain metal NPs deposited on 2D defective graphene sheets that constitute the graphitic carbon matrix.Figure 1 shows selected dark-field TEM images taken in three different areas for the samples under study, illustrating that the metal NPs are homogeneously distributed through the carbon matrix as a consequence of the preparation procedure.The white zones of the images indicate where the Cu and Zn metals are present, while the black background indicates the absence of these metals.These images show that the metals are spread out through the carbon matrix.The particle size distribution and the average dimension, ranging from 1.0 ± 0.2 to 1.4 ± 0.4 nm, were determined by measuring the size of a statistically relevant number of those metal NPs.Similar average particle size values for the series of samples are also collected in Table 1, while the corresponding particle size distribution histograms are inserted in the DF-TEM images presented in Figure 1.The absence of large NPs can be observed in the images by the absence of bright dots.This small dimension of the Cu-ZnO particles in spite of the relatively high loading of metal (up to 14 wt.%) reflects the occurrence of a strong interaction of the Cu NPs with the defective N-doped graphitic carbon that thwarts the growth of the Cu particle even though the pyrolysis is carried out at 900 • C. High-temperature annealing is known to cause agglomeration of small metal NPs, as has been observed in other cases [51].Thus, the presence of N on the graphitic carbon matrix is a prerequisite to obtain small metal NPs, since in the absence of N doping the resulting particle size is considerably larger, even in the range of 100 nm [54].There are precedents in the literature claiming the interaction of N atoms on graphene with supported metals, such as Pt resulting in the formation of single atoms or small clusters [55-57].High-resolution FESEM images at the 100-400 nm scale recorded for the Cu@(N)C and Cu-ZnO@(N)C materials reveal a fluffy, poorly packed, and highly porous morphology of the graphitic carbon matrix acting as a support for the metal-metal oxide NPs.This porous structure is inherited from chitosan aerogel beads dried in super-critical CO2 in which aggregation of chitosan fibrils by hydrogen bridges has been minimized [61].Figure 2 and Figure S2 show selected HR-FESEM images for Cu@(N)C (sample 1) and Cu-ZnO@(N)C samples 2-6, illustrating the porous, coral-like structure of the carbonaceous matrix, resulting from the graphitization of the polysaccharide fibrils of the precursor.As expected, no metal NPs could be observed in the HR-FESEM images due to their lower resolution, in agreement with the nanometric particle size of Cu-ZnO NPs measured by HR-TEM.To gain information on the interaction between the Cu-ZnO NPs and the (N)G support and to determine the oxidation state of the fresh Cu-ZnO@(N)C samples, XPS analysis of catalyst 4 was performed as a representative sample of the series.XPS analysis of sample 4 revealed the presence of the expected Cu, Zn, C, N, and O elements, but with remarkably different proportions on the surface compared to the analytical data of the bulk material.The elemental proportion based on XPS is provided in Table S2 of the supporting information, while the XPS peaks for the elements and their best deconvolutions Samples 1-6 were also characterized by XRD and Raman spectroscopy (Figure S1).In the XRD patterns, the expected diffraction peaks corresponding to metallic Cu (PDF No. 70-3038) [58] and ZnO (JCPDS No. 36-1451) [59] could be clearly identified for samples 1-6. Figure S1a shows the XRD patterns, indicating the assignment of the peaks, either to Cu metal or ZnO.The relative intensity of the peaks corresponding to ZnO was in accordance with the relative proportion of ZnO in the sample with respect to Cu.In addition, Figure S1 also shows the formation of graphitic carbon with a broad peak appearing around 25 • in the XRD pattern.The Scherrer equation was used to quantify the size of Cu and ZnO NPs that are based on the full width at half height of the most intense peaks in XRD resulting in values of 7.8 and 25.0 nm, respectively.These average sizes are much larger than those determined by TEM, appearing in Table 1.This discrepancy could be due to the fact that XRD measures the most crystalline particles in the Cu-ZnO(N)C samples and the small metal particles observed in TEM do not contribute much to the XRD.
The defective nature of N-doped graphene was established by Raman spectroscopy, where the characteristic G and D bands appearing at 1590 and 1350 cm −1 , together with resolved overtone 2D at 2700 cm −1 , were recorded [60].Figure S1b plots the representative Raman spectra recorded for each sample 1-6.The relative intensity ratio of the G vs. the D band was about 1.15 and their width at half height is in accordance with values for N-doped graphitic carbons previously reported from chitosan [60].
High-resolution FESEM images at the 100-400 nm scale recorded for the Cu@(N)C and Cu-ZnO@(N)C materials reveal a fluffy, poorly packed, and highly porous morphology of the graphitic carbon matrix acting as a support for the metal-metal oxide NPs.This porous structure is inherited from chitosan aerogel beads dried in super-critical CO 2 in which aggregation of chitosan fibrils by hydrogen bridges has been minimized [61].Figure 2 and Figure S2 show selected HR-FESEM images for Cu@(N)C (sample 1) and Cu-ZnO@(N)C samples 2-6, illustrating the porous, coral-like structure of the carbonaceous matrix, resulting from the graphitization of the polysaccharide fibrils of the precursor.As expected, no metal NPs could be observed in the HR-FESEM images due to their lower resolution, in agreement with the nanometric particle size of Cu-ZnO NPs measured by HR-TEM.
statistical particle size distribution of samples 1-6 according to the main frame image.
High-resolution FESEM images at the 100-400 nm scale recorded for the Cu@(N)C and Cu-ZnO@(N)C materials reveal a fluffy, poorly packed, and highly porous morphology of the graphitic carbon matrix acting as a support for the metal-metal oxide NPs.This porous structure is inherited from chitosan aerogel beads dried in super-critical CO2 in which aggregation of chitosan fibrils by hydrogen bridges has been minimized [61].Figure 2 and Figure S2 show selected HR-FESEM images for Cu@(N)C (sample 1) and Cu-ZnO@(N)C samples 2-6, illustrating the porous, coral-like structure of the carbonaceous matrix, resulting from the graphitization of the polysaccharide fibrils of the precursor.As expected, no metal NPs could be observed in the HR-FESEM images due to their lower resolution, in agreement with the nanometric particle size of Cu-ZnO NPs measured by HR-TEM.To gain information on the interaction between the Cu-ZnO NPs and the (N)G support and to determine the oxidation state of the fresh Cu-ZnO@(N)C samples, XPS analysis of catalyst 4 was performed as a representative sample of the series.XPS analysis of sample 4 revealed the presence of the expected Cu, Zn, C, N, and O elements, but with remarkably different proportions on the surface compared to the analytical data of the bulk material.The elemental proportion based on XPS is provided in Table S2 of the supporting information, while the XPS peaks for the elements and their best deconvolutions are presented in Figure S3.As it can be seen there, the percentages of the Cu and Zn elements on the surface are much below the values expected by the bulk analysis, with C, O, and N of the graphitic carbon being the prevalent surface elements.This result is in accordance with the carbon matrix wrapping the metal NPs.In addition, the surface Cu/Zn XPS ratio is about 0.5, far from the 4.2 ratio measured for the bulk sample.Since samples 2-5 are obtained by Zn impregnation on preformed Cu@(N)C, it appears that the Zn element is more external than the originally introduced Cu metal due to it being incorporated in the samples last.Besides elemental composition, analysis of the high resolution XPS Zn 2p core level peak indicates that it corresponds well to a single ZnO component with a To gain information on the interaction between the Cu-ZnO NPs and the (N)G support and to determine the oxidation state of the fresh Cu-ZnO@(N)C samples, XPS analysis of catalyst 4 was performed as a representative sample of the series.XPS analysis of sample 4 revealed the presence of the expected Cu, Zn, C, N, and O elements, but with remarkably different proportions on the surface compared to the analytical data of the bulk material.The elemental proportion based on XPS is provided in Table S2 of the supporting information, while the XPS peaks for the elements and their best deconvolutions are presented in Figure S3.As it can be seen there, the percentages of the Cu and Zn elements on the surface are much below the values expected by the bulk analysis, with C, O, and N of the graphitic carbon being the prevalent surface elements.This result is in accordance with the carbon matrix wrapping the metal NPs.In addition, the surface Cu/Zn XPS ratio is about 0.5, far from the 4.2 ratio measured for the bulk sample.Since samples 2-5 are obtained by Zn impregnation on preformed Cu@(N)C, it appears that the Zn element is more external than the originally introduced Cu metal due to it being incorporated in the samples last.Besides elemental composition, analysis of the high resolution XPS Zn 2p core level peak indicates that it corresponds well to a single ZnO component with a binding energy value of 1022 eV, in accordance with the literature [62].In contrast, the Cu 2p spectra indicate two components attributable to Cu 0 and Cu II oxidation states appearing at 932.3 and 934.4 eV, respectively.These values are downshifted by 0.7 eV in the case of Cu 0 and upshifted by 0.9 eV for Cu II , with respect to the reported literature values for these two oxidation states [63].These shifts in the binding energy support the occurrence of a strong Cu-(N)G interaction due to the preparation procedure based on high-temperature graphitization, as has been reported earlier [64].
For the sake of comparison and to put into context the catalytic activity of samples 1-6, three samples consisting of Cu-ZnO supported on Al 2 O 3 (Cu-ZnO/Al 2 O 3 ) were also prepared.A similar two-step impregnation procedure to that used for the preparation of Cu-ZnO@(N)C was followed for the preparation of Cu-ZnO/Al 2 O 3 , trying to reproduce the method employed in the preparation of Cu-ZnO@(N)C with Al 2 O 3 as a support.A second Cu-ZnO/Al 2 O 3 -wp (wp meaning without pyrolysis) was also prepared, consecutively adsorbing Cu and Zn, but without submitting the sample to pyrolysis.A third Cu-ZnO/Al 2 O 3 was prepared by simultaneous Cu and Zn salts impregnation and subsequent mild baking at 250 • C (Cu-ZnO/Al 2 O 3 -imp, imp meaning impregnation).

Catalytic Activity
The objective of the present study is to establish the performance of Cu-ZnO active sites supported on N-doped graphitic carbon as catalysts for the selective partial hydrogenation of CO 2 to methanol, following the lead of previously reported Fe-Co@(N)C catalysts that exhibit very high selectivity for CO 2 hydrogenation to methane [43], CO [44], or C 2+ [45], depending on the metal particle size and composition.The use of inorganic supports has been widely studied in the literature [65,66], but there is a limited effort made to gain information about the performance of graphitic carbons as a support in metal catalysts for methanol formation from CO 2 [67].
The catalytic experiments were carried out in a pressurized stainless-steel reactor operating at 40 bar under 4 mL of continuous flow and a H 2 to CO 2 ratio of 4. After considering the known thermodynamic limitations [24], the range of temperatures studied was between 150 and 300 • C in 50 • C steps that were maintained for 1 h before going to the next temperature increase.Systematic calculation of the Weisz-Prater number for each catalyst and conditions indicates that in none of the cases the reaction was under diffusion control (see Table S3 in supporting information).Previous controls operating at 300 • C in the absence of a catalyst showed that the austenite stainless steel reactor converts 1.1% of CO 2 with a selectivity to CO and CH 4 of 92.4 and 7.6%, respectively.This low CO 2 conversion and CO selectivity was maintained with the time of stream in a 3 h test and decreased with the reaction temperature.Table S4 in the supporting information lists the CO 2 conversion and product selectivity of these previous control experiments in the absence of a catalyst at different temperatures.Methanol was undetectable in these control experiments.
Catalyst 1 (Cu@(N)C) containing only Cu did not promote methanol formation and catalyzed mostly the RWGS (Equation ( 2)), accompanied by some formation of CH 4 and C 2+ products, with a combined selectivity of about 18%.This relatively high proportion of hydrocarbons in catalyst 1 could indicate the occurrence of some Fischer-Tropsch synthesis in which CO undergoes further hydrogenation of these hydrocarbons.
In contrast to the blank controls and the results with sample 1 lacking ZnO, methanol formation was observed in most of the reactions carried out in the presence of graphitic carbon-supported Cu-ZnO catalysts.The catalytic results achieved by sample 1 make clear that the presence of ZnO is required to drive the selectivity towards methanol.Analysis of the reaction products shows that besides the formation of methanol and CO as the major products, lesser amounts of CH 4 , C 2 H 6 , and C 3 H 8 (below a combined percentage of 5%) are also formed (see Tables S5-S12 and Figure S4).It should be commented that although alkali metal ions can be promoters of the catalytic activity of transition metals in CO 2 hydrogenation reactions and we have used NaOH during the preparation of the present Cu-ZnO@(N)G samples, the Na content in these samples is negligible.
To optimize the methanol formation, the performance of Cu-ZnO@(N)C samples was studied at different temperatures from 150 to 300 • C. Figure 3b shows the variation in CO 2 conversion and product selectivity in the case of sample 4, while the results of other samples are presented in Figure S4.As expected, CO 2 conversion increased for all the catalysts upon increasing the temperature in the 150-300 • C range.Methanol selectivity showed the opposite trend, being higher at lower CO 2 conversions and decreasing in favor of CO at higher temperatures.This general behavior agrees with the thermodynamics of the two main competing reactions presented in Equations ( 1) and ( 2), with methanol formation being exothermic and RWGS giving CO, being endothermic [24].At temperatures of 350 • C or higher, CO 2 conversions increase, but methanol selectivity becomes negligible or even methanol formation becomes undetectable.From the dependence of the CO 2 conversion on the temperature (Figure S6 in supporting information), the apparent activation energies (E a ) for CO 2 hydrogenation in the range of temperatures between 150 and 300 • C were obtained for each catalyst of the series (see Figure S4).Except for sample 3, the E a was in the range of 50 to 43 kJ mol −1 , with the lowest value being for sample 4.
sion and selectivity for samples 1-6 working at 300 °C, 40 bar, at a H2/CO2 ratio of 4. As can be seen, besides differences in CO2 conversion, CO was the main product for all catalysts at 300 °C, although the formation of methanol was observed for all the series of Cu-ZnO@(N)C samples 2-6. Figure S4 provides CO2 conversion and selectivity for samples 1-6 at 200 and 250 °C.The relative activity order of samples 1-6 found at 300 °C is maintained at 200 and 250 °C, with sample 4 being the best performing catalyst.From the screening of the Cu-ZnO@(N)C catalysts under study shown in Figure 3a, sample 4 was the best performing catalyst, reaching at 200 °C a conversion of 1.7% with the maximum methanol selectivity of 89.7%.The highest methanol productivity achieved for sample 4 as a catalyst was 2.57 molCH3OH kgcatalyst −1 h −1 or 82.24 gCH3OH kgcatalyst −1 h −1 .To put these data into context, Table S13 in the supporting information provides a comparison of reported data in the literature.Although a comparison of catalyst performance The series of catalysts using N-doped graphitic carbon as a support did not show a clear influence of the Cu/Zn ratio on methanol selectivity, which was mostly dependent on CO 2 conversion.For the same CO 2 conversion, similar methanol selectivity values were reached regardless of if the Cu/Zn ratio is high (sample 6, Cu/Zn 8.1) or low (sample 2, Cu/Zn 0.53).The main influence of the Cu/Zn ratio appears to be in CO 2 conversion, which was the highest for sample 4, the sample with the lowest E a .To illustrate differences in the catalytic performance depending on the Cu/Zn ratio, Figure 3a presents the CO 2 conversion and selectivity for samples 1-6 working at 300 • C, 40 bar, at a H 2 /CO 2 ratio of 4. As can be seen, besides differences in CO 2 conversion, CO was the main product for all catalysts at 300 • C, although the formation of methanol was observed for all the series of Cu-ZnO@(N)C samples 2-6.From the screening of the Cu-ZnO@(N)C catalysts under study shown in Figure 3a, sample 4 was the best performing catalyst, reaching at 200 • C a conversion of 1.7% with the maximum methanol selectivity of 89.7%.The highest methanol productivity achieved for sample 4 as a catalyst was 2.57 mol CH3OH kg catalyst −1 h −1 or 82.24 g CH3OH kg catalyst −1 h −1 .To put these data into context, Table S13 in the supporting information provides a comparison of reported data in the literature.Although a comparison of catalyst performance measured under different conditions must be taken always cautiously, Table S13 shows that the methanol productivity achieved in the present study is comparable or overcomes those previously reported for the best catalysts for methanol synthesis from CO 2 .
In fact, one key issue that is not reflected in Table S13 is catalyst stability.In the present case, sample 4 exhibits a remarkable catalytic stability, particularly in comparison to the benchmark Cu-ZnO/Al 2 O 3 .In a long run of 56 h under the reaction conditions, the catalytic activity of sample 4 did not decrease (Figure S7).Furthermore, after this long test at 300 • C, the 56 h used sample 4 exhibit, for two subsequent cycles of increased temperature from 150 to 300 • C in 1 h step of 50 • C temperature increase, identical catalytic data to those measured for the fresh sample presented in Figure S8 and Table S12.These data confirm again the catalyst's stability.This catalyst stability for long reaction runs was also observed for catalysts 2 and 3.It is of note that Cu-ZnO/Al 2 O 3 -imp did not produce methanol in the range from 200 to 300 • C, with CO and CH 4 being the only products formed (Table S12).Chart S1a,b present chromatograms of the reaction mixture to illustrate the quality of the analysis.
To determine how far from the equilibrium the catalytic values reached for catalysts 2-4, equilibrium data were calculated by minimization of the Gibbs free energy of the system using the RGIBBS module of the Aspenplus© program.For simplicity, these calculations consider the formation of methanol and CO as the only products, ignoring the small percentage of methane and other hydrocarbons detected in the product mixture in very low proportions (less than 5% according to Tables S2-S11).The composition of the equilibrium mixture was determined in the temperature range between 150 and 300 • C, introducing the operation conditions as the starting composition of the system.The results are presented in Table 2.As it can be seen there, a decrease in the CO 2 conversion from 250 to 350 • C is predicted thermodynamically due to the opposite enthalpy signs of the methanol synthesis (exoergonic, Equation ( 1)) and the RWGS (endoergonic, Equation ( 2)).These calculations indicate that the selectivity of methanol 350 • C should be very low.Comparison with the experimental data for catalyst 4 shows that while at 150 and 200 • C, CO 2 conversion is very far from the equilibrium values, the mixture composition approaches the expected equilibrium values in the temperature range from 250 to 300 • C, with theoretical conversions in the range of 25% and methanol selectivity between 60 and 12% under the operation conditions.This catalytic stability of the Cu-ZnO@(N)C samples contrasts with the performance of Cu-ZnO/Al 2 O 3 , obtained following a two-step impregnation procedure similar to that used for the preparation of Cu-ZnO@(N)C (see experimental section).Cu-ZnO/Al 2 O 3 can be taken as a benchmark catalyst with which one can compare the performance of Cu-ZnO@(N)C samples.The catalytic activity of Cu-ZnO/Al 2 O 3 is summarized in Table S5.As can be seen there, the CO 2 conversion with the fresh Cu-ZnO/Al 2 O 3 samples was somewhat lower than that achieved by the Cu-ZnO@(N)C samples, with a maximum 31% methanol selectivity at 300 • C at 12% CO 2 conversion.Therefore, the performance of the Cu-ZnO@(N)C samples compares well with that of fresh Cu-ZnO/Al 2 O 3 (Table S13), with selectivity to methanol decreasing with temperature and CO 2 conversion similarly in both cases.However, it was observed that Cu-ZnO/Al 2 O 3 undergoes a notable deactivation with time on stream, becoming severely deactivated in a few tens of hours.After 40 h reaction, the spent Cu-ZnO/Al 2 O 3 catalyst is black in color, suggesting coke deposition during the process (see Figure S5 in supporting information).Combustion chemical analysis of a black deactivated Cu-ZnO/Al 2 O 3 catalyst showed the presence of over 1% carbon in the material.In the literature, it has been reported that Cu-ZnO/Al 2 O 3 undergoes deactivation by sintering of the metal NPs [68] and by coke deposition [69]; not surprisingly, this deactivation happens here for the Cu-ZnO/Al 2 O 3 sample in the time scale of tens of hours.Cu-ZnO/Al 2 O 3 -wp not submitted to pyrolysis behaves similarly, in agreement with the XRD, which shows no difference in the Al 2 O 3 crystalline phase due to pyrolysis.Thus, activity and stability data confirm a similar performance and much better stability of Cu-ZnO@(N)C with respect to the Cu-ZnO/Al 2 O 3 reference catalyst.
The above results show the advantage of (N)C as a support in comparison to Al 2 O 3 .To understand the role of (N)C, the type of N atom having a stronger interaction with Cu-ZnO, the CO 2 adsorption of the Cu-Zn cluster, and the plausible reaction mechanism, DFT calculations were carried out.

Modelling and DFT Calculations
It is well known that the strong metal-support interaction between transition metal clusters and defective carbon (such as N-doped graphene) greatly contributes to the enhancement of the catalytic performance by tuning the electronic structures and improving the stability of, in particular, Cu and Zn-doped Cu clusters [70].The use of Cu-based catalysts for methanol synthesis via the CO 2 dehydrogenation route is as old as the process itself and Cu-ZnO/Al 2 O 3 catalysts are applicable on an industrial scale in the syngas route [66].On the other hand, CO 2 hydrogenation to methanol over Cu and Zn-doped Cu clusters supported on graphitic carbon has rarely been reported and its active site is still unclear [71][72][73][74].
To better understand the mechanism of methanol synthesis from CO 2 , periodic DFT calculations on ZnCu/N-C model catalysts were performed.A 5 × 5 super-cell of graphene with one pyridinic-N and one graphitic-N was built (Figure 4), with the two N-doped atoms being far enough away from each other to make their cross interaction negligible.A Cu cluster with 13 atoms (Cu 13 ) was simulated (Figure 4a), starting from the experimental geometry of these Cu clusters [75].
calculations were carried out.

Modelling and DFT Calculations
It is well known that the strong metal-support interaction between transition metal clusters and defective carbon (such as N-doped graphene) greatly contributes to the enhancement of the catalytic performance by tuning the electronic structures and improving the stability of, in particular, Cu and Zn-doped Cu clusters [70].The use of Cu-based catalysts for methanol synthesis via the CO2 dehydrogenation route is as old as the process itself and Cu-ZnO/Al2O3 catalysts are applicable on an industrial scale in the syngas route [66].On the other hand, CO2 hydrogenation to methanol over Cu and Zn-doped Cu clusters supported on graphitic carbon has rarely been reported and its active site is still unclear [71][72][73][74].
To be er understand the mechanism of methanol synthesis from CO2, periodic DFT calculations on ZnCu/N-C model catalysts were performed.A 5 × 5 super-cell of graphene with one pyridinic-N and one graphitic-N was built (Figure 4), with the two N-doped atoms being far enough away from each other to make their cross interaction negligible.A Cu cluster with 13 atoms (Cu13) was simulated (Figure 4a), starting from the experimental geometry of these Cu clusters [75].The binding strength of Cu13 on different sites of N-doped graphene was first examined by periodic DFT calculations.The binding energies of Cu13 on the pyridinic-N site, graphitic-N site, and graphitic-C site were calculated to be −5.70,−2.62, and −3.03 eV, respectively (Figure 4c-e).The results show that pyridinic N is the most favorable site for anchoring the Cu13 cluster (Figure S9) and the location of Cu13 on this site will be used here onwards.This proposal agrees with the Hirshfeld charge distribution (Table S14), showing a charge transfer from Cu to N-doped graphene of −0.75 e − , demonstrating that there is a strong interaction between the Cu cluster and N-doped graphene.Apparently, the The binding strength of Cu 13 on different sites of N-doped graphene was first examined by periodic DFT calculations.The binding energies of Cu 13 on the pyridinic-N site, graphitic-N site, and graphitic-C site were calculated to be −5.70,−2.62, and −3.03 eV, respectively (Figure 4c-e).The results show that pyridinic N is the most favorable site for anchoring the Cu 13 cluster (Figure S9) and the location of Cu 13 on this site will be used here onwards.This proposal agrees with the Hirshfeld charge distribution (Table S14), showing a charge transfer from Cu to N-doped graphene of −0.75 e − , demonstrating that there is a strong interaction between the Cu cluster and N-doped graphene.Apparently, the lone electron pair of pyridinic N interacts strongly with the Cu 13 cluster in comparison with graphitic N or C atoms.This strong catalyst-support interaction is believed to provide stability and support to the clusters, preventing their aggregation and maintaining their catalytic activity over multiple reaction cycles.
Additionally to the Cu 13 cluster, a Zn-doped Cu cluster (Zn 1 Cu 12 ) was also considered.The different substituted sites of the Zn 1 Cu 12 cluster on pyridinic-N of N-doped graphene were geometry optimized (Figures S10 and S11).Site 12 exhibits the minimum energy, indicating that it is the most stable structure.Thus, we used this model for subsequent calculations on the mechanism of CO 2 hydrogenation to methanol.In this model, the Zn atom is at the core surrounded by the twelve external Cu atoms.
To determine the active site of the reaction, different CO 2 adsorption configurations and sites were computed for the Zn 1 Cu 12 cluster on the pyridinic-N of N-doped graphene.Figure S12 shows that the line-to-line parallel adsorption has the strongest binding energy.Various non-equivalent line-to-line parallel adsorption sites were subsequently considered for CO 2 adsorption (Figure S13).The results show that one of the Cu-Cu sites (denoted as '4-9 ′ , see Figure S13) is the most active for CO 2 adsorption, compared not only among all Cu-Cu sites, but also among N-C, C-C, or Zn-Cu sites (Figure S14).In this absorption mode, two Cu atoms interact simultaneously with the C and the two O atoms of CO 2 that are quasi perpendicularly aligned with the graphene sheet.Therefore, subsequent intermediates and transition states of CO 2 conversion to methanol were studied on 4-9 Cu-Cu sites for the Zn 1 Cu 12 cluster on N-doped graphene.Based on reported studies [76,77], two main reaction pathways for CO 2 conversion to methanol were considered in our calculations: the (i) RWGS+CO+hydro pathway: RWGS reaction to produce a CO intermediate followed by its hydrogenation to methanol and (ii) formate pathway: initial hydrogenation of CO 2 to a *HCOO intermediate followed by its hydrogenation and dissociation to methanol.As will be commented below, and in agreement with previous studies [76,77], the present calculations also show the preferred formate pathway on ZnCu/C-N via *HCOOH, *H 2 COOH, and *CH 3 O intermediates over the RWGS+CO-hydro pathway for methanol synthesis (Figure 5 and Table S15).
Figure S12 shows that the line-to-line parallel adsorption has the strongest binding energy.Various non-equivalent line-to-line parallel adsorption sites were subsequently considered for CO2 adsorption (Figure S13).The results show that one of the Cu-Cu sites (denoted as '4-9′, see Figure S13) is the most active for CO2 adsorption, compared not only among all Cu-Cu sites, but also among N-C, C-C, or Zn-Cu sites (Figure S14).In this absorption mode, two Cu atoms interact simultaneously with the C and the two O atoms of CO2 that are quasi perpendicularly aligned with the graphene sheet.
Therefore, subsequent intermediates and transition states of CO2 conversion to methanol were studied on 4-9 Cu-Cu sites for the Zn1Cu12 cluster on N-doped graphene.Based on reported studies [76,77], two main reaction pathways for CO2 conversion to methanol were considered in our calculations: the (i) RWGS+CO+hydro pathway: RWGS reaction to produce a CO intermediate followed by its hydrogenation to methanol and (ii) formate pathway: initial hydrogenation of CO2 to a *HCOO intermediate followed by its hydrogenation and dissociation to methanol.As will be commented below, and in agreement with previous studies [76,77], the present calculations also show the preferred formate pathway on ZnCu/C-N via *HCOOH, *H2COOH, and *CH3O intermediates over the RWGS+CO-hydro pathway for methanol synthesis (Figure 5 and Table S15).[65].Consequently, the rate-determining step is believed to occur in the hydrogenation of HCOO*.Accordingly, our results reveal a low-   [65].Consequently, the rate-determining step is believed to occur in the hydrogenation of HCOO*.Accordingly, our results reveal a lowenergy barrier for the formation of HCOO* (0.35 eV, TS1, Figure S15a, Table S15).Then, the rate-determining step is the hydrogenation of *HCOO to form *HCOOH, with an activation energy of 2.12 eV (TS2, Figure S15b, Table S15).Subsequent hydrogenation leads to the formation of *H 2 COOH, with an energy barrier of 0.79 eV (TS3, Figure S15c, Table S15).This is followed by the dissociation of *H 2 COOH into *H 2 CO + *OH, with a barrier of 1.10 eV (TS4, Figure S15d, Table S15), and the hydrogenation of *OH to water, with a barrier of 1.77 eV (TS5, Figure S15e, Table S15).Finally, *H 2 CO is hydrogenated twice to form *H 3 CO (0.84 eV, TS6, Figure S15f, Table S15) and *CH 3 OH (1.56 eV, TS7, Figure S15g, Table S15).These results indicate that the reaction of *HCOO + *H→*HCOOH, with the largest energy barrier (2.12 eV), is the rate-determining step for the synthesis of methanol.Both the low barrier for HCOO* formation and the subsequent largest reaction barrier are lower than the maximum barrier in the RWGS+CO-hydro pathway discussed below, indicating that the formation of methanol should be preferred through the formate pathway.

Figure 3 .
Figure 3. CO2 conversion and selectivity for samples 1-6 having different molar ratios of Cu and Zn at 300 °C (a) under the same conditions for sample 4 (b) at temperatures from 150 to 300 °C.Reaction conditions: H2/CO2 ratio of 4, total flow 4 mL/min, 40 bar, 40 mg catalyst.Blank test: no catalyst.

Figure 3 .
Figure 3. CO 2 conversion and selectivity for samples 1-6 having different molar ratios of Cu and Zn at 300 • C (a) under the same conditions for sample 4 (b) at temperatures from 150 to 300 • C. Reaction conditions: H 2 /CO 2 ratio of 4, total flow 4 mL/min, 40 bar, 40 mg catalyst.Blank test: no catalyst.
Figure S4 provides CO 2 conversion and selectivity for samples 1-6 at 200 and 250 • C. The relative activity order of samples 1-6 found at 300 • C is maintained at 200 and 250 • C, with sample 4 being the best performing catalyst.

3. 3
.1.Formate Pathway HCOO* has been considered as the main intermediate for hydrogenating CO2 to CH3OH through the formate pathway

Figure 5 .
Figure 5. Potential energy diagram for the hydrogenation of CO 2 to CH 3 OH.Left: RWGS+CO+hydro pathway.Right: Formate pathway.

3. 3
.1.Formate Pathway HCOO* has been considered as the main intermediate for hydrogenating CO 2 to CH 3 OH through the formate pathway

Table 2 .
Calculated equilibrium compositions using Aspenplus as a function of the reaction temperature.Initial conditions: Pressure 40 bar, H 2 3.2 mL/min, CO 2 0.8 mL/min.