Fishbone derived-hydroxyapatite supported Ni-Zr nanocatalyst for CO2 methanation: Synergistic effects of support and zirconia
⁎Corresponding authors. ntri@ict.vast.vn (Nguyen Tri), hkpha@hcmut.edu.vn (Huynh Ky Phuong Ha)
-
Received: ,
Accepted: ,
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.
Abstract

Abstract
Herein, the effectiveness of nickel-zirconia nanocatalyst decorated on hydroxyapatite (HA) is examined as a potential candidate for the methanation of carbon dioxide (CO2). Co-impregnation was used to produce undoped nickel as well as zirconia-doped nickel catalysts that were supported on hydroxyapatite derived from fishbone. The prepared catalysts were characterized using several methods. With the molar ratio of H2/CO2 fixed at 4/1, the catalytic activity was evaluated through the CO2 methanation at low temperatures ranging of 250–––400 °C. A high dispersion of nickel particles over the catalyst surface with a size range of 30–––40 nm was provided as a function of HA and zirconia. This led to an improvement in the reducibility and CO2 adsorption capacity of the catalysts and an enhancement of the efficiency in CO2 methanation. At a reaction temperature of 400 °C, the data reveal that the 10Ni/HA sample doped with 4 wt% zirconia was the most effective. At this temperature, the CO2 conversion reached 92.9%, and the CH4 selectivity reached approximately 100%.
Keywords
CO2 methanation
Hydroxyapatite
Fishbone
Zirconia
Nickel nanocatalyst
1 Introduction
In the period 2015–––2040, the global energy demand is expected to increase by about 30% compared to the present time (Zhao et al., 2022), and to date, the main supply of energy still comes from fossil fuels. These natural resources need millions of years to form; however, the increasing demand of humanity is leading to the gradual exhaustion of this fuel. Besides, obtaining energy through burning fossil fuels also releases greenhouse gases, leading to global warming (Wai et al., 2021). In addition, the fossil fuel extraction also poses many threats to the environment, affecting vegetation and topsoil, increasing soil erosion, and losing habitat for many species. Under climate change and environmental pollution pressure, energy transition and carbon emission reduction is an inevitable orientation that deserves utmost attention. To meet the goal of declining the CO2 concentration in the second half of this century as the commitment to the “Net-Zero global CO2 emission” (Jiang et al., 2022); CO2 conversion technology is the center of interest. Besides the recent development of photocatalytic CO2 reduction which may adapt to the upcoming issues (Jiang et al., 2022; Muiruri et al., 2022), the conventional thermal CO2 methanation (Sabatier reaction) also seems appealing to investigators. Especially, when this reaction is a part of the “Power to Gas” technology, where hydrogen is achieved from water electrolysis splitting, and the optional excess hydrogen can be transformed into CH4 for storage via the Sabatier reaction. Then, the generated methane is stored or used for many different purposes (Mebrahtu et al., 2019). The advantages of the innovation in hydrogen production (electrolysis and other technologies (Ye et al., 2023), as well as renewable energy, are creating pathway for the CO2 methanation and the PtG process to be able to be applied in the future energy grid.
The process of CO2 methanation is an exothermic reaction, mainly occurring in the temperature range of 200–––400 °C, and in the reaction process, depending on different influencing factors, by-products such as carbon, CO, etc, can be produced. Besides tuning the reaction conditions, the catalysts also play an important role to obtain the efficient transforming CO2 into CH4. Many catalysts based on precious metals (Rd, Rh, Ru, Ir, etc) and transition metals (Ni, Co, Fe, Mn) with various supports (SiO2, Al2O3, TiO2, and SiO2-Al2O3) and medium pore materials have been studied to improve the efficiency of converting CO2 into methane gas (Lee et al., 2021; Komatsu and Furukawa, 2015; Namvar et al., 2023; Bakar et al., 2015). Due to the low cost and accessibility, Ni-metal-based catalysts are considered to be the most common, exhibiting high activity, good selectivity. The morphology, crystal lattice phase, size, dispersion and supporting interactions of the Ni catalyst, which can be tuned by altering the support, have a significant influence on the CO2 methanation (Kuznecova and Gusca, 2017; Van et al., 2020). Applying and studying the support are common orientations to improve the catalysts’ performance.
On the other hand, converting by-products or waste in agriculture or industry into applicable materials is a helpful strategy in the concept of sustainable development. One of the materials being studied and applied as support nowadays is hydroxyapatite (HA). HA has the molecular formula Ca10(PO4)6(OH)2, exists in crystalline form, hexagonal structure in different shapes depending on the synthesis method or other conditions in the synthesis process. The HA crystals are interconnected, creating gaps and porosity for the material. This material can be synthesized from organic materials such as animal bones, corals, seashells, eggshells, fish scales, shrimp shells, etc., or chemical methods using starting materials such as compounds containing calcium and phosphorus. Though HA originated from natural or chemical sources have the same nature and chemical composition, synthesizing from industrial by-products or waste has a more positive response to the environment. The FAO (Food and Agricultural Organization) estimates that only about 50–––60% of fish caught is for human consumption, the rest consists of by-products such as bones, skin, and offal. It is considered a very low value by-product. Today, this by-product is used in animal feed production; however, a large amount of the remaining waste is wasted and discharged directly into the environment, causing a great impact on the environment. Therefore, taking advantage of wasted by-products to synthesize new materials brings economic benefits and contributes to solving environmental problems.
With the thermal stable property, the application of HA as support for Ni-based catalyst has been conducted in our previous study in the CO2 methanation reaction. In order to increase the activity of the Ni/HA catalyst, it’s recommended to add a second metal (M) such as Fe, Zr, Co, La, Mg, or Ru to form an alloy system Ni–M or metal oxide to improve reducing properties, enhance mobility, improve metal dispersion and prevent metal sintering (Su et al., 2016). Among those promoters, ZrO2 is widely used due to its outstanding properties of enhanced activity, high thermal stability, and low carbon deposition rate due to its large number of active electron sites on the surface, base sites (OH groups, O2– adsorbed), acid-base pairs Zr4+–O2– and oxygen vacancies. ZrO2′s interactions with metals and other oxides can affect CO2 adsorption and activation, enhance H2 dissociation and atomic hydrogen diffusion, and change the reaction and binding pathways of the intermediate to react to further conversion (Li and Chen, 2019). This rare earth oxide has been studied on some Ni-based catalysts system, however, the addition and effect of this promoter on the Ni/HA catalyst haven’t been elucidated yet. Hence, in this work, the interaction of ZrO2 with Ni active sites and the HA support is explored in the thermal catalytic CO2 methanation.
2 Experimental
2.1 Fabrication of catalyst
The raw salmon bone was collected from seafood manufacturers in Vietnam. All the chemicals used in the experiment were purchased from Merck and used without further purification. The preparation of HA was taken place as described in our previous research (Tri et al., 2020). 1.0 g of the dried salmon bone was introduced to 50 mL of distilled water, then 14 mL H3PO4 1 wt% was added dropwise into the above mixture. The mixture was adjusted to pH = 10 by NH4OH 5 wt% solution and stirred vigorously for 2 h. After stirring, the mixture was transferred into Teflon-line auto-clave and heated to 120 °C for 7 h. The precipitate after the hydrolysis was isolated by centrifuging, washed with distilled water and ethanol, then dried at 80 °C for 12 h. The HA support was obtained after calcining the dry precipitate in the air at 800 °C for 1 h.
Ni-ZrO2/HA was prepared by the impregnation method. Calculated amounts of Ni(NO3)2·6H2O, ZrOCl2·8H2O, and HA support were mixed with a decent amount of distilled water. After mixing, the mixture was dried at 3 levels of temperature (80, 100, and 120 °C) for 2 h for each temperature. The calcination was carried out at 600 °C for 4 h (a heating rate of 10 °C/min) to obtain the NiO-ZrO2/HA. Before the catalytic activity evaluation, a pre-treatment of NiO-ZrO2/HA in a pure H2 flow (3 L/h) was taken place and the Ni-ZrO2/HA catalyst was officially obtained. The Ni content was fixed at 10 wt%. The catalysts were assigned 10NixZr/HA, in which × value is the wt.% of ZrO2 and × = 2, 4, and 6.
2.2 Characterization
The physical–chemical characteristics of the obtained samples were studied by several methods (Van et al., 2020), including X-ray diffraction on a Bruker D2 Phaser Xray diffractometer with Cu Kα radiation and recorded in 2θ = 10 − 80° with a step size of 0.03°; EDS spectrum on a JEOL JST-IT 200 instrument; ultraviolet–visible diffuse reflectance spectroscopy (UV–Vis DRS) on a V-770 UV–visible/NIR spectrophotometer; nitrogen adsorption − desorption isotherms at − 197 °C on a Nova 2200e instrument with the sample pretreated in a 1.8 L/h N2 flow at 300 °C for 2 h; field-emission scanning electron microscopy on Hitachi S4800 instrument; high-resolution transmission electron microscopy (HR-TEM) on JEOL JEM 1400 instrument; hydrogen temperature-programmed reduction (H2-TPR) carried out on a microreactor in a gas mixture of 10% H2/N2 at a flow rate of 1.8 L/h using a gas chromatograph GOW-MAC 69–––350 with a TCD detector and the used sample of 50 mg; and CO2-programmed temperature desorption (CO2-TPD) also carried out on a microreactor in He at a flow rate of 1.8 L/h from 50 to 800 °C at a ramping rate of 10 °C/min with 100 mg of catalyst adopted in a quartz reactor and reduced at 450 °C for 1 h in pure H2 at a flow rate of 1.8 L/h prior to being cooled down to 50 °C.
2.3 Catalytic activity
For catalytic activity testing, a quartz tube placed in a furnace was used as the reactor system. The catalyst bed is a mixture of 0.5 g quartz and 0.5 g catalyst in the same size of 0.25–––0.50 mm. The reduction process was carried out under a 3 L/h H2 flow at 400 °C for 2 h with a heating rate of 10 °C/min. After the pre-treatment, the catalytic activity of all samples was tested at the temperature range of 250–––400 °C with a total flow rate of 6 L/h consisting of a CO2/H2/N2 ratio of 3/12/85 mol%. To determine the CO2 conversion and CH4 selectivity in product flow, the reaction mixture was analyzed on a gas chromatograph (Agilent 5890 Series II) using a thermal conductivity detector with a capillary column (TG-BON Q, 30 m length, 12 μm inter diameter, 0.32 μm thickness) and a flame ionization detector with a capillary column (DB624, 30 m length, 0.32 mm outer diameter, 0.25 μm thickness). The gas mixtures of CH4 6 mol.% in N2 and CO 6 mol.% in N2 were used as a reference to monitor CH4 and CO produced in the output. CO2 conversion (
3 Results and discussions
3.1 Physico-chemical characteristics
The XRD patterns of catalysts with different Zr contents (Fig. 1a) shows that the presence of HA hexagonal phase in all samples is indicated by the diffraction peaks at 2θ = 10.9, 16.9, 21.6, 22.8, 26.1, 27.8, 28.7, 31.8, 32.0, 32.9, 34.0, 39.9, 46.8, 49.5, 50.5, 51.3, 52.0, 53.0, 60.2, 61.7, 64.0, 65.1, and 77.0° (JCPDS card No.09–432). Besides, the β-TCP Ca3(PO4)2 phase with low concentration was also detected with low intensive diffraction peaks at 2θ = 13.7, 31.2, 34.7, 35.6, 37.6, 57.1, 58.0, and 78.3° (JCPDS card No.09–169). The presence of β-TCP Ca3(PO4)2 relates to pH and sintering temperature (Le et al., 2012; Boukha et al., 2019). The same phenomenon was also observed in Taimai’s study (Tamai et al., 2007), in which a ratio of Ca/P < 1.67 induced the decomposition of HA, resulting in the formation of β-TCP Ca3(PO4)2. The characteristic diffraction peaks of Ni appear at 2θ = 44.7, 51.9, and 76.8° (JCPDS No.04–0850). According to Standard Card mp-20194, the characteristic peaks for ZrO2 appear at 2θ = 24.3, 29.7, 33.7, 34.8, 42.1, 49.3, 50,1, 57.9, 59.3, 61.7, 70.9, and 73.5°. However, these ZrO2 peaks are hardly observed. This phenomenon could account for the relatively low content of ZrO2, or maybe some of the Ni’s peaks have coincided with the HA peaks (2θ = 33.7, 49.3, 50.1, and 61.7°). The texture properties of HA support and undoped nickel as well as zirconia-doped nickel catalysts are given in the Table 1. The addition of zirconia to the 10Ni/HA catalyst shows the effect of decreasing the Ni crystalline size (15.3 nm down to 5.5 nm), which agrees with several pieces of literature investigating the impact of ZrO2 promoter in the Ni-based catalysts (Roh et al., 2004; Cai et al., 2011; Chong et al., 2019). The improved dispersal of active sites could result in the amelioration of the catalyst’s overall activity. The Ni crystalline size calculated by the Debye-Scherrer equation is 5.5 nm for the reduced 10Ni4Zr/HA sample, three times smaller than that of the pure 10Ni/HA sample.

- XRD patterns (a), H2-TPR (b) and CO2-TPD (c) profiles, the hydrogen consumption (H) from H2-TPR results and the desorbed CO2 amount (
Samples | dHAa, nm | dNia, nm | SBETb, m2.g−1 | dporeb, nm | Vporeb, cm3.g−1 | Tmaxc, °C |
---|---|---|---|---|---|---|
HA | 8.7 | – | 13.8 | 2.48 | 0.018 | 680 |
10Ni/HA | 8.7 | 15.3 | 15.6 | 2.36 | 0.019 | 320, 340, 370, 400 |
10Ni2Zr/HA | 8.7 | 9.2 | 17.3 | 2.45 | 0.023 | 370, 479, 642 |
10Ni4Zr/HA | 8.7 | 5.5 | 19.7 | 2.56 | 0.032 | 340, 370, 509, 671 |
10Ni6Zr/HA | 8.7 | 5.5 | 18.9 | 2.50 | 0.028 | 340, 370, 493, 671 |
The H2-TPR results (Fig. 1b and Table 1) were obtained to assess the reducibility of the supported Ni species (Boukha et al., 2022). Fig. 1b illustrates that the H2-TPR results of zirconia-doped catalysts exhibit typical peaks corresponding to the same temperature areas as the profiles of unmodified catalysts. The 10Ni/HA sample’s TPR profile exhibited four reduction peaks below 450 °C, which can respectively be classified as 4 types of NiO species (α, β1, β2, and γ). The first characteristic is ascribed to the reduction of surface oxygen species (α species), which induces the oxygen vacancies on the catalyst’s surface. The second and third reduction peaks can be attributed to the reduction of NiO that exhibits a weak and medium interaction with HA (β species), which further classified into β1 and β2 (Boukha et al., 2019). The last reduction peak at higher temperature, which assigned to the γ- NiO species, might be induced from the bulk NiO particles that strongly attached to the HA support surface. Intriguingly, the addition of Zr alters the peaks toward lower temperatures, demonstrating the promoter effect of zirconium on the reduction of Ni species. More specifically, the reduction temperatures of these catalysts are lower than those of the unmodified catalysts, suggesting improved metal phase dispersion since ZrO2 minimizes the weakening of the NiO–support contact and enhances the capacity to scatter NiO on the HA surface (Boukha et al., 2022). In addition, the Fig. 1b demonstrates that the zirconia-doped catalysts exhibit faint reduction signals at 600 °C; these reduction signals pertain to the reduction process of the NiO–ZrO2 interaction phases or surface ZrO2 in the doped catalysts. The ZrO2 doping level is modest (2–6 wt%); hence the decrease signal is not strong. Furthermore, regardless of the Zr-modified sample, the amounts of consumed H2 are considerably greater than those required for the reduction of a stoichiometric NiO phase. The observed tendency can be attributed to the presence of highly reducible Ni sites, most likely those located close to the NiO-ZrO2 interface. This is consistent with the observed increase in the density of defective metallic Ni crystallites with increasing Zr loading, as demonstrated by XRD data for reduced samples. Notably, our results show that Zr promotes the reducibility of 10NixZr/HA catalysts.
To evaluate the CO2 adsorption capacity of the ZrO2 modified catalyst for the methanation process, the CO2-TPD desorption method was performed and the results are shown in Fig. 1c. All catalysts show a wide and high intensity CO2 desorption region in the temperature range of 80–––300 °C, indicating that the catalyst surface contains a large number of medium and weak base centers. Low CO2 desorption peaks at temperatures from 500 °C onwards, especially at above 700 °C are possibly accounted for strong base sites. It was found that CO2 adsorbed on medium-intensity basic sites could be readily desorbed from the catalyst surface (Pan et al., 2014) to promote CO2 methanation. From Table 1, the CO2 desorption capacity of the 10NixZr/HA catalysts is higher than that of the 10Ni/HA catalyst. This could be associated to the more uniform Ni nanoparticles dispersing on the ZrO2-HA surface. According to researches (Diallo-Garcia et al., 2014; Mihet and Lazar, 2018), CO2 is primarily adsorbed on surface hydroxyl groups and metal nanoparticles affect CO2 adsorption capacity. The addition of ZrO2 also provided new active sites that could absorb CO2, since ZrO2 is a quantitative oxide and the hydroxyl provided by it can also increase the number of reactive sites and enhance CO2 adsorption. ZrO2 can inhibit carbon deposition and Ni sintering due to its strong oxygen storage capacity (Yang et al., 2020). The enhanced Ni dispersion and the slight decrease of primary sites could suggest a more harmonious more harmony combined adsorption of H2 and CO2, proposing a better catalytic activity in the hydrogenation of CO2 to methane.
The hydrogen consumption (in H2-TPR analysis) and CO2 desorption capacity (in CO2-TPD analysis) are illustrated in Fig. 1d. As being mentioned above, the addition of ZrO2 promotes the reduction of NiO species, corresponding to the increased amount of consumed hydrogen in H2-TPR analysis. Although the strong metal-support interaction can enhance the catalytic stability through preventing the Ni sites from agglomerating, weakening those interactions by adding ZrO2 may also improve the catalytic activity with more presenting Ni metal sites. The maximum hydrogen consumption is observed in the 10Ni4Zr/HA’s profile, indicating more ZrO2 loading might cover Ni sites and slightly decreases the Ni site density on the surface. The dispersion of Ni onto the HA surface didn’t alter the basic property of the catalyst, implied by the approximately same CO2 desorption capacity as seen in Fig. 1d. Adding ZrO2 produces more CO2 adsorption sites, and the CO2 adsorption capacity reaches its peak in the 4 wt% ZrO2 modification. The excess ZrO2 (higher than 4 wt%) does not induce more basic sites.
Fig. 1e indicates that the isotherm N2 adsorption–desorption diagrams of all the examined samples are of type III with hysteresis loops corresponding to the H3 class of the IUPAC classification (Li et al., 2023). The presence of a medium pore structure, as indicated by isotherms of type III, is advantageous for CO2 adsorption as well as the transport of reactants and intermediates during the reaction (Wang et al., 2019). The absence of capillary condensation in the isotherms of the catalyst sample at low pressures indicates that HA does not have a porous capillary structure (El Shafei et al., 2004). Fig. 1f reveals that the pore size distribution for all samples ranges between 10–––60 Å, which is typical for materials with medium pore sizes. The 10Ni4Zr/HA catalyst has a maximum mean pore diameter of 26 Å, with a slightly higher pore volume compared to other analyzed samples. This textural property is suitable for use for application in catalytic CO2 methanation (the kinetic diameters of CO2 and H2 correspond to 0.33 nm and 0.29 nm) (Li et al., 2006). The specific surface area of the modified catalyst is moderately more considerable than that of the unmodified one, which could be induced by the smaller nanoparticles, having a larger area/volume ratio compared to their bulk material.
The FESEM images (Fig. 2) show that all catalysts have rough and porous surfaces, and the particles have a spherical shape. Two particle sizes were observed on the catalyst surface: large HA particles and smaller particles of the active metal phase. The unmodified catalyst has a surface of highly concentrated and dense nanoparticles, stacked on top of each other, explaining why this catalyst has the lowest specific surface area (Fig. 2a). Meanwhile, on the surface of ZrO2 modified Ni/HA catalysts (Fig. 2b-d), it is observed that small spherical particles, with a size of approximate 10 nm, are sparsely distributed and relatively more evenly spread. The FESEM result is consistent with the increase in specific surface area of the modified catalyst compared with the unmodified one in BET result.

- SEM images of samples; a) 10Ni/HA, b) 10Ni2Zr/HA, c) 10Ni4Zr/HA, and d) 10Ni6Zr/HA.
The SAED and HR-TEM results (Fig. 3a-c) show the microstructural properties of the 10Ni4Zr/HA sample. The SAED image (Fig. 3a) shows that the nanoparticles are uniformly crystalline with several observed diffraction rings of HA (d = 0.34 and 0.28 nm assigned to (0 0 2) and (2 1 1) planes), and Ni (d = 0.20 and 0.17 nm associated to (1 1 1) and (2 0 0) plane). In addition, a faded ring with a respective d-spacing of 0.24 nm could be attributed to the (1 1 1) plane of NiO, indicating a part of the catalyst is not reduced into Ni0. The HRTEM images (Fig. 3b and c) show that the nanoparticles exist with a size range of 30–––40 nm, which is consistent with the captured SEM images and the crystal size calculated from the XRD result. The particles tend to clump together. The HRTEM image also shows bright patches corresponding to ZrO2 and Ni centers corresponding to dark patches. ZrO2 particles tend to agglomerate to form large masses with capillary pores and Ni evenly distributed on the surface and in the developed pores.

- SAED (a), HRTEM (b and c) images, EDS spectrum and component (d) and mapping (e) of elements, and the UV–Vis DRS spectrum (f) of 10Ni4Zr/HA sample.
Based on the EDS results (Fig. 3d), it was found that the mass composition of the elements on the catalyst surface is relatively close to the theoretical calculation, which are O (35.81%), P (15.93%), Ca (34.26%), Ni (10%), and Zr (4%). Slight deviations between experimental and theoretical values exist, possibly due to the impregnation method, since adequate contact cannot be guaranteed in all cases between the metal salt solution and the support during sample drying (Daza et al., 2009). The Ca/P molar ratio < 1.67 agrees with the appeared diffraction peaks of β-TCP Ca3(PO4)2 analyzed in the XRD spectrum. Moreover, the colourful dots on the EDS mapping results (Fig. 3e) demonstrated the uniform distribution of components on the material surface and the high dispersion of Ni on the support. The UV–Vis analysis results of the typical catalyst, 10Ni4Zr/HA (Fig. 3f), reveal three distinct absorbance regions. The first region lies within the wavelength range of 320 nm, characteristic of ZrO2 absorption (Soni et al., 2022). The second region falls within the wavelength range of 380 nm, indicative of Ni or NiO absorption (Srihasam et al., 2020). Finally, the third region is observed at a wavelength of 585 nm, representing HA absorption (Pal et al., 2020). Based on these findings, it can be deduced that the material 10Ni4Zr/HA has been successfully synthesized.
3.2 Activity of catalysts
The catalytic CO2 conversion of catalysts is given in Fig. 4a as temperature functions. Under the same reaction condition, in the temperature range of 250–––400 °C, the CO2 conversion of the ZrO2-modified catalyst samples was significantly higher than that of the unmodified catalyst sample. At low temperature, the CO2 conversion increased slowly, but from about 300 °C, the CO2 conversion increased rapidly and reached a maximum value at 400 °C for all modified samples. ZrO2 is a weak hydrophilic (Li and Chen, 2019), beneficial for water desorption and thermodynamically promoting CH4 formation. ZrO2 plays many roles in different catalysts, a reasonable load of ZrO2 can enhance the active metal's dispersion, reducing the crystal size of the catalyst, enlarging the surface area and increasing the Ni0 active site (from H2-TPR result), improving H2 adsorption and CO2 conversion.

- CO2 conversion (a) and CH4 selectivity (b) of 10NixZr/HA catalysts in CO2 methanation.
The effect of ZrO2 can be explained by the mechanism by which ZrO2 adsorbs CO2, while Ni adsorbs and dissociates H2. The adsorbed CO2 will interact with bridging OH groups on the surface of ZrO2 to form bicarbonate, an essential intermediate for methane synthesis. The role of ZrO2 is also involved in creating vacancy oxygen sites. This enhancement is due to many oxygen vacancies on the surface of m-ZrO2, which enhances the adsorption capacity. t-ZrO2 with abundant emptiness also affects the catalytic performance (Li and Chen, 2019).
As shown in SEM images, the compactness of the catalyst particles in NiZr/HA samples was less than that of Ni/HA. It is demonstrated that zirconia can limit agglomeration on the support, promoting Ni particles' dispersion and stabilization. Combined with the results of BET and XRD, 10Ni4Zr/HA has a specific surface area of 19.7 m2/g and the best Ni dispersion on the HA support, which could be one of the main reasons for the highest catalytic activity. At 400 °C, the catalyst achieved 92.9% CO2 conversion, 8.3% higher than the unmodified one. However, when the denaturation content increased from 4 to 6 wt%, the CO2 conversion decreased, although the ZrO2 addition reduced the crystal size. This could be due to decreases in the Ni0 active site (calculated by the hydrogen consumption (H) in Fig. 1d) as well as the dispersion of catalyst shown in the SEM images (Fig. 2) when increasing the ZrO2 content up to 6 wt%, suggesting that the ZrO2 content has an optimal range in affecting the catalyst performance. When the ZrO2 content exceeds the required value, the activity of the catalyst decreases due to the influence of the specific surface area of the catalyst, the dispersion of Ni decreases, or the size of NiO crystals increases too much higher than the value desired. In the case of the catalyst containing higher ZrO2 content, the lower amount of Ni can prevent the conversion of CO2 to CO via RWGS, and the right amount of ZrO2 will benefit methane formation. This fact proved that 10Ni4Zr/HA is one of the catalysts having high activity for CO2 methanation.
The enhanced interaction between Ni and Zr not only promotes the overflow of H spilling from Ni to the support, which is beneficial for CO2 hydrogenation and intermediates but also alleviates the clogging of the pore structure by restricting the average size of the Ni particles on the used catalysts to a size smaller than the average pore size. In addition, the increasing Ni particle size improved the low-temperature CO2 conversion by enhancing the CO2 adsorption of the catalysts to OCO*. The addition of Zr inhibited carbon deposition by reducing the number of active metal sites and increasing the surface oxygen, slightly improving the CH4 selectivity. However, adding too much Zr will reduce the catalyst's reactivity because the irregular pore structure increases the Ni particle size, leading to the ability to adsorb OCO* too strongly and reducing the activity of the catalyst. Another property leading to the distinction in product selectivity in the temperature range of 250–––300 °C (Fig. 4b) could be the CO2 adsorption capacity, relating to the catalysts’ CO2-TPD profiles. While suitable amount of ZrO2 with moderate basic sites could promote the formation of CH4 in the low temperatures, the excess loading of ZrO2 inducing a lower Ni density on the surface could favor the formation of nickel carbonyl, which later convert into CO in the product flow. The slow kinetic nature of the reaction in the low temperature combining with more CO2 adsorbed on the catalyst’s surface than H2 might take responsibility for the generation of the metallic carbonyl species (Aziz et al., 2015; Wen et al., 2021). Therefore, adding an appropriate amount of Zr can improve the reactivity of the catalysts. Within the scope of this work, the 10Ni4Zr/HA catalyst showed optimal CO2 conversion. A schematic of CO2 methanation on ZrO2-doped nickel nanocatalyst supported on hydroxyapatite from Salmon bone could be represented in Fig. 5.

- Schematic representation of CO2 methanation on NiZr/HA catalyst.
Reaction mechanism investigations were carried out to determine how methane is generated on the Ni/HA catalyst surface. In recent years, the CO2 methanation reaction mechanism has been intensively investigated in order to design a more active and selective catalyst for the synthesis of methane (Wai et al., 2020). Pan et al. proposed a reaction mechanism in which formate species is the primary reaction intermediate while carbonates and hydrogen carbonates are formed in CO2 methanation (Jia et al., 2020). As shown in [steps (1)–(11)], a reaction mechanism for CO2 methanation over Ni/HA catalyst is postulated on the basis of these previous studies. The H2 molecule splits into two hydrogen atoms (step 1) on the Ni surface (*). In steps (2)-(5), gaseous CO2 chemisorbs and reacts with surface hydroxyl (OH*) to form hydrogen-carbonates (HCO3*). The hydrogen-carbonate species becomes formate (HCOO*) by reacting with H atoms on the surface (step (6)). In step 4, surface H atoms react with surface O atoms to replenish the surface hydroxyl groups. Through additional hydrogenation (steps (7)-(8)), surface formate also transforms into OCH2OH* species and then into methoxy species (OCH3*). The OCH3* species is then hydrogenated (step (9)) to form CH4. Additionally, methoxy formation is also associated with the presence of Lewis acidic sites (Kim et al., 2019). Since HA catalytic support is utilized in this work for CO2 methanation, the Ca (II) sites in the HA unit cell structure are Lewis acidic (Yook et al., 2022), indicating that the HA support plays a role in stabilizing methoxy species generation and influencing methane synthesis.
To provide insight into the suitability of the 10Ni4Zr/HA catalyst, we compare its efficacy to that of various Ni catalysts at the same of reaction temperature and the CO2: H2 molar ratio reported in the literature in Table 2. Numerous supports, including SBA-15 (Van et al., 2020); MCM-41 (Yang et al., 2020), CNTs (Romero-Sáez et al., 2018); SiO2 (Guo and Lu, 2014), Al2O3 (Mutz et al., 2017), Pr2O3–CeO4, Sm2O3–CeO4 (Siakavelas et al., 2021), etc have been investigated. In the majority of cases, Ni content higher than 10% wt were supported on them. Despite their relatively modest Ni content, we may conclude that they can compete favourably with these reference formulations. Because HA can self-exchange cations and anions, this can be explained by the characteristic properties of HA. The active metal centers for the reaction can be introduced into the HA structure via ion exchange, and Ca2+ ions give the HA surface a weak base nature with the ability to dissociate CO2 and the OH group to stabilize the active metal, which promotes CO2 hydrogenation to CH4. In addition to the active metal's dispersion, the catalyst's fundamental properties and the amount of surface oxygen in the catalyst significantly influence its reactivity. Adding alkaline earth metals (Ca present in HA) can increase the number of active metal sites on the catalyst that influence CO2 adsorption and activation. Moreover, it is evident that the positive effect of Zr on the activity of the Ni catalyst is more pronounced in our HA-based materials, indicating that our green support is suitable for the CO2 methanation reaction. It should be emphasized that the enhanced structural, chemical, and catalytic properties of Ni-based catalysts supported on hydroxyapatite-derived Salmon bone waste represent an advantageous alternative to conventional basic catalysts.
Catalyst | Ni, wt% | WGHSV, mL/(g.h) | CCO2, mol% | XCO2, % | YCH4, % | H, molCO2/(molNi.h) | Y, molCH4/(molNi.h) | Refs. |
---|---|---|---|---|---|---|---|---|
50Ni4Ce/SBA | 50 | 12,000 | 20 | 90 | 99 | 12.5 | 12.4 | (Van et al., 2020) |
10Ni/@Zr(0.01)@MCM-41 | 10 | 6,420 | 10 | 88.6 | 96.5 | 16.5 | 15.9 | (Yang et al., 2020) |
5Ni/Ce0.85Zr0.15O2 | 5 | 10,000 | 16 | 60 | 100 | 55.7 | 55.7 | (Iglesias et al., 2019) |
5Ni/Zr/CNTs-COI | 5 | 75,000 | 10 | 20 | 70 | 87.1 | 60.9 | (Romero-Sáez et al., 2018) |
10Ni/Ce | 10 | 25,000 | 10 | 39.4 | 100 | 28.6 | 28.6 | (Siakavelas et al., 2021) |
10Ni/Sm-Ce | 10 | 25,000 | 10 | 44.9 | 100 | 32.6 | 32.6 | (Siakavelas et al., 2021) |
10Ni/Pr-Ce | 10 | 25,000 | 10 | 54.5 | 100 | 39.5 | 39.5 | (Siakavelas et al., 2021) |
10Ni/Mg-Ce | 10 | 25,000 | 10 | 43.2 | 100 | 31.3 | 31.3 | (Siakavelas et al., 2021) |
20Ni3Fe/Al2O3 | 20 | 134,000 | 10 | 82 | 99.0 | 15.9 | 15.8 | (Mutz et al., 2017) |
10Ni4Ca/SiO2 | 10 | 15,000 | 10 | 64.9 | 98.6 | 28.2 | 27.9 | (Guo and Lu, 2014) |
20NiAl-LDH | 20 | 2,400 | 18.75 | 75.1 | 95.0 | 4.9 | 4.7 | (He et al., 2013) |
10Ni4Zr/HA | 10 | 12,000 | 3 | 76.6 | 99.8 | 8.0 | 8.0 | This work |
H – CO2 conversion efficiency calculated per Ni mole, molCO2/(molNi.h); Y - CH4 yield calculated per Ni mole, molCH4/(molNi.h).
3.3 The stability of 10Ni4Zr/HA catalyst
The stability test, performed at 350 ◦C, reveals a robust performance of the optimized 10Ni4Zr/HA catalyst very notable stability during long time on stream (Fig. 6a). After an activation period of 50 h of reaction, reaching a CO2 conversion close to 76% and a high selectivity towards methane production (>95%). In addition, compared with the characterizations of the initial fresh 10Ni4Zr/HA catalyst (Fig. 2c and 3b), the characterizations of the spent one (Fig. 6b-d) revealed insignificant differences in their properties and the distribution of active phases on the HA surface based on the impregnation technique employed.

- The activity stability of the best catalyst (10Ni4Zr/HA) in CO2 methanation at 350 °C during 50 h on stream (a), XRD pattern (d), and SEM (b) and TEM (c) images of spent catalyst.
4 Conclusion
A simple and environmentally friendly approach was developed for synthesizing Ni/HA and ZrO2-doped Ni/HA catalysts for CO2 methanation. Ca2+ ions in the HA structure promoted CO2 adsorption, increasing catalytic activity. Adding the zirconia to the Ni/HA catalyst enhances NiO dispersion on the HA surface, leading to smaller NiO crystals, an increased NiO reduction and CO2 adsorption capacity, and better catalytic activity. At 400 °C, the 10Ni/HA catalyst-doped 4 wt% ZrO2 showed outstanding catalytic performance, with a CO2 conversion of 92.9% and a CH4 selectivity of almost 100%. As a result, one of the most promising techniques in developing effective catalyst designs toward the development of effective CO2 hydrogenation to CH4 fuel on an industrial scale is employing a zirconia-doped Ni/HA catalyst.
CRediT authorship contribution statement
Ngoc Doan Trang Tran: Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Thi Ngoc Han Che: Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Thi Thuy Van Nguyen: Conceptualization, Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Ba Long Do: Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Thanh Gia-Thien Ho: Conceptualization, Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Phung Anh Nguyen: Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Thi Thuy Phuong Pham: Conceptualization, Methodology, Resources, Writing – original draft, Validation, Investigation, Formal analysis. Nguyen Tri: Conceptualization, Methodology, Writing – review & editing, Supervision. Huynh Ky Phuong Ha: Conceptualization, Methodology, Writing – review & editing, Supervision.
Acknowledgments
This study is supported by the Vietnam Academy of Science and Technology under the NCXS02.02/22-23 project.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
References
- CO2 methanation over heterogeneous catalysts: Recent progress and future prospects. Green Chem.. 2015;17:2647-2663.
- [Google Scholar]
- The effect of noble metals on catalytic methanation reaction over supported Mn/Ni oxide based catalysts. Arab. J. Chem.. 2015;8:632-643.
- [Google Scholar]
- Influence of Ca/P ratio on the catalytic performance of Ni/hydroxyapatite samples in dry reforming of methane. Appl. Catal. A: Gen.. 2019;580:34-45.
- [Google Scholar]
- Influence of Ca/P ratio on the catalytic performance of Ni/hydroxyapatite samples in dry reforming of methane. Appl. Catal. A: Gen.. 2019;580:34-45.
- [Google Scholar]
- Study on the promotional effect of lanthana addition on the performance of hydroxyapatite-supported Ni catalysts for the CO2 methanation reaction. Appl. Catal. B: Environ.. 2022;314:121500
- [Google Scholar]
- Methanation of carbon dioxide on Ni/ZrO2-Al2O3 catalysts: Effects of ZrO2 promoter and preparation method of novel ZrO2-Al2O3 carrier. J. Nat. Gas Chem.. 2011;20:318-324.
- [Google Scholar]
- Syngas production via CO2 reforming of CH4 over Ni-based SBA-15: Promotional effect of promoters (Ce, Mg, and Zr) Mater. Today Energy. 2019;12:408-417.
- [Google Scholar]
- Dry reforming of methane using Ni–Ce catalysts supported on a modified mineral clay. Appl. Catal. A: Gen.. 2009;364:65-74.
- [Google Scholar]
- Identification of surface basic sites and acid–base pairs of hydroxyapatite. J. Phys. Chem. C. 2014;118:12744-12757.
- [Google Scholar]
- Fractal analysis of hydroxyapatite from nitrogen isotherms. J. Colloid Interface Sci.. 2004;277:410-416.
- [Google Scholar]
- The difference of roles of alkaline-earth metal oxides on silica-supported nickel catalysts for CO2 methanation. RSC Adv.. 2014;4:58171-58177.
- [Google Scholar]
- A surface defect-promoted Ni nanocatalyst with simultaneously enhanced activity and stability. Chem. Mater.. 2013;25:1040-1046.
- [Google Scholar]
- Zr promotion effect in CO2 methanation over ceria supported nickel catalysts. Int. J. Hydrog. Energy. 2019;44:1710-1719.
- [Google Scholar]
- Selective hydrogenation of CO2 to methanol over Ni/In2O3 catalyst. J. Energy Chem.. 2020;50:409-415.
- [Google Scholar]
- Role of oxygen vacancy in metal oxides for photocatalytic CO2 reduction. Appl. Catal. B: Environ. 2022122079
- [Google Scholar]
- Surface lewis acidity of periphery oxide species as a general kinetic descriptor for CO2 hydrogenation to methanol on supported copper nanoparticles. ACS Catal.. 2019;9:10409-10417.
- [Google Scholar]
- Intermetallic compound nanoparticles dispersed on the surface of oxide support as active and selective catalysts. Mater. Trans.. 2015;56:460-467.
- [Google Scholar]
- Property based ranking of CO and CO2 methanation catalysts. Energy Procedia. 2017;128:255-260.
- [Google Scholar]
- Effect of pH and temperature on the morphology and phases of co-precipitated hydroxyapatite. J. Solgel Sci. Technol.. 2012;61:592-599.
- [Google Scholar]
- Recent trend in thermal catalytic low temperature CO2 methanation: A critical review. Catal. Today. 2021;368:2-19.
- [Google Scholar]
- CO2 hydrogenation to methanol over ZrO2-containing catalysts: insights into ZrO2 induced synergy. ACS Catal.. 2019;9:7840-7861.
- [Google Scholar]
- Improved SAPO-34 membranes for CO2/CH4 separations. Adv. Mater.. 2006;18:2601-2603.
- [Google Scholar]
- Experimental and DFT study on single atom solution for carbon dioxide methanation. Fuel. 2023;351:128911
- [Google Scholar]
- CO2 methanation: principles and challenges. In: Studies in Surface Science and Catalysis. Elsevier; 2019. p. :85-103.
- [Google Scholar]
- Methanation of CO2 on Ni/γ-Al2O3: influence of Pt, Pd or Rh promotion. Catal. Today. 2018;306:294-299.
- [Google Scholar]
- Recent advance in nanostructured materials innovation towards photocatalytic CO2 reduction. Appl. Catal. A: Gen.. 2022;648:118927
- [Google Scholar]
- Potential of an alumina-supported Ni3Fe catalyst in the methanation of CO2: Impact of alloy formation on activity and stability. ACS Catal.. 2017;7:6802-6814.
- [Google Scholar]
- Multidisciplinary green approaches (ultrasonic, co-precipitation, hydrothermal, and microwave) for fabrication and characterization of Erbium-promoted Ni-Al2O3 catalyst for CO2 methanation. Arab. J. Chem.. 2023;16:104494
- [Google Scholar]
- Hydroxyapatite - A promising sunscreen filter. J. Aust. Ceram. Soc.. 2020;56:345-351.
- [Google Scholar]
- Insight into the reaction route of CO2 methanation: Promotion effect of medium basic sites. Catal. Commun.. 2014;45:74-78.
- [Google Scholar]
- Carbon dioxide reforming of methane over co-precipitated Ni-CeO2, Ni-ZrO2 and Ni-Ce-ZrO2 catalysts. Catal. Today. 2004;93:39-44.
- [Google Scholar]
- CO2 methanation over nickel-ZrO2 catalyst supported on carbon nanotubes: A comparison between two impregnation strategies. Appl. Catal. B: Environ.. 2018;237:817-825.
- [Google Scholar]
- 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]
- Synthesis and characterization of zirconia nanocrystalline powder by thermal treatment method. Mater. Today:. Proc.. 2022;54:908-911.
- [Google Scholar]
- Phytogenic generation of NiO nanoparticles using Stevia leaf extract and evaluation of their in-vitro antioxidant and antimicrobial properties. Biomolecules. 2020;10:89.
- [Google Scholar]
- Catalytic carbon dioxide hydrogenation to methane: A review of recent studies. J. Energy Chem.. 2016;25:553-565.
- [Google Scholar]
- Synthesis of a novel beta-tricalcium phosphate/hydroxyapatite biphasic calcium phosphate containing niobium ions and evaluation of its osteogenic properties. J. Artif. Organs. 2007;10:22-28.
- [Google Scholar]
- Hydrothermal and calcination regimes and characteristics of nanohydroxyapatite synthesized from Salmon bones. J. Biochem. Technol.. 2020;11:82-87.
- [Google Scholar]
- Positive effects of CeO2 promoter and co-reactant/CO on methanation of CO2-rich gas over Ni/SBA-15 catalyst. Mater. Trans.. 2020;61:1332-1338.
- [Google Scholar]
- Influence of surface formate species on methane selectivity for carbon dioxide methanation over nickel hydroxyapatite catalyst. ChemCatChem. 2020;12:6410-6419.
- [Google Scholar]
- Performance analysis of sabatier reaction on direct hydrogen inlet rates based on solar-to-gas conversion system. Int. J. Hydrog. Energy. 2021;46:26801-26808.
- [Google Scholar]
- Biomass derived N-doped biochar as efficient catalyst supports for CO2 methanation. J. CO2 Util.. 2019;34:733-741.
- [Google Scholar]
- Exploring the influence of nickel precursors on constructing efficient Ni-based CO2 methanation catalysts assisted with in-situ technologies. Appl. Catal. B: Environ.. 2021;297:120486
- [Google Scholar]
- CO2 methanation over nickel-based catalysts supported on MCM-41 with in situ doping of zirconium. J. CO2 Util.. 2020;42:101304
- [Google Scholar]
- Effect of controlling thiophene rings on DA polymer photocatalysts accessed via direct arylation for hydrogen production. Molecules. 2023;28:4507.
- [Google Scholar]
- Design strategies for hydroxyapatite-based materials to enhance their catalytic performance and applicability. Adv. Mater.. 2022;2204938
- [Google Scholar]
- The determinants of renewable energy sources for the fueling of green and sustainable economy. Energy. 2022;238:122029
- [Google Scholar]