Although Fischer-Tropsch (FT) synthesis was discovered by Franz Fischer and Hans Tropsch in the 1920s in Germany [1], it has gained immense attention in last few years due to depletion of non-renewable energy sources. FT synthesis is an environmental friendly route for alternative fuels and can produce liquid fuels from carbon sources by coal-to-liquid (CTL), natural gas-to-liquid (GTL) and biomass-to-liquid (BTL) [2] processes. Three types of reactors have been utilized commercially for FT synthesis: Fixed bed, fluidized bed, and slurry bubble column bed by leading GTL companies like Shell, Sasol, Exxon Mobil, and Energy Int. [3]. There is a minimum scale limit of this FT process to be economical; for example, the Pearl GTL, a collaboration between Shell and Qatar petroleum, producing 140,000 bpd (barrels per day) is considered as a profitable economic scale for the FT GTL mesoporous silica supports: SBA-15 and KIT-6. In order to translate new advancements in the laboratory as well as industry on both catalysis and microreactors for FT synthesis, chemical kinetics is a key issue in developing mathematical models for the reactors. However, to our knowledge, the kinetics of FT synthesis using mesoporous materials in a microreactor is relatively unknown in the literature. So, in order to understand more about the interaction between silica mesoporous materials and the metal, and especially kinetics, three different types of silica mesoporous materials (MCM-41, SBA-15, and KIT-6) containing cobalt and ruthenium metals were synthesized by one-pot hydrothermal method. To address the thermodynamic stability of the catalysts, CO-conversion using these three catalysts in the 3-D printed SS microchannel microreactor was also investigated.
The textural properties of the catalysts were evaluated using nitrogen Brunner-Emmett-Teller (BET) physisorption analysis. Table 1 shows the BET surface area, pore volume and the average pore diameter of all three catalysts. The surface areas of the catalysts are different depending upon the type of silica support. While the surface area of CoRu-MCM-41 was 1025 m 2 /g, that of CoRu-SBA-15 and CoRu-KIT-6 was around 691 m 2 /g and 690 m 2 /g, respectively. The general trend is consistent with that reported in the literature [48,52,53]. The pore diameter in the range of 3.2-5.3 nm was obtained from BJH desorption plot. The pore volume was in the range of 0.77-0.92 cm 3 /g. Figure 1a shows nitrogen adsorption-desorption isotherms for all the catalysts with pore size distribution of mesoporous materials. These isotherms represent the category of Type IV isotherms which is typical for mesoporous materials as mentioned in IUPAC classification [54]. These isotherms are classified into three different types of regions. The initial part of isotherm is a linear increment of nitrogen uptake at lower relative pressures (P/P 0 = 0-0.2) called Type II isotherm. This is due to the adsorption of N 2 on monolayer and multilayer within the pore walls of the catalyst. For relative pressure in the range of P/P 0 = 0.2-0.4, there is an exponential increment in the isotherms which indicates the ordered mesoporous structure of the catalysts. Especially, the steepness of CoRu-MCM-41 is sharp when compared to the other two samples which indicate that MCM-41 support is more ordered in the nature of all the catalysts. Finally, the third region, in the relative pressure range of P/P 0 = 0.4-0.95, has a long plateau for all the catalysts and it corresponds to the multilayer adsorption on the outer surface of the catalyst. The hysteresis loop for the samples is associated with condensation of N 2 uptake in the interstitial voids of mesopores of the support [55]. Figure 1b shows pore size distribution obtained from BJH desorption plots. A sharp single peak for pore size with narrow distribution is observed for all the catalysts covering uniformly the pores with sizes in the range of 3.2 to 5.3 nm as shown in Table 1. The pore sizes of KIT-6 support appear to be larger and wider when compared to that of MCM-41 and SBA-15 supports. Furthermore, the pore distribution of MCM-41 is bi-modal, having major pores distributed in the range of 2 to 3 nm and minimal pore distribution between 3-4 nm. a = Variation range ±2%, b = Variation range ±3%, c = Variation range ±5%.
The metal loadings in the catalysts (wt%) and the surface morphology were obtained by SEM-EDX analysis. Figure S1 shows the SEM-EDX images of a typical MCM-41 catalyst showing uniform metal distribution with porous morphology. The actual and intended metal loadings are quite similar (Table 1) and suggest that the one-pot hydrothermal synthesis is one of the best routes to prepare mesoporous materials with uniform metal distribution. This uniformity plays a key role in the activity of FT catalysts; it not only decreases sintering but also increases the thermal stability of the catalysts for long-term studies.
The size of the metal particles and the structure of the mesoporous support in all catalysts were obtained from TEM studies. The high magnification images in Figure 2 show the uniform ordered hexagonal pores present in the support. It is also worth noting that MCM-41 support has well defined hexagonal pores when compared to KIT-6 and SBA-15 and this is consistent with the BET surface area and the low angle XRD studies (discussed below). A uniform metal distribution with black dots, as shown in Figure S2, having almost circular in shape is clearly evident in the mesoporous silica matrix.
The metal loadings in the catalysts (wt%) and the surface morphology were obtained by SEM-EDX analysis. Figure S1 shows the SEM-EDX images of a typical MCM-41 catalyst showing uniform metal distribution with porous morphology. The actual and intended metal loadings are quite similar (Table 1) and suggest that the one-pot hydrothermal synthesis is one of the best routes to prepare mesoporous materials with uniform metal distribution. This uniformity plays a key role in the activity of FT catalysts; it not only decreases sintering but also increases the thermal stability of the catalysts for long-term studies.
The size of the metal particles and the structure of the mesoporous support in all catalysts were obtained from TEM studies. The high magnification images in Figure 2 show the uniform ordered hexagonal pores present in the support. It is also worth noting that MCM-41 support has well defined hexagonal pores when compared to KIT-6 and SBA-15 and this is consistent with the BET surface area and the low angle XRD studies (discussed below). A uniform metal distribution with black dots, as shown in Figure S2, having almost circular in shape is clearly evident in the mesoporous silica matrix. In order to obtain information about the structural phases of the catalysts, XRD studies were carried out. Figure 3 shows the small angle XRD diffraction patterns for different mesoporous silica supported catalysts. The variations of peaks are probably due to the presence of metal nanoparticles present in the catalyst. For CoRu-MCM-41 catalyst, a sharp intense peak between 2-theta values 2-3° and two broad peaks between 2-theta values 4-5.5° corresponds to (100), (110), and (200) reflections of hexagonal mesoporous structure. This confirms that these catalysts are highly ordered mesoporous in nature with no deformation of hexagonal framework even after the addition of metals and this is consistent with the observed TEM images. For CoRu-SBA-15 catalyst, the peak between 2-theta value 1-2° indicates the mesoporous structure with 2D hexagonal symmetry with p6mm space group and
In order to obtain information about the structural phases of the catalysts, XRD studies were carried out. Figure 3 shows the small angle XRD diffraction patterns for different mesoporous silica supported catalysts. The variations of peaks are probably due to the presence of metal nanoparticles present in the catalyst. For CoRu-MCM-41 catalyst, a sharp intense peak between 2-theta values 2-3 • and two broad peaks between 2-theta values 4-5.5 • corresponds to (100), (110), and (200) reflections of hexagonal mesoporous structure. This confirms that these catalysts are highly ordered mesoporous in nature with no deformation of hexagonal framework even after the addition of metals and this is consistent with the observed TEM images. For CoRu-SBA-15 catalyst, the peak between 2-theta value 1-2 • indicates the mesoporous structure with 2D hexagonal symmetry with p6mm space group and long range ordered mesoporous structure [56]. For CoRu-KIT-6 catalyst, the peak at 2-theta value 0.94 • corresponds to (211) plane and two low intensity peaks between 1.5-2 • ascribes to (420) and (332) diffraction planes. These planes confirmed the characteristic three-dimensional nature of mesoporous KIT-6 reported in the literature [57].
In order to obtain information about the structural phases of the catalysts, XRD studies were carried out. Figure 3 shows the small angle XRD diffraction patterns for different mesoporous silica supported catalysts. The variations of peaks are probably due to the presence of metal nanoparticles present in the catalyst. For CoRu-MCM-41 catalyst, a sharp intense peak between 2-theta values 2-3° and two broad peaks between 2-theta values 4-5.5° corresponds to (100), (110), and (200) reflections of hexagonal mesoporous structure. This confirms that these catalysts are highly ordered mesoporous in nature with no deformation of hexagonal framework even after the addition of metals and this is consistent with the observed TEM images. For CoRu-SBA-15 catalyst, the peak between 2-theta value 1-2° indicates the mesoporous structure with 2D hexagonal symmetry with p6mm space group and long range ordered mesoporous structure [56]. For CoRu-KIT-6 catalyst, the peak at 2-theta value 0.94° corresponds to (211) plane and two low intensity peaks between 1.5-2° ascribes to (420) and (332) diffraction planes. These planes confirmed the characteristic three-dimensional nature of mesoporous KIT-6 reported in the literature [57]. The wide-angle XRD (WAXRD) analysis was carried out to determine the crystallinity of metal oxides in different mesoporous supports. Figure 4 shows the WAXRD patterns of these samples. The observed 2θ angles are compared with the JCPDS (Joint Committee on Powder Diffraction Standards) database. For all catalysts, the peaks at 18.90
To determine the oxidation states of Co and Ru in MCM-41, KIT-6, and SBA-15, XPS studies were performed. Figure S3 shows the XPS spectra of Si 2p and O 1s containing a single spectrum which is centered at 104 eV and 532 eV, respectively, and confirms the presence of silicates in the sample. Figure 5a shows the Co 2p spectra for all the samples; the Co 2p 3/2, and Co 2p 1/2 peaks are clearly observed to indicate the presence of cobalt in two oxidation states in the silica matrix [60][61][62].
The peaks centered at 780.5 eV and 796.2 eV are associated with Co 2p 3/2 and Co 2p 1/2 , respectively in the MCM-41 matrix. Whereas, in the case of KIT-6, the peaks for Co 2p 3/2 and Co 2p 1/2 are observed at 779.6 eV and 795.8 eV, respectively. For SBA-15, the similar peaks are noticed at 779.7 eV and 794.8 eV, respectively. It is clear from these data that the binding energy for cobalt in the MCM-41 matrix is distinctly higher when compared to that of cobalt in KIT-6 and SBA-15. This suggests that cobalt in two oxidation states in MCM-41 is in a different environment from that of SBA-15 and KIT-6. This is also consistent with temperature programmed reduction (TPR) profile showing much higher reduction temperatures for CoRu-MCM-41 catalyst as discussed below. Similar XPS spectra for Co 2p were observed and analyzed by Bhoware et al., [63]. Figure 5b shows the XPS spectra for the ruthenium metal in the catalyst. The presence of Ru in the sample is confirmed by the Ru 3d spectra which is centered almost at 284.8 eV and it is associated with the Ru 3d 3/2 oxidation state [64]. However, in contrast to cobalt, there is no significant difference in the binding energy of the Ru metal ions in different mesoporous silica supports.
Catalysts 2019, 9, x FOR PEER REVIEW 6 of 29
The wide-angle XRD (WAXRD) analysis was carried out to determine the crystallinity of metal oxides in different mesoporous supports. Figure 4 shows the WAXRD patterns of these samples. The observed 2θ angles are compared with the JCPDS (Joint Committee on Powder Diffraction Standards) database. For all catalysts, the peaks at 18.90° (111), 31.09° (220), 36.74° (311), 38.36° (222), 44.72° (400), 59.25° (511), and 65.26° (440) correspond to the cubic structure of Co3O4 (JCPDS-42-1467) [58,59]. The orthorhombic structure of RuO2 (JCPDS-88-0323) is consistent with the observed peaks at 28.18° (110), 35.27° (101), and 54.56° (211) in all the catalysts.
To determine the oxidation states of Co and Ru in MCM-41, KIT-6, and SBA-15, XPS studies were performed. Figure S3 shows the XPS spectra of Si 2p and O 1s containing a single spectrum which is centered at 104 eV and 532 eV, respectively, and confirms the presence of silicates in the sample. Figure 5a shows the Co 2p spectra for all the samples; the Co 2p3/2, and Co 2p1/2 peaks are clearly observed to indicate the presence of cobalt in two oxidation states in the silica matrix [60][61][62]. The peaks centered at 780.5 eV and 796.2 eV are associated with Co 2p3/2 and Co 2p1/2, respectively in the MCM-41 matrix. Whereas, in the case of KIT-6, the peaks for Co 2p3/2 and Co 2p1/2 are observed at 779.6 eV and 795.8 eV, respectively. For SBA-15, the similar peaks are noticed at 779.7 eV and 794.8 eV, respectively. It is clear from these data that the binding energy for cobalt in the MCM-41 matrix is distinctly higher when compared to that of cobalt in KIT-6 and SBA-15. This suggests that cobalt in two oxidation states in MCM-41 is in a different environment from that of SBA-15 and KIT-6. This is also consistent with temperature programmed reduction (TPR) profile showing much higher reduction temperatures for CoRu-MCM-41 catalyst as discussed below. Similar XPS spectra for Co 2p were observed and analyzed by Bhoware et al., [63]. Figure 5b shows the XPS spectra for the ruthenium metal in the catalyst. The presence of Ru in the sample is confirmed by the Ru 3d spectra which is centered almost at 284.8 eV and it is associated with the Ru 3d3/2 oxidation state [64]. However, in contrast to cobalt, there is no significant difference in the binding energy of the Ru metal ions in different mesoporous silica supports. Temperature programmed reduction (TPR) is an ideal technique to analyze the reduction behavior of metal oxides in mesoporous silica. It helps to investigate the interaction between metal and the support by providing information on physiochemical properties of the material. All the calcined catalysts are treated with 10%H2 to record TPR profiles for Co and Ru metal oxides shown in Figure 6. All the samples contain well defined peaks for ruthenium at low reduction temperatures and cobalt at much higher reduction temperatures. The TPR profiles of all the catalysts show that the ruthenium oxide is reduced with H2 at a relatively lower temperature between 100 °C and 250 °C with main hydrogen consumption peaks for Ru 3+ to Ru 0 which is also reported by Panpranot et al.
[65]. However, the reduction behavior of Co3O4 (Co3O4 → CoO → Co 0 ) [66] in three samples to Co 0 is remarkably different depending on the type of support. For MCM-41, as reported by Lim et al., the Temperature programmed reduction (TPR) is an ideal technique to analyze the reduction behavior of metal oxides in mesoporous silica. It helps to investigate the interaction between metal and the support by providing information on physiochemical properties of the material. All the calcined catalysts are treated with 10%H 2 to record TPR profiles for Co and Ru metal oxides shown in Figure 6. All the samples contain well defined peaks for ruthenium at low reduction temperatures and cobalt at much higher reduction temperatures. The TPR profiles of all the catalysts show that the ruthenium oxide is reduced with H 2 at a relatively lower temperature between 100 • C and 250 • C with main hydrogen consumption peaks for Ru 3+ to Ru 0 which is also reported by Panpranot et al. [65]. However, the reduction behavior of Co 3 O 4 (Co 3 O 4 → CoO → Co 0 ) [66] in three samples to Co 0 is remarkably different depending on the type of support. For MCM-41, as reported by Lim et al., the small peak centered around 310 • C ascribes the reduction of cobalt to CoO, while the second main peak corresponds to the reduction of CoO to metal ions Co 2+ into the silica network [67]. The last hydrogen uptake has a peak centered almost 780 • C which suggests that the cobalt and the MCM-41 support have strong interaction which is also confirmed by the binding energy obtained from XPS in Figure 5a [48,68]. This could be due to the formation of a spinel structure as cobalt silicates [69] and consistent with the XPS and XRD data. Unlike MCM-41, the reduction temperatures of Co inSBA-15 and KIT-6 were quite low around 365 • C and 375 • C, respectively, confirming weaker metal interactions with SBA-15 and KIT-6 supports [53]. However, no separate three peaks were observed for the reduction of cobalt, this may be due to the absence of silicates in SBA-15 and KIT-6 samples. The shift in the reduction peaks to the lower temperatures can also be due to the incorporation of Ru metal in the support [70]. Although the 5% weight of Ru is maintained in the catalyst sample, there might be a slight difference in the actual loadings of the Ru metal as shown in EDX Table 2. Qin et al., have studied the effect of the Ru metal on Co-SBA-15 catalyst at different loading and found a remarkable effect on the activity of catalyst during the FT synthesis [70]. Thus, the overall interactions of metal-metal and metal-support have a strong influence on the reducibility and reactivity of the catalysts for FT synthesis. Since the operating temperature zone of the FT synthesis is less than 350 • C, the activity of the catalyst is more dependent on the ease of reducibility of the metal oxides to pure metals (active sites) in the support. Although MCM-41 exhibits higher hydrogen consumption than other catalysts, the amount of hydrogen adsorbed by CoRu-KIT-6 in the temperature range of 25-311 °C is higher when compared to that by CoRu-MCM-41 and CoRu-SBA-15 in the same range of temperature.
In order to have a better understanding of the effect of metal and support interaction on catalysts, kinetic studies were carried out in the SS microchannel microreactors. The main difficulty to describe the FT kinetics is the complexity of its mechanism and a larger number of possible chemical species involved. Kinetic models of FT synthesis using Co based catalysts are less abundant than Fe based catalysts in literature [71]. Most of the existing models are mainly based on power-law models where Langmuir-Hinselwood (LH) type equations have been used by different researchers [72][73][74][75][76]. Although the simple power-law expression is widely recognized in the field of catalysis, it was recognized to have limited application in FT synthesis due to the narrow range of reaction conditions [77,78]. However, LH type equations are widely used for prediction of rates over a wide range of reaction conditions. As an example, Yates and Satterfield [72] worked on Co-catalyst and fitted the rate data obtained at 220-240 °C. They found that the rate data were best fitted with simple LH expression. Rautavuoma and van dar Baan [79] reported the rate of reaction at 1 atm pressure and 250 °C . They observed that reaction proceeds through CO dissociation and formation of -CH2-surface intermediate. The amount of hydrogen consumed by CoRu-MCM-41, CoRu-SBA-15, and CoRu-KIT-6 catalysts was calculated and quantified to be 0.108, 0.05, and 0.032 mmol H 2 /gm, respectively in the temperature range of 25 to 1000 • C. The amount of hydrogen consumed by CoRu-MCM-41 was found almost two times more than that of CoRu-SBA-15 and 3 times more than by CoRu-KIT-6. Higher hydrogen uptake by CoRu-MCM-41 is most likely due to the reduction of Co-silicates ~750 • C.
Although MCM-41 exhibits higher hydrogen consumption than other catalysts, the amount of hydrogen adsorbed by CoRu-KIT-6 in the temperature range of 25-311 • C is higher when compared to that by CoRu-MCM-41 and CoRu-SBA-15 in the same range of temperature.
In order to have a better understanding of the effect of metal and support interaction on catalysts, kinetic studies were carried out in the SS microchannel microreactors. The main difficulty to describe the FT kinetics is the complexity of its mechanism and a larger number of possible chemical species involved. Kinetic models of FT synthesis using Co based catalysts are less abundant than Fe based catalysts in literature [71]. Most of the existing models are mainly based on power-law models where Langmuir-Hinselwood (LH) type equations have been used by different researchers [72][73][74][75][76]. Although the simple power-law expression is widely recognized in the field of catalysis, it was recognized to have limited application in FT synthesis due to the narrow range of reaction conditions [77,78]. However, LH type equations are widely used for prediction of rates over a wide range of reaction conditions. As an example, Yates and Satterfield [72] worked on Co-catalyst and fitted the rate data obtained at 220-240 • C. They found that the rate data were best fitted with simple LH expression. Rautavuoma and van dar Baan [79] reported the rate of reaction at 1 atm pressure and 250 • C. They observed that reaction proceeds through CO dissociation and formation of -CH 2 -surface intermediate.
In the microchannels of the microreactor, the flow of reactants is basically laminar. The complexity of the microchannel microreactor increases due to the parabolic type velocity profile; so, an average velocity profile is approximated during the development of the model for the microreactor system [80]. The outlet concentrations of the limiting reactant (CO), which was related to the rate of reaction, were calculated by an in-line GCMS. The following differential equation was used for a reactor model defined as Equation ( 1):
W cat = Wt. of the catalyst F in,CO = Molar feed rate of CO X CO = Conversion of CO -r CO = Disappearance rate of CO Equation (2), below, is used to calculate the disappearance rate of CO
The following boundary conditions (BC) were used:
The partial pressure of the compound was calculated using the following equations:
where p i is the partial pressure of the component, P T is the total pressure of the reactor at the inlet (1 atm) and N c is the total number of components. m i is the number of moles of component i.
In order to determine the most suitable kinetic model for a particular catalyst, all possible combinations of FT reactions were considered and rate equations were developed based on CO conversion. A number of Langmuir-Hinshelwood and Eley-Rideal models have been developed for kinetics of CO hydrogenation to hydrocarbons in different types of reactors over the past few years [80][81][82]. In this study, it was assumed that the FT reactions occur only at active sites and proposed six possible mechanisms for FT reactions to develop kinetic models in the microchannel microreactor as shown in Table 2. In order to derive an appropriate model that describes a suitable FT equation, we considered six cases with different elementary reaction steps for each case as shown in Table 2.
No Elementary Reaction
In FT-1, CO is adsorbed on active site (*) of catalyst to form a CO* intermediate. Then, CO* reacts with H 2 to give the products C (hydrocarbons) and D (H 2 O). Similarly, H 2 can be adsorbed on the catalyst site (*) to form H 2 * intermediate and this intermediate subsequently reacted with CO to yield products in the FT-2 mechanism. There are two steps of adsorption in the FT-3 mechanism. In 1 st and 2 nd steps, CO and H 2 both are adsorbed on catalyst active site (*) to form two intermediates (CO* and H 2 *). In the last step (surface reaction), these two intermediates react with each other to give products C and D. The FT-4 model consisted of three different steps. In the 1 st step (adsorption), CO is adsorbed on catalyst site (*) to form CO*. In 2 nd step (surface reaction), the intermediate (CO*) reacts with H 2 to form another intermediate (COH 2 *). In the last step (desorption), the final intermediate gave products (C and D). Like FT-4, FT-5 also consists of three different steps-adsorption, surface reaction, and desorption. In the 1 st step (adsorption), H 2 is adsorbed on catalyst site (*) to form the Using the six models described above, six rate equations can be deduced for FT reactions by considering surface reaction and rate-limiting desorption as shown in Table 3. (See Appendix A for the rate equation derived using the FT-3 kinetic model).
Rate Controlling Step (RCS) Kinetic Equation
All the models presented in Table 3 were verified against experimental data to obtain the best suitable mechanism with the best fit. The models FT-1, FT-2, FT-4, and FT-5 are based on Eley-Rideal-type mechanism. In this case, one reactant gets adsorbed and another reactant reacts directly from the gas phase to form intermediates that yield products. Other models FT-3 and FT-6 are based on the Langmuir-Hinshelwood mechanism, which means all reactants are adsorbed on the catalyst surface before the products are formed. The kinetic parameter (k) and equilibrium constants (K 1 , K 2 , K 3 ) at each temperature were evaluated by non-linear regression analysis based on Levenberg-Marquart algorithm in POLYMATH software by minimizing the sum of squared residuals of reaction rates [83,84]. The objective function is defined as:
where, N is the number of total observations, r cal i and r exp i are calculated from the model equation and experimental rates at the different reaction temperatures.
The rate constants and the equilibrium constants can be related to The Arrhenius equation and van't Hoff laws as shown below:
where k i and K i are reaction and equilibrium constants, respectively. E ai and ∆H i are the apparent activation energy and standard enthalpy change of i species.
The influence of the weight hourly space velocity (WHSV) on the CO conversion for three different catalysts at 1 atm and H 2 /CO molar ratio 2 with error bar is shown in Figure 7a-c. The reactions were carried out at three different temperatures (210 • C, 240 • C, and 270 • C). CO conversion increases with the increase of space velocity and temperature. While CO conversion increases quickly with the increase of space velocity at the beginning, it remains almost constant at higher space velocity as the reaction reaches the equilibrium state. The variation of CO conversion was within 10% as observed during these reactions.
reactions were carried out at three different temperatures (210 °C , 240 °C , and 270 °C). CO conversion increases with the increase of space velocity and temperature. While CO conversion increases quickly with the increase of space velocity at the beginning, it remains almost constant at higher space velocity as the reaction reaches the equilibrium state. The variation of CO conversion was within 10% as observed during these reactions.
The kinetic models derived from Langmuir-Hinshelwood and Eley-Rideal mechanisms consider elementary reactions consuming CO and H2 to produce hydrocarbons and water. Recently, the mechanistic aspects of FT synthesis were well investigated by computational catalysis studies
The kinetic models derived from Langmuir-Hinshelwood and Eley-Rideal mechanisms consider elementary reactions consuming CO and H 2 to produce hydrocarbons and water. Recently, the mechanistic aspects of FT synthesis were well investigated by computational catalysis studies using DFT-based quantum chemical models [85][86][87]. However, the use of Langmuir-Hinshelwood and Eley-Rideal models facilitates understanding of the FT mechanism more easily. In this work, all the kinetic models derived based on these mechanisms are investigated and fitted with experimental data to check the feasibility of the proposed mechanism. The objective function (F) in Equation (4) was utilized to measure the goodness of the model to select the best-fitted mechanism for FT synthesis. The kinetic parameters obtained for all the mechanisms for CoRu-MCM-41 are shown in Table 5. It can be inferred from data that the value of the rate constant (k) increases with increasing temperature with only one of the 6 mechanisms which is FT-3. Therefore, for CoRu-MCM-41, the model FT-3 is best fitted with the kinetic data. In order to determine the activation energy and the frequency factor from the Arrhenius equation, the logarithm of the rate constant was plotted against the inverse of reaction temperature as shown in Figure 8. The activation energy was determined to be 32.21 kJ/mol and the frequency factor was 4099 kmol/KgCat.hr. (atm) 2 for CoRu-MCM-41 catalyst. Table 6 shows that the rate constant increases with the increase of reaction temperature. However, the two adsorption equilibrium constants (K 1 and K 2 ) decrease with the increase of reaction temperature. Since adsorption is an exothermic process, the adsorption equilibrium constant decreases with rise in temperature.
In order to determine the activation energy and the frequency factor from the Arrhenius equation, the logarithm of the rate constant was plotted against the inverse of reaction temperature as shown in Figure 8. The activation energy was determined to be 32.21 kJ/mol and the frequency factor was 4099 kmol/KgCat.hr. (atm) 2 for CoRu-MCM-41 catalyst. Table 6 shows that the rate constant increases with the increase of reaction temperature. However, the two adsorption equilibrium constants (K1 and K2) decrease with the increase of reaction temperature. Since adsorption is an exothermic process, the adsorption equilibrium constant decreases with rise in temperature. Table 4 shows the experimental data of the kinetic runs for CoRu-SBA-15 at 210 °C, 240°C, and 270 °C. When the data are fit against all the kinetic models, FT-2 was the best-fitted model obtained for CoRu-SBA-15. The kinetic parameters for all of the proposed mechanisms are shown in Table 6. Table 4 shows the experimental data of the kinetic runs for CoRu-SBA-15 at 210 • C, 240 • C, and 270 • C. When the data are fit against all the kinetic models, FT-2 was the best-fitted model obtained for CoRu-SBA-15. The kinetic parameters for all of the proposed mechanisms are shown in Table 6.
The activation energy and frequency factor of this catalyst were evaluated from the Arrhenius equation by plotting the logarithm of the rate constant to the inverse of reaction temperature as shown in Figure 9. The activation energy was determined to be 13.39 kJ/mol and the frequency factor was 254 kmol/KgCat.hr. atm. Table 6 shows that for FT-3 mechanism, the reaction rate constant increases with reaction temperature and adsorption equilibrium constant (K 1 ) decreases with the increase of reaction temperature. Table 4 shows the experimental data for CoRu-KIT-6 catalyst. All the models were fitted with this experimental kinetic data and FT-6 was found to be the best fitted model for CoRu-KIT-6. The kinetic parameters for all the proposed mechanisms are shown in Table 7 and the activation energy from Figure 10 was determined to be 12.59 kJ/mol and the frequency factor was 39 kmol/KgCat.hr. atm. on our experimental and kinetic model data, it can be concluded that only one of the six mechanisms for each catalyst is statistically relevant for fitting the model. The rate constant (k), for some of the other five mechanisms, does not show an increasing trend with the increase in temperature or remains constant, while the equilibrium constants, K 1 and K 2 , did not show the decreasing trend with the increase in temperature. Thus, FT-3, FT-2, and FT-6 mechanisms were considered as kinetically relevant model equations for CoRu-MCM-41, CoRu-SBA-15, and CoRu-KIT-6, respectively.
The results from our studies in a microreactor are similar to those reported in literature. Mansouri et al., developed a similar mechanism to estimate kinetic parameters for FT synthesis using cobalt-based catalyst with silica support and found that the experimental data were best fitted with surface reaction mechanism proposed based on Langmuir-Hineshelwood model and the optimal activation for the proposed kinetic model was found to be 31.57 kJ/mol [88]. Very recently, Sonal et al., detailed mechanistic approach for FT synthesis based on Langmuir-Hinshelwood-Hougen-Watson (LHHW) and Eley-Rideal using Fe-Co based catalyst. They claimed that a mechanism based on the adsorption and desorption have a satisfactory fit to the experimental data with the activation energies for the formation of methane, paraffin and olefin to be around 70 kJ/mol, 113 kJ/mol, and 91 kJ/mol, respectively [89]. In order to improve the efficiency of FT synthesis significantly, detailed kinetic rate expressions were derived which is very similar to our work reported in literature for both fixed bed reactor as well as microreactor using iron or cobalt-based catalyst [80,81,90,91]. The elementary steps in the above studies were used to develop kinetic mechanisms considering FT reactions with and without water gas shift (WHS) reactions occurring on the surface of the catalysts forming intermediates with active sites. A similar approach was considered in this present study where CO*, COH 2 * are assumed to form as intermediates, where * is an active site of the catalyst. From the activation energy calculations, shown in Figures 8910, the FT activation energy is observed in the order, CoRu-MCM-41 > CoRu-SBA-15 > CoRu-KIT-6. Almost 20 kJ/mol less activation energy was obtained for SBA-15 supported catalyst than that of MCM-41 catalyst. The activation energy of KIT-6 supported catalyst is a bit less than that of SBA-15 supported catalyst. This reflects that activation energy depends on metal-support interactions in different mesoporous catalysts. The variation of activity in different mesoporous catalysts might arise due to experimental uncertainties and operating conditions [92]. In addition, the FT activation energy is sensitive to the reactor system. Sun et al., [93] reported that the activation energy in a microchannel microreactor is smaller than that observed in a fixed bed reactor (FBR).
Figure 11 shows the variation of the model predicted rate with the experimental rate for all catalysts. The best fitted models i.e., FT-3, FT-2, and FT-6 for CoRu-MCM-41, CoRu-SBA-15, CoRu-KIT-6, respectively, were chosen to plot the graph for predicted and experimental rates. The correlation coefficient (R 2 ) for all cases was above or equal to 0.95 with an error band of ±15% to ±30%. This indicates that the error between experimental and predicted values lies within the statistical permissible limits at all reaction temperatures for all the catalysts and consistent with the mechanistic models proposed in the literature. Moazami et al., conducted kinetic studies for FT synthesis in a fixed bed reactor with cobalt-based catalyst over silica support and found that 60% of the results were predicted with a relative error of less than 15%, while the rest of the proposed kinetic models has error less than 32% with confidence interval of 0.99 [94]. They also proposed a pseudo-homogenous one-dimensional model to evaluate the kinetic performance of the catalyst and achieved less than 8% error with the predicted data for kinetic experiments [95]. More recently, Marchese et al., performed kinetic studies with Co-Pt/γ-Al 2 O 3 catalyst in a lab-scale tubular reactor and reported an error band around ± 25% with a confidence level of 0.95 stating it lies in the suitable acceptable limits with many mechanistic models proposed in the literature [89,[96][97][98].
Sun et al., [93] reported that the activation energy in a microchannel microreactor is smaller than that observed in a fixed bed reactor (FBR).
Figure 11 shows the variation of the model predicted rate with the experimental rate for all catalysts. The best fitted models i.e., FT-3, FT-2, and FT-6 for CoRu-MCM-41, CoRu-SBA-15, CoRu-KIT-6, respectively, were chosen to plot the graph for predicted and experimental rates. The correlation coefficient (R 2 ) for all cases was above or equal to 0.95 with an error band of ±15% to ±30%. This indicates that the error between experimental and predicted values lies within the statistical permissible limits at all reaction temperatures for all the catalysts and consistent with the mechanistic models proposed in the literature. Moazami et al., conducted kinetic studies for FT synthesis in a fixed bed reactor with cobalt-based catalyst over silica support and found that 60% of the results were predicted with a relative error of less than 15%, while the rest of the proposed kinetic models has error less than 32% with confidence interval of 0.99 [94]. They also proposed a pseudo-homogenous one-dimensional model to evaluate the kinetic performance of the catalyst and achieved less than 8% error with the predicted data for kinetic experiments [95]. More recently, Marchese et al., performed kinetic studies with Co-Pt/γ-Al2O3 catalyst in a lab-scale tubular reactor and reported an error band around ± 25% with a confidence level of 0.95 stating it lies in the suitable acceptable limits with many mechanistic models proposed in the literature [89,[96][97][98].
In order to further understand how the interaction of Co and Ru metals with different mesoporous silica supports affects the stability of the FT catalysts, deactivation studies were performed. Figure 12 shows the deactivation rates of the catalysts tested continuously for 60 h at 240
In order to further understand how the interaction of Co and Ru metals with different mesoporous silica supports affects the stability of the FT catalysts, deactivation studies were performed. Figure 12 shows the deactivation rates of the catalysts tested continuously for 60 h at 240 • C, 1 atm, and H 2 :CO ratio of 2:1. All the catalysts maintained fairly consistent CO conversion with very little fluctuation during the first 10 h. More specifically, the catalysts maintained 65%-79% CO conversion in the first 10 h of the reaction with CoRu-KIT-6 exhibiting the highest conversion and CoRu-MCM-41, the lowest. The activity of all the catalysts dropped by 20% after 24 h and then started to decline further. At the end of 60 h, the activity of the CoRu-MCM-41 dropped by 70% whereas, in the case of CoRu-KIT-6 and CoRu-SBA-15, the CO conversion decreased by 64% and 84%, respectively. Our results suggest that in terms of stability, the support has a significant impact on FT performance. More significantly, MCM-41 and KIT-6 supports are more stable when compared to the FT stability studies with SBA-15.
The reagents used for catalysis synthesis were of analytical grade with no further purification. Tetramethyl orthosilicate, 99% (TMOS) and ammonium hydroxide, Tetraethyl orthosilicate reagent grade, 98% (TEOS), Pluronic acid (P123), Hydrochloric acid (HCl), Cetyltrimethylammoniumbromide (CTAB), Co(NO3)2.6H2O, RuCl3•xH2O were purchased from Sigma Aldrich. Ethanol (anhydrous), Butanol and acetone, ACS grade, were obtained from Fischer Scientific, Branchburg, New Jersey, USA.
The microchannel microreactor and the respective cover channel were fabricated using 3D printing technology. Typically, the microreactor and its cover channel are designed using AutoCAD software which is schematically shown in Figure 13. The design is based on the split and recombination principle which has 11 microchannels of 500 μm × 500 μm × 2.4 cm as reaction zone in between them. This stainless-steel 3D printed microreactor is assembled in a custom-built heating block with an inlet and outlet system which facilitates the flow of syngas through the channels. Many deactivation mechanisms have been proposed for FT studies that include catalysts poisoning, sintering, oxidation, the effect of water, carbon deposition and surface reconstruction [99]. The syngas used in our studies is a mixture of ultrahigh pure 5.0 CO and H 2 gases; therefore, there is very little or no chance of catalyst deactivation due to poisoning by the gas feed at the inlet to the reactor. It was also observed in our previous studies that the support (SiO 2 vs TiO 2 ) can enhance the stability of the catalyst to resist deactivation [38,40]. Iglesia et al., noticed that silica supported materials are less stable when compared to the other supports like Al 2 O 3 [100] in their FT studies with Co-based catalyst. They also reported that CoRu-TiO 2 and CoRu-SiO 2 were found to have almost the same activation energy upon the addition of Ru to the Co catalyst; however, there were strong differences in the deactivation rates of catalysts depending upon the support [101]. Based on our CO-conversion studies, the ability of catalysts to withstand or retard the FT deactivation rate was in the order of CoRu-KIT-6 > CoRu-MCM-41 > CoRu-SBA-15.
The reagents used for catalysis synthesis were of analytical grade with no further purification. Tetramethyl orthosilicate, 99% (TMOS) and ammonium hydroxide, Tetraethyl orthosilicate reagent grade, 98% (TEOS), Pluronic acid (P123), Hydrochloric acid (HCl), Cetyltrimethylammoniumbromide (CTAB), Co(NO 3 ) 2. 6H 2 O, RuCl 3 •xH 2 O were purchased from Sigma Aldrich. Ethanol (anhydrous), Butanol and acetone, ACS grade, were obtained from Fischer Scientific, Branchburg, New Jersey, USA.
The microchannel microreactor and the respective cover channel were fabricated using 3D printing technology. Typically, the microreactor and its cover channel are designed using AutoCAD software which is schematically shown in Figure 13. The design is based on the split and recombination principle which has 11 microchannels of 500 µm × 500 µm × 2.4 cm as reaction zone in between them. This stainless-steel 3D printed microreactor is assembled in a custom-built heating block with an inlet and outlet system which facilitates the flow of syngas through the channels.
For the synthesis of CoRu-SBA-15: TEOS, CTAB, water, ethanol, pluronic acid, and hydrochloric acid were mixed in molar ratios of 1:0.081:41:7.5:0.0168:5.981. In a typical synthesis procedure, P123 was dissolved in 2M HCl at 35 °C till a clear solution was obtained. Another solution was prepared by dissolving CTAB in DI water at 35 °C until a homogenous mixture was produced. These two solutions were mixed and stirred for 35 min. Ethanol containing metal precursors were added dropwise into the solution and stirred for 30 min. Afterwards, TEOS which was limiting reagent in this procedure was also added dropwise and stirred for 20 h at 35 °C. The aqueous solution was aged for 48 h at 98 °C followed by air drying for 24 h. The material was then oven-dried at 110 °C for 24 h. Finally, the dried material is then calcined in stepwise fashion with heating rate of 1 °C/min at 350 °C, 450 °C, and 550 °C for 8 h each, respectively to remove CTAB and pluronic acid.
In the case of CoRu-KIT-6, TEOS, P123, HCl, DI water, and butanol were mixed in a molar ratio of 1:0.017:1.83:195:1.31 [102]. For a typical procedure, P123 was added to HCl at 35 °C till a clear solution was obtained. A separate solution was prepared with butanol containing metal precursors and poured to the previous solution and stirred until a homogeneous solution was obtained. To this mixture, TEOS, which was the limiting reagent, was added dropwise and stirred at 500 rpm for 24 h. The final solution was aged for 24 h at 100 °C, followed by air drying for 24 h under the fume hood. The material is oven-dried at 110 °C for 24 h and then calcined at 550 °C for 4 h, to remove P123, the structure directing agent, SDA, with heating and cooling rates of 1 °C/min.
The catalyst is loaded into the microchannels of the microreactor using a PVA suspension containing the catalyst, DI water, binder PVA (polyvinyl alcohol 98%-99% hydrolyzed MW: 31000) and acetic acid of weight ratio 1:5:0.25:0.05. This suspension with well-dispersed catalyst was dipcoated and dried in air and then calcined in presence of air at 400 °C for 2 h with heating and cooling rates of 5 °C/min. Figure 13d shows the SEM image of the catalyst coated microreactor prior to the in-situ reduction.
Three types of Co-Ru based nanocatalysts supported by different mesoporous silica supports i.e., MCM-41, SBA-15 and KIT-6 supports are used in this study. A constant metal loading of 10%Co and ~5% Ru in weight was maintained in all preparations and this metal loading was also determined using the amount of the precursor. Three catalysts using different mesoporous support-10%Co5%Ru-MCM-41, 10%Co5%Ru-SBA-15 and 10%Co5%Ru-KIT-6-were synthesized using the one-pot hydrothermal procedure ( as shown below) [48]. The catalysts were labeled as CoRu-MCM-41, CoRu-SBA-15, and CoRu-KIT-6 in this manuscript.
For the synthesis of CoRu-MCM-41, TMOS, CTAB, DI-water, and ethanol were used in a molar ratio of 1:0.13:130:20 as described elsewhere [48]. In short, CTAB was dissolved in DI-water at 30 • C to produce a clear solution. The metal precursors were dissolved in ethanol in a separate beaker. The precursor, TMOS, which is a limiting agent for this synthesis, was added dropwise to the mixture of the two solutions prepared previously. Ammonium hydroxide was added dropwise to precipitate metal hydroxides in the solution, till the final pH was ~10. The precipitate was stirred for 3 h, followed by 18 h of aging at 65 • C. The precipitate was then washed with DI water till the filtrate reached a pH of 7, finally washed with ethanol and filtered. The filtered material is air-dried for a day and then oven-dried at 110 • C for 24 h. The dried catalyst is calcined at 550 • C for 16 h with a ramp rate of 2 • C/min to remove the CTAB template.
For the synthesis of CoRu-SBA-15: TEOS, CTAB, water, ethanol, pluronic acid, and hydrochloric acid were mixed in molar ratios of 1:0.081:41:7.5:0.0168:5.981. In a typical synthesis procedure, P123 was dissolved in 2M HCl at 35 • C till a clear solution was obtained. Another solution was prepared by dissolving CTAB in DI water at 35 • C until a homogenous mixture was produced. These two solutions were mixed and stirred for 35 min. Ethanol containing metal precursors were added dropwise into the solution and stirred for 30 min. Afterwards, TEOS which was limiting reagent in this procedure was also added dropwise and stirred for 20 h at 35 • C. The aqueous solution was aged for 48 h at 98 • C followed by air drying for 24 h. The material was then oven-dried at 110 • C for 24 h. Finally, the dried material is then calcined in stepwise fashion with heating rate of 1 • C/min at 350 • C, 450 • C, and 550 • C for 8 h each, respectively to remove CTAB and pluronic acid.
In the case of CoRu-KIT-6, TEOS, P123, HCl, DI water, and butanol were mixed in a molar ratio of 1:0.017:1.83:195:1.31 [102]. For a typical procedure, P123 was added to HCl at 35 • C till a clear solution was obtained. A separate solution was prepared with butanol containing metal precursors and poured to the previous solution and stirred until a homogeneous solution was obtained. To this mixture, TEOS, which was the limiting reagent, was added dropwise and stirred at 500 rpm for 24 h. The final solution was aged for 24 h at 100 • C, followed by air drying for 24 h under the fume hood. The material is oven-dried at 110 • C for 24 h and then calcined at 550 • C for 4 h, to remove P123, the structure directing agent, SDA, with heating and cooling rates of 1 • C/min.
The catalyst is loaded into the microchannels of the microreactor using a PVA suspension containing the catalyst, DI water, binder PVA alcohol 98%-99% hydrolyzed MW: 31000) and acetic acid of weight ratio 1:5:0.25:0.05. This suspension with well-dispersed catalyst was dip-coated and dried in air and then calcined in presence of air at 400 • C for 2 h with heating and cooling rates of 5 • C/min. Figure 13d shows the SEM image of the catalyst coated microreactor prior to the in-situ reduction.
Specific surface area, pore size, pore volume and TPR studies of the catalyst were carried out using Micromeritics, 3-Flex instrument. The Brunner-Emmett-Teller (BET) method was used to calculate the surface area of the catalyst where an equation was obtained from adsorption isotherm in the relative pressure range of 0.07-0.03. The surface area was calculated from adsorption isotherm in the relative pressure range of P/P 0 = 0.07-0.3 using the Brunner-Emmett-Teller (BET) equation. The total volume per gram of catalyst was determined from the amount of N 2 adsorbed at P/P 0 = 1. The N 2 desorption from the catalyst surface provides information about the pore size distribution using BJH (Barret-Joyner-Halenda) plots [103]. The H 2 temperature programmed reduction (TPR) analysis was also done with the same instrument which has a TCD detector to monitor the reduction signals of the catalyst. Around 50 mg of the catalyst was loaded into the quartz sample tube in which a stream of 10% H 2 /Ar at flowrate 110 mL/min was passed through and the temperature is increased to 1000 • C with 10 • C/min ramp rate. The small and wide-angle powder x-ray diffraction (XRD) were recorded using D8 Discover X-ray and Rigaku SmartLab X-ray diffractometers, respectively, with Cu K-alpha radiation (wavelength = 0.15418 nm) radiation generated at 40 mA and 40 kV. The step size and time per step used in these measurements are 0.05 • and 3 secs/step, respectively. The crystal sizes of the metal oxides were determined using the Scherrer equation. In the Scherrer equation below, τ stands for the crystal size, λ is the wavelength of the Cu Kα radiation, β is the full width half maximum and θ is the Bragg diffraction angle. τ = 0.9λ β * Cosθ (7) The morphology and the size of the catalysts were analyzed using transmission electron (TEM Carl Zeiss Libra 120) at 120 KeV and scanning electron microscopy (Zeiss Auriga FIB/FESEM). The sample for TEM was prepared by dispersing a small quantity of catalyst in 3 mL of ethanol followed by vortex dispersion and sonication for a few minutes. Then the suspension was drop coated on a carbon-coated copper grid of 300 µm mesh size, followed by drying in an oven at 100 • C for 12 h.
The elemental composition and oxidation states of the metals were analyzed using Energy Dispersive X-ray spectrometry (Zeiss Auriga FIB/FESEM obtained from Carl Zeiss, Oberkochen, Germany) and oxidation states by X-ray photoelectron spectroscopy (XPS-Escalab Xi+-Thermo Scientific obtained from Thermo Scientific, West Sussex, UK), respectively.
An in-house LabVIEW automated experimental setup was built to carry out the FT experiments for precise control over the operating conditions. The experimental setup is shown in Figure 14. The flowrates of the syngas mixture which is a mixture of hydrogen and carbon monoxide were controlled by precalibrated mass flow controllers obtained from cole parmer with flow rates ranging from 0-1 sccm. Nitrogen was used as a carrier gas into the system and was controlled by Aalborg mass flow controller with a maximum flow rate of 10 sccm. The upstream and downstream pressures were continuously monitored by pressure gauges obtained from Cole-Parmer and the data are fed to Aalborg solenoid valve from which the reaction pressure is controlled and kept constant throughout the reaction. All these controllers are operated by LabVIEW 2018 program. The product stream is directly fed to the GC-MS (Agilent Technologies 7890B GC and 5977 MSD). Prior to the start of the kinetic experiments, the microreactors were reduced ex-situ in Carbolite Gero tubular furnace with 10% H 2 Ar. To compensate the losses while transferring the microreactor to the heating block the microreactor containing the catalyst was reduced again in-situ for 6 h at 350 • C before the start of FT reaction. The kinetic studies were performed by varying the weight hourly space velocity (WHSV = Wcat/F CO,in , where W Cat = weight of the catalyst and F CO,in = molar flow rate of CO in feed) in the range of ~25.2-3.15 kgcat.h/kmol. The reactions were performed with syngas having a feed molar ratio (H 2 /CO) of 2:1 at 210 • C, 240 • C, and 270 • C while the reaction pressure was maintained at 1 atm. Based on our previous FT studies using this setup and preliminary runs, all reactions reached a steady state after an hour at each setpoint of WHSV. Deactivation studies were also performed for all three catalysts at 240 • C using syngas feed molar ration of 2:1. CO conversion was calculated based on following equation:
Three different types of mesoporous silica supported Co-Ru based catalysts were synthesized using the one-pot hydrothermal method and performance for FT was evaluated. These catalysts resulted in high surface area with hexagonal ordered mesoporous structure as supported by BET, low angle XRD and TEM studies. The interaction between the metal and different types of support has a significant effect on the kinetic and stability studies of FT synthesis. Six mechanistic models
Three different types of mesoporous silica supported Co-Ru based catalysts were synthesized using the one-pot hydrothermal method and performance for FT was evaluated. These catalysts resulted in high surface area with hexagonal ordered mesoporous structure as supported by BET, low angle XRD and TEM studies. The interaction between the metal and different types of support has a significant effect on the kinetic and stability studies of FT synthesis. Six mechanistic models were developed based on the Langmuir-Hinshelwood and Eley-Radiel mechanisms. The best fitted model for all catalysts was obtained on the basis of non-linear regression by comparing with the objective function which is equation-7 in this paper. The proposed model FT-3 was best fitted with the kinetic data for CoRu-MCM-41 catalyst. Whereas, FT-2 and FT-6 were well fitted with the kinetic data for CoRu-SBA-15 and CoRu-KIT-6, respectively. CoRu-KIT-6 was found to be more active than other catalysts with a low activation energy of 12.59 kJ/mol, whereas for CoRu-SBA-15 and CoRu-MCM-41 the activation energies are 13.39 and 32.21 kJ/mol, respectively. An average error of 5.65%, 1.76%, and 3.70% was obtained for catalysts CoRu-MCM-41, CoRu-SBA-15, CoRu-KIT-6, respectively considering the best fitted model explained above for FT synthesis. The predicted data provided by kinetic models were satisfactory with the experimental data. These results highlight the potential of the mechanistic FT models as well as reaction mechanisms to further improve the performance of FT synthesis. In addition, this information can help to design more active and selective catalysts for the optimized FT process. Furthermore, all catalysts exhibited significant resistance to the deactivation rate following the order CoRu-KIT-6> CoRu-MCM-41>CoRu-SBA-15. This study suggests that even if the support is of same type, the structure of the support plays a vital role in catalyst performance for FT synthesis.