Hydrogenation of Glucose on a Carbon-Supported Ru Catalyst: Optimization of the Reaction Conditions ()
1. Introduction
The use of renewable raw materials and the development of different strategies for the valorization of biomass, together with heterogeneous catalysis, are interesting alternatives for the substitution of fossil fuel resources in order to obtain intermediates and final chemical products of high added value [1] . Cellulose is the most abundant carbohydrate and it is inedible. For this reason, it is expected to be the first objective of the chemistry of biomass resources [2] . By degrading cellulosic biomass, mono-, di- and oligosaccharides can be obtained, which in turn can be transformed into valuable final chemical products or intermediates [3] . Among the most important chemical products derived from cellulosic biomass are sugar alcohols.
In 2004, the US Department of Energy [4] identified D-sorbitol as one of the 12 most important value-added chemicals that can be obtained from biomass. D-sorbitol is a potential source of alkanes for liquid biofuels [5] . Esters and other derivatives of D-sorbitol are also important additives and intermediates in the cosmetics and pharmaceutical industries, with the molecule being a key substitute for sugar in the food and beverage industry [6] .
Sorbitol is the main product of glucose hydrogenation [4] . However, while hydrogenation of D-glucose to D-sorbitol seems like a simple reaction, D-glucose is not only selectively converted to D-sorbitol, but it can also give many other byproducts (see Figure 1). Some of these by-products can even cause catalyst deactivation [7] . Sorbitol, like other hexitols, can also form chelate complexes with nickel (or Cu, Fe) readily in an aqueous solution. These complexes are stable under alkaline conditions, which results in nickel being left in the sorbitol solution. [8] [9] [10] .
Among the possible competitive reactions, the Lobry rearrangement of Bruyn-Alberda van Ekenstein of D-glucose results in the formation of D-mannose and D-fructose. Aldehydes (glyceraldehyde, formaldehyde) and
Figure 1. Reaction network of the hydrogenation of d-glucose.
ketones (dihydroxyacetone) could originate by cleavage in an alkaline medium of sugar molecules. This reaction occurs via an enediol intermediate via a proton-transfer mechanism [11] . The hydrogenation in the liquid phase of D-glucose generally leads as a main product to D-sorbitol but alkaline conditions could favor both the procedures of glucose isomerization and sorbitol hydrogenolysis, giving mannitol as main by-product and lower alcohols such as propanediol and ethylene glycol [12] . On the other hand, under acidic conditions, 5-hydroxymethylfurfural and levulinic acid can be formed by dehydration of glucose [13] . The way to adjust the selectivity to the desired product is through a correct selection of the active metal, the catalytic support, and the reaction conditions such as temperature and hydrogen pressure.
Historically, glucose has been hydrogenated to sorbitol over Ni catalysts. One of the advantages of the use of a Ni catalyst is its relatively low price compared to other more expensive catalytic systems. However, Ni presents the problem of leaching and as a consequence, purification of the obtained sorbitol is necessary [14] . It has been found that ruthenium-based catalysts outperform nickel-based catalysts due to their high activity, lower load requirement, and lower deactivation [8] [14] [15] [16] . Clearly, the price of metallic ruthenium as a raw material is much higher than that of nickel. However, other issues must be taken into account when comparing the final cost of each catalyst. Thus, the metallic load of the catalysts and the possibility of their reuse must be analyzed. Thus, although it is more expensive than nickel, it may be of interest for industrial uses.
The kinetics of glucose hydrogenation over different Ru catalysts varying reaction conditions have been recently reviewed by Ahmed and Hameed [17] . The authors described published results of kinetic data and the effect of the initial concentration of glucose, temperature and pressure on the rate constants and the selectivity to sorbitol. However, most of these results were obtained at high values of pressure and temperature.
As a contribution, in this work, the aqueous phase hydrogenation (APH) of glucose on a 3 wt% Ru/C catalyst is studied, establishing mild conditions of temperature, pressure, and agitation which are lower than the ones generally reported and optimize the obtention of the desired product. Kinetic modeling was developed using models based on Langmuir-Hinshelwood-Hougen-Watson (LHHW) kinetics.
2. Experimental
2.1. Catalyst Preparation and Characterization
A Ru catalyst, with a metal loading of 3 wt%, was prepared by impregnation with excess solution on a commercial activated carbon (C) (NORIT, ground, and sieved to 60 - 100 mesh) according to reference [18] .
Textural properties of the support, pore volume (VP), and specific surface area (SBET) were established by N2 physisorption at 77 K using a Micromeritics ASAP 2020 equipment.
The size distribution of metallic particles was determined by transmission electron microscopy (TEM) with a JEOL 100 CX II microscope. The mean particle size (dTEM) was calculated using Equation (1):
(1)
where di is the diameter of particle and ni is the number of particles with size di.
The reducibility of the catalyst was studied by Temperature-Programmed Reduction (TPR) experiments using a laboratory-constructed equipment and X-ray photoelectron spectroscopy (XPS) analysis were performed on a Specs Multi-technique system (SPECS) equipped with a dual Mg/Al X-ray source and a PHOIBOS 150 hemispherical analyzer operating in the Fixed Analyzer Transmission (FAT) mode. Both characterizations were conducted following previously published protocols [18] .
2.2. Catalytic Tests
The tests were carried out in a high-pressure reactor, using 0.25 g of catalyst and 50 mL of water as solvent. The reaction time was 7 h. The stirring speed, the temperature, the pressure, and the initial glucose concentration were varied between 250 - 700 rpm, 343 - 383 K, 0.5 - 2 MPa, and 0.033 - 0.133 M, respectively. The variables were modified one at a time, keeping the others constant in each test.
The reaction samples were analyzed offline by liquid chromatography in a UHPLC chromatograph equipped with a RI-detector (operating at 313 K). Separation of the components was achieved by a Phenomenex Rezex RCM Ca+2 Monosaccharide (300 × 7.8 mm) column operated at 353 K. Elution was performed at 0.6 mL/min flow rate of mobile phase (milli-Q-water).
The conversion of glucose (Glu) was calculated using the following expression:
(2)
Analytical criteria were verified to determine the absence of diffusional limitations and the experimental results were interpreted by kinetic modeling using heterogeneous models of Langmuir-Hinshelwood-Hougen-Watson (LHHW).
3. Results and Discussion
3.1. Catalyst Characterization
The activated carbon used as support had a pore volume of 0.53 cm3∙g−1 and a SBET of 1000 m2∙g−1 measured by N2 physisorption. The 3 wt% Ru/C catalyst has been deeply characterized, as previously published [18] . In order to better explain the kinetic results, the main characteristics of the catalyst are gathered in Table 1.
Briefly, it can be said that TEM analysis showed the presence of small and
Table 1. Ru 3 wt%/C catalyst preparation and characterization.
a.
[19] .
uniform-sized particles that are widely dispersed, and Ru is in the metallic form as was confirmed by XPS and TPR analysis.
3.2. Mass Transfer Limitation
Analytical criteria were used to compare the initial reaction rate with the maximum transport rates in the physical stages. In the case of liquid-phase hydrogenation reactions, the diffusion of reagents and/or products in the gas-liquid interface, the liquid-solid interface, and internal diffusion in the particles must be evaluated.
The reaction conditions used to determine the initial reaction rate were 363 K, 1.25 MPa of H2 pressure, and a stirring speed of 625 rpm, using 0.25 g of 3 wt% Ru/C catalyst and 0.90 g of glucose as substrate. From the conversion vs. time curve, the initial reaction rate was calculated at 10% conversion.
3.2.1. Gas-Liquid Mass Transfer Limitations
To ensure that there is no resistance to H2 transport at the gas-liquid interface, the criterion determined by the ratio between the observed reaction rate and the maximum transfer rate based on the liquid volume (
) was used.
(3)
The criterion states that, if the reaction rate observed is less than 10% of the maximum gas-liquid mass transfer rate, then the transport process is much faster than the observed chemical reaction, thus not limiting the overall process.
For the calculation of the
coefficient several correlations determined by different authors can be used. One of them is the proposed by Meille et al. [20] who studied the determination of the
coefficient in laboratory tank reactors, with magnetic stirring and without the use of baffles, with volumes between 25 and 300 cm3 (systems with similar properties to those used in this work). Considering a reactor of 150 cm3 of capacity, which is the most similar to the one used in this work, the correlation is:
(4)
The results gathered in Table 2 show that under the conditions employed in the present work, there are no limitations for the H2 transport at the gas-liquid interface and within the reaction volume.
3.2.2. Liquid-Solid Mass Transfer Limitations
In order to determine the absence of mass transfer limitations of the reagents in the liquid-solid interface, the selected criterion is similar to that seen for the gas-liquid interface. In this case (
) for both reagents (H2 and glucose), the ratio between the reaction rate and the mass transport rate must be less than 10% (Equation (5)).
(5)
As can be seen in Table 2, the criterion is satisfactorily fulfilled for the conditions tested, for both glucose and H2.
3.2.3. Internal Diffusion Limitation
To verify if there is resistance to mass transfer inside the catalyst particles, the Weisz-Prater criterion was used [21] [22] for both reactants (H2 and glucose). According to this criterion, the mass transfer limitations are negligible if it is satisfied that:
(6)
For both H2 and glucose, the Weisz-Prater criterion is fulfilled, since for both reactants the value of Φ is much smaller than 1 (see Table 2). Thus, it can be said that there are no limitations to the intraparticle mass transfer.
3.3. Influence of Stirring Speed
Agitation of the reaction mixture in a slurry reactor is important not only to keep the catalyst particles uniformly suspended, but also to facilitate the transport of the gaseous reactant through the gas-liquid interface and the dissolved reagent in the liquid medium towards the surface of the catalyst, there being a
Table 2. Parameters and correlations used in the verification of transport limitation.
relation between the increase of the stirring speed and the increase of the mass transfer coefficient.
In order to ensure the absence of mass transfer limitations in the gas-liquid interface, experiments were carried out varying the stirring speed between 250 and 700 rpm at 363 K and 1.25 MPa of H2 pressure, with 0.25 g of Ru/C catalyst and 0.90 g of glucose as substrate. Figure 2 shows the initial reaction rates (ri) for each stirring speed, calculated from the slopes of the conversion vs. time curves (at 10% conversion). With stirring speeds lower than approximately 600 rpm, (ri) decreases, indicating the presence of limitations to the transfer of hydrogen at the gas-liquid interface. For this reason, a speed of 625 rpm was adopted for subsequent tests.
3.4. Influence of Reaction Parameters
The influence of temperature on the catalytic activity at 1.25 MPa H2 and 0.1 M of glucose was investigated (Figure 3), varying its value from 343 K to 383 K.
Below 353 K, the reaction rate decreases significantly, and only 80% glucose conversion is reached at 343 K after 7 h of reaction. Above 373 K, it is observed that the time necessary to achieve complete conversion is reduced to 5 h. However,
Figure 2. Effect of the stirring speed on the initial D-glucose hydrogenation rate. [
,
,
,
;
].
Figure 3. Effect of temperature on glucose APH with Ru (3%)/C. In symbols, experimental results and solid lines the results using the best model (3a). [
,
,
;
].
in the total time studied (7 h), at 363 K, 99% conversion is achieved. Thus, although time is longer, high conversion is achieved in milder condition. In all the cases, the selectivity was 100% to sorbitol, proving that consecutive or parallel reactions such as isomerizations do not occur (see Figure 1).
The apparent activation energy (Ea) was determined by linear regression using an Arrhenius type function. A value of 30.8 kJ/mol was obtained. This value is similar to those reported in literature [17] . Thus, kinetic data using different ruthenium catalysts and various experimental conditions have been reported. In this sense, an activation energy of 32.9 kJ/mol and first-order kinetics with respect to glucose concentration has been reported by Mishra et al. [23] over Ru/HY zeolite. Crezee et al. [24] reported an activation energy of 55 kJ/mol at 4 MPa of H2 and 0.75 g/L of 5% Ru/C. Similar values have also been reported in studies carried out with other monosaccharides. Li et al. [25] reported that the apparent activation energy (Ea) of maltose hydrogenation using a Ru-P catalyst was 27 kJ/mol (6.45 kcal/mol) and for Ru-B was 32 kJ/mol (7.65 kcal/mol), using a maltose solution of 40% w/v, 363 K, 2 MPa pressure and stirring at 1200 rpm.
The effect of partial pressure of hydrogen (at 363 K and 0.1 M glucose) is presented in Figure 4.
It is observed that in every case an almost total conversion was obtained, increasing
with the increase of pressure. The variation of the pressure does not produce changes in the selectivity, obtaining sorbitol as the only product.
The reaction order with respect to hydrogen was calculated by considering for
a power law rate equation:
(7)
The linear regression of data obtained shows a reaction order with respect to the H2 of 0.586. This value was in reasonable agreement with that reported by Zhang et al. of 0.66 over silica-supported Pt nanoparticles [6] and also with that informed by Tukac of 0.65 [26] . This latter value was obtained in the catalytic hydrogenation of glucose, using a 40 wt% aqueous solution of D-glucose and a
Figure 4. Effect of H2 pressure on glucose APH with Ru (3%)/C. In symbols, experimental results and in solid lines the results using the best model (3a). [
,
,
;
].
kieselgur-supported nickel catalysts (12% NiO, 2% Cr2O3). The temperature range was 388 to 433 K and pressure range of 0.5 - 10 MPa [19] .
The results of the variation of the initial glucose concentration (at 363 K and 1.25 MPa) are given in Figure 5.
Although the final conversion reached after 7 h of reaction decreases, the activity of the catalyst increases with increasing glucose concentration. The reaction order with respect to glucose was determined by linear regression from the data with Equation (8) and a value of 0.406 was obtained. In the literature, it has been reported that glucose concentration presents a relatively complex influence upon the hydrogenation rate. Guo et al. [27] studying a Ru-B amorphous catalyst, observed that increasing glucose concentration the initial rate increased linearly up to 40 wt% and then presented a plateau at higher glucose concentration. This kinetic behavior could be explained by the adsorption strength of glucose and hydrogen on the catalyst. The glucose molecule adsorbs strongly on the catalyst reaching saturated adsorption even at low concentration. Only when the glucose concentration was very low in the liquid phase, and in consequence, the surface adsorption was unsaturated, the hydrogenation rate grew up with the glucose concentration.
The catalytic hydrogenation of D-glucose is generally described by means of Langmuir-Hinshelwood kinetics with a change from first-order dependency in D-glucose at low concentrations to zero-order behavior at high concentrations [8] .
3.5. Kinetic Modeling
Based on the above-discussed results, the following hypotheses for the formulation of the LHHW models were considered:
1) The superficial hydrogenation reaction is irreversible (total conversion of glucose to sorbitol was obtained in the experiments performed).
2) Hydrogen concentration in the liquid phase is constant, due to the constant pressure of hydrogen throughout the test, the high volume of solvent and the
Figure 5. Effect of initial glucose concentration on glucose APH with Ru (3%)/C. In symbols, experimental results and in solid lines the results using the best model (3a) [
,
,
,
].
efficient mixing.
Three general models were proposed: a first model of the Eley-Rideal type, with only glucose adsorption on the solid catalyst; a second model with adsorption of molecular hydrogen in addition to glucose adsorption and a third model considering the dissociative adsorption of hydrogen. Both models, 2 and 3, considered the competitive adsorption of H2 and glucose. By assuming different rate limiting steps (r.l.s.), (H2 adsorption, glucose adsorption or surface chemical reaction) and the possibility or not of sorbitol (the reaction product) adsorption, 19 different LHHW sub models were analyzed.
The equations were solved numerically using the algorithm of Runge-Kutta-Merson with the data obtained at 363 K. The estimation of model parameters was performed by non-linear regression, using a Levenberg-Marquardt algorithm that minimizes the objective function:
(8)
Table 3 presents the proposed LHHW models, with the estimated parameters
Table 3. Overview of model equations and kinetic modeling results using LHHW models.
and the correlation coefficients (r). The discrimination between the models was determined, after discarding those with parameters without physical sense, according to the highest r value and using the model selection criterion (MSC), according to the following equation:
(9)
The most significant model is the one that leads to the highest MSC value.
Models 2c, 2d, 2f, 3b, 3c, 3d and 3f, have some of their estimated parameters with negative and/or imaginary values, which made them lack of physical sense for the kinetic and equilibrium constants. Because of that, these models are rejected.
In addition, those models that consider that sorbitol is not adsorbed, or that the desorption of this product is the rate limiting step (1b, 1d, 1e, 2b, 2g, 3g) have a poor fit. The adsorption of sorbitol cannot be neglected, however its desorption does not constitute a limitation in the reaction rate (0 << 100 in the models with the best fit).
The best model among the ones here proposed to interpret the aqueous phase hydrogenation of glucose with 3 wt% Ru/C catalyst (Model 3a) considers that both reagents, glucose and H2 are adsorbed on the catalyst (competitive adsorption) to give sorbitol, which is then desorbed. The H2 adsorption is dissociative and the rate limiting step is the surface chemical reaction. Equations (i)-(iv) show the elementary steps describing this process, where (iii) is the r.l.s.:
(i)
(ii)
(iii)
(iv)
In Figure 6, the experimental results for the conversion of glucose as a function of time with the optimal conditions are presented in symbols, while the values
Figure 6. (Left) Aqueous phase hidrogenation (APH) of glucose with Ru (3%)/C in the optimal conditions determined. In symbols, experimental results and in solid lines the results using model 3a. (Right) Evolution of residuals. [
,
,
,
;
].
estimated by this model are represented in solid lines. As can be seen, a good agreement between the experimental and predicted data is obtained, which validates the kinetic model.
Figure 7 shows the model fit of the other conditions evaluated. Along with this, the distribution of residuals as a function of time is also shown and followed an acceptable random trend, being consistent with the random error hypothesis included in the regression and giving additional support to the adequacy of the model.
Subsequently, the dependence of the kinetic and equilibrium constants with temperature was established using the following expressions:
(10)
(11)
(12)
(13)
Figure 7. (Left) Aqueous phase hidrogenation (APH) of glucose with Ru (3%)/C. In symbols, experimental results and in solid lines the results using model 3a. (Right) Evolution of residuals. [
,
;
].
And the model 3a was adjusted with the data obtained for 343 - 383 K. The results are shown in Table 4, while the fit obtained for the different temperatures evaluated can be observed in Figure 8.
Table 4. Modeling results using all experiments and including the temperature dependency.
Figure 8. (left) Aqueous phase hidrogenation (APH) of glucose with Ru (3%)/C at different temperatures. In symbols, experimental results and in solid lines the results using model 3a. (Right) Evolution of residuals. [
,
,
;
].
The D-glucose adsorption enthalpy HA presents a positive value. This is indicating that the adsorption constant of D-glucose increases with temperature. This is not thermodynamically consistent with adsorption phenomena. However, this behavior has been described by other authors in terms of the mechanism of D-glucose adsorption [24] .
From the kinetic studies carried out, and considering the reaction mechanism proposed by model 3a, Figure 9 represents a surface of ruthenium atoms supported on carbon on which hydrogen is adsorbed dissociatively and glucose does so through its hemiacetal group, leading to the formation of sorbitol, following the elementary steps described from (i) to (iv).
Figure 9. Elementary steps of the chemical reaction for proposed model 3a to the hydrogenation of glucose on Ru (3%)/C.
4. Conclusions
The 3 wt% Ru/C catalyst studied was highly selective towards sorbitol under the conditions tested (363 K and 1.25 MPa). These experimental conditions are milder than the ones generally reported for the APH reaction.
The kinetic study of D-glucose hydrogenation was carried out. The experimental data could be interpreted through LHHW heterogeneous models. The best fit was obtained when the competitive adsorption of both glucose and hydrogen was considered. Hydrogen is dissociative adsorbed and the rate-limiting step is the surface chemical reaction.
Acknowledgements
The authors would like to acknowledge the following institutions of Argentina for the financial support, UNNOBA (project SIB 2022, Exp. 1664/2022), CONICET (project PIP 0134/2021) and ANPCyT FONCyT (PICT 2019/1962).