Hydrolysis–Dehydration of Cellulose: Efﬁciency of NbZr Catalysts under Batch and Flow Conditions

: Niobium oxide supported on ZrO 2 and mixed oxide of NbO x -ZrO 2 was prepared and characterized. Mechanical treatment was followed by the microwave heating procedure of catalysts with more advanced textural parameters. The amount of Lewis (LAS) and Brønsted (BAS) acid sites rose with the increasing Nb content in the catalysts. The catalytic properties of samples of niobia-zirconia (NbZr samples, NbZr catalysts) were studied in a cellulose hydrolysis–dehydration reaction at 453 K under an inert Ar atmosphere in a batch reactor. Glucose and 5-hydroxumethylfurfural (5-HMF) were the major products. The initial reaction rate could be tuned by the density of acid sites on the surface of solid. At a low density of acid sites (0.1–0.3 µ mol · m − 2 ), the initial reaction rate had a pronounced inverse correlation. Increasing the LAS/BAS from 0.3 to 2.5 slightly stimulated the formation of the target products. The catalytic properties of NbZr catalysts prepared by microwave treatment were studied in cellulose transformation in a ﬂow set-up. Glucose was found to be the major product. The maximum yield of glucose was observed in the presence of the sample of 17%Nb/ZrO 2 . Increasing Nb content resulted in the formation of Nb-associated acid centers and, in turn, increasing catalyst acidity and activity.


Introduction
A promising alternative to limited fossils is a plant biomass which could provide a sustainable development based on renewable biobased raw materials [1,2].Cellulose, being the most abundant inedible plant source, is a perspective renewable for the production of energy and chemicals [3,4].The main product of cellulose hydrolysis is glucose, which is a substrate for the production of valuable chemicals such as formic acid [5][6][7][8], 5-hydroxymathylfurfural (5-HMF) [9], and sorbitol [10].Notably, 5-HMF seems to be a perspective sleeping-platform bioorganic molecule and could be applied as a substrate for the synthesis of chemicals, biofuels, pharmaceuticals, plastics, etc. [11].Thus, 2,5dimethylfuran and 2-methylfuran, which are promising biofuels and octane booster additives, could be synthesized starting with 5-HMF.Compounds of 2,5-furandicarboxylic and 5-hydroxymethylfuranoic acids are also furan derivatives of 5-HMF with high potential in the fuel and polymer industries [12].Moreover, 5-hydroxymethylfurfural is a precursor of non-furan derivatives, namely levulinic and adipic acids, 1,6-hexanediol, caprolactam, caprolactone, etc. [9].
It is known that the synthesis of 5-HMF from cellulose proceeds through three successive stages: (1) the hydrolysis of cellulose to glucose, (2) the isomerization of glucose to fructose, and then (3) the dehydration of fructose to 5-HMF.All the stages can be produced It is known that the synthesis of 5-HMF from cellulose proceeds through three successive stages: (1) the hydrolysis of cellulose to glucose, (2) the isomerization of glucose to fructose, and then (3) the dehydration of fructose to 5-HMF.All the stages can be produced in a one vessel (one-pot process) without the isolation and purification of the reaction intermediates, glucose and fructose (Scheme 1) [9,13,14].Scheme 1. One-pot synthesis of 5-HMF from cellulose.
However, direct processing the polysaccharide to 5-HMF is quite challenging as cellulose is insoluble in water and has a crystal structure.According to the literature, solid catalytic systems of heterogeneous catalysis present a perspective approach for the transformation of bioresources to valuable chemicals [1].Heteropoly compounds [15], carbon materials [16], polymer resins, noble metal nanoparticles, and zeolites [17] can be used for the processing cellulose into 5-HMF.Water-insoluble oxides (ZrO2, TiO2, SiO2, Nb2O5, etc.) are of particular interest for cellulose processing [18].The high potential of zirconia catalysts in hydrolysis of cellulose was reported previously [19,20].Recently, Yang et al. reported a 30% yield of reducing sugars of cellulose over zirconium oxide impregnated by SO4 2− and supported on montmorillonite [21].Catalysts based on niobium oxide can effectively dehydrate monosaccharaides to 5-HMF.Thus, Ngee et al. reported an 88% yield of 5-HMF from fructose over sulfated niobium oxide at 393 K in the presence of a DMSO cosolvent [22].
Materials based on niobium and zirconium oxides are of considerable interest in the processing of cellulose into 5-HMF due to the existing acid-base sites [19] and their high stability in hot water [23][24][25][26].We previously investigated [27,28]one-pot hydrolysis-dehydration of mechanically activated cellulose in pure hot water to produce glucose and 5hydroxymethylfurfural in the presence of amorphous niobium oxide supported on a zirconium dioxide (NbOx/ZrO2).It was found that activity of NbOx/ZrO2 systems depended on the Nb2O5 content and decreased as follows: 2.8%Nb/ZrO2~1.1%Nb/ZrO2≫ ZrO2 > 0.5%Nb/ZrO2.The authors explained the trend by the changing aggregation state of niobium oxide species that led to the variation of the acid-base properties of the Nb/ZrO2 composite.The highest glucose and 5-HMF yields (22 and 16 mol%, respectively) were obtained over 2.8%Nb/ZrO2.
The nature of acidic sites is an important parameter for one-pot cellulose transformation processes.As reported previously, the isomerization of glucose to fructose proceeds on Lewis acid sites (LAS), while Brønsted acid sites (BAS) are responsible for hydrolysis and dehydration [26,29].According to Nijhuis et al. [30], a maximum 55% yield of 5-HMF in the dehydration of glucose over titanium, zirconium, and niobium phosphate catalysts was observed at a BAS/LAS ratio of about 1.The catalytic properties of Nb/WO3 also can be tuned by the variation of the BAS/LAS ratio.According to Yue et al. [24], a maximum 38% yield of 5-HMF was observed at a BAS/LAS ratio of 1.5.Thus, the design of catalytic materials with a high water tolerance and a suitable ratio of BAS/LAS is a challenge.Note that to our knowledge, an influence of the BAS/LAS ratio of NbOx/ZrO2 catalysts in cellulose hydrolysis-dehydration has not been reported previously.
A type of catalytic reactor can also affect the one-pot transformation of polysaccharides.Cellulose transformation is usually realized in batch reactors (autoclaves) in pure water.However, a flow reactor also can be used.Thus, fructose, sucrose, and inulin were dehydrated in the presence of phosphate-containing niobium catalysts under both batch and flow reactors [31].It was demonstrated that a 51.2% conversion of fructose with a 59.4% selectivity of 5-HMF were achieved in a batch reactor at 373 K after 3 h, while a However, direct processing the polysaccharide to 5-HMF is quite challenging as cellulose is insoluble in water and has a crystal structure.According to the literature, solid catalytic systems of heterogeneous catalysis present a perspective approach for the transformation of bioresources to valuable chemicals [1].Heteropoly compounds [15], carbon materials [16], polymer resins, noble metal nanoparticles, and zeolites [17] can be used for the processing cellulose into 5-HMF.Water-insoluble oxides (ZrO 2 , TiO 2 , SiO 2 , Nb 2 O 5 , etc.) are of particular interest for cellulose processing [18].The high potential of zirconia catalysts in hydrolysis of cellulose was reported previously [19,20].Recently, Yang et al. reported a 30% yield of reducing sugars of cellulose over zirconium oxide impregnated by SO 4 2− and supported on montmorillonite [21].Catalysts based on niobium oxide can effectively dehydrate monosaccharaides to 5-HMF.Thus, Ngee et al. reported an 88% yield of 5-HMF from fructose over sulfated niobium oxide at 393 K in the presence of a DMSO co-solvent [22].
Materials based on niobium and zirconium oxides are of considerable interest in the processing of cellulose into 5-HMF due to the existing acid-base sites [19] and their high stability in hot water [23][24][25][26].We previously investigated [27,28] one-pot hydrolysisdehydration of mechanically activated cellulose in pure hot water to produce glucose and 5-hydroxymethylfurfural in the presence of amorphous niobium oxide supported on a zirconium dioxide (NbO x /ZrO 2 ).It was found that activity of NbO x /ZrO 2 systems depended on the Nb 2 O 5 content and decreased as follows: 2.8%Nb/ZrO 2 ~1.1%Nb/ZrO 2 ZrO 2 > 0.5%Nb/ZrO 2 .The authors explained the trend by the changing aggregation state of niobium oxide species that led to the variation of the acid-base properties of the Nb/ZrO 2 composite.The highest glucose and 5-HMF yields (22 and 16 mol%, respectively) were obtained over 2.8%Nb/ZrO 2 .
The nature of acidic sites is an important parameter for one-pot cellulose transformation processes.As reported previously, the isomerization of glucose to fructose proceeds on Lewis acid sites (LAS), while Brønsted acid sites (BAS) are responsible for hydrolysis and dehydration [26,29].According to Nijhuis et al. [30], a maximum 55% yield of 5-HMF in the dehydration of glucose over titanium, zirconium, and niobium phosphate catalysts was observed at a BAS/LAS ratio of about 1.The catalytic properties of Nb/WO 3 also can be tuned by the variation of the BAS/LAS ratio.According to Yue et al. [24], a maximum 38% yield of 5-HMF was observed at a BAS/LAS ratio of 1.5.Thus, the design of catalytic materials with a high water tolerance and a suitable ratio of BAS/LAS is a challenge.Note that to our knowledge, an influence of the BAS/LAS ratio of NbO x /ZrO 2 catalysts in cellulose hydrolysis-dehydration has not been reported previously.
A type of catalytic reactor can also affect the one-pot transformation of polysaccharides.Cellulose transformation is usually realized in batch reactors (autoclaves) in pure water.However, a flow reactor also can be used.Thus, fructose, sucrose, and inulin were dehydrated in the presence of phosphate-containing niobium catalysts under both batch and flow reactors [31].It was demonstrated that a 51.2% conversion of fructose with a 59.4% selectivity of 5-HMF were achieved in a batch reactor at 373 K after 3 h, while a 45.5% conversion of fructose and a 58.0%selectivity of 5-HMF were achieved in a flow reactor at 358 K.The conversion of mono-and disaccharides to 5-HMF in the presence of catalysts based on niobium under flow conditions was investigated by Carniti et al. [26,[32][33][34].They demonstrated that the conversion of fructose and selectivity of 5-HMF were 68% and 34%, respectively, in the presence of niobium phosphate at 373 K [34].At the same, the transformation of cellobiose (disaccharide) at 403 K made it possible to achieve 5-HMF in 21% yield (31% selectivity, 70% substrate conversion) [26].To our knowledge, the conversion of cellulose in the presence of solid oxide NbO x and/or ZrO 2 catalytic systems in a flow reactor has not been carried out.
This work presents data on the evaluation of the influence of the BAS/LAS ratio and acid site density of Nb/ZrO 2 catalysts and on the possibility of using the catalysts under study to obtain glucose and 5-HMF under hydrothermal batch and flow conditions.We also aimed to reveal that an aggregation state of niobium oxide or acid properties of catalysts would be a key parameter which affects the catalyst activity.

Synthesis and Characterization of Catalysts
Catalysts of niobium and zirconium oxides (also referred to as NbZr catalysts, NbZr samples, NbZr materials) were prepared two ways (Table 1).Such approach allowed us to reveal the effects of the agglomeration of the niobia particles, surface acidity, and ratio of Brønsted acid sites (BAS) and Lewis acid sites (LAS) (BAS/LAS ratio) on the reaction rate and selectivity toward 5-HMF in the one-pot hydrolysis-dehydration of cellulose.

Structural properties.
First of all, we investigated structure of NbZr catalysts.X-ray diffraction patterns (XRD) of the samples are shown in Figure 1.It was found that the method of catalyst preparation affected the structure of ZrO 2 .In XRD patterns of NbZr-MAW samples, we revealed characteristic diffraction peaks of monoclinic phase of ZrO 2 at diffraction angles (2Θ) 24.4  (m-ZrO 2 , JCPDS cards no.37-1484) [35].Moreover, the four peaks at 30.2, 35.2, 50.4,and 60.0 • belonged to the tetragonal phase of ZrO 2 (t-ZrO 2 , JCPDS cards no.17-0923) [35].Analysis of the XRD patterns of NbZr-W pointed to the amorphous phase for both ZrO 2 and Nb 2 O 5 (Table 1).These results agree with the textural properties of NbZr materials.
Textural properties.The specific surface areas (S BET ) of the NbZr materials were calculated from the nitrogen adsorption isotherms using the Brunauer-Emmett-Teller (BET) method.The main results are shown in Table 1.As can be observed, the Nb content affected the textural properties of the NbZr-MAW materials.The increase in Nb content from 0 to 6 wt.% led to an increase in the specific surface area from 10 to 60 m 2 •g −1 .At the same, the mean pore size changed in the reverse order.Note that NbZr-MAW0 and NbZr-MAW1 were macroporous solids, while NbZr-MAW2 and NbZr-MAW3 were mesoporous solids.Textural properties.The specific surface areas (SBET) of the NbZr materials were calculated from the nitrogen adsorption isotherms using the Brunauer-Emmett-Teller (BET) method.The main results are shown in Table 1.As can be observed, the Nb content affected the textural properties of the NbZr-МАW materials.The increase in Nb content from 0 to 6 wt.% led to an increase in the specific surface area from 10 to 60 m 2 •g −1 .At the same, the mean pore size changed in the reverse order.Note that NbZr-МАW0 and NbZr-МАW1 were macroporous solids, while NbZr-МАW2 and NbZr-МАW3 were mesoporous solids.
The specific surface area of the NbZr-W materials was higher in comparison with that of the NbZr-MAW samples; this can be explained by the difference in phase state (Table 1).The specific surface area also depends on the Nb content.The specific surface area rose with the increase in Nb loading.Moreover, all NbZr-W samples were mesoporous solids.The total pore volume does not change substantially when the Nb content increased from 0 to 4 wt.%, and the total pore content increased 1.8 times at 17 wt.% of Nb content in the sample.
It can be seen that synthesis via mechanical treatment and microwave heating makes it possible to obtain macro-/mesoporous materials.At the same time, mesoporous materials with high specific surface area can be obtained by the method of the joint microwave heating of niobia and zirconia precursors.Previously, we reported textural data on pure ZrO2 samples [36].It was found that oxide catalysts prepared by microwave heating, similar to the present work, had a higher specific surface area compared to catalysts produced by a combination of mechanical activation and microwave.
Electronic properties.The effect of method preparation and Nb content on the electronic properties of NbZr samples was investigated by Ultraviolet-visible diffuse reflectance spectroscopy (DR-UV-Vis).DR-UV-Vis spectra of NbZr samples are shown in Figures S1 and S2 (Supplementary Materials).The broad absorption in the 250-400 nm region associated with O 2− → Nb 5+ charge transfer transition was observed in the spectrum of Nb2O5 [37].The DR-UV-Vis spectra of all NbZr samples also contained this band.This band shifted to lower energy with the increase in the niobium content, indicating that Nb 5+ and Zr 4+ ions were electronically dependent.
The optical band edge can be discussed in the Tauc's equation related to the determination of band energy (Eg) (Equation (1)): where α is the absorbance, h is photon energy, B is an independent parameter of the photon's energy for the respective transitions, and Eg is the band gap.The plot of (αh)  The specific surface area of the NbZr-W materials was higher in comparison with that of the NbZr-MAW samples; this can be explained by the difference in phase state (Table 1).The specific surface area also depends on the Nb content.The specific surface area rose with the increase in Nb loading.Moreover, all NbZr-W samples were mesoporous solids.The total pore volume does not change substantially when the Nb content increased from 0 to 4 wt.%, and the total pore content increased 1.8 times at 17 wt.% of Nb content in the sample.
It can be seen that synthesis via mechanical treatment and microwave heating makes it possible to obtain macro-/mesoporous materials.At the same time, mesoporous materials with high specific surface area can be obtained by the method of the joint microwave heating of niobia and zirconia precursors.Previously, we reported textural data on pure ZrO 2 samples [36].It was found that oxide catalysts prepared by microwave heating, similar to the present work, had a higher specific surface area compared to catalysts produced by a combination of mechanical activation and microwave.
Electronic properties.The effect of method preparation and Nb content on the electronic properties of NbZr samples was investigated by Ultraviolet-visible diffuse reflectance spectroscopy (DR-UV-Vis).DR-UV-Vis spectra of NbZr samples are shown in Figures S1 and S2 (Supplementary Materials).The broad absorption in the 250-400 nm region associated with O 2− → Nb 5+ charge transfer transition was observed in the spectrum of Nb 2 O 5 [37].The DR-UV-Vis spectra of all NbZr samples also contained this band.This band shifted to lower energy with the increase in the niobium content, indicating that Nb 5+ and Zr 4+ ions were electronically dependent.
The optical band edge can be discussed in the Tauc's equation related to the determination of band energy (E g ) (Equation ( 1)): where α is the absorbance, hν is photon energy, B is an independent parameter of the photon's energy for the respective transitions, and E g is the band gap.The plot of (αhν) 1/2 versus hν are displayed in Figures S1 and S2 (Supplementary Materials) for NbZr samples.The E g values were determined from the steep absorption edge.The findings show that the E g value for the studied NbZr samples depended on the Zr/Nb ratio and decreased with the increasing Nb content (Figure 2).A similar trend was demonstrated for a ZrO 2 /Nb 2 O 5 composite [38].E g values of the samples with ZrO 2 :Nb 2 O 5 mole ratios from 150:1 to 4:1 (from 1.3 to 33.5 mol % of Nb) decreased linearly with the Nb content [38].This phenomenon was explained by the formation of the niobium oxide phase and changing oligomeric state of its particles.
with the increasing Nb content (Figure 2).A similar trend was demonstrated for ZrO2/Nb2O5 composite [38].Eg values of the samples with ZrO2:Nb2O5 mole ratios fro 150:1 to 4:1 (from 1.3 to 33.5 mol % of Nb) decreased linearly with the Nb content [38 This phenomenon was explained by the formation of the niobium oxide phase and chan ing oligomeric state of its particles.The Eg values of 4.45 eV for niobium oxalate [39] and Nb-MCM-41 samples [40] co responded to the predominant formation of isolated NbO4 structures on the surface.values equal to 3.9 eV were assigned to polymeric NbOx structures, for example, to niob acid Nb2O5 nH2O [39].According to that, NbOx polymer structures on the surface of NbZ MAW catalysts were formed in all synthesized samples containing 0.6, 1.2, and 6 wt.% niobium (Figure 2).The catalysts of mixed oxides NbZr-W with a niobium content of and 4 wt.% had isolated NbO4 on the surface.The NbOx polymer structures (Eg = 4.0 eV were formed only on the surface of the NbZr-W3 sample which contained 17 wt % of N Although increasing the red shift corresponded to rising the degree of polymerization the NbOx structures, the formation of bulk niobium oxide Nb2O5 (Eg = 3.42-3.2eV [41 was not observed for all the catalysts. Nature of acid sites.The nature of active sites was investigated by IR spectroscop IR spectra of NbZr-MAW2 and NbZr-MAW3 samples preliminarily dehydroxylated 673 K are shown in Figure S3 (Supplementary Materials).The spectrum of the NbZ MAW2 sample contained a broad band with a maximum of 3600 cm −1 and, on its shoulde a weak band around 3660 cm −1 , which can be assigned to the terminal and bridging -O groups, respectively.At the same time, the intensity of the band at 3660 cm −1 was high in comparison with that of the band at 3600 cm −1 in the spectrum of NbZr-MAW3.We ca assume that this difference arose from the higher Nb content in NbZr-MAW3 (6 wt.% than in NbZr-MAW2 (1.2 wt.%) that may have favored the appearance of hydroge bonded -OH groups.
The number and strength of acid sites in NbZr samples was analyzed by IR spectro copy using pyridine as the probe molecule Figure S4 (Supplementary Materials).The ma results are presented in Table 2. From, these data, the NbZr-W and NbZr-MAW sampl contained both LAS and BAS on the surface.The number of acid sites also depended o the Nb content.Thus, the amount of LAS and BAS rose with the increasing Nb content the NbZr-МАW samples prepared via mechanical treatment/microwave heating.Th change of the molar ratio of LAS/BAS had a similar trend.However, for NbZr-W sampl prepared via the joint microwave heating of niobia and zirconia precursors, the quanti The E g values of 4.45 eV for niobium oxalate [39] and Nb-MCM-41 samples [40] corresponded to the predominant formation of isolated NbO 4 structures on the surface.E g values equal to 3.9 eV were assigned to polymeric NbO x structures, for example, to niobic acid Nb 2 O 5 nH 2 O [39].According to that, NbO x polymer structures on the surface of NbZr-MAW catalysts were formed in all synthesized samples containing 0.6, 1.2, and 6 wt.% of niobium (Figure 2).The catalysts of mixed oxides NbZr-W with a niobium content of 2 and 4 wt.% had isolated NbO 4 on the surface.The NbO x polymer structures (E g = 4.0 eV) were formed only on the surface of the NbZr-W3 sample which contained 17 wt % of Nb.Although increasing the red shift corresponded to rising the degree of polymerization of the NbO x structures, the formation of bulk niobium oxide Nb 2 O 5 (E g = 3.42-3.2eV [41]) was not observed for all the catalysts.
Nature of acid sites.The nature of active sites was investigated by IR spectroscopy.IR spectra of NbZr-MAW2 and NbZr-MAW3 samples preliminarily dehydroxylated at 673 K are shown in Figure S3 (Supplementary Materials).The spectrum of the NbZr-MAW2 sample contained a broad band with a maximum of 3600 cm −1 and, on its shoulder, a weak band around 3660 cm −1 , which can be assigned to the terminal and bridging -OH groups, respectively.At the same time, the intensity of the band at 3660 cm −1 was higher in comparison with that of the band at 3600 cm −1 in the spectrum of NbZr-MAW3.We can assume that this difference arose from the higher Nb content in NbZr-MAW3 (6 wt.%) than in NbZr-MAW2 (1.2 wt.%) that may have favored the appearance of hydrogen-bonded -OH groups.
The number and strength of acid sites in NbZr samples was analyzed by IR spectroscopy using pyridine as the probe molecule Figure S4 (Supplementary Materials).The main results are presented in Table 2. From, these data, the NbZr-W and NbZr-MAW samples contained both LAS and BAS on the surface.The number of acid sites also depended on the Nb content.Thus, the amount of LAS and BAS rose with the increasing Nb content in the NbZr-MAW samples prepared via mechanical treatment/microwave heating.The change of the molar ratio of LAS/BAS had a similar trend.However, for NbZr-W samples prepared via the joint microwave heating of niobia and zirconia precursors, the quantity of niobium had a small effect on the number of LAS and BAS and LAS/BAS molar ratio (Table 2).
The method of preparation of the NbZr samples affected the strength of the acid sites.The strength of the acid sites of the NbZr-W samples was higher (pH susp 2.1-2.4)than that of NbZr-MAW (pH susp 5.7-8.2) (Table 3).We suggest that this was due to the different ordering of the materials, as indicated by the XRD data. 1 Experimental conditions: 10 g L −1 of cellulose in 45 mL of reaction mixture, 1 g L −1 of catalyst for NbZr-W samples and 0.5 g L −1 catalyst for NbZr-MAW samples, 1 mPa of Ar, 453 K, 7 h. 2 Initial rate of formation of target products (glucose + 5-HMF). 3Cellobiose, fructose, furfural, acetic, formic, malic, levulinic, lactic, glycolic, and tartaric acids were detected.

Catalytic Properties of NbZr Catalysts in a Batch Reactor
The catalytic properties of NbZr catalysts were investigated in a cellulose hydrolysisdehydration reaction at 10 g L −1 of cellulose in 45 mL of reaction mixture under an Ar atmosphere (1 mPa), 453 K, for 7 h.First of all, it was found that reaction in the presence of all studied NbZr catalysts was heterogeneous, following the hydrothermal stability tests of the NbZr catalysts.According to ICP-AES analysis, only 10 −4 -10 −3 wt.% of Nb and Zr were leached form solids.
Second, the analysis of reaction mixtures by HPLC shows that glucose and 5-HMF were the major products.Moreover, cellobiose, fructose, furfural, acetic, formic, malic, levulinic, lactic, glycolic, and tartaric acids were registered.The main results are shown in Table 4.
Effect of catalyst loading.The effect of catalyst content in the reaction mixture was investigated in the presence of NbZr-W3 and NbZr-MAW3.Correlations between the weight ratio of catalyst/cellulose and the yield of products are shown in Figure 3.One can see that the yield of glucose and 5-HMF depended on the cellulose/catalyst ratio.Thus, the total yield of glucose and 5-HMF did not exceed 40% of all products produced at an NbZr-W3/cellulose ratio of 1. Their total yield rose to 70% when the NbZr-W3/cellulose was increased to 5 (Figure 3A).Note that the glucose/5-HMF ratio increased with the decreasing amount of catalyst.
see that the yield of glucose and 5-HMF depended on the cellulose/catalyst ratio.Thus, the total yield of glucose and 5-HMF did not exceed 40% of all products produced at an NbZr-W3/cellulose ratio of 1. Their total yield rose to 70% when the NbZr-W3/cellulose was increased to 5 (Figure 3A).Note that the glucose/5-HMF ratio increased with the decreasing amount of catalyst.
In the presence of NbZr-МАW3, the yield of the major products similarly depended on the catalyst loading (Figure 3В).Thus, when the NbZr-МАW3/cellulose ratio was 1/1, the yields of glucose and 5-HMF were 2 and 10%, respectively.The yield of glucose and 5-HMF increased to 14 and 13%, respectively, when decreasing the catalyst loading to 1 g L −1 (NbZr-МАW3/cellulose ratio = 1/10).The changes in the total yield of major products and glucose/5-HMF ratio were caused by a decrease in the number of active centers available for the reaction.It should be noted that the subsequent transformations of the formed target reaction products to the side products and humins were catalyzed by the same active sites as the desired reactions.The optimum ratio of catalyst/cellulose made it possible to find a balance between the formation of glucose and 5-HMF and their degradation.Note that the efficiencies of NbZr-W3 and NbZr-MAW3 were different due to the different specific surface area and concentrations of active centers of NbZr-W3 and NbZr-MAW3 (Tables 1 and 2).Thus, the NbZr-W3 sample had an SBET of 154 m 2 g −1 and a total number of acid sites of 33.4 µmol g −1 , while the specific surface area of NbZr-МАW3 was 60 m 2 g −1 with an acid site concentration of 51,2 µmol g −1 (Table 2).According to the findings, the high yields of glucose and 5-HMF in the presence of NbZr-W and NbZr-MAW can be obtained at catalyst/cellulose weight ratio of 1/20 and 1/10, respectively.

Effect of acidity.
The acidity of catalysts is a key parameter which determines catalyst activity [26,29,42,43].A thorough analysis of the data showed no direct relationship between the catalytic activity of NbZr catalysts and the total amount of LAS and BAS (Figure S5, Supplementary Materials).Thus, we compared NbZr-МАW2 and NbZr-МАW3 (Table 1: In the presence of NbZr-MAW3, the yield of the major products similarly depended on the catalyst loading (Figure 3B).Thus, when the NbZr-MAW3/cellulose ratio was 1/1, the yields of glucose and 5-HMF were 2 and 10%, respectively.The yield of glucose and 5-HMF increased to 14 and 13%, respectively, when decreasing the catalyst loading to 1 g L −1 (NbZr-MAW3/cellulose ratio = 1/10).
The changes in the total yield of major products and glucose/5-HMF ratio were caused by a decrease in the number of active centers available for the reaction.It should be noted that the subsequent transformations of the formed target reaction products to the side products and humins were catalyzed by the same active sites as the desired reactions.The optimum ratio of catalyst/cellulose made it possible to find a balance between the formation of glucose and 5-HMF and their degradation.Note that the efficiencies of NbZr-W3 and NbZr-MAW3 were different due to the different specific surface area and concentrations of active centers of NbZr-W3 and NbZr-MAW3 (Tables 1 and 2).Thus, the NbZr-W3 sample had an S BET of 154 m 2 g −1 and a total number of acid sites of 33.4 µmol g −1 , while the specific surface area of NbZr-MAW3 was 60 m 2 g −1 with an acid site concentration of 51.2 µmol g −1 (Table 2).According to the findings, the high yields of glucose and 5-HMF in the presence of NbZr-W and NbZr-MAW can be obtained at catalyst/cellulose weight ratio of 1/20 and 1/10, respectively.
Effect of acidity.The acidity of catalysts is a key parameter which determines catalyst activity [26,29,42,43].A thorough analysis of the data showed no direct relationship between the catalytic activity of NbZr catalysts and the total amount of LAS and BAS (Figure S5, Supplementary Materials).Thus, we compared NbZr-MAW2 and NbZr-MAW3 (Table 4, lines 3 and 4).The catalysts had similar amounts of BAS and S BET (51-60 m 2 g −1 ) and differed from each other in the LAS content by two times (NbZr-MAW3 > NbZr-MAW2).In the presence of both catalysts, the glucose yield was equals to 14%, while the 5-HMF yield and their initial reaction rate) were close.A similar result was found for the NbZr-W series.
It was reported previously that BAS centers are responsible for the hydrolysis of cellulose to glucose and dehydration of sugars to 5-HMF.On the other hand, LAS plays an important role in the isomerization of glucose to fructose (Scheme 1) [44].Thus, it could be proposed that an appropriate combination of BAS and LAS could stimulate the reaction.However, during the investigation, the effect of the LAS/BAS ratio on the yields of glucose and 5-HMF was negligible.As can be seen from Figure 4A, the increase in the LAS/BAS ratio yield of major products did not increase significantly.However, according to [44], increasing of LAS ratio favors the formation of the target products.Since catalysts have different textural characteristics, it can be assumed that their activity should be determined by the density of acid sites on the surface of NbZr samples.The correlation between the initial reaction rate (R) and the density of acid sites (ρ AS ) is shown in Figure 4B.
It was reported previously that BAS centers are responsible for the hydrolysis of cellulose to glucose and dehydration of sugars to 5-HMF.On the other hand, LAS plays an important role in the isomerization of glucose to fructose (Scheme 1) [44].Thus, it could be proposed that an appropriate combination of BAS and LAS could stimulate the reaction.However, during the investigation, the effect of the LAS/BAS ratio on the yields of glucose and 5-HMF was negligible.As can be seen from Figure 4A, the increase in the LAS/BAS ratio yield of major products did not increase significantly.However, according to [44], increasing of LAS ratio favors the formation of the target products.Since catalysts have different textural characteristics, it can be assumed that their activity should be determined by the density of acid sites on the surface of NbZr samples.The correlation between the initial reaction rate (R) and the density of acid sites (AS) is shown in Figure 4B.It is clearly seen that the density of acid sites on the surface of the solid is an important parameter affecting the activity of NbZr samples.Two concentration areas can be distinguished: a low density of acid sites (0.1-0.3 µmol•m −2 ) (region I) and a high density of acid sites (>0.3 µmol•m −2 ) (region II).In the first region, the initial reaction rate decreased with the increasing density of acid sites, while in the second region, the effect of ρ AS was negligible.
Since the samples had a different amount of acid sites, we estimated the efficiency of the NbZr catalyst using Equation (2): where R-initial reaction rate (mmol L −1 h −1 ); m cat -amount of catalyst loaded (g); and [AS]-concentration of acid sites in a catalyst (mol g −1 ).The data shown in Figure 4C suggest that increasing the density of acid sites led to the decreasing activity of one active site.This trend can be the result of a decrease in the accessibility of the active site for reagents.Interestingly, the density of acid sites, not the aggregation of niobium, is a key parameter that significantly adjusted the catalyst activity.Thus, NbZr-W2 was more active compared to NbZr-W1 (Figure 4C), although both catalysts contained isolated Nb structures (Figure 2).

Catalytic Properties of NbZr Catalysts in a Flow Reactor
We were also interested in understanding the one-pot hydrolysis-dehydration of cellulose in the presence of NbZr catalysts in a flow set-up.This is because the majority of studies have focused on the studying the process in batch reactors.Thus, Sasaki et al. [45] investigated the transformation of microcrystalline cellulose under flow conditions at 563-673 K and 25 MPa.The authors reported 75% yields of different products (oligosaccharides, erythrose, dihydroxyacetone, fructose, glucose, glyceraldehyde).
Rogalinski et al. [46,47] studied cellulose transformation under both batch and flow modes and lower temperature of 483-583 K.The yields of glucose and its destruction products reached 5 and 10%, respectively, while the cellulose conversion was 25%.It was demonstrated that application of solid catalysts made it possible to decrease the reaction temperature.Overall, 75 and 10% yields of glucose and 5-HMF, respectively, were achieved when transforming microcrystalline cellulose at 473 K and 25 MPa pressure in the presence of Sibunit carbon materials treated to increase its acidity prior to the reaction [16].
Herein, we investigated the catalytic behavior of NbZr catalysts prepared by microwave treatment in a flow set-up.The series was chosen as it demonstrated high activity under batch conditions.According to the analysis of reaction mixtures by HPLC, glucose was the major product of the reaction.Moreover, reaction mixture contained water soluble oligosaccharides, cellobiose, glucose, fructose, 5-HMF, and furfural.
Effect of catalyst loading.The effect of catalyst loading on the reaction rate and yield of glucose was studied in the presence of NbZr-W2.The weight ratios of the NbZr-W2/cellulose were 1:1, 1:2.5, 1:5, and 1:10.Kinetic curves of glucose formation are shown in Figure 5A.These data reveal that the highest yield of glucose of 21% was observed when the catalyst/substrate ratio was 1/5.The decreasing of catalyst loading caused a simultaneous decrease in the glucose yield due to an insufficient number of active centers in the reactor.On the other hand, when the cellulose/catalyst ratio was 1/1(2.5), a decrease in the glucose yield was observed due to the significant transformation of glucose to side products.
Effect of Nb content.The effect of Nb content in NbZr catalyst on the yield of glucose was investigated in the presence of NbZr-W samples.The experiments were carried out at 473 K, a pressure of 25 MPa, a cellulose/catalyst ratio of 1/5, and a stop of water flow of 20 min.Kinetic curves of glucose accumulation are shown in Figure 5B.In the presence of zirconium oxide not promoted by NbO x , a 20% yield of glucose was detected.Interestingly, the major product yield decreased to 13% when 2% of niobium was in the NbZr-W1 catalyst.Subsequently, glucose yields rose to 21 and 26% when Nb amount was 4 and 17%, respectively.(Figure 5B).
A similar trend of product yields on the niobium content we reported previously for cellulose processing in a batch reactor [27].Decreasing activity of the catalyst, which had a minor Nb amount, accounted for the formation of isolated structures of Nb.Increasing the Nb content resulted in the formation of Nb-associated acid centers and, in turn, increased catalyst acidity.Thus, among the catalysts under study, NbZr-W3 catalyst was the most perspective, and a glucose yield of 26% was reached.
On the other hand, when the cellulose/catalyst ratio was 1/1(2.5), a decrease in the glucose yield was observed due to the significant transformation of glucose to side products.Effect of Nb content.The effect of Nb content in NbZr catalyst on the yield of glucose was investigated in the presence of NbZr-W samples.The experiments were carried out at 473 K, a pressure of 25 MPa, a cellulose/catalyst ratio of 1/5, and a stop of water flow of 20 min.Kinetic curves of glucose accumulation are shown in Figure 5B.In the presence of zirconium oxide not promoted by NbOx, a 20% yield of glucose was detected.Interestingly, the major product yield decreased to 13% when 2% of niobium was in the NbZr-W1 catalyst.Subsequently, glucose yields rose to 21 and 26% when Nb amount was 4 and 17%, respectively.(Figure 5В).
A similar trend of product yields on the niobium content we reported previously for cellulose processing in a batch reactor [27].Decreasing activity of the catalyst, which had a minor Nb amount, accounted for the formation of isolated structures of Nb.Increasing the Nb content resulted in the formation of Nb-associated acid centers and, in turn, increased catalyst acidity.Thus, among the catalysts under study, NbZr-W3 catalyst was the most perspective, and a glucose yield of 26% was reached,

Catalyst Characterization
Analyses of the Nb and Zr content in samples were carried out by inductively coupled plasma-atomic emission spectrometry (ICP-AES) using a Perkin-Elmer instrument OPTIMA 4300.
The porous structures of samples were determined by the adsorption of N 2 at 77 K using an ASAP 2020 Plus the Quantachrome Autosorb-iQ instrument proposed by Micromeritics.The specific surface area (S BET ) was calculated from the adsorption data using the maxima arrangement modification of the BET method over the relative pressure range between 0.05 and 0.20.The total pore volume (V Σ ) was evaluated at a p/p 0 of 0.99.
The diffraction images (XRD) were obtained with the use of a Bruker D8 Advabced diffractometer (Manheim, Germany) with a non-monochromate Cu Kα radiation (λ = 1.5418A), focusing geometry θ-θ in the scanning mode within the range of angles from 3 to 75 • , with a step of 0.05 • .
Infrared (IR) spectra were recorded by Shimadzu FTIR-8300S spectrometer (Kyoto, Japan) with a DRS-8000 diffusion reflectance cell in the range of 400 and 6000 cm −1 .Spectra resolution was 4 cm −1 .All spectra are presented in the F(R) Kubelka-Munk scale (Equation ( 4)): where R is the reflection coefficient.IR spectroscopy using pyridine as probe molecules was used for the analysis of surface acidity (Supplementary Materials, Section S2.3).
Ultraviolet-visible diffuse reflectance spectroscopy (DR-UV-Vis) spectra were recorded on a UV-2501 PC Shimadzu spectrometer with an IRS-250A accessory in the 190-900 nm range with a resolution of 2 nm on samples in powder form placed into a special cell for DR-UV-Vis measurement.BaSO 4 was used as the standard for measurements.
The water pH suspension (pH susp ) of the catalysts was determined as follows.A 10 mL vial with water was purged with Ar for 10 min, and then 100 mg of the catalyst was added to water and the vial was hermetically sealed.The suspension was stirred on a shaker.The pH of mixture was measured in certain time intervals until equilibrium was reached.The pH of the reaction solutions was measured with an Anion 4100 pH tester (TD Anion, Novosibirsk, Russia).
The total organic carbon balance (TOC) was measured by Multi N/C 2100S TOC Analyzer (Analytik Jena, Jena, Germany).

Catalytic Test in a Batch Reactor
Hydrothermal stability tests were carried out in a batch reactor (Autoclave Engineers, Erie, PA, USA) at 453 K, Ar pressure 1 MPa, and vigorous stirring 1000 rpm.To investigate catalyst stability, a catalyst (450 mg) and water (45 mL) were placed into the autoclave.The stability test time was 3 h.Leaching of Zr and Nb after hydrothermal treatment was carried out by ICP-AES analysis.
Hydrolysis-dehydration of cellulose was investigated in the batch reactor first at 453 K under an argon atmosphere.High-pressure autoclave produced by Autoclave Engineers (USA) was used for the batch experiments.Initial cellulose loading was 10 g•L −1 .The amount of catalyst varied in the range of 0.25-10 g•L −1 .The volume of the reaction mixture was 45 mL.
Analysis of the reaction mixtures was carried out by HPLC.Shimadzu Prominence LC-20 apparatus was used.The chromatograph was equipped by a RI detector and Rezex ROA-Organic Acids column (Phenomenex, 300 × 5.0 mm, Torrance, CA, USA).Additional information on catalyst testing techniques can be found in our previous papers [27,36,55].

Catalytic Test in a Flow Reactor
Hydrolysis-dehydration of cellulose under flow conditions was carried out in a set-up equipped with a pump, back pressure regulator, and temperature controller [16].The scheme of the apparatus of hydrolysis-dehydration of cellulose under flow conditions is presented in Figure S6 (Supplementary Materials).The reactor was a stainless-steel tube (length 200 mm, inner diameter 2.7 mm).The experiments were carried out at a temperature of 473 K and 25 MPa of water.The experiment time was 4 h.A typical experiment was carried out as follows.A mechanical mixture of a catalyst (25-250 mg) and cellulose (250 mg) was placed in the reactor.Pure water was pumped through a thermostated reactor at a rate of 0.5 mL•min −1 .The installation worked in the stop-flow mode.Water was pumped through the reactor for 40 min, then the flow was stopped for 20 min [16].At the outlet of the reactor, the reaction mixture was cooled to room temperature and collected in separate flasks every hour for analysis of the reaction products.
Analysis of reaction mixtures to determine the content of glucose and other products of cellulose depolymerization was carried out by HPLC on an Ultimate 3000 Thermo Scientific ultrahigh pressure liquid chromatograph (UHPLC) (Waltham, MA, USA) equipped with a Shodex SH1011 (Showa Denco, Tokyo, Japan) column thermostated at 333 K and a RI detector.An aqueous solution of 5 mM H 2 SO 4 was used as the eluent.The eluent flow rate was 0.6 mL•min −1 .
The yield of the products was calculated as follows (Equation ( 5)) [55]: where Y-a product yield (%), C product -a product concentration (mol L −1 ), V-reaction mixture volume (L), m cellulose -cellulose loading (g), Nc-stoichiometric coefficient between a glucan unit and a product (1 for glucose and 5-HMF, 6 for formic acid etc.), and m glucanmolar weight of glucan unit in cellulose (162 g mol −1 ).The initial reaction rate was calculated according to Equation ( 6): where R-initial reaction rate (mmol L −1 h −1 ); C Glu and C 5-HMF -concentration of glucose and 5-HMF (mmol L −1 ), respectively, achieved at the second kinetic point (1h); C 0 Glu and C 0 5-HMF -concentration of glucose and 5-HMF (mmol L −1 ), respectively, the first kinetic point; and t reaction time (h).

Figure 2 .
Figure 2. Influence of Nb content on a band gap (Eg).

5 Figure 2 .
Figure 2. Influence of Nb content on a band gap (E g ).

Figure 4 .
Figure 4. (A) Correlation between the yields of glucose and 5-HMF and the LAS/BAS ratio.(B) Correlation between the initial reaction rate and the density of acid sites.(C) Correlation between the TOFAS and the density of acid sites.

3 Figure 4 .
Figure 4. (A) Correlation between the yields of glucose and 5-HMF and the LAS/BAS ratio.(B) Correlation between the initial reaction rate and the density of acid sites.(C) Correlation between the TOF AS and the density of acid sites.

Figure 5 .
Figure 5. (A) Kinetic curves of glucose yield at different NbZr-W2/cellulose ratios in the cellulose hydrolysis-dehydration under flow conditions.(B) Kinetic curves of glucose yield at different Nb contents in the NbZr-W catalysts in the cellulose hydrolysis-dehydration under flow conditions.(C) Correlations between Nb content in the NbZr-W samples and the yield of glucose for 4 h and the density of acid sites on the surface of NbZr-W catalysts (AS).(Experimental conditions: 473 К, water pressure pump 20 MPa, stop flow, 20 min.).

Figure 5 .
Figure 5. (A) Kinetic curves of glucose yield at different NbZr-W2/cellulose ratios in the cellulose hydrolysis-dehydration under flow conditions.(B) Kinetic curves of glucose yield at different Nb contents in the NbZr-W catalysts in the cellulose hydrolysis-dehydration under flow conditions.(C) Correlations between Nb content in the NbZr-W samples and the yield of glucose for 4 h and the density of acid sites on the surface of NbZr-W catalysts (ρ AS ).(Experimental conditions: 473 K, water pressure pump 20 MPa, stop flow, 20 min.).

Table 1 .
Textural characteristics of the NbZr catalysts.

Table 2 .
The BAS and LAS concentrations and BAS strength in the NbZr-W, NbZr-MAW, and Nb 2 O 5 samples.Amount of BAS; 2 Amount of LAS; 3 Density of acid sites on the surface of solid (BAS + LAS)/S BET , µmol•m −2 . 1

Table 3 .
Cellulose conversion, yields of glucose and 5-HMF and initial reaction rate of cellulose hydrolysis-dehydration 1 .

Table 4 .
Catalytic properties of the heterogeneous catalysts in the transformation of cellulose via hydrolysis-dehydration in pure water.

Table 5 .
Conditions for the synthesis of Nb-Zr catalysts.