Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 21027 (2023 ) Cite this article Reforming Catalyst
This study aims to investigate the role of hydrogen on CNTs synthesis and kinetics of CNTs formation. The CNTs were synthesized by catalytic chemical vapor deposition of methane over FeMo/MgO catalyst. The experimental results revealed that hydrogen plays an important role in the structural changes of catalyst during the pre-reduction process. The catalyst structure fully transformed into metallic FeMo phases, resulting in an increased yield of 5 folds higher than those of the non-reduced catalyst. However, the slightly larger diameter and lower crystallinity ratio of CNTs was obtained. The hydrogen co-feeding during the synthesis can slightly increase the CNTs yield. After achieving the optimum amount of hydrogen addition, further increase in hydrogen would inhibit the methane decomposition, resulting in lower product yield. The hydrogenation of carbon to methane was proceeded in hydrogen co-feed process. However, the hydrogenation was non-selective to allotropes of carbon. Therefore, the addition of hydrogen would not benefit neither maintaining the catalyst stability nor improving the crystallinity of the CNT products. The kinetic model of CNTs formation, derived from the two types of active site of dissociative adsorption of methane, corresponded well to the experimental results. The rate of CNTs formation greatly increases with the partial pressure of methane but decreases when saturation is exceeded. The activation energy was found to be 13.22 kJ mol−1, showing the rate controlling step to be in the process of mass transfer.
Carbon nanotubes (CNTs) were discovered over three decades ago in 19911. The unique properties of CNTs have been extensively researched and applied in several applications such as semi-conductors, batteries, reinforced materials and universal composite materials2. The synthesis of CNTs can be simply achieved by a transformation of carbon-containing substances under suitable conditions to carbon atoms, recombination of six-member ring graphene structure. The catalytic decomposition of hydrocarbons such as methane can produce CNTs grown on the catalyst surface, while hydrogen (H2) as gas product is produced3,4,5,6. From above, H2 gas generally found during CNTs synthesis. There are many researched show that all of these process success to produce a CNTs, which are H2 free process7, H2 was used only to pre-reduction of catalyst process but not to CNTs formation process8, or H2 was used in the entire process from reduction of catalyst to synthesis of CNTs process9,10 as shown in Table 1.
The roles of H2 on CNTs synthesis has been extensively studied. Piedigrosso et al.11 reported that addition of H2 would reduce carbon yield by elimination of both CNTs and amorphous carbon through hydrogenation. In contrast, some research works found that the co-feeding of H2 during CNTs synthesis can increase CNT yield12,13,14,15,16,17,18 by selective elimination of amorphous carbon. The other effects of H2 were widely reported; the H2 concentration made different CNTs diameter14,15,18, included effect to formation of SWCNTs some researcher found H2 hinder formation of SWCNTs19,20. On the other hand, the H2-free or less concentration of H2 process could not synthesis a SWCNTs21, including to MWCNTs22. The reason of this effect was reported in the same way that is H2 has a role to control the size, phase, and morphology of catalysts. For CNTs synthesis, catalyst is the one necessary factor to control CNTs formation, for example the same Fe-based catalyst but one is α-Fe phase and other one is Fe3C phase, the formation of CNTs occur with a huge different in yield, diameter, and CNT characteristics. He et al.23 and Torres et al.6 reported Fe3C phase resulted in the formation of bamboo-like CNTs, additionally Fe3C forming smaller diameter CNTs but less yield than α-Fe12,24.
Research has shown that the addition of Mo can increase the efficiency of CNTs synthesis by increasing yield and enhancing better stability at high temperatures8,25. Furthermore, mechanisms of the effects of hydrogen on the synthesis of CNTs using FeMo/MgO as a catalyst have not been reported.
Recently, the CNTs manufactures have been developed largely that directly affecting to the continued decline in sale prices and expected that the price will further go down globally26, even so CNTs are still expensive, approximately price of SWCNTs and MWCNTs are 38,000฿27 and 3400฿28 per gram, respectively. The market for MWCNTs has a decline in large-scale production; however, there still remains global demand of > 2000–2500 tons per annum with increased demand in composites, automotive and aerospace applications and especially as battery additives in Asia29. In contrast, H2 demand is increasing consistently. Therefore, to compete in the CNTs market, manufacturers should balance costs to suit the CNTs properties.
This work focuses on studying the role of H2 to whole process of CNTs synthesis by using FeMo/MgO as a catalyst and finding out in-depth behavior of each process to guide further improvements in CNTs manufacturing process as well as its reaction kinetic.
FeMo/MgO catalyst with 30%wt metals loading was prepared by impregnation method. The iron (III) nitrate nonahydrate (Fe(NO3)3.9H2O, 98%, AR, Loba chem.) and ammonium heptamolybdate ((NH4)6Mo7O24, AR, KEMAUS ) solution were prepared as precursors of Fe and Mo with mass ratio of Fe:Mo equal to 2:1, then slowly drop of Mo solution into Fe solution to avoid precipitation. The FeMo solution was dropped on magnesium oxide fine particle (MgO, pharma, Applichem Panreac.). The mixture was stirred and dried on hot plate stirrer at 90 °C until forming of paste, then dried in oven with the same temperature for 4 h and calcination at 500 °C for 3 h. The FeMo/MgO in oxide form is ready to use for the reaction. The images of the experimental process are shown in Fig. 1.
Preparation of FeMo/MgO by impregnation method.
CNTs were synthesized by the catalytic chemical vapor deposition (CCVD) method using methane as a carbon source. The schematic diagram was shown in Fig. 2. Firstly, 0.5 g of FeMo/MgO catalyst was put in a quartz boat, then placed in the quartz tube reactor center. This CCVD process started with N2 fed into the reactor and heated up from room temperature until it reached desired temperature by linearly increasing at a rate of 10 °C min−1.
Schematic diagram of CCVD method.
The total gas flow rate was controlled with following each condition, which calibrated by soap film meter. The effluent gas was collected and analyzed with Gas Chromatography (GC-TCD 8A, Shimadzu). Conditions for GC analysis were set as follows (3 mm ∅ × 2 m column, INJ/DET: 120 °C, COL: 100 °C, active carbon packing: 0.2–0.25 mm, and He carrier: 40 ml min−1). Until the end of the reaction (180 min), the system was cooled to room temperature in N2 gas flow. The CNTs were formed at the surface of the catalyst in the quartz boat.
The catalytic performance in terms of CH4 conversion, Carbon yield, g-CNTs/g-catalysts, and calculated H2 flow rate were calculated according to Eqs. (1)–(4), respectively.
The kinetics investigation in “Kinetic study of CNTs synthesis” section will perform experiment with H2 pre-reduction process (R-woH2). The product yield used to calculate yield rate as an Eq. (5). The reduced catalyst weight after the reduction process will be 0.75 times to the initial catalyst weight. Rate of reaction will be calculated by Eq. (6)
The study of morphology, internal, and external structure was examined using imaging process. The external structure of as-prepared and purified product was observed by Field Emission Scanning Electron Microscopes (FE-SEM, HITACHI SU-8010). The internal structure, type of CNTs, number of walled, CNTs growth mechanism (Tip, based growth), CNTs and catalyst diameter size were observed by Transmission Electron Microscopes (TEM, JEOL JEM-2100 Plus). Sample preparation for TEM analyzed was prepared by following step; the product was dispersed in ethanol (99.5v/v%, Sigma-Aldrich) and then dropped onto a 300 mesh of copper grid coated with carbon film. As-prepared CNTs and catalysts diameter size were measured from SEM, and TEM images by using ImageJ software. Crystallinity and structure of CNTs and catalysts were analyzed for study of changing in the sample, because of the controlled processes. X-ray diffraction (XRD) with the angle range scan 5°–80° and TEM selected area electron diffraction (SAED) were used to specify crystallinity, phase change, and crystal size. The crystallinity of the synthesized carbon products, which validated with qualitative analysis of the proportion between crystalline carbon and amorphous carbon was carried out using 532 nm laser light source Raman spectroscopy (Raman NT-MDT model: NTEGR-SPECTRA). The purity of as-prepared product was measured from the thermo-gravimetric analysis (TGA), by observing the weight loss in expected temperature range.
In the study roles of H2, there are two variables that need to be studied which are, the pre-reduction of catalyst and the H2 feeding during CNTs formation. Four processes: H2 free process (no pre-reduction process and without H2 feed during CNTs formation, nR-woH2), H2 co-feed process (no pre-reduced process and with H2 feed during CNTs formation, nR-wH2), H2 pre-reduction process (pre-reduction process and without H2 feed during CNTs formation, R-woH2), and H2 combined process (pre-reduction process with H2 feed during CNTs formation, R-wH2) which concluded in Table 2, were used to investigate H2 roles through catalytic activity and product properties with following results.
The catalytic activity was examined over CCVD for synthesis of CNTs process. From the experiment, FeMo/MgO catalyst is capable of converting CH4 in all processes as shown by methane conversion in Fig. 3.
Catalytic activity through (solid line) methane conversion and (dotted line) product weight.
Methane conversions decreased over time because of catalyst deactivation. The formation of product CNTs obstructed the active surface of catalysts. The catalyst performance depends on methane conversion, deactivated rate, and quantity of final product. It was found that R-woH2 was the most effective process with a 3.41 g final product followed with R-wH2 (3.33 g), nR-wH2 (1.19 g), and nR-woH2 (0.93 g), respectively. The thermal stability of the product was observed with thermo-gravimetric analysis (TGA) in air atmosphere, as shown in Fig. 4.
TGA of as-prepared CNTs by comparing each H2 processes.
From Fig. 4, it was found that product begins to decompose at temperatures of 600 °C and above. Decomposition curves showed similar decomposed characteristics or referring to the same group of products. Residuals from decomposition under air atmosphere are catalyst and supported materials which are highly thermal stability than carbon product, therefore we can get carbon content from TGA. The percentage carbon content corresponds to the yield. Products from pre-reduction of catalyst process (R-woH2, R-wH2) have higher carbon content than non-reduced processes.
The morphology of as-prepared products from various processes was investigated by image processing through scanning electron microscope (SEM) as shown in Fig. 5. Dense CNTs were found in all SEM images included with some amorphous carbon. Considering the CNTs products, the diameter of MWCNTs was measured from over 200 samples at 10,000×. magnification combined with over 100 samples at 50,000×. magnification by using ImageJ software to produce a diameter size distribution of as-prepared CNTs from each process.
SEM images of as-prepared CNTs with magnification ×10,000, ×50,000, and CNTs diameter size distribution from (a) nR-woH2, (b) nR-wH2, (c) R-woH2, and (d) R-wH2, respectively.
The diameter size distribution of CNTs represents non-uniform product, ranging from tiny CNTs (smaller than 5 nm) to huge CNTs (larger than 100 nm), as a result of improper catalyst preparation method, thus affecting the active size which is expressed as the diameter of CNTs36. Considering the mean diameter, CNTs from non-reduced processes: nR-woH2 and nR-wH2 were approximately the same sizes (29.29 and 32.72 nm) meanwhile, the reduced process: R-woH2 and R-wH2 have similar in diameter (43.25 and 45.41 nm) but were larger than CNTs-nR processes.
Qualitative analysis of the proportion between crystalline carbon and amorphous carbon was carried out using 532 nm laser light source Raman spectroscopy. Raman is used to determine proportions between graphitize carbon and disorder or amorphous carbon by the relative height of G band (~ 1580 cm−1) and D band (~ 1350 cm−1) that represent graphitic and disorder carbon, respectively. The G/D peak intensity ratios (IG/ID) of a carbonaceous product were investigated more than 5 points for each process to confirm the results of the analysis and prevent false sampling. Examples of the Raman analysis result are shown in Fig. 6. Average IG/ID ratio and other detailed data for the study of overview H2 roles are shown in Table 3.
The representive Raman shifts of CNTs prepared by different conditions.
From the Raman analysis, it was found that the products formed by non-reduced processes have higher IG/ID ratio than reduced processes which approximately are 3.1 and 2.5, respectively. Observed that the addition of H2 during the formation of CNTs did not cause a difference in the crystallinity ratio (more detailed in “Hydrogenation reaction” section). Therefore, there was a possibility that the crystallinity ratio was also controlled by the catalyst structure.
In the study of H2 roles, the samples were clearly divided into two groups which are non-reduced processes and reduced process. Non-reduced process (nR-woH2 and nR-wH2) produces products with smaller diameters and higher crystallinity ratios than reduced processes (R-woH2 and R-wH2), but reduced processes have a higher yield than 5 times (~ 1 to ~ 5) of non-reduced process.
The structure of the catalyst and as-prepared product in each process were investigated by XRD method to describe the mechanisms that cause changes in CNTs. Firstly, the H2-temperature programmed reduction (H2-TPR) analysis was used to investigate H2 consumption during temperature changes, which is a similar experimental to pre-reduction of catalyst process. In the reduced process, it was confirmed by H2-TPR that pre-reduction of catalyst from room temperature to 900 °C would allow the catalyst to remain structurally stable or without structural changes in the presence of hydrogen in the system3, as shown in Figure S1, that H2 is not consumed at 900 °C and above.
XRD patterns of catalysts consisted of calcined MgO, calcined FeMo/MgO, and reduced FeMo/MgO are shown in Fig. 7. The catalyst phases were identified by the diffraction peaks. The diffraction peaks of magnesium oxide (MgO, PDF 01-071-1176) were observed in Fig. 7a which were cubic structures. Figure 7b exhibited the diffraction peaks of MgO impregnated with iron (Fe) and molybdenum (Mo), where only monoclinic MgMo4 at 23.2° (021), and 26.3° (220) 37 was found. Xu et al.38 reported an XRD analysis of calcined FeMo/MgO, revealing the structure of MgO, Fe3O4, and MgMo4. The results of the xrd diffraction were similarly blunt same as this work. XRD patterns of reduced FeMo/MgO was show in Fig. 7c, tetragonal FeMo exhibited the diffraction peak at 37.3° (002), 44.6° (411), and 64.7° (223)39.
XRD patterns of (a) calcined MgO (b) calicned FeMo/MgO, (c) reduced FeMo/MgO.
Figure 7 found that the calcination process did not detect Fe-Mo interaction, FeMo was found after the catalytic reduction process. Therefore, the reduction of catalyst process changes the catalyst structure from the oxide form to metallic structure and establishes an interaction between iron and molybdenum.
Comparison of product structure between reduced FeMo/MgO, non-reduced CNTs (nR-woH2), and reduced CNTs (R-woH2) is shown in Fig. 8. FeMo undetectable after CNTs growth from both of non-reduced and reduced processes, which may cause by carbon obscuring or active metal was encapsulated within CNTs. Diffraction patterns of hexagonal carbon or carbon nanotubes was found at 25.9° (002)40, which shown in Fig. 8b and c.
Comparing XRD patterns of (a) reduced FeMo/MgO (b) CNTs-nR-woH2, (c) CNTs-nR-wH2.
The XRD patterns of as-prepared CNTs are shown in Fig. 9. It was found that the product characteristics are divided into two groups in the same way as the properties of CNTs, which are non-reduced processes and reduced processes.
XRD patterns of as-prepared CNTs from (a) nR-woH2, (b) nR-wH2, (c) R-woH2, (d) R-wH2.
The details of the XRD analysis results in the 30°–60° and 68°–80° ranges are shown in Fig. 10A and B, respectively. The diffraction pattern exhibited iron carbide (Fe3C) that can be found in all processes. Meanwhile, diffraction peak of molybdenum carbide (Mo2C) is discovered in non-reduced processes. Therefore, in a non-reduced process, the catalyst is converted from oxide to carbide form. In addition, iron, and molybdenum interacted poorly, resulting in the formation of molybdenum carbide structures. These results show why synthesized CNTs by non-reduced and reduced processes differ.
XRD patterns of as-prepared CNTs in ranges (A) 30°–60° (B) 68°–80° from (a) nR-woH2, (b) nR-wH2, (c) R-woH2, (d) R-wH2.
Several previous studies have mentioned similar behaviors with Fe catalysts. He et al.23 reported that the carbide catalyst (Fe3C) will make CNTs with a smaller diameter but lower yield than the metallic catalyst (α-Fe). In this work, we found that non-reduced processes will be made the carbide catalyst, but reduced processes will be made metallic catalyst (FeMo). This is why products from the non-reduced and reduced processes are different.
The carbide formation of catalyst is described by Behr et al.15. They studied the H2-to-CH4 ratio in the feed gas. Carbide formation will occur when lower H2– to –CH4 ratio, while the metallic formation will occur when a higher H2-to-CH4 ratio. Reduced processes in this work occur in H2 and N2 atmosphere, in the meaning there is no carbon (CH4) in the system, that made a FeMo metallic, but for non-reduced processes, catalyst will be reduced by H2 gas dissociated from CH4, because CH4 can be decomposed by the oxide form of the catalyst41. The reaction product hydrogen gas may make the reduction form of the catalyst with presence of carbon atoms on catalyst surface that caused a carbide catalyst.
Therefore, from the study of H2 roles, it was found that the characteristics of CNTs, whether yield, diameter, and proportion of crystallinity, are determined by the structure of catalyst that was controlled by the pre-reduction process.
Effects of H2 concentration were studied with H2 co-feed or nR-wH2 process. The catalytic performance was investigated at varied CH4:H2 ratio which are 1:0 (50:0 ml min−1, nR-woH2), 2:1 (50:25 ml min−1), 1:1 (50:50 ml min−1), 1:2 (50:100 ml min−1, nR-wH2), and 1:3 (50:150 ml min−1). N2 was used to adjust total flow rate to 200 ml min−1. Methane conversion was investigated and H2 flow rate was calculated by Eq. 4 then plotted against reaction time as shown in Fig. 11. As-prepared CNTs were measured for average diameter and summarizing the results with other information in Fig. 12 and Table 4.
Effect of H2 concentration on (a) CH4 conversion (b) calculated H2 flow rate over FeMo/MgO catalyst at 900 °C.
g-CNTs/g-catalyst, and average CNTs diameter as a function of CH4:H2 flow rate ratio. Error bars are the standard deviation for each sample and is a measure of the diameter distribution.
Methane conversion in each H2 concentration was slightly different as shown in Fig. 11a, where the CH4:H2 ratio was 1:1 with the highest methane conversion. In Fig. 11b, exhibited H2 flow rate that is higher than amount of H2 fed into the system. Since H2 is one of the products resulting from the decomposition of methane, it was found that H2 flow rate decreases over time corresponding to the reduced methane conversion.
Figure 12 shows the yield of the produce (g-CNTs/g-catalyst) with an optimum point of 1:1, corresponding to the conversion of methane. When the H2 content is increased above 1:1, the yield is less. Observing the flow rate of H2, found that there is high amount of H2 in system. The study of Kashiwaya Y., and Watanabe M. reported that methane decomposition is reversible reaction in gas phase before CNTs growth42. Therefore, according to Le Chatelier’s principle, when the product quantity (H2) increases, will prevent the dissociation of methane, resulting in fewer CNTs43,44. Additionally, the addition of the proper amount of H2 allows catalyst to restructure that ready for the formation of higher CNTs. Observed from H2-free process (nR-woH2) with the lowest methane conversion (~ 8%) at 5 min, then the conversion increased. Because initially the catalyst was in an oxide form with lower methane conversion performance, and when partially H2 was formed, the catalyst changed its structure and thus the conversion increased.
The results of the study on adjustment concentration of H2 gas, when adding H2, yield was higher but decreased when excess H2. Therefore, adding H2 to CNTs product has an optimal point, that consistent with various research12,13,14,15,16,17,18.
Measure and establish the CNTs diameter size distribution from the study of effects of H2 concentration as shown in Fig. 12 and Table 4, found that as-prepared CNTs from this work are highly distributed. As a result, it was unable to confirm the results of the study of the diameter changes from the adjustment of the H2 concentration.
However, several studies indicate that H2 addition helps to eliminate amorphous carbon as well as maintain catalyst stability. To study about this issue, the study of catalytic hydrogenation reaction was conducted in the next section.
This section examined the catalytic hydrogenation reaction that was thought to play an important role in the synthesis of CNTs, by investigating whether this reaction can help stabilize the catalyst during the synthesis of CNTs, included with an improvement of CNTs quality.
Piedigrosso et al.11 reported that, CNTs are capable of catalytic hydrogenation, but there is no selective eliminate of solid carbon, either crystalline carbon or amorphous carbon. Catalytic hydrogenation may enhance the catalyst performance by removing the carbon covering the catalyst surface, which caused the deactivation of catalyst. In this work, an experiment was conducted to study the synthesis of CNTs in reduced-process, by comparing between processes with and without H2 feed (R-woH2 and R-wH2) with a reaction time of 5 h (300 min) to study catalyst stability of both processes (Figure S2). It was found that the conversion of methane continues to decrease at the same rate in both processes. Comparison the product weight difference between two processes, at 3 h, the difference was about 0.07 g, but at 5 h, the difference was about 0.26 g. Therefore, it seems that addition of H2 does not keep catalyst stability, but also reduces the weight of the product.
The catalytic hydrogenation reaction was tested to prove these issues. As-prepared CNTs from R-woH2 and R-wH2 process were tested for catalytic hydrogenation reaction at 900 °C with flow rates of 100 and 50 ml min−1 of H2 and N2, respectively. Effluent gas, product weight, and crystalline ratio were analyzed to identify the behavior that occurred.
Figure 13 exhibits flow rate of CH4 and CO released from the system during the catalytic hydrogenation reaction test. It also shows the calculated quantity of CNTs consumed as a %yield. From the experiment, it was found that CO was initially formed. Methane occurred throughout the experiment, with the flow rate decreasing slightly over time. After the experiment, approximately 3% of CNTs were consumed. In addition, CH4 and CO generated rate have no different for each process, although these two samples are expected to have different incidences, because R-wH2 has been occurred in catalytic hydrogenation (during the formation of CNTs) for longer than R-woH2. If catalytic hydrogenation is a selective reaction, a difference will occur in some way.
Flow rate of CH4 and CO, and %yield as a funtion of catalytic hydrogenation reaction time.
The thermal stability of the product comparing between As-received CNTs (from process R-woH2 and R-wH2) and CNTs after process Catalytic hydrogenation reaction (CHT) is shown in Fig. 14, and the average crystallinity ratio by Raman characteristic is shown in Table 5.
TGA of CNTs for study of catalytic hydrogenation reaction.
From the thermal decomposition in Fig. 14, it was found that all samples were not different, with similar carbon content and decomposed at the same temperature range. Consider the proportion of crystallinity. Raman analysis has shown no difference in IG/ID ratio with an average of 2.6.
Therefore, from the study of the catalytic hydrogenation reaction. It can be concluded that the catalytic hydrogenation reaction occurs during the synthesis of CNTs, whereas CNTs are used as a precursor to form methane and carbon monoxide. This reaction was not selective to eliminate either crystalline or amorphous carbon. As a result, it does not improve the stability of the catalyst, by etching carbon from surface of catalyst. The addition of H2 to the system will not enhance the properties of CNTs and produce less product, because of the catalytic hydrogenation reaction and the shifted of reaction equilibrium.
From the study of the role of H2 in “Effect of the presence of H2” section, it was concluded that an important role of H2 is to determine the catalyst structure by H2 pre-reduction of catalyst before CNTs formation. Whereas the addition of H2 during the formation of CNTs was not found to have a beneficial effect on the synthesis of CNTs and resulted in less product yield.
The aim of the kinetic study of CNTs synthesis was to investigate the mechanisms occurring during the synthesis of CNTs, with the hope that the data could be used for further process scale-up. Data collection for use in kinetic studies in CNTs synthesis is quite difficult because solids form while reaction that made an unsteady state reaction. Research works have chosen a variety methods to investigated reaction rate, such as measuring the height of VACNTs45, analysis of hydrogen output gas46, or batch analysis by collecting data on changing weight of the catalyst bed47,48. In this work, batch analysis was used. The reaction rate in each process was collected as shown in Fig. 15a. The rate of product yield was calculated according to the Eq. (6), with units of gram product divided by grams of catalyst and reaction time in second. A dramatic decrease in yield rate can be observed in all processes. From the study in “Roles of hydrogen” section it was found that R-woH2 process had the highest yield and was able to reduce the complexity of analysis by removing H2 feeding. Therefore, this process was chosen for further kinetic studies with 30 min reaction time. The partial pressure of CH4 in feed gas was varied to identify the reaction order and effects of CH4 concentration. The experimental data for kinetic study was shown in Table S1. Partial pressure of CH4 was conducted by varying of CH4 flow rate.
(a) The yield rate of CNTs synthesis at 900oC, (b) Effect of partial pressure of methane on the yield rate of CNTs formation at 900oC with line indicating the calculated rates by kinetic model, (c) Effect of temperature on the CNTs yield rate, and (d) The Arrhenius plot for estimation of the activation energy.
Yadav et al.49 derived model 1 and 2 by assuming dissociative adsorption of methane followed by removal of hydrogen from the adsorbed methyl group as the rate determining step. The reaction step was considered as irreversible reaction and the difference between model 1 and 2 is the one type of active site (x) or two types (x and y). Pseudo first order reaction kinetic, kinetic model 1 and model 2 were used to compare for this reaction, leading to the following rate Eqs. (7–9). The rate constants are given by k’, k1 and k2, respectively. To estimate the parameters k, Eqs. (7–9) can be linearized and parameters were estimated from the slope and intercepts.
The experimental and the simulated values from models were plotted as shown in Fig. 15b. The partial pressure of methane equal to 1 does not consider as it is a pure methane system which is inconsistent with the other proportions. The rate of formation of CNTs increased with an increase in the partial pressure of methane, then reaches as saturation level followed by decrease in the rate at higher partial pressures of methane. The value of the rate constant from each model was shown in Table 6. Yadav et al. resulted as well. It can be seen that model 1 predicts the reaction rate same as model 2 which is better than 1st order reaction over a covered condition in this work. These results show that the decomposition of methane into CNTs and hydrogen is might be 1st order reaction but, model 1 and model 2 can predict more accurately.
The effect of temperature was studied in the range of 850–1000 °C with partial pressure of methane equal to 0.24 atm. (50 ml min−1). It was observed that the rate of CNTs formation increased with temperature until higher than 950 °C because of catalyst begin to be agglomerated that decreased catalyst performance, as shown in Fig. 15c.
To estimate the activation energy (Ea), Arrhenius plot is needed. Then the rate of reaction (lnRCNTs) was plotted versus T−1 as shown in Fig. 15d. The rate of CNTs formation at 1000 °C was not considered for Ea estimation, as agglomeration of active phase occurred.
The activation energy was found to be 13.22 kJ mol−1, compared to the other tasks in summarize table of Yadav et al.49 work, our Ea was much lower. It may be in a range of mass transfer limit50, consistent with an assumption of pseudo first order reaction that used to describe mass transfer limit the reaction rate. Additionally, our experiment found the behavior of the CNTs formation rapidly occurs, and a compaction of the bed layer will occur, which may hinder the formation of CNTs, showing that this experimental process undesignated for study mass transfer limit process. Therefore, it is possible that Ea has a lower value compared to other works. Moreover, Ea determinations should be carried out in reactions with low conversion and reaction temperatures to get accurate and reasonable values 51,52.
In this study, the role of H2 on CNTs synthesis was investigated by means of pre-reduction of catalyst and co-feeding during CNTs synthesis. For FeMo/MgO catalyst, the reduced process could produce CNTs with the largest yield ascribed to the presence of highly active FeMo metallic phases. Meanwhile, the Mo2C phase was formed in the non-reduced process, resulting in the formation of CNTs with smaller diameter and higher crystallinity. The co-feeding of H2 during CNTs synthesis in the non-reduced process can increase the yield of CNTs formation. However, the yield became lower when H2 exceeded the equilibrium of reaction. In addition, the presence of H2 during the CNTs synthesis would promote the catalytic hydrogenation of carbon to methane. This reaction was not selective for elimination of amorphous carbon. As a result, the co-feeding of H2 would not improve neither the CNTs crystallinity nor the stability of the catalyst. The kinetic study of CNTs synthesis shows that methane decomposition to CNTs and hydrogen is a pseudo first order reaction with mass transfer as rate-controlling step. The rate of CNTs formation increases with an increased methane partial pressure. The activation energy derived from the experimental data was found to be relatively low at 13.22 kJ mol−1. In sum, the synthesis of MWCNTs by CCVD of methane using FeMo/MgO catalyst can be achieved under H2 free condition. The characteristics of as-grown CNTs would greatly depend on the structure of the catalyst. The adding of hydrogen shows an important role in changing the catalyst structure during the pre-reduction of catalyst.
All data related to the finding of this study are accessible upon request from the corresponding author Apinan Soottitantawat.
Iijima, S. Helical microtubules of graphitic carbon. Nature 354, 56–58 (1991).
Article ADS CAS Google Scholar
Cheptubes. Carbon Nanotubes Properties and Applications. https://www.cheaptubes.com/carbon-nanotubes-properties-and-applications/
Torres, D., Pinilla, J., Lázaro, M., Moliner, R. & Suelves, I. Hydrogen and multiwall carbon nanotubes production by catalytic decomposition of methane: Thermogravimetric analysis and scaling-up of Fe–Mo catalysts. Int. J. Hydrog. Energy 39, 3698–3709 (2014).
Łamacz, A. & Łabojko, G. CNT and H2 production during CH4 decomposition over Ni/CeZrO2. II. Catalyst performance and its regeneration in a fluidized bed. ChemEngineering 3, 25 (2019).
Fan, Z., Weng, W., Zhou, J., Gu, D. & Xiao, W. Catalytic decomposition of methane to produce hydrogen: A review. J. Energy Chem. 58, 415–430 (2020).
Torres, D., Pinilla, J. L. & Suelves, I. Cobalt doping of α-Fe/Al2O3 catalysts for the production of hydrogen and high-quality carbon nanotubes by thermal decomposition of methane. Int. J. Hydrog. Energy 45, 19313–19323. https://doi.org/10.1016/j.ijhydene.2020.05.104 (2020).
Wang, H. et al. Selective synthesis of single walled carbon nanotubes on metal (iron, nickel or cobalt) sulfate-based catalysts. Carbon 129, 128–136 (2018).
Abdullahi, I., Sakulchaicharoen, N. & Herrera, J. E. Selective synthesis of single-walled carbon nanotubes on Fe–MgO catalyst by chemical vapor deposition of methane. Diam. Relat. Mater. 41, 84–93 (2014).
Article ADS CAS Google Scholar
Chang, K.-H., Kim, J.-D. & Yoon, K.-W. Google Patents (2018).
Inoue, S., Lojindarat, S., Kawamoto, T., Matsumura, Y. & Charinpanitkul, T. Spontaneous and controlled-diameter synthesis of single-walled and few-walled carbon nanotubes. Chem. Phys. Lett. 699, 88–92 (2018).
Article ADS CAS Google Scholar
Piedigrosso , P. , Colomer , JF , Fonseca , A. & Nagy , JB in AIP Conference Proceedings 16–19 (1998).
Pinheiro, P., Schouler, M., Gadelle, P., Mermoux, M. & Dooryhee, E. Effect of hydrogen on the orientation of carbon layers in deposits from the carbon monoxide disproportionation reaction over Co/Al2O3 catalysts. Carbon 38, 1469–1479 (2000).
Xiong, G. Y., Suda, Y., Wang, D. Z., Huang, J. Y. & Ren, Z. F. Effect of temperature, pressure, and gas ratio of methane to hydrogen on the synthesis of double-walled carbon nanotubes by chemical vapour deposition. Nanotechnology 16, 532–535. https://doi.org/10.1088/0957-4484/16/4/033 (2005).
Article ADS CAS Google Scholar
Biris, A. R. et al. Effect of hydrogen on the growth and morphology of single wall carbon nanotubes synthesized on a FeMo/MgO catalytic system. Phys. Lett. A 372, 3051–3057. https://doi.org/10.1016/j.physleta.2008.01.023 (2008).
Article ADS CAS Google Scholar
Behr, M. J., Gaulding, E. A., Mkhoyan, K. A. & Aydil, E. S. Effect of hydrogen on catalyst nanoparticles in carbon nanotube growth. J. Appl. Phys. https://doi.org/10.1063/1.3467971 (2010).
Reynolds, C., Duong, B. & Seraphin, S. Effects of hydrogen flow rate on carbon nanotube growth. J. Undergrad. Res. Phys. (2010).
Schünemann, C. et al. Catalyst poisoning by amorphous carbon during carbon nanotube growth: Fact or fiction?. ACS Nano 5, 8928–8934 (2011).
Li, Y., Ji, K., Duan, Y., Meng, G. & Dai, Z. Effect of hydrogen concentration on the growth of carbon nanotube arrays for gecko-inspired adhesive applications. Coatings https://doi.org/10.3390/coatings7120221 (2017).
Hao, Y., Qingwen, L., Jin, Z. & Zhongfan, L. The effect of hydrogen on the formation of nitrogen-doped carbon nanotubes via catalytic pyrolysis of acetonitrile. Chem. Phys. Lett. 380, 347–351. https://doi.org/10.1016/j.cplett.2003.09.031 (2003).
Article ADS CAS Google Scholar
Lamouroux, E., Serp, P., Kihn, Y. & Kalck, P. Identification of key parameters for the selective growth of single or double wall carbon nanotubes on FeMo/Al2O3 CVD catalysts. Appl. Catal. A Gen. 323, 162–173. https://doi.org/10.1016/j.apcata.2007.02.019 (2007).
Hart, A. J., Slocum, A. H. & Royer, L. Growth of conformal single-walled carbon nanotube films from Mo/Fe/Al2O3 deposited by electron beam evaporation. Carbon 44, 348–359. https://doi.org/10.1016/j.carbon.2005.07.008 (2006).
Chung, U. C., Lee, D. B., Jeong, Y. U., Ha, M. J. & Chung, W. S. Effect of H2 gas on carbon nanotubes synthesis. Mater. Sci. Forum 475–479, 3559–3562. https://doi.org/10.4028/www.scientific.net/MSF.475-479.3559 (2005).
He, Z. et al. Iron catalysts for the growth of carbon nanofibers: Fe, Fe3C or both?. Chem. Mater. 23, 5379–5387. https://doi.org/10.1021/cm202315j (2011).
Emmenegger, C. et al. Synthesis of carbon nanotubes over Fe catalyst on aluminium and suggested growth mechanism. Carbon 41, 539–547 (2003).
Kalaichelvan, K. in IEEE-International Conference on Advances In Engineering, Science And Management (ICAESM-2012). 429–433 (IEEE).
Markets, R.A. The Global Carbon Nanotubes (CNT) Market (2018–2023) is Projected to Grow at a CAGR of 16.7%—Technological Advancements and Decreasing Production Cost is Driving Growth. https://www.prnewswire.com/news-releases/the-global-carbon-nanotubes-cnt-market-2018-2023-is-projected-to-grow-at-a-cagr-of-16-7---technological-advancements-and-decreasing-production-cost-is-driving-growth-300752102.html (2018).
Sigma-Aldrich. 755710 Sigma-Aldrich carbon nanotube, single-walled >70% (TGA). https://www.sigmaaldrich.com/catalog/product/aldrich/755710?lang=en®ion=TH
Sigma-Aldrich. 755133 Sigma-Aldrich Carbon nanotube, multi-walled thin, <5% Metal Oxide(TGA). https://www.sigmaaldrich.com/catalog/product/aldrich/755133?lang=en®ion=TH#productDetailSafetyRelatedDocs
ResearchAndMarkets. Global Carbon Nanotubes Market Report 2020: Production Capacities for MWCNTS and SWCNTs, Historical and Forecast to 2030 - ResearchAndMarkets.com. https://www.businesswire.com/news/home/20201111005365/en/Global-Carbon-Nanotubes-Market-Report-2020-Production-Capacities-for-MWCNTS-and-SWCNTs-Historical-and-Forecast-to-2030---ResearchAndMarkets.com (2020).
Alijani, H. & Shariatinia, Z. Synthesis of high growth rate SWCNTs and their magnetite cobalt sulfide nanohybrid as super-adsorbent for mercury removal. Chem. Eng. Res. Des. 129, 132–149 (2018).
Chiang, W.-H. & Sankaran, R. M. The influence of bimetallic catalyst composition on single-walled carbon nanotube yield. Carbon 50, 1044–1050 (2012).
Liu, W.-W., Aziz, A., Chai, S.-P., Mohamed, A. R. & Hashim, U. Synthesis of single-walled carbon nanotubes: Effects of active metals, catalyst supports, and metal loading percentage. J. Nanomater. 2013 (2013).
Yadav, M. D., Patwardhan, A. W., Joshi, J. B. & Dasgupta, K. Selective synthesis of metallic and semi-conducting single-walled carbon nanotube by floating catalyst chemical vapour deposition. Diam. Relat. Mater. 97, 107432 (2019).
Chang, K.-H., Kim, J.-D. & Yoon, K.-W. Google Patents (2017).
Wang, H., Gu, G., Zhu, S., Yang, L. & Chen, Y. Synthesis of (9, 8) single-walled carbon nanotubes on CoSO4/SiO2 catalysts: The effect of Co mass loadings. Carbon 169, 288–296 (2020).
Cheung, C. L., Kurtz, A., Park, H. & Lieber, C. M. Diameter-controlled synthesis of carbon nanotubes. J. Phys. Chem. B 106, 2429–2433 (2002).
Tan, S., Ojovan, M. I., Hyatt, N. C. & Hand, R. J. MoO3 incorporation in magnesium aluminosilicate glasses. J. Nuclear Mater. 458, 335–342 (2015).
Article ADS CAS Google Scholar
Xu, X. et al. Controllable synthesis of carbon nanotubes by changing the Mo content in bimetallic Fe–Mo/MgO catalyst. Mater. Chem. Phys. 127, 379–384 (2011).
Fellicia, D. et al. Study of sigma phase in duplex SAF 2507. IOP Conf. Ser. Mater. Sci. Eng. 202, 012039. https://doi.org/10.1088/1757-899X/202/1/012039 (2017).
Das, R., Bee Abd Hamid, S., Eaqub Ali, M., Ramakrishna, S. & Yongzhi, W. Carbon nanotubes characterization by X-ray powder diffraction–a review. Curr. Nanosci. 11, 23–35 (2015).
Duc Vu Quyen, N., Quang Khieu, D., Tuyen, T. N., Xuan Tin, D. & Thi Hoang Diem, B. Carbon nanotubes: Synthesis via chemical vapour deposition without hydrogen, surface modification, and application. J. Chem. 2019 (2019).
Kashiwaya, Y. & Watanabe, M. Kinetic analysis of the decomposition reaction of CH4 injecting into molten slag. ISIJ Int. 52, 1394–1403 (2012).
Liu, Q., Wu, P., He, J., Liu, C. & Jiang, W. Catalytic decomposition of methane by two-step cascade catalytic process: Simultaneous production of hydrogen and carbon nanotubes. Chem. Eng. Res. Des. 163, 96–106 (2020).
Wang, J., Jin, L., Yu, G. & Hu, H. Effect of hydrogen addition on formation of hydrogen and carbon from methane decomposition over Ni/Al2O3. Can. J. Chem. Eng. 98, 536–543 (2020).
Wirth, C. T., Zhang, C., Zhong, G., Hofmann, S. & Robertson, J. Diffusion-and reaction-limited growth of carbon nanotube forests. ACS Nano 3, 3560–3566 (2009).
Article CAS PubMed Google Scholar
Pirard, S. L., Douven, S., Bossuot, C., Heyen, G. & Pirard, J.-P. A kinetic study of multi-walled carbon nanotube synthesis by catalytic chemical vapor deposition using a Fe–Co/Al2O3 catalyst. Carbon 45, 1167–1175 (2007).
Ni, L. et al. Kinetic study of carbon nanotube synthesis over Mo/Co/MgO catalysts. Carbon 44, 2265–2272 (2006).
Hsieh, C.-T., Lin, Y.-T., Chen, W.-Y. & Wei, J.-L. Parameter setting on growth of carbon nanotubes over transition metal/alumina catalysts in a fluidized bed reactor. Powder Technol. 192, 16–22 (2009).
Yadav, M. D., Dasgupta, K., Patwardhan, A. W., Kaushal, A. & Joshi, J. B. Kinetic study of single-walled carbon nanotube synthesis by thermocatalytic decomposition of methane using floating catalyst chemical vapour deposition. Chem. Eng. Sci. 196, 91–103 (2019).
Fogler, H. S. Essentials of chemical reaction engineering (Pearson Education, 2010).
Bartholomew, C. H. & Farrauto, R. J. Fundamentals of industrial catalytic processes. (John Wiley & Sons, 2011).
Hartmann, M., Machoke, A. G. & Schwieger, W. Catalytic test reactions for the evaluation of hierarchical zeolites. Chem. Soc. Rev. 45, 3313–3330 (2016).
Article CAS PubMed Google Scholar
This research work is financially supported by PTT Exploration and Production Public Company Limited (PTTEP). The authors are thankful for equipment support (Raman spectroscopy) from Mahidol University—Frontier Research Facility (MU-FRF). Authors would like to thank Nawapol Udpuay, Chawalit Takoon and Dr. Suwilai Chaveanghong, scientists at MU-FRF for their kind assistance in instrumental operation and technical support. In addition, scientists of MU-FRF, Nawapol Udpuay, Chawalit Takoon, and Dr. Suwilai Chaveanghong for their kind assistance in instrumental operation and technical supports.
Center of Excellence in Particle and Material Processing Technology, Department of Chemical Engineering, Faculty of Engineering, Chulalongkorn University, Bangkok, 10330, Thailand
Soottitantawats of Chawalkul Chotmunkhongs, Tawatchai Charinpanitkul & Apina
Department of Chemical Engineering, Faculty of Engineering, Mahidol University, Nakhon Pathom, 73170, Thailand
Sakhon Ratchahat & Weerawut Chaiwat
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
C.C.: Conceptualization, Investigation, Methodology, Data curation, Validation, and Writing—Original draft preparation. S.R.: Data curation, Validation, Supervisor, Resources, Funding acquisition, and Writing—Reviewing and Editing. W.C.: Resources, Validation, Funding acquisition, and Reviewing. T.C.: Resources, Funding acquisition and Validation. A.S.: Methodology, Validation, Visualization, Data curation, Supervision, Funding acquisition, and Writing—Reviewing and Editing.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Chotmunkhongsin, C., Ratchahat, S., Chaiwat, W. et al. Synthesis of MWCNTs by chemical vapor deposition of methane using FeMo/MgO catalyst: role of hydrogen and kinetic study. Sci Rep 13, 21027 (2023). https://doi.org/10.1038/s41598-023-48456-z
DOI: https://doi.org/10.1038/s41598-023-48456-z
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Zsm 5 Catalyst Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.