Abstract

In a transition to a circular economy, second-generation biomass energy has come to the forefront. The present study is aimed at characterizing biochar and byproducts of the pyrolysis of star anise residue (ANI) in the N2 and CO2 atmospheres as well as the kinetics and optimal reaction mechanisms based on the Flynn–Wall–Ozawa and Coats-Redfern methods. The ANI pyrolysis involved three stages, with the first one (161.5–559.1°C) as the main phase. The activation energy was lower in the N2 atmosphere than in the CO2 atmosphere (179.44–190.17 kJ/mol). The primary volatile products generated during the ANI pyrolysis were small molecule products (H2O, CO2, CO, and CH4), organic acids, alcohols, and ketones. The atmosphere type exerted a minimal impact on the types of gases released, with the CO2 atmosphere increasing CO and CH4 emissions. The pyrolytic oil of ANI contained a variety of organic compounds, including alcohols, phenols, ketones, acids, sugars, and other nitrogen- and oxygen-containing cyclic compounds, with its predominant compounds being acids, esters, ketones, and sugars. The elevated temperature range of 300–700°C enhanced the charring degree of the ANI biochar. The biochar showed stronger aromaticity in the CO2 atmosphere but better granularity in the N2 atmosphere. This study introduced an innovative perspective by showcasing the potential of ANI as a promising biomass source for energy generation and underscored its abundance, sustainability, and applicability as a raw material in fragrance production. It also emphasized the significance of CO2-reuse technology as a means to mitigate CO2 emissions. The findings of this work offer a theoretical and practical basis for the comprehensive utilization and efficient disposal of star anise residues.

1. Introduction

The utilization of nonrenewable energy sources has resulted in adverse environmental consequences, such as climate change and acid deposition [13]. The locally widespread adoption and multiobjective optimization of renewable energy provide one of the most critical leverages toward the global-scale realization of carbon peaking and carbon neutrality. In a transition to a circular economy, second-generation biomass energy has become an increasingly popular substitute for petrochemicals across many countries owing to its sustainability, spatiotemporal abundance, diversity, noncompetition with limited productive lands, and versatility of byproducts [4]. Huang et al. [5] reported that the Chinese star anise (Illicium verum) covers 4,727 km2, more than 85% of its globally total land area, and the dry weight aboveground biomass (e.g., branches, leaves, and fruits) of star anise (ANI) in China was approximately 206.9 thousand tons (kt) in 2012 and 135 kt in 2015. In China, ANI is mainly used as a raw material for spice production, accounting for more than 90% of its use, with its total oil output of approximately 4.1 kt in 2015, mainly extracted via water vapor distillation. According to the extraction rate of 7.2% for its dried fruit, the distillation residues of ANI in Guangxi annually amount to 53.5 kt [6], thus indicating its vast amount and bioenergy potential.

The thermochemical conversion methods play a pivotal role in enhancing our understanding of the biomass gasification processes. They enable us to delve into the underlying reaction mechanisms, choose suitable gasifying agents, pinpoint optimal reaction temperatures, and consider the unique characteristics of biomass [7]. Valizadeh et al. [7] conducted the process optimization of gasification under diverse conditions toward the enhancement of hydrogen production and found that the incorporation of catalysts significantly boosted the gasification efficiency while curbing the formation of undesirable tar byproducts. As a crucial avenue for transforming biomass into valuable biofuels and bioenergy, studies on gasification and pyrolysis technologies carry substantial importance in overcoming their technical and economic obstacles as well as in facilitating their advancement, commercialization, and industrial implementation. Seo et al. [8] stated that the conversion of crop residues to biochar offers a dual benefit of serving as a direct solid fuel and a catalyst in the production of various biofuels. Also, the in situ or postmodification of biochar can enhance its catalytic performance by improving its pore structure, specific surface area, and surface functionality [8].

Currently, the comprehensive kinetics and gas-to-biochar byproducts of the ANI pyrolysis and their atmosphere dependency remain lacking, although many studies have investigated the biomass pyrolysis process. For example, Fu et al. [9] used thermogravimetric-Fourier transform infrared spectroscopy (TG-FTIR) analysis to reveal the temperature-dependent selectivity and stability of volatile products of the pyrolysis of coffee residues. Zhang et al. [10] found that the pyrolysis of cow dung between 600 and 1100°C occurred more vigorously in CO2 than in N2. Liu et al. [11] utilized the TG-derivative TG (DTG) data to elucidate the extent to which the atmosphere type affected the copyrolysis of durian shell and dyeing sludge. Wen et al. [12] determined via the Horowitz–Metzger method that the three-dimensional diffusion model best described the pyrolysis mechanism of incense sticks in N2 or CO2. By combining TG-DTG-FTIR analysis and the isoconversional methods, such as Starink, Flynn–Wall–Ozawa (FWO), and distributed activation energy model (DAEM), Huang et al. [13] pointed out that the activation energy was lower for the pyrolysis of aboveground water hyacinth (172.09–173.09 kJ/mol) than for that of belowground water hyacinth (230.11–232.06 kJ/mol). Torres-Sciancalepore et al. [14] found that the pyrolytic gaseous products of Rosa rugosa seed (RSW) were mainly CH4 and H2, while the phenolic compounds in its pyrolytic oil are of significance for certain industries.

The FWO method is a widely employed technique in thermodynamic and kinetic studies, entailing measuring changes in sample weight or specific properties at varying heating rates. By illustrating the relationship between sample mass loss and temperature, this method can extract essential kinetic parameters (e.g., activation energy) and reaction rates of pyrolytic degradation and their temperature dependency according to the Arrhenius formula. On the other hand, the Coats-Redfern (CR) method is another valuable technique extensively used in thermodynamic and kinetic studies, serving to characterize the kinetic parameters and thermal behavior and stability of materials. Fernandez et al. [15] found that the FWO method resulted in a more accurate estimation of activation energy than did the Starink, Friedman, and Kissinger-Akahira-Sunose (KAS) methods. They also pinpointed the best-fit pyrolytic mechanisms by employing the CR method [15]. In the present study, the FWO and CR methods were adopted to estimate the kinetic parameters and models of the ANI pyrolysis and validate the accuracy of the collected data.

Kim et al. [16] pointed out that the mixed H2-CO2 atmosphere utilized during the pyrolysis of oak trees expedited the thermal lysis of volatile organic compounds (VOCs) and also triggered reactions with VOCs. Lee et al. [17] stated that the lignin content of biomass was the main driver of its CO2 sensitivity. Cho et al. [18] emphasized that adjusting the atmospheric CO2 concentration could reduce the production of tar and control the ratio of CO to H2 in pyrolytic gases during the copyrolysis of subbituminous coal and lignocellulosic biomass. However, significant knowledge gaps still exist regarding the extent to which the CO2 atmosphere and its concentration influence the behavior of the biomass pyrolysis and its biochar characteristics. The CO2 atmosphere was previously reported to promote the formation of pyrolytic carbon in a certain temperature range but exert an opposite effect above that range [19]. Duan et al. [20] found that the coal pyrolysis in the CO2 atmosphere increased both volatile yield and rate of production. As the main product of industrial activities, such as thermal power generation, food fermentation, and cement manufacturing, CO2 that can be captured and reused can effectively reduce CO2 emission and help to achieve global carbon peaking and carbon neutrality. Currently, carbon capture and separation technology further facilitates the reuse of CO2 [21]. Therefore, this study selected N2 and CO2 to reveal the atmosphere dependency of the ANI pyrolysis and its dynamics.

Thus, this study is aimed at (1) investigating the performances and gas-to-biochar byproducts of the ANI pyrolysis in the N2 and CO2 atmospheres via TG-FTIR analysis and (2) quantifying their activation energy and optimal mechanism functions via the FWO and CR methods. The structure and key components of the present study are summarized in Figure 1. Figure 1 shows the experiment of anise in this paper and the analysis method of activation energy value and pyrolysis mechanism function. The novelties of this study lie in the elucidation of the atmosphere dependency of the ANI pyrolysis behavior and its effect on its biochar.

2. Materials and Methods

2.1. Preparation and Characterization of Star Anise Samples

The experiment was conducted with the commercially available star anise (ANI). Prior to the experiments, the ANI samples were dried at 105°C for 24 h, loaded into a sample bag, sealed, and stored in a dryer for later use. The ash (), moisture (), and volatile () contents of ANI were determined in accordance with China Solid Biofuels (GB/T28731-2012). Fixed carbon (FC) content was calculated using the following equation: . Sulfur (S), nitrogen (N), carbon (C), and hydrogen (H) were obtained by using an elemental analyzer (Vario EL cube by Elementar, Germany).

2.2. TG Analysis

Before the experiments, approximately 4.00 g samples were oven-dried at 65°C for 48 h. The TG and DTG curves were obtained by performing the TG experiments using a TG analyzer (NETZSCH STA 409 PC, Germany). In each TG experiment, a clean crucible was loaded with  mg of the sample and placed into the TG reactor, where N2 and CO2 were employed as carrier gases. Heating of the reactor was carried out at rates of 5, 10, 20, and 40°C/min from 30°C to 1000°C, with a 50 mL/min gas flow rate. To avoid nonrandom error, a blank baseline was acquired at each heating rate by running an empty alumina crucible under the same conditions.

2.3. TG-FTIR Analysis

The gas evolution during the biomass pyrolysis was analyzed by coupling a TG analyzer (NETZSCH STA 409 PC, Germany) and a Fourier transform infrared spectrometer (IS50 FTIR, Thermo, USA). In each experiment,  mg of sample was loaded into a TG reactor and heated to 1000°C, while the heating rate was 20°C/min, with a gas flow rate of 50 mL/min under a carrier gas of N2 or CO2.

2.4. Py-GC/MS Analysis

The pyrolysis experiments were conducted using a vertical microfurnace pyrolyzer (PY-3030D/7890B-5977A). Approximately 6 mg of ANI was placed into a quartz tube, which was then positioned on a sample tray. Once the temperature inside the reactor reached 600°C, the sample was introduced into the reactor for pyrolytic cracking and maintained at this temperature for 38 min in the helium atmosphere.

2.5. Preparation of ANI Biochar

Approximately  g of the sample tiled in the crucible was placed in a tubular furnace with N2 or CO2 as a carrier gas, heated at 10°C/min from 30°C to 300°C, 500°C, and 700°C for 1 h, taken out of the tubular furnace after being cooled down to room temperature, weighed, and used to estimate the harvest rate as follows: . The resulting samples were named according to the treatments of the atmosphere type and preparation temperature, namely, N2-300, N2-500, N2-700, CO2-300, CO2-500, and CO2-700.

2.6. Characteristic Performance Parameters

The pyrolysis characteristics of ANI were derived from the following six indices: initial volatilization temperature (), peak temperature (), average weight loss rate (), maximum pyrolysis rate (), total weight loss rate (), and comprehensive pyrolysis index (CPI). The CPI value was calculated as follows: where is the total weight loss and is the temperature difference of . A positive correlation exists between the pyrolysis performance and CPI value.

2.7. Thermodynamic and Kinetic Analyses

The main purpose of quantifying the pyrolysis kinetics is to explore the pyrolysis reaction mechanism and the relationship between the reactant structure and reaction capacity so as to control the reaction through the management of the operational settings. The thermal decomposition rate of solids can be expressed as follows:

The value of can be expressed using the Arrhenius equation as , which leads to the following:

For the nonisothermal heterogeneous systems, substituting the heating rate, , into Eq. (2) yields the following: where is the conversion degree (%) expressed as , where is the initial raw material mass (mg), is the raw material mass (mg) at a given time in the reaction, is the remaining solid mass (mg), time () in min, is the preexponential factor (s-1), (°C), is the gas constant (J/(mol·K), and is the apparent activation energy (J/mol).

Assuming that , the integral of Eq. (4) is as follows:

2.7.1. Flynn–Wall–Ozawa (FWO) Method

The FWO method can be used to calculate the reaction activation energy when the mechanism function is unknown, which can avoid the errors caused by the assumed mechanism function. Using the integral method: , the following can be written:

Equation (6) can be simplified as follows:

Since , Eq. (7) can be transformed into the following:

By measuring temperature, , and conversion degree, , at varying heating rates, , based on the isconversional methods, the activation energy value, , was derived from the slope of the best-fit regression line between versus .

2.7.2. The Coats-Redfern (CR) Method

The CR method is used to analyze the first-order kinetic parameters of thermal degradation and identify the most suitable reaction mechanism. As the most common model, the CR method was selected to compare with the activation energy value estimated via the FWO method to find the best reaction mechanism of the ANI pyrolysis. Assuming , its substitution in Eq. (5) yields the following:

The model that had the highest coefficient of determination () was chosen as the best-fit model of the pyrolysis mechanism, with the values of and corresponding to the slope and intercept of the best-fit line for the plot versus . The mechanism function model used in this study is shown in Table 1.

3. Results and Discussion

3.1. Physicochemical Driving Forces

Table 2 presents the results of the elemental and industrial analyses of the ANI pyrolysis. The elemental analysis adhered to standard procedures, involving the measurement of specific elements within the sample and the calculation of their respective percentages. It should be noted that while the reported elemental content may not sum up to 100% due to the potential presence of other unidentified elements or compounds in the sample, the analysis was conducted in line with well-established protocols and methodologies. There are technical challenges that can make the precise measurement of certain elements difficult, thus posing inherent limitations and uncertainties associated with the analytical process. Such factors as the sensitivity of the analysis method, losses during sample preparation, and operational errors can all influence the accuracy of the measurements. Also, the sample itself may contain trace elements or other components that are either not accounted for in the analysis or are challenging to detect. Thus, the total elemental content determined through analysis typically falls slightly below 100%.

The moisture content of ANI (4.74%) indicated its suitability for pyrolysis since biomass with moisture content below 10 wt% minimizes waste heat. Its low sulfur (S) content (0.16 wt%) also suggested a low risk of sulfur oxide (SOx) emission, whereas its relatively high nitrogen (N) content (1.22 wt%) highlighted the need to prevent nitrogen oxide (NOx) emission during the actual industrial production. The C/H ratio of 0.62 indicated a reduced likelihood of benzene derivative formation during the ANI pyrolysis as the higher ratios are typically more conducive to the production of such derivatives [25]. ANI displayed a high carbon (C) content of 44.32 wt% and oxygen (O) content of 40.83 wt%, signifying the presence of abundant organic compounds and enhanced thermal reactivity. Overall, the ANI pyrolysis exhibits substantial potential with the advantage of low SOx emission but a possibility of NOx emission, which should be carefully managed during industrial-scale production. Despite the uncertainties associated with them, these findings provide valuable insights into the pyrolysis process of ANI.

3.2. Analysis of the TG-DTG Curves
3.2.1. Atmosphere Dependency of the Pyrolysis of Star Anise

According to Figure 2, the pyrolysis of ANI in both atmospheres was observed to occur in the following three distinct stages: stage I (40–161.5°C): the water evaporation [2629]; stage II (161.5–559.1°C): the decomposition of volatiles; and stage III (559.1–1,000°C): the decomposition of residues, such as charcoal, minerals, and ash. In the N2 and CO2 atmospheres, the weight loss curves peaked at 318.5°C and 321°C, respectively. Stage I of the pyrolysis process corresponded to the initial loss of moisture in the sample. This moisture included both free water and bound water. Free water is present on the solid surface and within the large pores of the biomass material, whereas bound water is trapped within the solid matrix through processes such as adsorption. Free water tends to vaporize more readily than bound water [14]. However, both water types can completely evaporate at below 160°C [14]. The decompositions of hemicellulose, cellulose, and lignin were the main causes of mass loss at this stage. The hemicellulose structure is the simplest one of the three organic components, and hence, is first degraded, followed by cellulose and lignin [12]. Lignin, hemicellulose, and cellulose were previously reported to typically experience a significant decomposition between 160 and 900°C, 220 and 315°C, and 315 and 400°C, respectively [30]. The trends of the curves of phases I and II showed similar reactions between 30 and 559.1°C. In the N2 atmosphere in stage III, the DTG curve flattened out, and the TG curve exhibited a small mass loss probably due to the slow decomposition of residues, such as charcoal, minerals, and ash [10]. In contrast to the curve in the N2, the DTG curve in this stage showed one more peak in the CO2 atmosphere. The decomposition of CaCO3 was inhibited due to excessive CO2 [10]. Therefore, the DTG peak in the CO2 atmosphere was generated by carbon residue via the CO2 gasification, whose reaction type was thus , the Boudouard reaction, easy to occur at a high temperature, exhibiting an exponential relationship with temperature [31].

Overall, the curve showed an earlier decline in CO2 than in N2. Thereafter, between 307.7 and 397.7°C, the DTG curve of this stage was lower in the CO2 atmosphere than in the N2 atmosphere. This range was where the TG curve in the two atmospheres began to separate and dropped faster in the N2 atmosphere than in the CO2 atmosphere. This difference indicated that the ANI pyrolysis encountered greater resistance in the CO2 atmosphere, potentially due to the involvement of CO2 in the cross-linking reactions on the C surface, thus impeding the release of volatiles, a phenomenon cited by Sher et al. [32]. Between 397.7 and 690.3°C, the DTG curves for the two atmospheres were similar, while their TG curves showed similar trends.

Overall, the ANI pyrolysis increased the reaction rate slightly more in the CO2 atmosphere than in the N2 atmosphere in the initial stage (161.5–307.7°C). However, between 307.7 and 397.7°C, the reaction rate was higher in the N2 atmosphere than in the CO2 atmosphere. In other words, the CO2 atmosphere inhibited the precipitation of ANI in the pyrolysis.

3.2.2. Heating Rate Dependency of the ANI Pyrolysis

The ANI pyrolysis exhibited similar patterns in the two atmospheres; thus, the TG-DTG curves in the N2 atmosphere (Figure 3(a)) were discussed hereafter. The TG-DTG curves show that the mass loss regime remained consistent across all the heating rates, whereas the peak area changed. The four heating rates in the N2 atmosphere corresponded to the peak temperatures of 302.4, 311.9, 324.6, and 330.2°C and to the maximum degradation rates of 4.44, 8.44, 16.75, and 24.67 (%/min), respectively. The DTG curve moved to the right with the increasing heating rate due to the heat transfer lag between the internal and external particles of the sample. With the low heating rate, the small temperature difference caused slow but full volatilization [33]. However, the heat outside the sample could not be transferred to its internal part as the heating rate was high due to the action of thermal inertia [12]. Thus, the large temperature difference simultaneously released volatiles and degradation materials and led to the DTG peak with the increased heating rate in the high-temperature zone. Collard and Blin [34] pointed out that a higher heating rate promoted the formation of volatile compounds, while a lower heating rate was more conducive to the generation of carbon black, as consistent with the curves and trends in Figure 3.

The comparison of the TG-DTG curves of the N2 and CO2 atmospheres (Figures 3(a) and 3(b)) showed a similar change between 30 and 559.1°C; however, the maximum pyrolysis rate was lower in the CO2 atmosphere than in the N2 atmosphere (Table 4), as discussed in Sher et al. [32], possibly due to CO2 hindering the release of volatiles via cross-linking on the char surface. Between 559.1 and 1000°C, the ANI pyrolysis completely differed in the two atmospheres, with the N2 atmosphere exhibiting a slow and steady decline due to the decomposition of minerals and CaCO3 [35]. A large peak of weight loss occurred at 713.4°C in the CO2 atmosphere due to the Boudouard reaction . According to Lahijani et al. [36], the Boudouard reaction shows a thermodynamically enhanced effect at above 720°C which might result in the weight loss observed in the ANI pyrolysis.

3.3. Gas Products of the ANI Pyrolysis

This study focuses on discussing the infrared spectrogram at a heating rate of 20°C/min as the types of gas released during the biomass pyrolysis were reported to remain unaffected by the heating rate [37]. Table 3 shows the pyrolysis gases and their characteristic frequencies which can be used to detect the atmosphere dependency of CO and CH4 emissions.

As shown in Table 3, the CO2 atmosphere exhibited two characteristic frequencies since its infrared activity captured as a linear molecule involved asymmetric expansion, surface bending, and bending vibrations out of its four forms (symmetric and asymmetric expansion, surface bending, and bending vibrations). In Figure 4(a), the two peaks of in-plane and out-of-plane bending vibrations of CO2 (2500–2250 cm-1) were close together in the infrared spectrogram.

Figures 4(a) and 4(b) show the IR spectrum at the peak temperatures (326 and 320°C) in the two atmospheres. The data within the ranges of 2500–2250 cm-1 and 780–500 cm-1 were excluded from Figure 4(b) as they represented the IR spectrum recorded in the CO2 atmosphere. As shown in Figure 4(a), the peaks between 4000 and 3400 cm-1 indicated the existence of H2O, due to the transfer of crystallization and bound water in the gaseous state and the reorganization of the functional groups of alcohols and acids in the gaseous organic products due to the breakage of their chemical bonds. The peak between 3250 and 2850 cm-1 indicated the expansion vibration of C-H, and thus, the discharge of CH4 during the pyrolysis. The variation in CH4 with the temperature in the N2 atmosphere in Figure 4(d) showed three peaks, with the two peaks at 230.8 and 339.7°C. Ma et al. [40] showed that the pathway to the CH4 formation involved free radical reaction at below 500°C, whereas at >500°C, the degradations of methyl, benzyl, methylene, and methoxy groups were involved. Huang et al. [13] stated that CH4 produced at the low temperature (200–650°C) was attributed to the cleavage of the side chain of lignin. According to Chen et al. [41], a high amount of CH4 can be produced because of the presence of aromatic rings and numerous methoxy in the lignin structure. In other words, the aromatic ring and methoxy cleavage in lignin were the main drivers of the CH4 production by the ANI pyrolysis at the low temperature (200–500°C). The third peak at 690.3°C was low relative to the first two peaks. Wang et al. [42] demonstrated that the CH4 production at the high temperature (>500°C) is due to the rupture of the aromatic ring at >500°C. The remaining material contained a small amount of lignin between 600 and 800°C. The peak of CO2 between 2500 and 2250 cm-1 indicated the production of CO2 during the pyrolysis. The unique double absorption peak between 2250 and 2150 cm-1 conformed to the IR spectrum peak characteristic of CO released from the ANI pyrolysis. In Figure 4(c), the CO curve of the ANI pyrolysis in the N2 atmosphere showed only one peak at 326.5°C. At below 400°C, CO2 and CO emissions primarily stem from the breakdown and transformation of thermally unstable carboxyl, carbon, and ether groups in the phenylpropane side chain [42]. The telescopic vibration of C-O and C-H bonds induced the peak absorption between 1900 and 950 cm-1. Zhan et al. [43] found that hemicellulose and cellulose pyrolysis results in the production of alcohols, aldehydes, acids, and ketones. The CH4 peak in the high-temperature region in Figure 4(d) pointed to the presence of an aromatic ring in ANI and suggested the production of phenol, alcohol, aldehyde, and acids during the ANI pyrolysis when combined with Figure 4(a).

Figure 4 shows that the species of the pyrolysis gas in the two atmospheres were consistent; however, their release intensity differed. As far as CO and CH4 were concerned as an example in Figure 4(c), the temperature-dependent release of CO in the CO2 atmosphere showed a weak peak at 326.2°C due to the degradation of octagonal ether and carbon group. At above 600°C, the curve of CO began to rise and peaked at 925.4°C because of the reaction between carbon residue and CO2. As shown in Figure 4(d), CH4 showed a similar trend. Based on the Lambert-Beer law [44], the absorption strength of the gas at a specific wavelength is correlated with the concentration of the gas. Given the absorption strength of the gases in Figures 4(c) and 4(d), more CO and CH4 were released in the CO2 atmosphere than in the N2 atmosphere. Overall, the atmosphere type exerted a weak impact on the type of pyrolytic gas release at the peak temperature but a significant effect on its release concentration. The CO2 atmosphere increased CO and CH4 emissions from the ANI pyrolysis.

3.4. Compounds of the Pyrolytic ANI Tar

Table 4 provides a comprehensive overview of the pyrolytic ANI tar at 600°C along with their retention times and corresponding peak area ratios, as analyzed via Py-GC/MS. This analysis revealed a diverse array of organic compounds within the material (49 compounds), encompassing alcohols, phenols, ketones, acids, sugars, and various cyclic compounds containing N and O. The following compounds stood out due to their substantial peak area ratios, each exceeding 5%: acetic acid, oxo- (10.2%); acetic acid, oxo-, methyl ester (8.94%); acetic acid (7.29%); 2-propanone, 1-hydroxy- (5.18%); and β-D-glucopyranose, 1,6-anh (19.28%). Among these, acetic acid and acetic acid, oxo-, methyl ester are versatile compounds with applications as reactants or solvents in various industrial and biological processes. In particular, acetic acid is a widely used solvent and raw material in numerous industrial applications. β-D-Glucopyranose, 1,6-anh holds significant importance in the food, pharmaceutical, and bioengineering sectors and is used in the preparation of food additives, pharmaceuticals, dietary supplements, and biofuels. Overall, the pyrolytic ANI tar obtained demonstrated a considerable industrial and economic value, given the presence of these valuable compounds with diverse applications.

3.5. Comprehensive Pyrolysis Performance

The thermodynamic parameters facilitate the optimal design of pyrolysis reactors and operational conditions and their scale-up to large-scale real-time applications. The results of the comprehensive pyrolysis index for ANI in the two atmospheres are shown in Table 5. Due to the phenomenon of thermal hysteresis, as the heating rate increased, there was an increase in the initial volatilization temperature () and peak temperature () of the ANI pyrolysis. At the varying heating rates between 266.9 and 293.0 °C in the N2 atmosphere, the corresponding initial volatilization temperature () was consistent with that of the hemicellulose decomposition. The increased value indicated that the increased heating rate was conducive to the ANI pyrolysis. With the low heating rate, the volatiles can be precipitated slowly and fully [33], whereas with the high heating rate, the heat outside the sample could not be timely transferred to the inner parts [12]. The large temperature difference caused the simultaneous volatilization of various volatiles but not sufficiently.

Except for the total mass loss (), the trend of each parameter with the increasing heating rate was consistent in both atmospheres. As shown in Table 5, both and were lower in the CO2 atmosphere than in the N2 atmosphere at the same heating rates. In other words, the CO2 atmosphere weakly suppressed the emission of volatiles. At the same heating rates, the total mass loss was higher in the CO2 atmosphere than in the N2 atmosphere. As also evidenced in Figure 3, this behavior was due to the Boudouard reaction between 559.1 and 1000°C in the CO2 atmosphere, whereas the slow decompositions of charcoal, minerals, and ash residues mainly occurred in the N2 atmosphere. The increased total mass loss in the CO2 atmosphere with the rising heating rate was also due to the Boudouard reaction.

Overall, the increased heating rate exerted a limited promotion effect on the ANI pyrolysis. The CO2 atmosphere advanced the main reaction stage of the pyrolysis and inhibited the release of volatiles, consistent with the conclusion that the CO2 atmosphere inhibited the precipitation of volatiles in the ANI pyrolysis. Both the higher total mass loss and the comprehensive pyrolysis index were caused by the Boudouard reaction of the carbonaceous residue in the CO2 atmosphere.

3.6. Thermodynamic and Kinetic Systems

The activation energy () is the basic parameter to evaluate the difficulty of the reaction. The smaller the value is, the smaller the threshold value of energy for the reaction is. This study focused on the second stage according to the FWO and CR methods at 20°C/min.

3.6.1. Apparent Activation Energy

As mentioned earlier, this study utilized the FWO method to determine the kinetic parameters for the pyrolysis of ANI, based on insights of Fernandez et al. [45] and Torres-Sciancalepore et al. [14]. The FWO-estimated apparent activation energy for the ANI pyrolysis in both atmospheres is illustrated in Figure 5. The kinetic parameters of the ANI pyrolysis were quantified between 244.4 and 427.2 °C in the N2 atmosphere and 243 and 425°C in the CO2 atmosphere, with the conversion () degree of 0.2–0.8 and four heating rates (). In Figure 5, a highly linear correlation existed between and in the two atmospheres, with the values of 0.994–0.997 in the N2 atmosphere and 0.994–0.999 in the CO2 atmosphere (Table 6). The average activation energy when was estimated at 175.9 kJ/mol (N2) and 182.3 kJ/mol (CO2).

Figure 6 illustrates the changes in the values of the pyrolysis in both atmospheres as a function of the value of . As shown in Figure 6, the degradations of hemicellulose, cellulose, and lignin occurred in this stage most likely due to their different pyrolysis characteristics [46]. In the N2 atmosphere, the activation energy rose with the pyrolysis. The pyrolysis involved the breaking, recombination, and transformation of reactant molecules. The change in the activation energy during this process can be attributed to two main reasons. First, at the lower conversion degree, when the reaction began, the interactions between reactant molecules were relatively relaxed. This behavior lowered energy barriers to initiate and drive the reaction. As the conversion degree rose, the reaction progressed to a certain extent, and the interactions between reactant molecules became stronger. Consequently, more energy was needed to overcome these interactions and sustain the reaction. Therefore, as the conversion degree rose, the reaction grew more challenging, necessitating a corresponding increase in activation energy. Second, the increased conversion degree implied that more reactant molecules were already transformed into products, reducing the quantity of remaining reactant molecules. This reduction in reactant molecules decreased the reaction rate. To maintain the progress of the reaction, higher activation energy was required to accelerate the reaction rate. Therefore, at the higher conversion degree, the demand for activation energy further rose.

When the TG-DTG curves of the two atmospheres were compared, during the pyrolysis of octane in the CO2 atmosphere between 559.1 and 1000°C, the Boudouard reaction () occurred. Consequently, the quantity of reactant molecules became a significant factor influencing the change in the activation energy, as shown in Figure 6. Given in Figures 2, 3(b), and 6, the transition from hemicellulose to cellulose and to resistant lignin [47] in the N2 atmosphere was consistent with the trend of values estimated by Zhang et al. [10] through the FWO method.

Between 284 and 333°C when , the value was higher in the CO2 atmosphere than in the N2 atmosphere. The CO2 atmosphere exerted a certain negative impact on the main pyrolysis stage, as also confirmed by Zhu et al. [19]. The pyrolysis-starting temperature shifted to the left in the CO2 atmosphere compared with that in the N2 atmosphere, which was not caused by CO2 lowering the energy base. Overall, the decompositions of lignin, hemicellulose, and cellulose primarily drove the mass loss during the main pyrolysis stage in the N2 atmosphere, in line with the conclusion drawn in Section 3.2.1. The release of volatiles was hampered in the CO2 atmosphere. The lower initial temperature of the pyrolysis in the CO2 atmosphere than in the N2 atmosphere at the same heating rates in Table 5 did not result from the lowering of the pyrolysis energy barrier by CO2.

3.6.2. Reaction Mechanism

The CR method was used to fit the data between 244.4 and 427.2°C in the N2 atmosphere and 243 and 425°C in the CO2 atmosphere. By comparing the 18 solid reaction mechanisms, the optimal mechanism function of the ANI pyrolysis was selected. Figure 7 and Table 7 show the best-fit models for both atmospheres. The estimated via the F2.5, 2.5-order reaction was similar to the average value estimated via the FWO method, and they exhibited a significant linear relationship in both atmospheres, with an value of 0.997 in the N2 atmosphere and 0.998 in the CO2 atmosphere with their values of 183.27 kJ/mol and 179.59 kJ/mol, respectively.

The CR method is used to estimate apparent activation energy as an average value of a particular reaction phase, which may not accurately portray the entire complexity of the pyrolysis process. When used to estimate the kinetic parameters, the CR method may result in larger errors than the FWO method. To affirm the reliability of the reaction mechanism, the comparison of the measured and predicted data from the two atmospheres (Figure 8) led to , thus accurately predicting the ANI pyrolysis. The mechanism function of the second pyrolysis stage in both atmospheres was best described by the F2.5 model with the following form: .

3.7. Characterization of Pyrolytic ANI Biochar
3.7.1. Internal Microstructure

SEM analysis of the surface morphology of the biochar can provide a basis for a better understanding of the pyrolysis process and its biochar characteristics. Based on SEM and biofilm electrostatic test (BET) analyses, Kalin et al. [48] found that agriculturally stable biochar showed increased organic matter, alkaline response, and well-developed aromatic porous structure. Jia et al. [49] investigated the extent to which biochar influenced crop growth based on SEM, BET, and FTIR. The microstructures of the pyrolytic biochar as a function of the temperature in both atmospheres are illustrated in Figure 9. In the N2 atmosphere, the biochar showed a relatively smooth granular structure at 300°C (Figure 9(a)), angular and sharp edges, increased surface-attached particles at 500°C (Figure 9(b)), and a severely broken structure and increased roughness, with its internal cavities opening up in a bursting pattern at 700°C (Figure 9(c)). These results indicated that the biochar structure was strongly deformed between 300 and 500°C. The pore-like appearance of the biochar surface produced at 700°C (Figure 9(c)) was due to the intensified volatilization [50]. Figure 9(c-1) shows fine particles on the surface and increased surface heterogeneity. The surface micromorphology of the biochar did not significantly differ between the two atmospheres (Figures 9(d), 9(d-1), 9(e-1), 9(f), and 9(f-1)). At 300°C, 500°C, and 700°C, the C content of the biochar generated by the N2 atmosphere accounted for 52.39%, 70.88%, and 73.95%, respectively, while that of the biochar pyrolyzed in the CO2 atmosphere accounted for 63.99%, 72.76%, and 73.26%. The elevated pyrolysis temperature between 300 and 700°C increased the roughness of the biochar surface and its carbonization degree.

Compared with that at 300°C in the N2 atmosphere, the biochar produced in the CO2 atmosphere showed severe surface fragmentation, increased surface roughness, and smooth edges (Figures 9(a), 9(a-1), 9(d), and 9(d-1)). This result indicated the crosslinking effect of CO2 and the carbon surface [32]. The biochar (Figures 9(e) and 9(e-1)) was more dispersed from the adherent whole into several large particles at 500°C than at 300°C (Figures 9(d) and 9(d-1)). Its surface began to show fine particles and a few obvious burst-like cavities compared with that at 500°C in the N2 atmosphere (Figures 9(b) and 9(b-1)). Its structural characteristics were closer to those of the biochar at 700°C in the N2 atmosphere (Figures 9(c) and 9(c-1)). At 700°C, the biochar in both atmospheres was relatively similar. In other words, between 300 and 700°C, its surface structure at the same temperature changed more strongly in the CO2 atmosphere than in the N2 atmosphere; however, this difference weakened gradually. The biochar showed a higher degree of charring in the CO2 gradually than in the N2 gradually between 300 and 500°C but the opposite trend at 700°C. With the increased temperature, the C content of the biochar was lower in the CO2 atmosphere than in the N2 atmosphere. The particle size of the biochar was larger in the CO2 atmosphere than in the N2 atmosphere, also indicating the cross-linking effect of CO2 and the carbon surface [32]. Overall, between 300 and 700°C, the CO2 atmosphere promoted the fragmentation of the biochar structure and charring more than did the N2 atmosphere. However, the increased temperature affected the carbonization degree of ANI less in the CO2 atmosphere than in the N2 atmosphere. The biochar was more granular in the N2 atmosphere than in the CO2 atmosphere.

As shown in Table 8, increasing the preparation temperature of the ANI biochar positively influenced its BET surface area and pore volume between 300 and 500°C. At the same preparation temperature, the ANI biochar produced in the CO2 atmosphere exhibited a higher BET surface area and pore volume than that in the N2 atmosphere. This finding aligned with the observations made through the SEM images, providing strong evidence for the consistency and reliability of the experimental results. Overall, to create the ANI biochar with a high surface area and substantial pore volume, substituting CO2 for N2 as the pyrolysis atmosphere effectively reduced the energy input required by the pyrolysis temperature. This process suggested a potential energy-saving approach for the production of high-quality biochar.

3.7.2. Functional Groups

The FTIR characterization of the ANI samples prepared at 300, 500, and 700°C in the two atmospheres is shown in Figure 10. The spectral curve of N2-300 was significantly different from that of N2-500 and N2-700 but showed a high similarity with that of the original sample. Therefore, the pyrolysis degree of the biochar at 300°C was less in the N2 atmosphere than in the CO2 atmosphere. The chemical bonds of the peaks in response to the pyrolysis temperature and atmosphere for the functional characteristics of the biochar are presented in Figure 10. The curve in Figure 10 was analyzed by comparing the FTIR spectra in Table 9. The peak near 3400 cm-1 fell with the temperature, where the corresponding functional group was -OH group, indicating that the increased temperature raised the degree of the dehydration reaction during the pyrolysis. Between 3000 and 2700 cm-1, the corresponding stretching vibration of aliphatic C-H, CH2, and CH3 occurred. With the increased temperature, the peak intensity gradually rose, including the release of CH4 during the pyrolysis. The small peaks in this range may be caused by the degradation of cellulose [51]. The peaks between 1750 and 1700 cm-1 weakened with the increased temperature, which corresponded to the stretching vibration of the C=O bond, and almost disappeared at 500°C. The bands due to the aromatic C=C bonding vibrations were clearly observed by the peaks between 1670 and 1600 cm-1, indicating that the biochar in the two atmospheres contained the aromatic C=C bond. Finally, the peak between 1000 and 650 cm-1 was aromatic C–H with deformed vibration. Given the spectral curves of N2-300, N2-500, and N2-700, the peak of C=C and the aromatic property of the biochar gradually weakened with the increased temperature. The spectral curves of N2-300, N2-500, and N2-700 flattened with the increased temperature, hydroxyl (-OH) and base (C=O) disappeared, and the biochar was graphitized at 700°C in the N2 atmosphere [52].

Except for the aromatic C=C and C-H bonds, with the increasing temperature, the functional groups showed consistent trends between both atmospheres. The peak intensity of the aromatic C=C bond of the biochar in the CO2 atmosphere did not change with the increased temperature but was higher than that in the N2 atmosphere at the same temperature. The peak of the aromatic C-H bond showed an upward trend with the temperature rise, indicating that the biochar was more aromatic in the CO2 atmosphere than in the N2 atmosphere. The comparison of the curves of N2-500 and CO2-500 showed that the influence of the atmosphere on the biochar was mainly reflected in its aromatic properties. Overall, all the functional groups of the ANI biochar decreased with the increased temperature in both atmospheres, except for the aromatic C=C bond, which did not change with the increased temperature, and the aromatic C-H bond, which increased with the increased temperature.

4. Conclusion

The star anise pyrolysis was triphasic, with the primary decomposition occurring between 161.5 and 559.1°C. The increasing heating rate resulted in a thermal hysteresis effect on the mass loss curve regardless of the atmosphere type. The lower value of in the CO2 atmosphere than in the N2 atmosphere indicated that CO2 was highly conducive to the response of the biochar pyrolysis. The higher value of in the N2 atmosphere than in the CO2 atmosphere showed that N2 was highly conducive to the release of volatiles from the ANI pyrolysis. The two atmospheres affected the release intensity of the pyrolysis gas, with the CO2 atmosphere increasing CO and CH4 emissions. The values of the main reaction phase were estimated at 175.93 kJ/mol in the N2 atmosphere and 182.39 kJ/mol in the CO2 atmosphere, whose mechanism functions belonged to the F2.5 model in the form of . The temperature rises between 300 and 700°C and increases the charcoal degree of the biochar and the roughness of its surface. The biochar showed strong aromatic property and integrity in the CO2 atmosphere, whereas the biochar exhibited enhanced particle property in the N2 atmosphere. The pyrolysis oil derived from ANI was composed of a diverse array of organic compounds, encompassing alcohols, phenols, ketones, acids, sugars, and other cyclic compounds containing N and O. Among these compounds, the dominant categories found in the oil were acids, esters, ketones, and sugars.

Nomenclature

:Temperature difference of
:Preexponential factor
ANI:Anise
ANIB:Anise biochar
:Ash
BET:Brunauer-Emmett-Teller
C:Carbon
CPI:Comprehensive pyrolysis index
CR:Coats-Redfern
DAEM:Distributed activation energy model
DTG:Differential thermogravimetry
:Apparent activation energy
FC:Fixed carbon
FTIR:Fourier transform infrared spectroscopy
FWO:Flynn–Wall–Ozawa
H:Hydrogen
HHV:Higher heating value
IR:Infrared
:
KAS:Kissinger-Akahira-Sunose
:Moisture
:Total weight loss rate
N:Nitrogen
NOx:Nitrogen oxide
O:Oxygen
Py-GCMS:Pyrolysis-gas chromatography-mass spectrometry
:Gas constant
:Coefficient of determination
:Maximum pyrolysis rate
RSW:Rosa rugosa seed
:Average weight loss rate
S:Sulfur
SEM:Scanning electron microscope
SOx:Sulfur oxide
:Time
:, thermodynamic temperature
:Celsius temperature
TG:Thermogravimetric analysis
TG-DTG:Thermogravimetric differential thermal analysis
TG-FTIR:Thermogravimetric-Fourier transform infrared spectroscopy
:Initial volatilization temperature
:Peak temperature
:Volatile
VOCs:Volatile organic compounds
:Initial material mass
:Remaining solid mass
:Raw material mass at a given time in reaction
:Conversion degree
:Heating rate.

Data Availability

Data is available on request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Authors’ Contributions

Hanlin Cao was responsible for the investigation, software, validation, formal analysis, and writing. Qingbao Luo was responsible for the methodology, software, formal analysis, writing—original draft, and writing—review and editing. Yao He was responsible for the software and data curation. Jiayu Liang was responsible for the software, formal analysis, and data curation. Sen Lin was responsible for the conceptualization, data curation, and writing—review and editing. Fatih Evrendilek was responsible for the writing—review and editing. Jianying Qi was responsible for the software and data curation. Zuoyi Yang was responsible for the formal analysis and funding acquisition. Sheng Zhong was responsible for the formal analysis and project administration. Jingyong Liu was responsible for the conceptualization, resources, methodology, project administration, funding acquisition, and writing review. Chunxiao Yang was responsible for the experiment. Wenxiao Huang was responsible for the formal analysis and funding acquisition.

Acknowledgments

The authors acknowledge Ms. Yang from the Analysis and Testing Center of Guangdong University of Technology for her help with TG-FTIR analysis. This work was supported by the Open Foundation of Key Laboratory of Radioactive and Rare Scattered Minerals, Ministry of Natural Resources (2022-RRSM-01).

Supplementary Materials

Graphical Abstract: the research content and the research scheme used in this paper are expressed in the form of pictures. (Supplementary Materials)