J. M.
Jimenez-Martin
a,
M.
El Tawil-Lucas
a,
C.
García-Jerez
a,
J.
Moreno
a,
A.
García
a,
B.
Hernández
a and
J.
Iglesias
*ab
aChemical & Environmental Engineering Group, Universidad Rey Juan Carlos, C/Tulipán s/n. Móstoles, 28933, Madrid, Spain. E-mail: jose.iglesias@urjc.es; Tel: +34 914888565
bInstituto de Tecnologías para la Sostenibilidad, Universidad Rey Juan Carlos, C/Tulipán s/n. Móstoles, 28933, Madrid, Spain
First published on 9th October 2024
Catalytic performance of potassium-exchanged Sn-containing USY zeolite ([K]Sn-USY) for the transformation of hexoses into methyl lactate in methanol has been studied. This work explores the effects of temperature and catalyst loading on this transformation process, shedding light on the kinetic aspects of this process in which the starting sugars undergo a complex network of different transformations to yield different biobased products with commercial interest. The proposed kinetic model integrates several key reactions: isomerization of the starting sugars, glycosidation with the alcohol solvent, retro-aldol cleavage of carbohydrates to minor sugars, and hydrolytic side transformations. By considering these steps, insights into the determining factors that influence overall transformation and high selectivity of the [K]Sn-USY zeolite for methyl lactate are clarified. Notably, this work highlights the formation of side products, including methoxy methyl 2-hydroxy-butanoate (MMHB) and glycolaldehyde dimethyl acetal (GADMA), which have been previously overlooked in other kinetic modelling and derived from the retroaldol cleavage of aldoses. The [K]Sn-USY catalysts exhibit promising activity and selectivity, making them attractive candidates for bio-based plastic production. A comprehensive understanding of the catalytic behavior and side product formation contributes to the optimization of these catalysts for sustainable carbohydrate valorization.
Production of alkyl lactates from carbohydrates in an alcohol medium and in the presence of Lewis acid catalysts, such as Sn-containing zeolites, consists of a cascade reaction with several steps. Initially, glucose is isomerized to fructose (a ketose isomer of glucose), enabling the retroaldol cleavage of the C–C bonds between carbons 3 and 4, thereby leading to the formation of glyceraldehyde and dihydroxyacetone. Both trioses can be converted to methyl lactate (MLa); this undergoes dehydration to 2-hydroxyacrolein, which undergoes keto–enol tautomerism to produce pyruvaldehyde. This intermediate is then acetalized with an alcohol solvent to form a hemiacetal, which is readily isomerized through 1,2 hydride shift to yield the corresponding alkyl lactate.17 Along these transformations, the production of a wide variety of side products is feasible, diminishing the yield towards the target alkyl lactate. Such side transformations include glycosidation or dehydration of the starting sugar substrate,4,18 retroaldol cleavage of sugars in aldose form—yielding a tetrose and a diose—19 or the acetalization of the hemiacetal.17
Identifying the optimal conditions of MLA production is challenging since the reaction involves numerous steps and employs highly versatile tin-containing zeolites. One approach lies in the development of kinetic models for Sn-Beta zeolites. Rajabbeigi et al.20 covered the kinetic study of the isomerization between glucose, fructose and mannose in water between 70 and 130 °C, taking into account the deactivation of the catalyst in the kinetic model. Tosi et al.21 conducted the evaluation of the conversion of glucose and fructose to methyl lactate at high temperature (160 °C), introducing the formation of methyl glycosides as side products at the beginning of the reaction. This work demonstrated sugar glycosidation is reversible under the tested conditions, and thus, glycosides act as masking agents of the substrate. Nevertheless, tin-containing zeolites with other structures are demonstrated to have different catalytic performances in the production of methyl lactate.22,23 In those cases, faujasite structures remain scarcely studied despite providing outstanding catalytic activity and higher yields towards alkyl lactates than Sn-β. Under these circumstances, further research is needed to unravel the kinetics and reaction rates in each of the steps for this particular zeolite.
This research focuses on the study of the catalytic performance of a potassium-exchanged Sn-containing USY zeolite in the transformation of glucose, fructose and mannose in methanol media. The study explores the influence of the temperature conditions between 100 and 160 °C and the catalyst loading, covering both isomerization and retro-aldol cleavage promotion ranges, as well as different side transformations at low reaction times (below 120 min). To uncover the determining steps of this transformation, a kinetic model of the detected transformations is proposed. The model integrates for the first time the isomerization of the three sugars and the equilibria with their methyl glycosides, and the formation of side compounds, MLa promoted by Sn catalysts and compounds generated by retro-aldol reactions. Integration of all the reactions helps to determine the influence of the starting substrate, the operating temperature, and the catalyst loading under the operating space. Understanding the effects of these variables on the yields of each product is critical for advancing the development of K-exchanged Sn-USY catalysts.
(1) |
(2) |
Fig. 1 Reaction scheme describing the different reaction routes present when treating C6 sugar monosaccharides in the presence of [K]Sn-USY in methanol medium. |
The reaction system comprises the Lewis-acid catalyzed isomerization between the three hexoses (glucose, mannose and fructose) together with their equilibria with methyl glycosides (Brønsted acid catalyzed), which act as substrate-masking agents during lactate synthesis.21 Apart from the equilibrium transformations between sugars and their corresponding methyl glycosides, the following irreversible transformations take place: Retro aldol condensation reactions (Lewis acid driven) occur from all the free sugars. In the case of aldoses (mannose and glucose), this splitting occurs in an asymmetric way, conducting a C2-fragment (glycolaldehyde) and a C4-sugar (tetrose). The former, undergoes a double acetalization in the methanol medium, finally providing glycolaldehyde dimethyl acetal (GADMA), whereas the latter evolves through a cascade of transformations to provide two different products through parallel transformations: methyl vinyl glycolate (MVG) and 4-methoxy, 2-hydroxy methyl butanoate (MMHB). Regarding the retro-aldol cleavage of fructose, the symmetric splitting of the molecule provides 2 C3 carbon backbone moieties: dihydroxyacetone (DHA) and glyceraldehyde (GLA), which are in equilibrium. GLA undergoes a cascade of transformations involving dehydration, tautomerization and 1,2-H shift processes that lead to MLa. Finally, a last reaction route is also present, in which fructose evolves through a Brønsted-acid catalyzed pathway to yield 5-hydroxymethyl furfural (HMF), its methyl ether (MMF), and the result from the hydrolytic transformation to methyl levulinate (MLe).
Fig. 2 Simplified reaction scheme with its A) reactions and B) their kinetic and equilibrium constants in the transformation of hexoses with [K]Sn-USY in methanol. |
The system proposed simplifies the whole mechanism to match the conclusions arising from the preliminary analysis of the experimental results. The epimerization of aldoses has been described to be more difficult than their isomerization into ketoses but occurs under the same conditions as sugar isomerization. Thus, epimerization of aldoses has been considered within the kinetic model. Regarding the complex series of transformations conducted to the different α-hydroxy esters (MLa, MVG, MMHB), most of the reaction intermediates have not been detected in quantifiable amounts, probably as a consequence of a low lifetime due to their large reactivity under the tested reaction conditions. In this way, these series of transformations have been lumped into a single transformation, with a kinetic rate corresponding to one of the limiting stages. This has been applied to GADMA and MMHB, which both are assumed to be directly produced from glucose. The equilibrium between glucose and mannose favors the formation of glucose so these products are only assumed to be generated from glucose. The formation of MVG is too low as compared to that of MMHB, so it is assumed negligible, and it was removed from the reaction scheme shown in Fig. 2. In a similar way, the cascade of transformations conducting to methyl lactate is lumped into a series of two reactions, one conducting to dihydroxyacetone from fructose, and a second one from DHA to methyl lactate. Finally, Brønsted acid-promoted reactions such as the dehydration of fructose into 5-hydroxymethyl furfural, and its subsequent hydrolytic transformation into methyl levulinate and methoxy methyl furfural are not included since the production rate is very small.
Taking into account these assumptions, a pseudo-homogeneous first-order kinetic model has been proposed with the equations summarized in Table 1. Eqn (3) to (5) reproduce the decomposition and formation of sugars. In the isomerization equilibria the equilibrium constants, defined as given in Table 2, must meet KGFKGMKFM = 1. Thus, KGM can be substituted in the eqn (3) and (5) by KGM = 1/(KFGKMF). The resultant equilibrium constants, KFG and KMF, are determined fitting the van't Hoff equation, eqn (13), to the reported thermodynamic data,25,26 see more details in the ESI.† Eqn (4), (6) and (8) define the formation and decomposition of methyl glycosides. The formation of these methyl glycosides is defined by the kinetic parameter of the reactions and the equilibrium with their corresponding sugar. For these reactions, no thermodynamic data is available so the equilibrium constants are determined by fitting them using the experimental data. The definitions of all the thermodynamic constants with the kinetic ones are provided in Table 2. Even though these equilibrium constants are mostly driven by the reaction's thermodynamics, it is important to note that they are also affected by the adsorption of each component on the zeolite.
(13) |
Apart from the equilibria, the formation of MLa follows a chain reaction mechanism with DHA as an intermediate as described in eqn (9) and (10). Two side products are included in the model: GADMA, which is generated from fructose as shown in eqn (11); and MMHB, which is generated as presented in eqn (12). The kinetic constants obtained for MLa, DHA, GADMA and MMHB include the effect of mass transfer, but this is negligible since the zeolites have particle sizes in the order of micrometers and the reactor is magnetically agitated at 500 rpm.
Kinetic parameters are related to the Arrhenius equation, eqn (20). The kinetic model was implemented in Python utilizing the optimization framework Pyomo. In the model, the differential equations are implemented following an orthogonal collocation method with the Pyomo.dae library. The objective function was defined as the minimization of the sum of residual squares (SSE), see eqn (21). The non-linear programming problem that defines the model was solved with Ipopt v.3.11 as the solver.
(20) |
(21) |
Fig. 3 A) XRD patterns and B) argon adsorption–desorption isotherms, recorded for the [K]Sn-USY catalyst. |
The total metal loading found in the functionalized zeolite by ICP-OES analysis accounted for 1.91 wt%. This result evidences the high efficiency of the metalation step, since almost quantitative incorporation of the tin added in this synthesis stage was achieved, as the theoretical tin loading to be achieved was 2.00 wt% Sn in the final material. This high efficiency is ascribed to a combination of the efficient dealumination step, which creates plenty of framework vacancies, and the high effectiveness of the triethylamine-driven metal grafting. Contacting the tin-functionalized zeolite with a KCl aqueous solution led to the incorporation of 0.4 wt% of potassium into the catalyst. This process exerts an important effect on the reduction of the total acidity of the catalyst, especially of strong Brønsted acid sites with ion exchange capacity. The result is a catalyst with a reduced Brønsted acid capacity, but with noticeable Lewis acidity, coming from the attached Sn species. However, the obtained material still presents a B/L = 9. This indicates that, despite the selective reduction of Brønsted acidity, mainly attributed to the remaining acid protons acting as tetrahedral aluminum compensating ions, which are not fully neutralized. The main reason seems that the remaining Brønsted acidity can be ascribed to the presence of unoccupied silanol nests, with weaker Brønsted acidity. Nevertheless, the obtained catalyst presents a Lewis-enhanced acidity and an expected improved performance on hexose transformation.
The environment of the zeolite-attached metal sites was evaluated through 119Sn solid-state MAS NMR and DR UV-vis analysis, the results are depicted in Fig. 4. 119Sn NMR spectra recorded for Sn-USY and its K-exchanged counterpart are depicted in Fig. 4A. These spectra feature main signals located at −660 ppm in the case of Sn-USY and at −670 ppm for [K]Sn-USY. Both signals are located in the region commonly ascribed to the presence of hydrated tin sites in tetrahedral conformation. Besides, the lack of signals around −600 ppm points to the absence of large SnO2 domains. The DR UV-vis spectra recorded for all the prepared materials, depicted in Fig. 4B, support the above conclusion on the absence of large domains of tin oxide. DR UV-vis spectrum recorded for the [K]Sn-USY catalyst presents a maximum absorption signal at 225 nm, together with a shoulder at 270 nm, both signals being conventionally ascribed to the presence of penta- and hexa-coordinated tin species. This result might be considered evidence for the presence of tin oxide species, but the removal of these signals when drying the sample under a vacuum at 300 °C points to an opposite direction. The spectra recorded for the dried sample are featured by a predominant signal located at 190 nm, corresponding to the ligand-to-metal charge transfer in tin species with tetrahedral coordination incorporated into the zeolite framework.27 In this way, the signals detected in the [K]Sn-USY exposed to the atmosphere humidity can be easily ascribed to hydrated tin sites, in which water acts as ligands, expanding the coordination sphere of the metal centers.28
From the previous results it can be concluded that the prepared catalysts display the right combination of textural and structural physicochemical properties together with the right functionalization degree with tin species in tetrahedral coordination, which provide Lewis acidity to the final zeolitic catalysts, and thus, these are expected to show a high catalytic performance in the transformation of C6 sugars into methyl lactate, as reported elsewhere.29
Reaction profiles suggest that both the isomerization of the different substrates and the etherification of free sugars to the corresponding methyl glycosides proceed at fast reaction rates, as it can be inferred from the high yields obtained towards these products during the early beginning of the reaction tests. Obviously, the temperature conditions strongly influenced the formation of these initial sugar derivatives, but at low-temperature conditions, the formation of methyl fructoside becomes the dominant reaction pathway, regardless of the starting substrate. As a consequence, maximum yields of MF (ranging from 60–80% of the starting free sugar) were produced after 120 min operating at 120 °C. Lower temperature conditions (100 °C) also produced methyl fructoside as the main product, though at slower reaction rates. Surprisingly, the formation of other glycosides (methyl glucoside, methyl mannoside) was quite reduced, indicating that the formation of fructose, and its etherification with the alcohol solvent is prevalent. This result suggests that the isomerization of the starting monosaccharide occurred at a substantially faster rate than glycosidation reactions, as it can be easily estimated from the abrupt consumption of all the sugar substrates during the first stages of the reaction tests. Among the glycosidation reactions, the higher productivity of methyl fructoside as compared to mannoside and glucoside suggests that the formation of the former is dominant. This different reactivity might be caused by differences in the ring closure step of the Fischer–Helferich glycosidation mechanism between aldoses and ketoses, which could be easier for the latter, thus enhancing their reactivity in this transformation.30 Finally, regarding the transformation of aldoses, the low productivity of glucose and mannose, or their methyl glycosides, when starting from mannose and glucose, respectively, suggest that epimerization reactions occur with very low extension, at least under the tested reaction conditions.
When considering higher temperature conditions, a substantial increase in the extension of parallel transformations to those previously described is observed. These other reactions include both the hydrolytic transformation of sugars to provide 5-HMF and its derivatives and the retro-aldol cleavage of sugars. Under these conditions, the production of methyl lactate becomes more dominant insofar as the reaction temperature increases, being produced in a higher extension at 160 °C and long reaction times (2 h) and reaching product yields in the range of 55–70%, depending on the starting substrate. Together with MLa, other products are also obtained, such as GADMA, MVG and MMHB, though in a limited extension – up to a yield of 2.3% for GADMA, 9.4% for MVG and 4.8% for MMHV, when treating fructose at 160 °C. On the other hand, hydrolytic pathways conducting to HMF, MMF and MLE, provide yields below 10% of the starting sugar in every case, because of the depressed Brønsted activity of [K]Sn-USY catalysts through ion exchange prevents these reactions.22 No great differences can be found in terms of product distributions when using the different hexose substrates, though in the case of fructose, faster production of MLa, even at lower temperature conditions was observed. This is an expected result, as the cleavage of fructose into trioses, finally yielding MLa, does not require the isomerisation stage, as it occurs for glucose or mannose.
The evaluation of the influence of the catalyst loading was performed by studying this parameter at three different levels. For this purpose, the reaction temperature was fixed at 140 °C, and the catalyst loading was varied between 0.2 and 0.8 g (substrate to catalyst mass ratio 2.5–9.5 wt/wt). The results from these experiments using glucose, mannose and fructose as starting sugar substrates are depicted in Fig. 6.
The main products obtained in the presence of different catalyst loadings were all the same as already described. The results evidenced the expected trend, that the higher the catalyst loadings the faster the transformation of sugars and the production of final compounds such as MLa. This can be easily inferred from the fast consumption of the starting monosaccharides during the early stages of the tests and the increasing MLa yields, and related products, obtained at longer reaction times as the catalyst loading increased. From the kinetic profiles, it is observed that the addition of increasing quantities of the catalyst produced no great differences in the product distributions, further than the transformation of the starting sugars and intermediates progressed towards the formation of methyl lactate and similar products. Besides, the relationship between the catalyst loading and the extension of the different chemical transformations seems to be linear, prompting the absence of external mass transfer control of the transformation. Finally, the slight differences in terms of reactivity between the different substrates were the same as already described when evaluating the influence of temperature. Here again, fructose was converted in a faster way as compared to glucose and mannose, whereas this last (mannose) evidenced a less reactive behavior.
Parameter | Pre-exponential parameter, ko (L kgcat−1 min−1) | Activation energy, Ea (kJ mol−1) |
---|---|---|
k 1 | 7511 ± 6.22 × 108 | 20.285 ± 54.212 |
k 4 | 2.7 × 10−12 ± 26651 | 38.963 ± 52.878 |
k 6 | 8884 ± 7.01 × 104 | 22.020 ± 20.02 |
k 7 | 2.73 × 10−9 ± 28.08 | 7.681 ± 10.08 |
k 9 | 2.582 × 10−16 ± 1.7 × 105 | 2.726 ± 45.03 |
k 11 | 264.232 ± 4.1 × 106 | 27.360 ± 35.3 |
k DHA | 3.382 × 105 ± 2.8 × 109 | 78.184 ± 7.194 |
k MLa | 5.289 × 108 ± 3.1 × 1012 | 103.877 ± 5.801 |
k GADMA | 15.822 ± 6.7 × 106 | 40.208 ± 44.58 |
k MMHB | 7.761 × 1011 ± 2.18 × 1019 | 142.031 ± 189.98 |
Equilibrium constant, Keq | ΔS (J mol−1 K−1) | ΔH (kJ K−1) |
---|---|---|
K F–MF | −297.63 | −118.04 |
K G–MG | 52.38 | 26.169 |
K M–MM | 306.06 | 111.28 |
This can also be observed in the evolution of the concentration of different substrates (mannose, fructose, and glucose) reported in Fig. 7 to 9. From the equilibrium between the sugars, it is observed that the degradation of mannose is promoted, which is in line with the observation of Tosi et al.21 In the equilibria of the sugars with the methyl glycosides, it can be seen that the isolation affects them with an increase of the concentration followed by a later drop. The comparison of the experimental values with the values reported for the equilibrium constants shows the expected behavior for MM and MG, but for MF the value obtained for KF–MF is around 1, and the concentration of MF is between 3 and 5 times higher than the one of fructose. This lower value obtained differs from the experimental observation since the thermodynamic equilibrium constant introduced for the sugar's equilibria has been modeled as in other previous works and it is only based on reaction thermodynamics. Other phenomena such as the adsorption of the sugars and methyl glycosides on the catalyst modify the equilibrium and the partition coefficients related to that adsorption are included within the obtained equilibrium constant. Future studies focusing on interphase characterization can provide further knowledge to decouple such phenomena. Similarly, future work can be carried out comparing other alcohol moieties through molecular design and molecular simulation could provide an initial guess of the entropy value to be used in the van't Hoff equation for the methyl glycosides. The value obtained for the entropy through correlations is higher than expected. On the contrary, the value for enthalpy is one order below that obtained for the reactions in the glycosides equilibria, something that can be expected.
Apart from the comparison of the equilibrium constants, the kinetic parameters were obtained. The activation energies estimated from the model for the methyl glycosides are around one order below those involved in the reactions corresponding to the equilibria between the sugars.
These great values are due to the high mixing rates of the reactor. At larger scales where mixing is not that efficient, the kinetic rates in the formation of methyl glycosides are expected to drop due to the mass transfer limitations. Apart from these limitations in the application of the model, it is also worth noting the errors are involved in the estimation of the kinetic parameters. Only three sets of experiments, each of them also with different initial substrates, have been employed in the estimation of the kinetic parameters, which implies the involvement of low information and subsequently large errors. Furthermore, the model proposed includes more terms than other approaches, which results in some parameters being cross-linked. As a result, there are some case studies where the activation energy obtained is slightly higher than double the one obtained by fitting all the experimental results together, which is the one reported as the average value.
The analysis of the kinetics for the two side products, MMHB and GADMA, show in both cases higher activation energies than the ones obtained for the methyl glycosides and the sugars equilibria. From both compounds, the one with the smaller activation energy is GADMA, but this energy is still double the one for the formation of mannose and fructose. This remarks that the preferred reaction paths for glucose are the transformation to fructose and methyl glucoside. Furthermore, at low temperatures, the formation of these reactions is preferred as it is shown in the concentration profiles reported in the ESI.† Thus, at these low temperatures, it is important to have methanol as an isolating agent since it can significantly help more on avoiding the formation of undesired compounds such as MMHB and DAGMA. In the case of glucose being employed as a substrate, it is partially stored in the form of MG and then it is released back to glucose for its transformation into fructose; see the profiles of MG in Fig. 6 to 8. Despite the prediction of the model, it does not ensure a good accuracy for these compounds at 140 °C and 160 °C, it can be observed that at low temperatures, 100 °C and 120 °C, the model perfectly fits the concentration, see figures reported in the ESI.† At low temperatures, the concentration of these compounds is higher, so it finally has a higher weight on the minimization of the square errors. Furthermore, it is also important to note that in the model, all retro-aldol reactions are included as a single one.
The kinetic parameters of MLa show that its formation is faster than the other two side products of the reaction system, GADMA and MMHB. The estimated activation energies are higher than those for the methyl glycosides and sugars equilibria, but the pre-exponential parameter is several orders higher. This indicates that the formation of MLa is more sensitive to the temperature and can be promoted by increasing it. The estimated kinetic parameter for MLa is also around 5 times higher than the one for the formation of the two other side products, GADMA and MMHB.
However, its generation is mostly promoted thanks to the formation of methyl glycosides as the isolating agents. This enhancement makes the system achieve a yield in the range of 0.25–0.3 mol L−1, see results at 160 °C in the ESI.† This maximum yield is lower when mannose is the original substrate, ∼0.2 mol L−1, than for fructose and glucose. The transformation of mannose derives from other side products (MMHB and GADMA) and favors the accumulation of fructose as MF. Apart from the analysis of the parameters in MLa production, the comparison of the model with the experimental data shows that the two steps defined in the model are not always observed in the experimental data. At low temperatures, the experimental data in the formation of MLa follows the profile of a first-order mechanism, but the model is still in the first part of the profile, see ESI.† This leads to an error in predicting the concentration of MLa at low temperatures, but at high temperatures, the error is reduced. However, it is important to note the concentration achieved in each case. At low temperatures, the yield of MLa is below 5 mol L−1; meanwhile, at high temperatures, the yield is around 35 mol L−1. Thus, the cases at high temperatures have more impact on the minimization of the SSEs.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4re00307a |
This journal is © The Royal Society of Chemistry 2025 |