Next Article in Journal
Formation of Pre-PCTA/DT Intermediates from 2-Chlorothiophenol on Silica Clusters: A Quantum Mechanical Study
Previous Article in Journal
Investigating Vα7.2+/CD161 T Cell and MAIT Cell Profiles Using Flow Cytometry in Healthy Subjects and Subjects with Atopic Dermatitis
Previous Article in Special Issue
Nb-Based Catalysts for the Valorization of Furfural into Valuable Product through in One-Pot Reaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogenolysis of Glycerol over NiCeZr Catalyst Modified with Mg, Cu, and Sn at the Surface Level

by
Norberto Vera-Hincapie
1,
Unai Iriarte-Velasco
2,
Jose Luis Ayastuy
1,* and
Miguel Ángel Gutiérrez-Ortiz
1
1
Department of Chemical Engineering, Faculty of Science and Technology, University of the Basque Country UPV/EHU, Sarriena S/N, 48940 Leioa, Spain
2
Department of Chemical Engineering, Faculty of Pharmacy, University of the Basque Country UPV/EHU, Paseo de la Universidad, 7, 01006 Vitoria, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(6), 3484; https://doi.org/10.3390/ijms25063484
Submission received: 26 February 2024 / Revised: 13 March 2024 / Accepted: 15 March 2024 / Published: 20 March 2024
(This article belongs to the Special Issue Heterogeneous Catalysis for Biomass Valorization)

Abstract

:
Biomass valorization is an essential strategy for converting organic resources into valuable energy and chemicals, contributing to the circular economy, and reducing carbon footprints. Glycerol, a byproduct of biodiesel production, can be used as a feedstock for a variety of high-value products and can contribute to reducing the carbon footprint. This study examines the impact of surface-level modifications of Mg, Cu, and Sn on Ni-Ce-Zr catalysts for the hydrogenolysis of glycerol, with in situ generated hydrogen. The aim of this approach is to enhance the efficiency and sustainability of the biomass valorization process. However, the surface modification resulted in a decrease in the global conversion of glycerol due to the reduced availability of metal sites. The study found that valuable products, such as H2 and CH4 in the gas phase, and 1,2-PG in the liquid phase, were obtained. The majority of the liquid fraction was observed, particularly for Cu- and Sn-doped catalysts, which was attributed to their increased acidity. The primary selectivity was towards the cleavage of the C–O bond. Post-reaction characterizations revealed that the primary causes of deactivation was leaching, which was reduced by the inclusion of Cu and Sn. These findings demonstrate the potential of Cu- and Sn-modified Ni-Ce-Zr catalysts to provide a sustainable pathway for converting glycerol into value-added chemicals.

1. Introduction

The demand for sustainable solutions to address environmental challenges posed by growing human populations, rapid industrialization, and excessive fossil fuel consumption has increased. Biomass valorization is a crucial strategy for converting organic resources into valuable energy and chemicals, contributing to the circular economy, and reducing carbon footprints.
Biomass valorization is the process of converting biomass into high-value energy and chemicals. This can be achieved through various methods, including thermochemical processes like gasification and pyrolysis [1,2], chemical processes such as hydrolysis [3], and biological methods like anaerobic digestion [4]. Each method has its own advantages and challenges. For example, thermochemical processes can be expensive, while biological methods may have lower yields. The integration of thermochemical and biological processes offers a promising solution to the challenges mentioned [5]. This enhances resource efficiency and reduces environmental impact. Integrated biorefineries, which combine various conversion technologies, mark a crucial step towards sustainable bioenergy and bioproduct production, including biofuels such as biodiesel [6,7].
Biodiesel production is a sustainable alternative to fossil fuels that also generates glycerol, a byproduct with significant potential in sustainable chemistry due to its versatile properties that makes it an ideal raw material for a wide range of products. The valorization routes for glycerol depend on the activation and reactivity of specific bonds such as C–C, C–O, C–H, and O–H. The preferential cleavage of these bonds determines the catalyst properties and process conditions [8,9]. Bifunctional catalysts containing both metal and acid/basic sites are suggested as the most suitable for cleaving C–C and C–O bonds [10].
Several methods have been investigated for converting glycerol into fuel alternatives using technologies developed for glycerol valorization. Glycerol steam reforming, which uses nickel-based and noble metal catalysts, has attracted significant interest due to the high concentration of hydrogen it produces [11,12,13]. At the same time, liquid-phase hydrogenolysis has emerged as another important approach [14]. This method uses hydrogen, which is often produced via reforming processes, to convert glycerol into a variety of valuable chemicals. Hydrogenolysis is particularly attractive because it operates under milder conditions (around 200 °C and 30 bar) [15]. This makes it an excellent option for glycerol valorization [16]. Coupling the hydrogenolysis with the aqueous phase reforming (APR) of glycerol produces H2 in situ, making it a more sustainable and safer alternative to traditional methods that require the external addition of H2. The APR and hydrogenolysis can be easily integrated into a one-pot reaction due to their comparable operational conditions. During this process, notable reactions occur, including dehydration–hydrogenation, which facilitates hydroxyacetone formation at Lewis acid sites, and the dehydrogenation–dehydration–hydrogenation sequence, which is particularly effective under basic conditions [17,18]. These pathways are crucial for the production of propylene glycols (PG), particularly 1,2-PG, due to its widespread use as a platform molecule in the industry.
Bimetallic catalysts have demonstrated promising outcomes in glycerol hydrogenolysis [19,20,21]. The addition of promoters such as Sn, Mg, Co, and Cu enhances selectivity through various mechanisms, such as blocking surface defects and reducing particle size [22]. Ni-Mg catalysts exhibit increased H2 selectivity by enhancing nickel’s electronegativity, reducing the dehydration pathway [23]. Ni-Cu catalysts were found to have better Ni dispersion, reducibility, and WGS, while also resisting oxidation, sintering, and coke formation [24]. On the other hand, Ni-Sn catalysts were found to suppress CO methanation by modifying nickel’s defective sites, offering increased sintering resistance [25].
This study investigates the effect of Mg, Cu, and Sn on the performance of NiCeZr catalysts in glycerol hydrogenolysis. The undoped NiCeZr catalyst has been proven effective for producing valuable liquids through glycerol hydrogenolysis [26]. In this study, we analyze the structural and surface changes that occur upon the addition of dopants and evaluate their effect on reaction efficiency and selectivity. The catalysts synthesized were characterized in both their fresh and reduced forms to establish a correlation between their physicochemical properties and catalytic performance. Additionally, the spent catalysts were characterized to identify the primary causes of deactivation.

2. Results and Discussion

2.1. Textural Properties and Bulk Chemical Composition

The NiCeZr base contained 9.2 wt.% of Ni and a Ce/Zr atom ratio of 85.9/14.1, which is close to the targeted values. The promoters were loaded according to the targeted values, except for Sn, which had a 20% deficit (Table 1). All the solids exhibited a type IV (a) isotherm with an H2 (b) hysteresis cycle (Figure S1, Supplementary Materials); these are characteristics of mesoporous solids with a bottleneck pore structure. Cu and Sn impregnation blocked the smallest pores, resulting in a slight increase in the average pore size compared to NiCeZr (Table S1, Supplementary Materials). However, the impregnation of Mg had little effect on the pore sizes.
The calcined NiCeZr assay had a SBET of 235 m2/g and decreased by 4–15% (depending on the modifier) after impregnation (Table 1). After impregnation with Mg and Cu, the pore volume tended to increase, likely due to the formation of larger pores, while impregnation with Sn reduced the Vpore. The surface tension of the water and isopropanol mixture used to impregnate Sn was lower than that of water alone. This could facilitate greater wetting and allow Sn to penetrate all pores more effectively, thereby decreasing pore volume.
The pore structure of the reduced solids remained unchanged since the isotherms showed little variation. The average pore size of Mg/NiCeZr and Cu/NiCeZr solids increased slightly compared to their calcined counterparts. The catalyst impregnated with Sn exhibited bimodal PSD curves. The SBET of the reduced solids decreased if compared to their calcined counterparts, (by 10–19%) except for Sn/NiCeZr, which was fully retained. Interestingly, all the modified solids exhibited a greater SBET than the base NiCeZr, suggesting that promoters enhanced the thermal stability of the solids.

2.2. Structural Characterization

The XRD diffractograms of the calcined solids (Figure 1A) showed broad diffraction peaks, indicative of a small crystallite size. The most intense signals were attributed to the NiCeZr solid solution that crystallized in the tetragonal phase (ICDD 01-088-2398). Additionally, a barely discernible peak at 44.3°, attributed to NiO (ICDD 01-078-0643), was observed, suggesting that some of the nickel had segregated as NiO. After impregnating the promoters, no distinct peaks for MgO, CuO, or SnOX were observed, indicating their high dispersion. These results were expected since the nominal loading of the promoters was about 1/20 of the monolayer.
Upon reduction, the dominant peaks observed in their calcined forms were preserved, with the absence of other peaks from the Ce-containing or Zr-containing phases, suggesting the CeZr mixed oxide remained. For the base NiCeZr assay, the peak from NiO disappeared and a new peak emerged at around 44.5°, which was attributable to the metallic Ni (ICDD 01-087-0712). This suggests that the nickel was reduced and remained as metallic clusters on the catalyst’s surface. A similar outcome was obtained for the Mg/NiCeZr material. The spectrum of the Cu/NiCeZr sample did not show the NiO signal. Instead, a new peak appeared at 43.7°, which was attributed to an intermetallic Ni-Cu alloy [27], together with a peak from metallic nickel (inset in the figure). The diffractogram of the reduced Sn/NiCeZr solid showed peaks from metallic nickel and NiO, the presence of the latter suggesting a lower reducibility of this solid. Thus, it could be concluded that, after reduction, nickel in different oxidation states coexisted within the Sn/NiCeZr solid, which was a similar outcome to previous studies [28].
The crystallographic parameters in Table 1 showed that the calcined materials had crystallite dimensions ranging from 2.3 to 2.8 nm. This small size was consistent with the broad diffraction peaks of low intensity in the XRD diffractograms and the high SBET. The solid particles were bigger than crystallites, indicating that particles formed by an aggregation of up to two crystallites (Table 1). The lattice length of the modified solids increased along the a-axis while decreasing along the c-axis, indicating structural changes upon the addition of the modifier.
After reduction, all samples showed an increase in both the dimensions of the NiCeZr crystallites and the lattice length along the a-axis. This expansion of the crystalline network could be attributed to the release of oxygen from the structure and the decrease in the oxidation state of specific cations. The repositioning of the cations and oxygen vacancies facilitated these structural changes [29]. The increase in the a-axis was found to be correlated with the atomic size of the modifier, following the sequence NiCeZr < Mg/NiCeZr ≈ Cu/NiCeZr < Sn/NiCeZr.
The size of the metallic nickel crystallites in the NiCeZr base catalyst was 11.4 nm and varied unevenly with the modifier. The most significant variation occurred in the Cu/NiCeZr catalyst (decreased to 5.0 nm) and the Sn/NiCeZr catalyst (increased to 20.1 nm), while the Mg/NiCeZr catalyst preserved the size of the Ni0 crystallites. These results suggest that a strong interaction between Ni and Cu led to an alloy with a reduced crystallite size.
Further characterization of the solid structures was carried out through Raman spectroscopy (Figure 2). This technique, complementing the XRD technique, is particularly useful in the analysis of NiCeZr mixed oxides [30]. The spectra of the calcined solids presented broad bands, which posed a challenge for their assignment. The broadness of the bands was caused by both the small grain size, as observed from the XRD results, and the presence of ordered defects such as anion vacancies and substituted cations [30]. For the base NiCeZr calcined solid, two primary signals were observed: one at approximately 525 cm−1, attributed to the Ni-O bonds, and another at approximately 615 cm−1, named the D band of defects, which was related to external defects resulting from the addition of Ni in CeZr solids [31]. In the surface-modified solids, a redshift of the Ni-O band was observed, suggesting alterations in the Ni neighborhood. This phenomenon was noted across all materials, indicating modifications in the electronic structure and potentially in the arrangement of atoms adjacent to Ni.
After reduction, the Raman profiles changed for all the catalysts. For example, for the base NiCeZr and Mg/NiCeZr catalysts, the D band vanished and the band from Ni-O bonds shifted to 550 cm−1. For the Cu/NiCeZr catalyst, the Ni-O band notably broadened, while it almost vanished for the Sn/NiCeZr catalyst. These solids were the only solids where the D band remained.

2.3. Redox Properties from the H2-TPR Analysis and Oxygen Storage Capacity

The reduction profiles of the base NiCeZr solid (Figure 3) contained three primary contributions, each corresponding to the reduction of different chemical species [32,33]. The contribution below 200 °C (named α-contribution) was associated with the reduction of easily reducible oxygen at oxygen vacancy sites (OV). The intermediate temperature contribution (the β-contribution) in the range between 200 and 360 °C involved the reduction of free NiO entities on the surface. Lastly, the γ-contribution, encompassing the range from the end of the β-contribution up to 800 °C, included the reduction of Ni2+ species in the solid solution, as well as of Ce4+ ions on the surface and in bulk. The reduction profile of the Mg/NiCeZr solid was similar to that of NiCeZr, the γ-peak getting broader and upshifting by 7 °C. These results suggested that the presence of Mg somewhat hindered the reduction of Ce4+ and/or Ni2+ in the ternary solution [34,35]. The reduction profile of the Cu/NiCeZr solid differed from that of the base NiCeZr in two main ways. Firstly, there were two well-defined peaks in the β-region, corresponding to the separated reduction of highly CuO to Cu0 (peaking at 245 °C) and NiO to Ni0 (peaking at 342 °C) [36,37].
Secondly, the γ-contribution was shifted downwards from 531 °C to 505 °C, suggesting that the interactions between the copper oxide and NiCeZr support enhances the reducibility of ceria and Ni2+ in the solid solution [37], due to the spillover effect on metallic Cu [38]. Finally, the reduction profile of Sn/NiCeZr was similar to that of the base NiCeZr. The profile showed a general upshift across the entire range, specifically by 30–45 °C for the β- and γ-contributions. This suggests that the presence of surface SnO2 hindered the reduction of surface NiO, and Ni2+ and Ce4+ in the solid solution [39]. The reduction profile of the supported SnO2 typically exhibits two contributions in the 400–600 °C range. These contributions are attributed to a sequential Sn4+ → Sn2+ → Sn0 reduction process [40]. In our sample, they were masked by the intense peak from the γ-contribution.
Table 2 shows that the reduction percentage for H2-TPR was slightly above 100% for all catalysts except for the Mg/NiCeZr assay, which was 78%. This might be due to the decrease in available NiO. The interaction between Ni2+ and Mg2+ ions can affect the distribution of Ni2+ ions on the surface and subsurface of the material, with the insertion of Ni2+ in the MgO lattice, making it much less reducible [41].
The Mg/NiCeZr and Sn/NiCeZr catalysts showed a decrease in hydrogen consumption related to both OV (α-contribution) and β-contribution. The former may be attributed to the fact that OV sites could serve as anchor points for these metal oxides [42]. However, the α-contribution in the Cu/NiCeZr catalyst was similar to that of the base catalyst, as new oxygen vacancies could be formed upon Cu impregnation [43]. The catalyst showed an increase in hydrogen consumption in the beta zone by 0.18 mmolH2/g compared to the base NiCeZr material. This increase was attributed to the concomitant reduction of CuO entities in this temperature region. In the Sn/NiCeZr solid, there was a slight increase in hydrogen uptake in the β- and γ-peaks, which was attributed to the reduction of SnO2 to metallic Sn.
In summary, the addition of Cu promoted the reducibility of Ni in NiCeZr catalysts, resulting in a lower temperature of the γ-contribution. However, Mg seems to hinder the reduction of both Ce4+ and Ni2+, possibly by altering the Ni2+ distribution and reducibility. The presence of Sn increases the reduction temperatures, suggesting a possible hindrance in the reduction pathways of surface NiO, Ni2+ in the lattice, and Ce4+. Each promoter clearly affects the reduction profile of the catalyst.
The OSC of the base NiCeZr was assessed by the O2/CO pulse injection method at 235 °C and accounted for 63 mmolO/gcat, as shown in Table 2. This value was attributed to the Ce3+/Ce4+ and oxygen vacancies of the catalyst [44]. The OSC varied little upon the addition of Mg and Cu, but notably decreased by 32% upon impregnating Sn, in accordance with the hindered reducibility of this sample.

2.4. Characterization of Acid and Metallic Sites

An in-depth analysis was conducted to characterize both metallic and acid sites in the Me/NiCeZr catalysts. The density of the acid sites was significantly influenced by the promoter (Table 2). Specifically, the acid sites’ density of the catalyst modified with MgO decreased by 17% with respect to the undoped solid, underscoring the basic nature of alkaline earth metals [35]. Conversely, the catalysts modified with Cu and Sn notably increased their surface acid sites’ density (increases by 17% and 42% for Cu/NiCeZr and Sn/NiCeZr, respectively), in agreement with the acidic characteristics of these modifiers [45,46].
The desorption profiles of ammonia in Figure S2 (Supplementary Materials) indicate a strong effect of the modifier. The profiles consist of various peaks at different temperatures, which are indicative of acid sites of different strengths. The assignment of the peaks to their respective strengths was performed according to [26]. The contribution of each acid-strength site is summarized in Table S2 (Supplementary Materials). The impregnation of Me primarily resulted in a slight increase in the contribution of intermediate strength acid sites and a concomitant decrease in strong sites. The type of dopant, however, does not significantly affect the percentage contribution of each center. In these types of solids, acidity originates from the Lewis type, with coordinatively unsaturated cations on their surfaces [47].
The density of metallic sites was determined through hydrogen chemisorption. The results presented in Table 2 show a decreasing trend in metal site density when NiCeZr was doped with various elements. The undoped NiCeZr catalyst had 0.22 metallic sites per nm2. After being doped with Mg, the density of metal sites decreased by 32%, which is consistent with the H2-TPR results and in agreement with other studies [48]. Doping with Cu and Sn resulted in a significant decrease in metallic site density (around 90%). In the case of the Cu/NiCeZr assay, the decrease was attributed to the formation of the NiCu alloy on the catalyst surface, as observed from the XRD results [49,50]. For Sn/NiCeZr, the observed effect could be attributed to the growth of the metallic nickel crystallite (doubled with respect NiCeZr) and the presence of SnOX, which is inert for H2 chemisorption, covering the metallic nickel, as observed through H2-TPR.

2.5. Surface Characterization Using XPS Analysis

Distinctive features of both Ce3+ and Ce4+ were identified in the Ce 3d core-level spectra of all calcined samples (Figure S3A, Supplementary Materials). Two doublets labeled as v′/u′ and v0/u0 were observed for Ce3+, while three doublets named v/u, v″/u″, and v′′′/u′′′ were detected for Ce4+, corresponding to the Ce 3d5/2 and Ce 3d3/2 levels, with an energy separation of 18.3–18.6 eV between both signals, in agreement with the literature [51]. Ce4+ was found to be the predominant form of cerium for all the catalysts (Table 3).
The analysis of the Ni 2p3/2 core-level spectra (Figure 4) revealed a consistent pattern, consisting of a main peak in the range of 851–859 eV and a satellite peak at around 862 eV, both indicative of the Ni2+ species [52]. These results confirmed that the calcined solids contained only Ni2+. The spectrum of Mg/NiCeZr was found to be identical to that of NiCeZr. However, the main peak shifted slightly towards a lower binding energy (−0.1 eV) for the Cu/NiCeZr catalyst. This shift could be associated with changes in the unfilled electronic holes of the D band, resulting from the charge transfer from Cu to the adjacent Ni [53]. On the other hand, the main peak shifted up by 0.2 eV for the Sn/NiCeZr catalyst, which was likely due to the higher electronegativity of Sn (1.96) compared to that of nickel (1.91) [54].
The Zr 3d core-level spectra for all calcined solids (Figure S4 in the Supplementary Materials) presented a doublet at approximately 183.05 eV and 185.5 eV and corresponded to the Zr 3d spin-orbit splitting components (Zr 3d5/2 and Zr 3d3/2) of Zr4+ [55,56]. Upon the impregnation of Me, the Zr 3d peaks upshifted in their BE, indicating that the chemical nature of the Zr4+ species in these materials differed from that in the parent NiCeZr [57]. Among the three dopants used, Sn caused the largest shift, suggesting that the strongest interaction occurs between Zr and Me.
The O 1s spectrum of the calcined catalysts exhibited three distinct peaks (Figure S5, Supplementary Materials). The peak located at around 529.5 eV corresponded to the oxygen in the lattice (OL); the peak at 531.0 eV was linked to oxygen species adsorbed on the vacancies (OV); and the peak at 533 eV could be attributed to chemisorbed species, such as water, carbonate, or OH compounds (Oads) [58]. Ce-containing materials often have Ce3+ ions associated with surface oxygen vacancies, which affect the oxygen storage capacity (OSC) [59]. Our catalysts followed this linear correlation (Figure S6, Supplementary Materials).
Figure 5 displays the Mg 2s, Cu 2p3/2, and Sn 3d5/2 features of calcined Me/NiCeZr catalysts. The Mg 2s core level peaked at 90.2 eV, which is slightly higher than that of MgO [60]. This could be attributed to Mg2+ species interacting with their surroundings, which is in agreement with H2-TPR. The Cu 2p3/2 spectrum of Cu/NiCeZr consisted of a main peak at 934.2 eV, accompanied by intense shake-up satellite peaks in the range of 939–946 eV. Both features show Cu2+ species [61]. The BE of the main peak in the Cu/NiCeZr catalyst was slightly higher than that of standard CuO [62], indicating a strong interaction between Cu and its surroundings. The absence of a secondary peak near 932.8 eV, typically associated with low-valence copper species such as Cu1+, suggests that the copper species on the surface were only present as Cu2+ [63]. The Sn 3d5/2 spectrum of the calcined Sn/NiCeZr catalyst showed a major contribution at 487.4 eV, corresponding to Sn4+ species, and an additional weak peak at 485.2 eV, corresponding to less oxidized tin species. The vacuum in the XPS chamber likely generated the formation of partially reduced tin. The Sn 3d5/2 peak of the Sn/NiCeZr catalyst shifted slightly upwards compared to standard SnO2 [64]. This indicates a strong interaction between Sn and its surroundings, which is consistent with the XRD and H2-TPR results.
For most of the elements, the core-level profiles of the reduced catalysts were notably different from those of the calcined catalysts. For example, the Ce 3d spectra of all catalysts (Figure S3B, Supplementary Materials) only showed Ce3+ features, indicating a complete reduction of Ce4+ ions to Ce3+, consistent with the H2-TPR results. Upon reduction, the Ni 2p spectra of the NiCeZr catalyst (Figure 4B) showed a major peak at 852.5 eV, indicating the presence of metallic Ni [65]. The presence of remnant Ni2+ in the form of the NiCeZr solid solution, with hindered reducibility, is indicated by the weak peak at around 857 eV. The intensity of the Ni0 peak was considerably higher than that of Ni2+, indicating a significant surface nickel reduction. Similar profiles were obtained for the reduced Me/NiCeZr catalysts after reduction. The intensity of the Ni2+ peak varied among the catalysts. It increased for the Mg and Sn-doped catalysts and decreased for the Cu/NiCeZr catalyst, which is consistent with the H2-TPR results. The Zr 3d signal in the reduced samples remained largely unchanged, reflecting its hard reducibility. The O 1s spectra (Figure S5B, Supplementary Materials) exhibited two changes upon reduction. The oxygen binding energy experienced a downshift of 0.2–0.4 eV, indicating the polarization of oxygen–metal bonds. Furthermore, the peak area of Oads decreased, which may be attributed to the desorption of chemisorbed species resulting from the reduction treatment.
The spectrum of Mg 2s in the calcined and reduced Mg/NiCeZr catalysts was identical, indicating that MgO was not reduced under the experimental conditions. The Cu 2p3/2 spectrum of the reduced Cu/NiCeZr catalyst indicates the disappearance of signals corresponding to Cu2+ and the appearance of a new signal at 932 eV, which is associated with Cu0 [66]. The Sn 3d5/2 spectrum of the Sn/NiCeZr catalyst after reduction showed two distinct peaks at approximately 484.2 and 487.1 eV. These peaks were identified as metallic Sn0 and Snδ+, respectively [64]. Based on the larger area of metallic Sn compared to oxidized Sn, it can be inferred that the majority of Sn was in a metallic state. The persistent presence of the Snδ+ signal was attributed to incomplete reduction of SnO2 at 600 °C.
In addition to oxidation state information obtained from BEs, the relative surface concentration of elements can be calculated using their signal intensities. Table 3 shows that the surface of the reference NiCeZr catalyst was enriched in Zr in the calcined forms (bulk Ce/Zr: 0.164), which is consistent with the literature [67]. The study found that the surface enrichment in Zr was lower for the doped solids. The results suggest that doped metals substituted Zr on the surface, which is consistent with the increase in the a-axis lattice length observed from the XRD results. The surface atomic ratio of Ni/(Ce + Zr) in the calcined NiCeZr solid was lower than that in the bulk (0.232). However, the addition of a promoter (Mg, Cu, or Sn) increased the ratio, surpassing the bulk composition. The surface concentration of Ce3+ increased with the addition of a promoter. A charge imbalance occurred, leading to the formation of oxygen vacancies, structural defects, and unsaturated chemical bonds on the catalyst’s surface [68]. The surface impregnation method used for its addition resulted in the atomic content of each promoter on the surface being two to three times that of the bulk, as expected.
Upon reduction, all surface cerium was reduced in accordance with the literature [69,70]. Additionally, the Ce/Zr ratio on the surface increased, which is consistent with the literature [71]. However, all catalysts exhibited a decrease in the Ni/(Ce + Zr) ratio, which was attributed to the decoration of Ni particles by cerium [47]. In the reduced materials, the atomic content of each promoter on the surface was preserved.

2.6. Hydrogenolysis of Glycerol

The performance of the synthesized Me/NiCeZr catalysts was evaluated for glycerol hydrogenolysis in the aqueous phase without the addition of external hydrogen. The presence of gas-phase H2 or hydrogenolysis products would imply the production of H2 by sacrificing part of the fed glycerol. The NiCeZr catalyst achieved the highest glycerol conversion rate (Figure 6), with most of the reacted carbon obtained in the liquid phase (Yliq > Ygas). The addition of any promoter resulted in lower glycerol conversion and reduced liquid yield. The type of promoter had a significant but distinct impact on both outputs. For instance, the Xgly decreased as follows: NiCeZr (75.3%) > Mg/NiCeZr (71.3%) > Cu/NiCeZr (64.7%) > Sn/NiCeZr (54.9%). However, Yliq followed the trend: NiCeZr (56.3%) > Cu/NiCeZr (51.6%) > Sn/NiCeZr (44.1%) > Mg/NiCeZr (37.6%). There was a clear correlation between the density of the acid sites (Table 2) and the percentage of Yliq in relation to the total carbon yield (Yliq + Ygas), irrespective of the glycerol conversion, in line with others [72].
Mg increased the gas yield (30.8%), though with a corresponding decrease in glycerol conversion and liquid yield. It was found that the addition of MgO, by providing basic sites and altering the electronegativity of Ni in the catalyst, favors hydrogen formation and reduces the likelihood of feed dehydration, translating into improved selectivity towards gas production rather than liquid products [73,74,75].
The reaction scheme for glycerol hydrogenolysis using metal–acid bifunctional catalysts involves three primary pathways: dehydration, dehydrogenation, and hydrogenation reactions [76] (Scheme 1). Hydroxyacetone is formed when the primary OH is cleaved, leading to the cleavage of the C–O bond and the creation of a C=O bond (Path A). Alternatively, 3-hydroxypropanal is formed when the secondary OH is cleaved (Path B). Alternatively, glyceraldehyde is formed by breaking the O–H bond (Path C). Path A involves dehydrating the terminal hydroxyl on acid sites to hydroxyacetone (HA). Subsequent hydrogenation on metal sites produces 1,2-propylene glycol (1,2-PG). Further dehydration and hydrogenation steps can produce mono-alcohols such as 1-propanol and 2-propanol. Path B starts with the dehydration of the secondary hydroxyl to 3-hydroxypropanal, which then undergoes hydrogenation to form 1,3-PG. This pathway is preferred on Brönsted acid sites, while Lewis acid sites tend to lead to the path A products. Finally, Path C involves secondary C-C breaking reactions that are responsible for the formation of C2 and C1 products. Path C consists of sequential dehydrogenation and decarbonylation (C–C scission) on metal sites, which produces CO and ethylene glycol (EG) [77]. The pathway could also proceed by converting EG into other products, such as ethanol, methanol, or acetaldehyde, through the dehydration and hydrogenation of ethylene glycol. The intermediates from all paths can be further hydrogenated to produce a diverse range of products. The Water–Gas Shift (WGS) and hydrogenation can convert the carbon monoxide into carbon dioxide or methane.

2.6.1. Gas Products

The analysis of gas product yields shows that the catalyst used has a significant impact (Table 4). The reference NiCeZr catalyst produced the highest yield of CO2 at 11.15%, followed by CH4 at 4.3% and H2 at 2.3%. CO and C2+ alkanes were less abundant products, with yields of 0.65% and 2.15%, respectively. The addition of Mg not only enhanced the yields of all gases but also resulted in the highest total Ygas recorded for this catalyst. This was particularly evident in the performance of YCH4, which experienced a significant increase to 9%. However, the modification with Cu reduces the yields of both H2 and CH4, slightly increasing those of CO and C2+. When using the Sn/NiCeZr catalyst, there was a significant reduction in the yield of all gases, with CH4 dropping to below 1%.
For all catalysts, the molar ratio of H2/CO2 was significantly below 1.0, which is the theoretical stoichiometric ratio for glycerol reforming (2.33). This suggests that the hydrogen produced was consumed in secondary reactions, such as the hydrogenation of liquid intermediates or the methanation of CO and CO2.
The hydrogen selectivity of the modified Me/NiCeZr catalysts was lower than that of the reference catalyst, despite the anticipated synergistic effects often associated with bimetallic catalyst systems [25]. It is noteworthy that the addition of Mg resulted in unexpectedly low hydrogen selectivity, despite the catalysts’ basicity promoting H2 production [78]. The catalyst’s high selectivity for methane suggests that the produced H2 was consumed in the reactor for CO/CO2 hydrogenation. In contrast, the catalysts modified with Cu and Sn exhibited similar hydrogen selectivity, which was associated with lower methane production. It is important to note that selectivity towards H2 is calculated based on hydrogenated products, while selectivity towards C compounds is calculated based on carbon. As a result, the sum of selectivities for H2 and all C compounds exceeds 100%.
Previous studies have reported varying effects on methane production when using different modified NiCeZr catalysts. The literature suggests that bimetallic interactions in Cu catalysts can decrease CH4 selectivity [79], while Sn modification may impact selectivity due to changes in Ni defect sites [80]. Additionally, Cu/NiCeZr and Sn/NiCeZr catalysts showed an increase in CO selectivity, which could indicate a subdued WGS reaction. The results are intriguing as Cu-modified catalysts are typically known to promote the WGS reaction [37]. It is possible that the high density of acid sites on both catalysts could hinder their WGS activity [74]. This is supported by the increase in selectivity for C2+ alkanes for both catalysts.

2.6.2. Liquid Products

As illustrated in Figure 6, a significant portion of the carbon feed is acquired in liquid form. Figure 7 displays the selectivity towards primary liquid products. The primary liquid products obtained were 1,2-PG, HA, ethanol, and EG, which accounted for over 95% of the total. Other minor products included methanol, acetone, and 1-propanol. The highest selectivity was found for 1,2-PG, with no notable differences among them (around 50%), except for Sn/NiCeZr catalysts. However, the selectivity to 1,2-PG of Sn/NiCeZr decreased to 37.9%, at the expense of an increase in selectivity for HA. The formation of 1,2-PG (though the hydrogenation of HA) may be suppressed due to the lower density of metal sites in the catalyst. The selectivity towards HA was increased by the bimetallic Cu/NiCeZr and Sn/NiCeZr catalysts, likely due to their higher acidity that enhances HA formation by dehydrating glycerol to HA. However, the Mg/NiCeZr catalyst, which potentially has a more electronegative Ni, showed lower dehydration activity, resulting in reduced HA selectivity (21.7%). The dehydration of glycerol produces HA, which is then hydrogenated to form 1,2-PG (Scheme 1). The difference in product distribution suggests that the metallic sites in NiCeZr and Mg/NiCeZr had a higher hydrogenation capacity.
There was a balance between the production of EG and ethanol, both products from the dehydrogenation path in Scheme 1. The formation of ethanol through the dehydration and subsequent hydrogenation of EG suggests enhanced hydrogenation activity at the metallic sites of NiCeZr and NiCeZr-Mg catalysts. This highlights the bifunctional nature of these catalysts, where a synergistic interplay between acid and metallic sites plays a crucial role in substrate transformation.
The liquid products showed selectivity towards the HA and 1,2-PG pathway, rather than the EG and ethanol pathway (Table 5). The promotion with Cu and Sn further enhanced this preference, which is consistent with their higher acid sites’ density.
For all catalysts studied, primary products (single C–O or C–C bond scission) were more prevalent than secondary products, with selectivity for primary products exceeding 75%. The undoped catalyst exhibited the highest selectivity for secondary products, while for the Cu- and Sn-doped catalysts, it decreased significantly, due to their very low metal availability.
The analysis of selectivity focused on the breakage of C–C and C–O bonds and revealed a predominance of C–O scission, with selectivity surpassing 92% (Table 5). The Cu/NiCeZr and Sn/NiCeZr catalysts exhibited the highest selectivity for C–O bond cleavage. The acidic sites had a significant impact on the competitive nature of C–C and C–O bond cleavages. In summary, the liquid product had a lower O/C ratio than the feedstock, indicating a lower oxygen content in the products. The catalysts that showed higher selectivity to C–O scission also demonstrated greater hydrodeoxygenation capacity.
Table S3 (Supplementary Materials) shows a comparison between the results of our catalysts with other works in the literature regarding the production of value-added products from glycerol without external hydrogen. Our catalysts showed high glycerol conversion and selectivity towards 1,2-propylene glycol, comparable or superior to previous studies. This achievement is particularly noteworthy as it was accomplished without the addition of external hydrogen, indicating a significant advantage in terms of process efficiency and sustainability.
The undoped NiCeZr catalyst stability studies (Figure S7, Supplementary Materials) demonstrate a gradual decrease in glycerol conversion over 30 h TOS, indicating a decline in catalytic efficiency for glycerol processing. The yield for liquid products also decreased, while the yield for gas remained constant. The selectivity for 1,2-PG may vary due to the dynamic nature of the catalytic sites during the reaction or a balance between the production of gaseous and liquid products. Nickel leaching (see below) is the most likely cause of the decrease in glycerol conversion and liquid product yield, resulting in the loss of active metal sites. Overall, the evidence points to nickel leaching as the primary factor in the decrease in glycerol conversion and liquid product yield. This loss may alter the acidity of the catalyst surface, potentially affecting the selectivity towards 1,2-PG. However, the steady production of gaseous byproducts suggests that these may arise from reaction pathways that are less affected by the nickel content and acidity of the catalyst. The breakdown of stable intermediates or byproducts that are not primarily formed through nickel’s catalytic action may be the cause.

2.7. Characterization of Spent Catalysts

Designing a catalyst for aqueous-phase biomass treatment applications presents numerous challenges, including deactivation due to metal leaching, phase transformations, or metal sintering, among other factors [17]. To explore the changes experienced by the catalysts, spent catalysts were examined.
The N2 adsorption–desorption isotherms of all the spent catalysts (Figure S8, Supplementary Materials) were of type IV, as for the reduced assays. However, the isotherms of the base NiCeZr and Mg/NiCeZr spent catalysts showed a hysteresis loop closure at slightly higher relative pressures than in the reduced materials (Figure S8 in Supplementary Materials), while for the spent Cu/NiCeZr this occurred at lower relative pressures. These results suggest that the desorption of the adsorbate from the pores is more challenging, possibly because of surface modifications or pore blockages that occurred during use. The PSD for the spent reference catalysts widened, while the PSD for the Mg/NiCeZr catalysts narrowed and shifted towards lower pores. Minor changes were observed for Cu/NiCeZr and Sn/NiCeZr catalysts. The SBET increased by 15–25% for all spent catalysts except for the Mg/NiCeZr catalyst, which showed a decrease. It seems that the BET surface area variation was significantly affected by the leaching of Ni and Me (see below). In summary, the texture of all catalysts changed under hydrothermal conditions. Cu and Sn were found to be the promoters that best preserved the textural properties.
The XRD profiles of the spent catalysts (Figure 8A) were similar to their respective reduced forms. However, upon closer inspection, it was discovered that there were changes in the peaks’ intensity and width when compared to their reduced forms. For example, the intensity of the Ni0 peak noticeably increased for the reference NiCeZr after usage, and a new peak at 51.8° emerged. Both features indicate an enlargement in the crystal size of metallic nickel. The spent Mg/NiCeZr catalyst showed a new peak from MgO (ICDD #00-045-0946), accompanied by a decrease in the intensity of the Ni0 peak, indicating a reduction in the crystal size of metallic nickel. The XRD profile of the spent Cu/NiCeZr catalyst revealed the absence of the Ni0 peak, leaving only the peak corresponding to the Ni-Cu alloy. In the case of Sn/NiCeZr, the intensity of the Ni0 peak slightly decreased. The study showed that the dopants were effective in preserving the metallic nickel at the core of the catalyst, despite potential surface oxidation [76]. The size of the nickel crystallites doubled in the spent NiCeZr catalyst compared to its reduced form and showed little variation for the Mg- and Sn-doped catalysts (Table 6). For Cu/NiCeZr catalysts, the size increased by 50%. These results confirm that the impregnation of Mg and Sn effectively prevented nickel coalescence during the reaction. Despite these changes in the size of the metallic nickel, the main XRD feature for the ceria–zirconia solid solution remained consistent, indicating no phase variation.
The Raman analysis conducted on the carbon region of the catalysts after the reaction (Figure 8B) revealed the presence of D and G bands, specifically for the Sn/NiCeZr catalyst. These bands indicate the formation of carbon-based residues that were not detected with the XRD analysis. The D band signals amorphous carbon types, whereas the G band points to graphitic forms [81]. The low intensity of these bands suggests a minimal presence of carbon deposits, except for the Sn-doped assay, which has the lowest oxygen vacancies and highest acidity. This reduced accumulation of carbon was attributed to the role of oxygen vacancies within the catalyst, which helps to mitigate carbon deposition [82].
The leaching of the active phase can cause irreversible catalyst deactivation, which is a common problem in hydrothermal environments. The NiCeZr catalyst showed the highest Ni and Ce leaching, which notably decreased upon doping, specifically with Cu and Sn. Among the dopants, only Mg experienced leaching, at values higher than 40%. The detachment of Mg could be promoted by the formation and subsequent elution of Mg organometallic complexes. Zr was revealed to be the most resistant to leaching. These results show that the dopants can help stabilize nickel particles under harsh hydrothermal conditions.
Accordingly, an in situ treatment with H2 could remove the coke and reduce the Ni surface that has been oxidized in the operation. However, the reduction treatment could only partially regenerate the catalyst due to the irreversible catalyst deactivation caused by metal leaching. Alternative strategies that aim to stabilize the metal nanoparticles against leaching could hinder catalyst deactivation.

3. Materials and Methods

3.1. Materials

3.1.1. Chemicals

The catalysts were synthesized using the following reagents: cerium(IV) nitrate hexahydrate (Ce(NO3)3·6H2O, Mn = 434.22, 99%), zirconium(IV) nitrate hydrate (ZrO(NO3)2·xH2O, Mn = 231.23 anhydrous basis, 99%), nickel nitrate hexahydrate (Ni(NO3)2·6H2O, Mn = 290.79, 99.99%), cetyltrimethylammonium bromide (CTAB) (C16H33N(CH3)3Br, Mn = 364.45, 99%), sodium hydroxide (NaOH, 98%), magnesium nitrate hexahydrate (Mg(NO3)2·6H2O, Mn = 256.4, 99.999%), copper(II) nitrate trihydrate (Cu(NO3)2·3H2O, Mn = 187.6 anhydrous basis, 99-104%), purchased from Sigma-Aldrich (Saint Louis, MO, USA), tributyltin acetate (C14H30O2Sn, Mn = 349.1, 97.0%), purchased from Merck (Darmstadt, Germany), isopropyl alcohol ((CH3)2CHOH, Mn = 60.10, 99.8%), purchased from PanReac (Castellar del Vallés, Spain), and deionized water.

3.1.2. Synthesis of Catalysts

The procedure for synthesizing the base NiCeZr catalyst (10 wt.%Ni, Ce/Zr = 15/85 at./at.) can be found elsewhere [26]. In brief, it was synthesized using a one-pot co-precipitation method with cetyltrimethylammonium bromide (CTAB) as a surfactant, and calcined for 4 h at 500 °C (heating ramp of 1 °C/min). The modifiers were added to NiCeZr using a wet impregnation method with aqueous solutions (50:50 isopropanol:water mixture for Sn) of the corresponding metal salt (at 1 wt.% concentration). The solid product was dried at 100 °C overnight and calcined at 500 °C for 4 h using a heating ramp of 1 °C/min. The modified solids were named Me/NiCeZr (Me = Mg, Cu or Sn).

3.2. Characterization of the Catalysts

The chemical composition of the catalysts was analyzed using inductively coupled plasma atomic emission spectroscopy (ICP-AES) (7700, Agilent, Santa Clara, CA, USA), following a standard acid digestion procedure. To quantify the metals leached into the final liquid product, inductively coupled plasma mass spectrometry (ICP-MS) (XSeries 2, Thermo Scientific, Waltham, MA, USA) was employed. The textural properties of the solids were evaluated by N2 adsorption–desorption isotherms at 77 K (TRISTAR II 3020, Micromeritics, Norcross, GA, USA). Before adsorption measurements, the samples were outgassed at 400 °C for 10 h to eliminate moisture and adsorbed gases. The BET method was used to determine the specific surface area, while the BJH (adsorption-branch) method was used to obtain the average pore size.
The powder X-ray diffraction (XRD) patterns of the solids were obtained using a PANalytical X’Pert-Pro instrument, covering a range of 5 to 90° 2θ with a step size of 0.026° and a duration of 598 s, and employing monochromatized CuKα radiation (λ = 1.5418 Å). The crystallite size was determined by applying the Scherrer equation. The identification of crystalline phases was conducted by comparing them with the ICDD database. The degree of particle agglomeration was estimated by calculating the ratio of the average particle volume, assumed to be spherical based on the SBET, to the crystallite volume obtained from the XRD results.
The Raman spectra were recorded using a Renishaw InVia Raman spectrometer coupled with a Leica DMLM microscope. The analysis employed a 514 nm ion argon laser (Modu-Laser) and a holographic grating of 1800 lines/mm. Each spectrum was acquired over a 20 s period with 10 accumulations, ensuring the laser power remained consistently below 2 mW. The measurements were conducted within a spectral window of 150–1500 cm−1.
The reducibility of the calcined solids was evaluated through a temperature-programmed reduction (H2-TPR) technique using a Micromeritics AutoChem 2920 apparatus. The samples were initially cleaned by purging with He at a flow rate of 50 mL/min at 500 °C for 30 min, followed by cooling to 40 °C. The H2-TPR analysis involved flowing a 5% H2/Ar gas mixture over the samples and incrementally increasing the temperature from room temperature to 900 °C at a rate of 10 °C/min. The Thermal Conductivity Detector (TCD) continuously monitored the outflow from the reactor, while the resulting water was trapped.
The oxygen storage capacity (OSC) was evaluated using a Micromeritics AutoChem 2920 instrument connected to a Pfeiffer Vacuum OmniStar mass spectrometer. After surface cleaning by heating at 500 °C in He flow for 30 min, the catalyst was reduced for an hour at 600 °C in a 5% H2/Ar flow. Then, the gas flow was switched to He and the sample cooled down to 235 °C. At this temperature, the samples experienced ten pulses of oxygen 5% O2/He, followed by ten pulses of 5% CO/He. The OSC was determined as the sum of the released oxygen moles after each CO pulse.
The metallic surface area was further characterized using static H2 chemisorption and analyzed using a Micromeritics ASAP 2020 instrument. Initially, the sample was subjected to surface cleaning under He flow at 350 °C to remove any adsorbed species. Subsequently, the temperature was reduced to 50 °C, followed by H2 reduction at 600 °C with a flow rate of 50 mL/min for 1 h. To eliminate any residual H2, the sample was then degassed for 1 h. The chemisorption process involved a double H2 adsorption isotherm conducted at 5 °C, incorporating an intermediate degassing step. This temperature was specifically chosen to reduce the occurrence of H2 spillover. The quantitative assessment of the chemisorbed hydrogen was determined by the difference between the first and second isotherms, providing insights into the metallic surface, metal sites, and dispersion.
The acidity of the catalysts was assessed by ammonia chemisorption followed by temperature-programmed desorption (NH3-TPD) (AutoChem 2920, Micromeritics, Norcross, GA, USA). The catalyst’s surface was initially cleaned by heating it to 500 °C under a stream of He for 30 min. Following this, the catalyst was reduced at 600 °C for an hour in a 5% H2/Ar stream, and then cooled to 90 °C under a flow of He. Then, six pulses of 10% NH3/He (loop volume 0.5312 mL) were introduced to quantify the chemisorbed ammonia. Subsequently, He was flowed for one hour to evacuate the system. Then, the sample’s temperature was gradually increased to 900 °C at a rate of 10 °C/min, following the reactor exhaust using the TCD. The number of acid sites was determined by integrating the areas of the NH3 pulses, while the acid strength was inferred from the desorption temperatures.
The oxidation states of the surface metals in both the calcined and reduced forms of the catalysts were investigated using X-ray photoelectron spectroscopy (XPS). The analysis was conducted on a SPECS spectrometer, which is equipped with a Phoibos 150 1DDLD analyzer and utilizes monochromatic Al Kα X-ray sources (1486.7 eV). The pass energy was finely tuned in increments of 0.05 eV to 30 eV. For samples requiring in situ reduction, the procedure was carried out at 600 °C for 1 h. Spectrometer calibration was achieved using the Ag 3d5/2 peak at 368.26 eV, and the binding energy was referenced to the C 1s peak of adventitious carbon at 284.6 eV. Peak deconvolution was performed after Shirley background subtraction, employing a combined Gaussian–Lorentzian function. This analysis was facilitated by the use of CASA XPS software, ensuring accurate interpretation of the data.

3.3. Catalytic Tests

Catalytic testing for glycerol hydrogenolysis was conducted in a bench-scale fixed-bed upflow reactor (Microactivity Effi, PID Eng&Tech, Alcobendas, Spain) using a 10 wt.% glycerol (99.5%, PanReac, Castellar del Vallés, Spain) aqueous solution. The Electronic Supplementary Information includes a detailed schematic representation of the experimental setup (Figure S9). The reaction conditions were maintained at 235 °C and 35 bar, with a WHSV of 12 h−1. The catalyst underwent an in situ reduction at 600 °C for 1 h under a 10% H2/He flow before the experiments. The total on-stream time (TOS) for the reaction was 3 h. Product separation into gaseous and liquid phases was achieved using a Peltier device. The gaseous output was continuously swept with a He flow of 40 mL/min after the backpressure regulator and analyzed online via microGC-TCD (Agilent 490), equipped with Al2O3-KCl, PPQ, and MS5A columns using He as a carrier, and MS5A column utilizing Ar. The liquid samples, collected hourly, were analyzed off-line using a GC-FID system (Agilent 6890N with a DB-Heavy Wax column). The carbon content in the liquid phase was determined using a Shimadzu TOC-L apparatus, ensuring a carbon balance exceeding 90% in all experiments.

3.4. Equations Used for Calculations for Catalytic Results

Glycerol conversion (Xgly) was determined using Equation (1):
X g l y ( % ) = 100 × F g l y i n F g l y o u t F g l y i n
where F g l y i n and F g l y o u t are the molar flow of glycerol at the reactor inlet and outlet, respectively. The following methodology was used to calculate gas and liquid product yields based on carbon content (Equations (2) and (3)):
Y g a s ( % ) = 100 × F C , g a s o u t 3 × F g l y i n
Y l i q ( % ) = 100 × F C , l i q o u t 3 × F g l y i n
F C , g a s o u t and F C , l i q o u t are the total molar flow of carbon in the outlet gas and liquid stream, respectively (excluding glycerol). These metrics range between 0 and 100%.
For the carbon-containing products, selectivity (Si) and yield (Yi) for product i were determined based on its carbon content, as follows (Equations (4) and (5)):
S i % = 100 × F C . i o u t i = 1 n F C , i o u t
Y i % = 100 × F C . i o u t 3 × F g l y i n
where F C , i o u t is the carbon flow in the product i. These metrics range from 0% (no production of C-containing i) to 100% (only i is obtained as a C-containing product).
The selectivity to hydrogen was calculated as the percentage of H2 in relation to the total hydrogen released in all gaseous products (Equation (6)).
S H 2 % = 100 × 2 × F H 2 o u t F H , g a s o u t
The hydrogen yield (YH2) represented the ratio of produced H2 to the ideally produced by glycerol APR (Equation (7)):
Y H 2 % = 100 × 1 7 × F H 2 o u t F g l y i n
where F H 2 o u t and where F H , g a s o u t are the H2 and total H flow in the outlet gas, respectively. These metrics range from 0% (no H2 exhaust) to 100% (H2 being the only H-containing molecules in the gas phase).

4. Conclusions

The investigation aimed to improve the applicability of NiCeZr catalyst doping at the surface level with Mg, Cu, and Sn, for the sustainable valorization of biomass through glycerol hydrogenolysis. Modifications slightly expanded pore sizes for Cu- and Sn-doped catalysts, but all doped variants experienced a decrease in surface area, particularly for Cu/NiCeZr. The redox behavior and oxygen storage analyses demonstrated the distinct effects of each dopant. Mg had a slight impeding effect on the reduction temperatures of the catalysts, while Cu facilitated and Sn increased the temperatures. The structural analyses showed the impact of dopants on catalyst properties. This also influenced the redox properties and oxygen vacancies of the catalysts.
The efficiency and selectivity of glycerol conversion were significantly impacted by the modifications. It was observed that dopants had a significant effect on the efficiency and selectivity of glycerol conversion in NiCeZr catalysts. The undoped catalyst had a conversion rate of 75.3%, with a selectivity rate of 93.5% for the C–O bond cleavage. However, after doping, the conversion rate slightly decreased, and the selectivity for cleaving C–O bonds increased for the Cu and Sn variants. Furthermore, the introduction of dopants led to an increased selectivity for cleaving C–C bonds. These results can be correlated with the density of acid and metal sites, as well as oxygen vacancies.
In post-reaction assessments, the spent catalysts were found to be structurally and functionally stable. The dopants were effective in preserving the metallic nickel core and mitigating sintering effects. However, the Mg/NiCeZr catalyst exhibits significant Mg leaching, which needs to be addressed for further optimization. The findings show the complex connection between changes in structure and catalytic performance, providing useful insights for creating effective and durable catalysts for biomass valorization.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25063484/s1. References [83,84,85,86,87,88,89,90] are cited in the Supplementary Materials.

Author Contributions

N.V.-H.: investigation, writing—original draft, and data curation; U.I.-V.: formal analysis and writing—review and editing; J.L.A.: funding acquisition, conceptualization, writing—review and editing, data curation, and project administration; M.Á.G.-O.: resources, funding acquisition, supervision, and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grant PID2019-106692RB-I00 funded by MCIN/AEI/10.13039/501100011033. N.V.H. would like to thank MICINN for the pre-doctoral grant (BES-2017-081580). Likewise, the authors thank for technical support provided by SGIker of UPV/EHU and European funding (ERDF and ESF).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. de Wild, P.J.; den Uil, H.; Reith, J.H.; Kiel, J.H.A.; Heeres, H.J. Biomass Valorisation by Staged Degasification. A New Pyrolysis-Based Thermochemical Conversion Option to Produce Value-Added Chemicals from Lignocellulosic Biomass. J. Anal. Appl. Pyrolysis 2009, 85, 124–133. [Google Scholar] [CrossRef]
  2. Mishra, S.; Upadhyay, R.K. Review on Biomass Gasification: Gasifiers, Gasifying Mediums, and Operational Parameters. Mater. Sci. Energy Technol. 2021, 4, 329–340. [Google Scholar] [CrossRef]
  3. Gao, Y.; Remón, J.; Matharu, A.S. Microwave-Assisted Hydrothermal Treatments for Biomass Valorisation: A Critical Review. Green Chem. 2021, 23, 3502–3525. [Google Scholar] [CrossRef]
  4. Nzeteu, C.; Coelho, F.; Davis, E.; Trego, A.; O’Flaherty, V. Current Trends in Biological Valorization of Waste-Derived Biomass: The Critical Role of VFAs to Fuel a Biorefinery. Fermentation 2022, 8, 445. [Google Scholar] [CrossRef]
  5. Seo, M.W.; Lee, S.H.; Nam, H.; Lee, D.; Tokmurzin, D.; Wang, S.; Park, Y.K. Recent Advances of Thermochemical Conversion Processes for Biorefinery. Bioresour. Technol. 2022, 343, 126109. [Google Scholar] [CrossRef]
  6. Awasthi, M.K.; Sarsaiya, S.; Patel, A.; Juneja, A.; Singh, R.P.; Yan, B.; Awasthi, S.K.; Jain, A.; Liu, T.; Duan, Y.; et al. Refining Biomass Residues for Sustainable Energy and Bio-Products: An Assessment of Technology, Its Importance, and Strategic Applications in Circular Bio-Economy. Renew. Sustain. Energy Rev. 2020, 127, 109876. [Google Scholar] [CrossRef]
  7. Calvo-Flores, F.G.; Martin-Martinez, F.J. Biorefineries: Achievements and Challenges for a Bio-Based Economy. Front. Chem. 2022, 10, 973417. [Google Scholar] [CrossRef] [PubMed]
  8. Checa, M.; Nogales-Delgado, S.; Montes, V.; Encinar, J.M. Recent Advances in Glycerol Catalytic Valorization: A Review. Catalysts 2020, 10, 1279. [Google Scholar] [CrossRef]
  9. Wan, W.; Ammal, S.C.; Lin, Z.; You, K.E.; Heyden, A.; Chen, J.G. Controlling Reaction Pathways of Selective C–O Bond Cleavage of Glycerol. Nat. Commun. 2018, 9, 4612. [Google Scholar] [CrossRef]
  10. Wawrzetz, A.; Peng, B.; Hrabar, A.; Jentys, A.; Lemonidou, A.A.; Lercher, J.A. Towards Understanding the Bifunctional Hydrodeoxygenation and Aqueous Phase Reforming of Glycerol. J. Catal. 2010, 269, 411–420. [Google Scholar] [CrossRef]
  11. Charisiou, N.D.; Siakavelas, G.; Tzounis, L.; Dou, B.; Sebastian, V.; Hinder, S.J.; Baker, M.A.; Polychronopoulou, K.; Goula, M.A. Ni/Y2O3–ZrO2 Catalyst for Hydrogen Production through the Glycerol Steam Reforming Reaction. Int. J. Hydrogen Energy 2020, 45, 10442–10460. [Google Scholar] [CrossRef]
  12. Dahdah, E.; Estephane, J.; Gennequin, C.; Aboukaïs, A.; Aouad, S.; Abi-Aad, E. Effect of La Promotion on Ni/Mg-Al Hydrotalcite Derived Catalysts for Glycerol Steam Reforming. J. Environ. Chem. Eng. 2020, 8, 104228. [Google Scholar] [CrossRef]
  13. Charisiou, N.D.; Siakavelas, G.; Papageridis, K.N.; Baklavaridis, A.; Tzounis, L.; Polychronopoulou, K.; Goula, M.A. Hydrogen Production via the Glycerol Steam Reforming Reaction over Nickel Supported on Alumina and Lanthana-Alumina Catalysts. Int. J. Hydrogen Energy 2017, 42, 13039–13060. [Google Scholar] [CrossRef]
  14. Gatti, M.N.; Nichio, N.N.; Pompeo, F. Advances for Biorefineries: Glycerol Hydrogenolysis to 1,3-Propylene Glycol. Reactions 2022, 3, 451–498. [Google Scholar] [CrossRef]
  15. Martin, A.; Armbruster, U.; Gandarias, I.; Arias, P.L. Glycerol Hydrogenolysis into Propanediols Using in Situ Generated Hydrogen—A Critical Review. Eur. J. Lipid Sci. Technol. 2013, 115, 9–27. [Google Scholar] [CrossRef]
  16. Nanda, M.R.; Yuan, Z.; Qin, W.; Xu, C. (Charles) Recent Advancements in Catalytic Conversion of Glycerol into Propylene Glycol: A Review. Catal. Rev. 2016, 58, 309–336. [Google Scholar] [CrossRef]
  17. Nakagawa, Y.; Tomishige, K. Heterogeneous Catalysis of the Glycerol Hydrogenolysis. Catal. Sci. Technol. 2011, 1, 179–190. [Google Scholar] [CrossRef]
  18. Vasiliadou, E.S.; Lemonidou, A.A. Kinetic Study of Liquid-Phase Glycerol Hydrogenolysis over Cu/SiO2 Catalyst. Chem. Eng. J. 2013, 231, 103–112. [Google Scholar] [CrossRef]
  19. Wang, Y.; Zhou, J.; Guo, X. Catalytic Hydrogenolysis of Glycerol to Propanediols: A Review. RSC Adv. 2015, 5, 74611–74628. [Google Scholar] [CrossRef]
  20. Zhou, J.; Guo, L.; Guo, X.; Mao, J.; Zhang, S. Selective Hydrogenolysis of Glycerol to Propanediols on Supported Cu-Containing Bimetallic Catalysts. Green Chem. 2010, 12, 1835–1843. [Google Scholar] [CrossRef]
  21. Maris, E.P.; Ketchie, W.C.; Murayama, M.; Davis, R.J. Glycerol Hydrogenolysis on Carbon-Supported PtRu and AuRu Bimetallic Catalysts. J. Catal. 2007, 251, 281–294. [Google Scholar] [CrossRef]
  22. Pipitone, G.; Zoppi, G.; Pirone, R.; Bensaid, S. A Critical Review on Catalyst Design for Aqueous Phase Reforming. Int. J. Hydrogen Energy 2022, 47, 151–180. [Google Scholar] [CrossRef]
  23. Dieuzeide, M.L.; Laborde, M.; Amadeo, N.; Cannilla, C.; Bonura, G.; Frusteri, F. Hydrogen Production by Glycerol Steam Reforming: How Mg Doping Affects the Catalytic Behaviour of Ni/Al2O3 Catalysts. Int. J. Hydrogen Energy 2016, 41, 157–166. [Google Scholar] [CrossRef]
  24. Wu, K.; Dou, B.; Zhang, H.; Liu, D.; Chen, H.; Xu, Y. Aqueous Phase Reforming of Biodiesel Byproduct Glycerol over Mesoporous Ni-Cu/CeO2 for Renewable Hydrogen Production. Fuel 2022, 308, 122014. [Google Scholar] [CrossRef]
  25. Shabaker, J.W.; Huber, G.W.; Dumesic, J.A. Aqueous-Phase Reforming of Oxygenated Hydrocarbons over Sn-Modified Ni Catalysts. J. Catal. 2004, 222, 180–191. [Google Scholar] [CrossRef]
  26. Vera-Hincapie, N.; Iriarte-Velasco, U.; Ayastuy, J.L.; Gutiérrez-Ortiz, M.A. Effect of the Synthesis Method on the Performance of Ni-CeO2-ZrO2 Catalysts for the Hydrogenolysis of Glycerol with Endogenous Hydrogen. J. Environ. Chem. Eng. 2023, 11, 111072. [Google Scholar] [CrossRef]
  27. Cangiano, M.D.L.A.; Ojeda, M.W.; Carreras, A.C.; González, J.A.; Ruiz, M.D.C. A Study of the Composition and Microstructure of Nanodispersed Cu-Ni Alloys Obtained by Different Routes from Copper and Nickel Oxides. Mater. Charact. 2010, 61, 1135–1146. [Google Scholar] [CrossRef]
  28. Dhanapal, K.; Narayanan, V.; Stephen, A. Influence of Sn on the Magnetic Ordering of Ni-Sn Alloy Synthesized Using Chemical Reduction Method. J. Magn. Magn. Mater. 2016, 406, 103–109. [Google Scholar] [CrossRef]
  29. Berwal, U.; Singh, V.; Sharma, R. Effect of Ce4+→Ce3+ Conversion on the Structural and Luminescence Properties of Ce4+ Doped Gd2Ti2O7 Pyrochlore Oxide. J. Lumin. 2023, 257, 119687. [Google Scholar] [CrossRef]
  30. Loridant, S. Raman Spectroscopy as a Powerful Tool to Characterize Ceria-Based Catalysts. Catal. Today 2021, 373, 98–111. [Google Scholar] [CrossRef]
  31. Sartoretti, E.; Novara, C.; Giorgis, F.; Piumetti, M.; Bensaid, S.; Russo, N.; Fino, D. In Situ Raman Analyses of the Soot Oxidation Reaction over Nanostructured Ceria-Based Catalysts. Sci. Rep. 2019, 9, 3875. [Google Scholar] [CrossRef]
  32. Shan, W.; Luo, M.; Ying, P.; Shen, W.; Li, C. Reduction Property and Catalytic Activity of Ce1-XNiXO2 Mixed Oxide Catalysts for CH4 Oxidation. Appl. Catal. A Gen. 2003, 246, 1–9. [Google Scholar] [CrossRef]
  33. Yisup, N.; Cao, Y.; Feng, W.L.; Dai, W.L.; Fan, K.N. Catalytic Oxidation of Methane over Novel Ce-Ni-O Mixed Oxide Catalysts Prepared by Oxalate Gel-Coprecipitation. Catal. Lett. 2005, 99, 207–213. [Google Scholar] [CrossRef]
  34. Shi, Q.; Liu, C.; Chen, W. Hydrogen Production from Steam Reforming of Ethanol over Ni/MgO-CeO2 Catalyst at Low Temperature. J. Rare Earths 2009, 27, 948–954. [Google Scholar] [CrossRef]
  35. Siakavelas, G.I.; Charisiou, N.D.; AlKhoori, S.; AlKhoori, A.A.; Sebastian, V.; Hinder, S.J.; Baker, M.A.; Yentekakis, I.V.; Polychronopoulou, K.; Goula, M.A. Highly Selective and Stable Nickel Catalysts Supported on Ceria Promoted with Sm2O3, Pr2O3 and MgO for the CO2 Methanation Reaction. Appl. Catal. B Environ. 2021, 282, 119562. [Google Scholar] [CrossRef]
  36. Saw, E.T.; Oemar, U.; Ang, M.L.; Kus, H.; Kawi, S. High-Temperature Water Gas Shift Reaction on Ni-Cu/CeO2 Catalysts: Effect of Ceria Nanocrystal Size on Carboxylate Formation. Catal. Sci. Technol. 2016, 6, 5336–5349. [Google Scholar] [CrossRef]
  37. Saw, E.T.; Oemar, U.; Tan, X.R.; Du, Y.; Borgna, A.; Hidajat, K.; Kawi, S. Bimetallic Ni-Cu Catalyst Supported on CeO2 for High-Temperature Water-Gas Shift Reaction: Methane Suppression via Enhanced CO Adsorption. J. Catal. 2014, 314, 32–46. [Google Scholar] [CrossRef]
  38. Courtois, X.; Perrichon, V.; Courtois, X. Distinct Roles of Copper in Bimetallic Copper-Rhodium Three-Way-Catalysts Deposited on Redox Supports. Distinct Roles of Copper in Bimetallic Copper-Rhodium Three-Way-Catalysts Deposited on Redox Supports. Appl. Catal. B Environ. 2005, 57, 63–72. [Google Scholar] [CrossRef]
  39. Solsona, B.; Concepción, P.; Demicol, B.; Hernández, S.; Delgado, J.J.; Calvino, J.J.; López Nieto, J.M. Selective Oxidative Dehydrogenation of Ethane over SnO2-Promoted NiO Catalysts. J. Catal. 2012, 295, 104–114. [Google Scholar] [CrossRef]
  40. Lee, J.S.; Kim, J.; Kang, M. Hydrogen Production from Ethanol Steam Reforming over SnO2-K2O/Zeolite y Catalyst. Bull. Korean Chem. Soc. 2011, 32, 1912–1920. [Google Scholar] [CrossRef]
  41. Parmaliana, A.; Arena, F.; Frusteri, F.; Giordano, N. Temperature-Programmed Reduction Study of NiO–MgO Interactions in Magnesia-Supported Ni Catalysts and NiO–MgO Physical Mixture. J. Chem. Soc. Faraday Trans. 1990, 86, 2663–2669. [Google Scholar] [CrossRef]
  42. Lin, F.; Chen, Z.; Gong, H.; Wang, X.; Chen, L.; Yu, H. Oxygen Vacancy Induced Strong Metal-Support Interactions on Ni/Ce0.8Zr0.2O2 Nanorod Catalysts for Promoting Steam Reforming of Toluene: Experimental and Computational Studies. Langmuir 2023, 39, 4495–4506. [Google Scholar] [CrossRef]
  43. Hu, Y.; Wang, N.; Zhou, Z. Synergetic Effect of Cu Active Sites and Oxygen Vacancies in Cu/CeO2-ZrO2 for the Water-Gas Shift Reaction. Catal. Sci. Technol. 2021, 11, 2518–2528. [Google Scholar] [CrossRef]
  44. Li, P.; Chen, X.; Li, Y.; Schwank, J.W. A Review on Oxygen Storage Capacity of CeO2 -Based Materials: Influence Factors, Measurement Techniques, and Applications in Reactions Related to Catalytic Automotive Emissions Control. Catal. Today 2019, 327, 90–115. [Google Scholar] [CrossRef]
  45. Kaminski, P.; Ziolek, M. Mobility of Gold, Copper and Cerium Species in Au, Cu/Ce, Zr-Oxides and Its Impact on Total Oxidation of Methanol. Appl. Catal. B Environ. 2016, 187, 328–341. [Google Scholar] [CrossRef]
  46. Fang, C.; Shi, L.; Li, H.; Huang, L.; Zhang, J.; Zhang, D. Creating Hierarchically Macro-/Mesoporous Sn/CeO2 for the Selective Catalytic Reduction of NO with NH3. RSC Adv. 2016, 6, 78727–78736. [Google Scholar] [CrossRef]
  47. Caballero, A.; Holgado, J.P.; Gonzalez-Delacruz, V.M.; Habas, S.E.; Herranz, T.; Salmeron, M. In Situ Spectroscopic Detection of SMSI Effect in a Ni/CeO2 System: Hydrogen-Induced Burial and Dig out of Metallic Nickel. Chem. Commun. 2010, 46, 1097–1099. [Google Scholar] [CrossRef]
  48. Jin, B.; Li, S.; Liang, X. Enhanced Activity and Stability of MgO-Promoted Ni/Al2O3 Catalyst for Dry Reforming of Methane: Role of MgO. Fuel 2021, 284, 119082. [Google Scholar] [CrossRef]
  49. Gutta, N.; Velisoju, V.K.; Chatla, A.; Boosa, V.; Tardio, J.; Patel, J.; Akula, V. Promotional Effect of Cu and Influence of Surface Ni-Cu Alloy for Enhanced H2 Yields from CH4 Decomposition over Cu-Modified Ni Supported on MCM-41 Catalyst. Energy Fuels 2018, 32, 4008–4015. [Google Scholar] [CrossRef]
  50. Miranda, B.C.; Chimentão, R.J.; Szanyi, J.; Braga, A.H.; Santos, J.B.O.; Gispert-Guirado, F.; Llorca, J.; Medina, F. Influence of Copper on Nickel-Based Catalysts in the Conversion of Glycerol. Appl. Catal. B Environ. 2015, 166–167, 166–180. [Google Scholar] [CrossRef]
  51. Khan, M.E.; Khan, M.M.; Cho, M.H. Ce3+-Ion, Surface Oxygen Vacancy, and Visible Light-Induced Photocatalytic Dye Degradation and Photocapacitive Performance of CeO2-Graphene Nanostructures. Sci. Rep. 2017, 7, 5928. [Google Scholar] [CrossRef]
  52. Weidler, N.; Schuch, J.; Knaus, F.; Stenner, P.; Hoch, S.; Maljusch, A.; Schäfer, R.; Kaiser, B.; Jaegermann, W. X-Ray Photoelectron Spectroscopic Investigation of Plasma-Enhanced Chemical Vapor Deposited NiOx, NiOx(OH)y, and CoNiOx(OH)y: Influence of the Chemical Composition on the Catalytic Activity for the Oxygen Evolution Reaction. J. Phys. Chem. C 2017, 121, 6455–6463. [Google Scholar] [CrossRef]
  53. Shen, Y.; Lua, A.C. Sol-Gel Synthesis of Ni and Ni Supported Catalysts for Hydrogen Production by Methane Decomposition. RSC Adv. 2014, 4, 42159–42167. [Google Scholar] [CrossRef]
  54. Casillas-Trujillo, L.; Osinger, B.; Lindblad, R.; Karlsson, D.; Abrikosov, A.I.; Fritze, S.; Von Fieandt, K.; Alling, B.; Hotz, I.; Jansson, U.; et al. Experimental and Theoretical Evidence of Charge Transfer in Multi-Component Alloys-How Chemical Interactions Reduce Atomic Size Mismatch. Mater. Chem. Front. 2021, 5, 5746–5759. [Google Scholar] [CrossRef]
  55. Azdad, Z.; Marot, L.; Moser, L.; Steiner, R.; Meyer, E. Valence Band Behaviour of Zirconium Oxide, Photoelectron and Auger Spectroscopy Study. Sci. Rep. 2018, 8, 16251. [Google Scholar] [CrossRef] [PubMed]
  56. Kim, B.J.; Park, H.R.; Lee, Y.L.; Ahn, S.Y.; Kim, K.J.; Hong, G.R.; Roh, H.S. Customized Ni-MgO-ZrO2catalysts for the Dry Reforming of Methane Using Coke Oven Gas: Optimizing the MgO Content. J. CO2 Util. 2023, 68, 102379. [Google Scholar] [CrossRef]
  57. Asencios, Y.J.O.; Yigit, N.; Wicht, T.; Stöger-Pollach, M.; Lucrédio, A.F.; Marcos, F.C.F.; Assaf, E.M.; Rupprechter, G. Partial Oxidation of Bio-Methane over Nickel Supported on MgO–ZrO2 Solid Solutions. Top. Catal. 2023, 66, 1539–1552. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, Y.; Aghamohammadi, S.; Li, D.; Li, K.; Farrauto, R. Structure Dependence of Nb2O5−X Supported Manganese Oxide for Catalytic Oxidation of Propane: Enhanced Oxidation Activity for MnOx on a Low Surface Area Nb2O5−X. Appl. Catal. B Environ. 2019, 244, 438–447. [Google Scholar] [CrossRef]
  59. Deng, J.; Zhou, Y.; Cui, Y.; Lan, L.; Wang, J.; Yuan, S.; Chen, Y. The Influence of H2O2 on the Properties of CeO2-ZrO2 Mixed Oxides. J. Mater. Sci. 2017, 52, 5242–5255. [Google Scholar] [CrossRef]
  60. Sádaba, I.; Ojeda, M.; Mariscal, R.; Fierro, J.L.G.; Granados, M.L. Catalytic and Structural Properties of Co-Precipitated Mg-Zr Mixed Oxides for Furfural Valorization via Aqueous Aldol Condensation with Acetone. Appl. Catal. B Environ. 2011, 101, 638–648. [Google Scholar] [CrossRef]
  61. Zhuang, H.D.; Bai, S.F.; Liu, X.M.; Yan, Z.F. Structure and Performance of Cu/ZrO2 Catalyst for the Synthesis of Methanol from CO2 Hydrogenation. Ranliao Huaxue Xuebao/J. Fuel Chem. Technol. 2010, 38, 462–467. [Google Scholar] [CrossRef]
  62. Vasquez, R.P. CuO by XPS. Surf. Sci. Spectra 1998, 5, 262–266. [Google Scholar] [CrossRef]
  63. Ghodselahi, T.; Vesaghi, M.A.; Shafiekhani, A.; Baghizadeh, A.; Lameii, M. XPS Study of the Cu@Cu2O Core-Shell Nanoparticles. Appl. Surf. Sci. 2008, 255, 2730–2734. [Google Scholar] [CrossRef]
  64. Hengne, A.M.; Samal, A.K.; Enakonda, L.R.; Harb, M.; Gevers, L.E.; Anjum, D.H.; Hedhili, M.N.; Saih, Y.; Huang, K.W.; Basset, J.M. Ni-Sn-Supported ZrO2 Catalysts Modified by Indium for Selective CO2 Hydrogenation to Methanol. ACS Omega 2018, 3, 3688–3701. [Google Scholar] [CrossRef] [PubMed]
  65. Degerman, D.; Lömker, P.; Goodwin, C.M.; Shipilin, M.; García-Martínez, F.; Schlueter, C.; Nilsson, A.; Amann, P. State of the Surface During CO Hydrogenation over Ni(111) and Ni(211) Probed by Operando X-Ray Photoelectron Spectroscopy. J. Phys. Chem. C 2023, 127, 4021–4032. [Google Scholar] [CrossRef]
  66. Concepción, P.; Boronat, M.; García-García, S.; Fernández, E.; Corma, A. Enhanced Stability of Cu Clusters of Low Atomicity against Oxidation. Effect on the Catalytic Redox Process. ACS Catal. 2017, 7, 3560–3568. [Google Scholar] [CrossRef]
  67. Fan, J.; Weng, D.; Wu, X.; Wu, X.; Ran, R. Modification of CeO2-ZrO2 Mixed Oxides by Coprecipitated/Impregnated Sr: Effect on the Microstructure and Oxygen Storage Capacity. J. Catal. 2008, 258, 177–186. [Google Scholar] [CrossRef]
  68. Serrano-Lotina, A.; Cruz, K.; Bañares, M.A.; Daturi, M.; Ávila, P. Bimetallic MnO2-Supported Catalysts for Selective Reduction of NO with NH3. Operando IR Studies. Appl. Surf. Sci. 2023, 610, 155550. [Google Scholar] [CrossRef]
  69. Popescu, I.; Piumetti, M.; Bensaid, S.; Marcu, I.C. Study of Ce-Cu Mixed Oxide Catalysts by: In Situ Electrical Conductivity Measurements. Phys. Chem. Chem. Phys. 2017, 19, 31929–31939. [Google Scholar] [CrossRef]
  70. Liang, P.; Zhang, Y.; Dong, X.; Zhang, J.; Li, B.; Zhang, Y.; Jiao, T.; Zhang, W.; Li, X. Effect of Sn(II)/Sn(IV) Additive on the Catalytic Combustion of Diethylamine and N2 Selectivity over CeZrOx Catalyst. Fuel 2024, 358, 130164. [Google Scholar] [CrossRef]
  71. Li, S.; Deng, J.; Wang, J.; Chen, Y.; Li, Y. Effects of Thermal Treatment Conditions on Redox Properties of Ceria-Zirconia Materials. J. Rare Earths 2023, 41, 1969–1975. [Google Scholar] [CrossRef]
  72. Coronado, I.; Stekrova, M.; Reinikainen, M.; Simell, P.; Lefferts, L.; Lehtonen, J. A Review of Catalytic Aqueous-Phase Reforming of Oxygenated Hydrocarbons Derived from Biorefinery Water Fractions. Int. J. Hydrogen Energy 2016, 41, 11003–11032. [Google Scholar] [CrossRef]
  73. Larimi, A.S.; Kazemeini, M.; Khorasheh, F. Highly Selective Doped Pt-MgO Nano-Sheets for Renewable Hydrogen Production from APR of Glycerol. Int. J. Hydrogen Energy 2016, 41, 17390–17398. [Google Scholar] [CrossRef]
  74. Guo, Y.; Azmat, M.U.; Liu, X.; Wang, Y.; Lu, G. Effect of Support’s Basic Properties on Hydrogen Production in Aqueous-Phase Reforming of Glycerol and Correlation between WGS and APR. Appl. Energy 2012, 92, 218–223. [Google Scholar] [CrossRef]
  75. Wen, G.; Xu, Y.; Ma, H.; Xu, Z.; Tian, Z. Production of Hydrogen by Aqueous-Phase Reforming of Glycerol. Int. J. Hydrogen Energy 2008, 33, 6657–6666. [Google Scholar] [CrossRef]
  76. Robinson, A.M.; Hensley, J.E.; Will Medlin, J. Bifunctional Catalysts for Upgrading of Biomass-Derived Oxygenates: A Review. ACS Catal. 2016, 6, 5026–5043. [Google Scholar] [CrossRef]
  77. Zhang, D.; Zhang, Q.; Zhou, Z.; Li, Z.; Meng, K.; Fang, T.; You, Z.; Zhang, G.; Yin, B.; Shen, J.; et al. Hydrogenolysis of Glycerol to 1,3-Propanediol: Are Spatial and Electronic Configuration of “Metal-Solid Acid” Interface Key for Active and Durable Catalysts? ChemCatChem 2022, 14, 1–21. [Google Scholar] [CrossRef]
  78. Pendem, C.; Sarkar, B.; Siddiqui, N.; Konathala, L.N.S.; Baskar, C.; Bal, R. K-Promoted Pt-Hydrotalcite Catalyst for Production of H2 by Aqueous Phase Reforming of Glycerol. ACS Sustain. Chem. Eng. 2018, 6, 2122–2131. [Google Scholar] [CrossRef]
  79. Manfro, R.L.; Pires, T.P.; Ribeiro, N.F.P.; Souza, M.M. Aqueous-Phase Reforming of Glycerol Using Ni-Cu Catalysts Prepared from Hydrotalcite-like Precursors. Catal. Sci. Technol. 2013, 3, 1278–1287. [Google Scholar] [CrossRef]
  80. Shabaker, J.W.; Simonetti, D.A.; Cortright, R.D.; Dumesic, J.A. Sn-Modified Ni Catalysts for Aqueous-Phase Reforming: Characterization and Deactivation Studies. J. Catal. 2005, 231, 67–76. [Google Scholar] [CrossRef]
  81. Zhang, C.C.; Hartlaub, S.; Petrovic, I.; Yilmaz, B. Raman Spectroscopy Characterization of Amorphous Coke Generated in Industrial Processes. ACS Omega 2022, 7, 2565–2570. [Google Scholar] [CrossRef]
  82. Rosli, S.N.A.; Abidin, S.Z.; Osazuwa, O.U.; Fan, X.; Jiao, Y. The Effect of Oxygen Mobility/Vacancy on Carbon Gasification in Nano Catalytic Dry Reforming of Methane: A Review. J. CO2 Util. 2022, 63, 102109. [Google Scholar] [CrossRef]
  83. Perosa, A.; Tundo, P. Selective Hydrogenolysis of Glycerol with Raney Nickel. Ind. Eng. Chem. Res. 2005, 44, 8535–8537. [Google Scholar] [CrossRef]
  84. Yu, W.; Xu, J.; Ma, H.; Chen, C.; Zhao, J.; Miao, H.; Song, Q. A Remarkable Enhancement of Catalytic Activity for KBH4 Treating the Carbothermal Reduced Ni/AC Catalyst in Glycerol Hydrogenolysis. Catal. Commun. 2010, 11, 493–497. [Google Scholar] [CrossRef]
  85. Marinoiu, A.; Ionita, G.; Gáspár, C.L.; Cobzaru, C.; Oprea, S. Glycerol Hydrogenolysis to Propylene Glycol. React. Kinet. Catal. Lett. 2009, 97, 315–320. [Google Scholar] [CrossRef]
  86. Syuhada, A.; Ameen, M.; Azizan, M.T.; Aqsha, A.; Yusoff, M.H.M.; Ramli, A.; Alnarabiji, M.S.; Sher, F. In-Situ Hydrogenolysis of Glycerol Using Hydrogen Produced via Aqueous Phase Reforming of Glycerol over Sonochemically Synthesized Nickel-Based Nano-Catalyst. Mol. Catal. 2021, 514, 111860. [Google Scholar] [CrossRef]
  87. Freitas, I.C.; Manfro, R.L.; Souza, M.M.V.M. Hydrogenolysis of Glycerol to Propylene Glycol in Continuous System without Hydrogen Addition over Cu-Ni Catalysts. Appl. Catal. B Environ. 2018, 220, 31–41. [Google Scholar] [CrossRef]
  88. Raso, R.; García, L.; Ruiz, J.; Oliva, M.; Arauzo, J. Aqueous Phase Hydrogenolysis of Glycerol over Ni/Al-Fe Catalysts without External Hydrogen Addition. Appl. Catal. B Environ. 2021, 283, 119598. [Google Scholar] [CrossRef]
  89. Cai, F.; Pan, D.; Ibrahim, J.J.; Zhang, J.; Xiao, G. Hydrogenolysis of Glycerol over Supported Bimetallic Ni/Cu Catalysts with and without External Hydrogen Addition in a Fixed-Bed Flow Reactor. Appl. Catal. A Gen. 2018, 564, 172–182. [Google Scholar] [CrossRef]
  90. Gallego-García, D.; Iriarte-Velasco, U.; Gutiérrez-Ortiz, M.A.; Ayastuy, J.L. Nickel Aluminate Spinel-Derived Catalysts for Aqueous-Phase Hydrogenolysis of Glycerol with in-Situ Hydrogen Production: Effect of Molybdenum Doping. Appl. Catal. B Environ. 2024, 344, 123671. [Google Scholar] [CrossRef]
Figure 1. XRD profiles of the (A) calcined, (B) reduced solids.
Figure 1. XRD profiles of the (A) calcined, (B) reduced solids.
Ijms 25 03484 g001
Figure 2. Raman spectra for (A) calcined and (B) reduced forms of the Me/NiCeZr materials.
Figure 2. Raman spectra for (A) calcined and (B) reduced forms of the Me/NiCeZr materials.
Ijms 25 03484 g002
Figure 3. Reduction profiles of Me/NiCeZr solids.
Figure 3. Reduction profiles of Me/NiCeZr solids.
Ijms 25 03484 g003
Figure 4. Detailed Ni 2p spectra for calcined (A) and reduced (B) catalysts.
Figure 4. Detailed Ni 2p spectra for calcined (A) and reduced (B) catalysts.
Ijms 25 03484 g004
Figure 5. XPS spectra of Mg 2s (A), Cu 2p3/2 (B), and Sn 3d5/2 (C) in the calcined and reduced forms of Me/NiCeZr catalysts.
Figure 5. XPS spectra of Mg 2s (A), Cu 2p3/2 (B), and Sn 3d5/2 (C) in the calcined and reduced forms of Me/NiCeZr catalysts.
Ijms 25 03484 g005
Figure 6. Effect of the composition of the catalyst on overall glycerol conversion and yields of liquid and gas. Reaction conditions: T = 235 °C, P = 35 bar, WHSV = 12 h−1, 10 wt.% glycerol, data at 3 h TOS.
Figure 6. Effect of the composition of the catalyst on overall glycerol conversion and yields of liquid and gas. Reaction conditions: T = 235 °C, P = 35 bar, WHSV = 12 h−1, 10 wt.% glycerol, data at 3 h TOS.
Ijms 25 03484 g006
Scheme 1. Reaction routes in the hydrogenolysis of glycerol.
Scheme 1. Reaction routes in the hydrogenolysis of glycerol.
Ijms 25 03484 sch001
Figure 7. Selectivity towards primary liquid products. * Others include methanol, acetone, and 1-propanol.
Figure 7. Selectivity towards primary liquid products. * Others include methanol, acetone, and 1-propanol.
Ijms 25 03484 g007
Figure 8. (A) XRD patterns and (B) Raman spectra of used catalysts.
Figure 8. (A) XRD patterns and (B) Raman spectra of used catalysts.
Ijms 25 03484 g008
Table 1. Chemical composition of the bulk, crystallite sizes, lattice parameter, and available metal sites.
Table 1. Chemical composition of the bulk, crystallite sizes, lattice parameter, and available metal sites.
CatalystMe
(wt.%)
Ni
(wt.%)
Ce/Zr
(at./at.)
SBET
(m2/g)
dNiCeZr.XRD
(nm)
DA aLattice *
(Å)
dNi0
(nm)
Ni0 Sites (×10−19)
(Sites/gcat) b
NiCeZr09.214.1/85.9235 (190)2.7 (2.9)1.6a = 3.6094 (3.6559)
c = 5.2539 (5.2705)
11.44.2
Mg/NiCeZr1.059.213.6/86.4226 (197)2.8 (3.3)1.6a = 3.6241 (3.6615)
c = 5.2315 (5.2181)
10.92.9
Cu/NiCeZr1.069.215.7/84.3224 (202)2.3 (2.9)1.9a = 3.6452 (3.6698)
c = 5.2014 (5.2232)
5.00.7
Sn/NiCeZr0.829.214.6/85.4200 (199)2.5 (3.0)2.0a = 3.6366 (3.6720)
b = 5.2244 (5.2239)
20.10.5
* Calculated as tetragonal: 1 d 2 = ( h 2 + k 2 ) a 2 + l 2 c 2 (a = b ≠ c). Values in parenthesis correspond to the reduced forms. a DA: degree of agglomeration, calculated as dBET/dXRD. b calculated from H2 chemisorption.
Table 2. Redox and surface characteristics of the materials.
Table 2. Redox and surface characteristics of the materials.
CatalystH2 Uptake (mmolH2/g)OSC
(mmolo/gcat)
Sites Density (Sites/nm2)Metal/Acid
(Sites/Sites)
Overall aPeakMetal bAcid c
αβγ
NiCeZr2.20 (106)0.110.441.65630.220.360.62
Mg/NiCeZr1.62 (78)0.030.251.34600.150.300.49
Cu/NiCeZr2.28 (101)0.100.621.56650.030.420.08
Sn/NiCeZr2.28 (103)0.040.461.78430.020.510.05
In parenthesis, the percentage of reduction with respect to the theoretical ratio for the complete reduction of all nickel, cerium, and Me ions, assuming they are present as Ni2+, Ce4+, Cu2+, and Sn4+, respectively. a Uptake up to 900 °C. b from H2 chemisorption. c from ammonia chemisorption.
Table 3. Surface chemical composition from the XPS analysis.
Table 3. Surface chemical composition from the XPS analysis.
CatalystCe/ZrNi/(Ce + Zr)Ce3+/(Ce4+ + Ce3+)Bulk Me/
(Ni + Ce + Zr)
Me/(Ni + Ce + Zr)Ni0/NiTotal
NiCeZr0.13 (0.17)0.17 (0.08)0.16 (1.0)00 (0)0.0 (0.70)
Mg/NiCeZr0.18 (0.21)0.33 (0.10)0.27 (1.0)0.050.11 (0.10)0.0 (0.79)
Cu/NiCeZr0.17 (0.22)0.39 (0.10) 0.19 (1.0)0.020.06 (0.05)0.0 (0.89)
Sn/NiCeZr0.16 (0.23)0.31 (0.10)0.19 (1.0)0.010.03 (0.02)0.0 (0.86)
In parenthesis, values for solids in situ reduced at 600 °C for 1 h in 10% H2/He flow.
Table 4. Gas product yields and selectivities.
Table 4. Gas product yields and selectivities.
CatalystYield (%)Selectivity (%)H2/CO2
H2CH4COCO2C2+H2CH4COCO2C2+
NiCeZr2.34.300.6511.152.1532.223.53.661.111.80.48
Mg/NiCeZr3.29.210.7617.553.3024.829.92.556.910.70.43
Cu/NiCeZr1.82.460.998.073.0132.516.96.855.520.80.51
Sn/NiCeZr0.60.690.431.601.6129.516.010.136.937.00.87
Table 5. Selectivities to primary and secondary products, type of cleaved bond, and yield to 1,2-PG (in %).
Table 5. Selectivities to primary and secondary products, type of cleaved bond, and yield to 1,2-PG (in %).
CatalystPrimary ProductsSecondary ProductsC–CC–OO/CY1,2-PG
NiCeZr75.524.56.693.40.6528.5
Mg/NiCeZr79.620.47.992.10.6619.3
Cu/NiCeZr87.912.15.194.90.6625.0
Sn/NiCeZr92.67.43.996.10.6716.7
Primary products: single C–O or C–C bond scission (EG, 1,2-PG, HA). Secondary products: more than a single C–O or C–C bond scission (ethanol, methanol, acetone, 1-propanol).
Table 6. Textural properties, crystallite size, and concentration of the leached metals of Me/NiCeZr catalysts after 3 h of use.
Table 6. Textural properties, crystallite size, and concentration of the leached metals of Me/NiCeZr catalysts after 3 h of use.
CatalystSBET
(m2/g)
ΔVpore
(%)
Δdpore
(%)
dNiº (nm)Leaching (wt.%)
NiMeCeZr
NiCeZr218 (+15%)+21+721.03.90.00.90.0
Mg/NiCeZr211 (−7%)−21−158.92.140.30.10.0
Cu/NiCeZr279 (+25%)+28+167.80.90.00.10.0
Sn/NiCeZr246 (+23%)+35−718.20.60.00.00.0
In parenthesis, percentage variation with respect to freshly reduced forms.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vera-Hincapie, N.; Iriarte-Velasco, U.; Ayastuy, J.L.; Gutiérrez-Ortiz, M.Á. Hydrogenolysis of Glycerol over NiCeZr Catalyst Modified with Mg, Cu, and Sn at the Surface Level. Int. J. Mol. Sci. 2024, 25, 3484. https://doi.org/10.3390/ijms25063484

AMA Style

Vera-Hincapie N, Iriarte-Velasco U, Ayastuy JL, Gutiérrez-Ortiz MÁ. Hydrogenolysis of Glycerol over NiCeZr Catalyst Modified with Mg, Cu, and Sn at the Surface Level. International Journal of Molecular Sciences. 2024; 25(6):3484. https://doi.org/10.3390/ijms25063484

Chicago/Turabian Style

Vera-Hincapie, Norberto, Unai Iriarte-Velasco, Jose Luis Ayastuy, and Miguel Ángel Gutiérrez-Ortiz. 2024. "Hydrogenolysis of Glycerol over NiCeZr Catalyst Modified with Mg, Cu, and Sn at the Surface Level" International Journal of Molecular Sciences 25, no. 6: 3484. https://doi.org/10.3390/ijms25063484

APA Style

Vera-Hincapie, N., Iriarte-Velasco, U., Ayastuy, J. L., & Gutiérrez-Ortiz, M. Á. (2024). Hydrogenolysis of Glycerol over NiCeZr Catalyst Modified with Mg, Cu, and Sn at the Surface Level. International Journal of Molecular Sciences, 25(6), 3484. https://doi.org/10.3390/ijms25063484

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop