2.6. Thermal Properties of C. aesculifolia: TGA-DTG
The thermal behavior generated by the mass degradation of the lignocellulosic material under study was characterized using a thermogravimetric analyzer. In
Figure 3, the two vertical axes show the loss of mass (left axis) at a certain velocity (right axis). The latter represents the first derivative (DTG) of the thermogravimetric curve (TGA) throughout the temperature range. The thermogram reveals a behavior similar to that of diverse lignocellulosic materials [
54].
Three regions at distinct heating speeds are discernible in
Figure 3 (β = 10–30 °C/min; 283–303 K/min). The first, located at temperatures between 300 and 415 K, can be identified with the loss of surface water from the biomass. As the graph advances, it presents the greatest loss of mass (≈60 %) in region 2 (415–680 K). For lignocellulosic materials, this zone is distinguished by the elimination of volatile materials, especially the two main polysaccharides: hemicellulose and cellulose. This is normally the area where the chemical reactions that decompose these compounds and the greatest generation of gases (methane, hydrogen, nitrogen, carbon monoxide and dioxide) occur, so it is called the active stage of pyrolysis. Finally, in the range of 680–1100 K, the degradation of lignin—the main cementing compound in these materials—occurs. Region 3 also reveals the formation of fixed carbon. In this stage, both curves (TGA and DTG) are practically flat.
In addition to the mass degradation curves (TGA),
Figure 3 displays the characteristic peaks of the derivative of TGA—that is, the DTG that represents the ratio of change in weight loss with respect to temperature throughout the thermal process. Like TGA, DTG provides valuable information on key changes in the biomass, such as the maximum reaction temperatures (peaks) and decomposition velocities at any given point [
55].
Parallel to TGA, the two main peaks in DTG represent the decomposition velocities of water, polysaccharides (hemicellulose and cellulose), and lignin. It is interesting to note that although this biomass is made up of these three compounds, as well as water, only two principal peaks encompass them. As mentioned, the polysaccharides degrade at distinct temperatures. Reports generally indicate that, due to heterogeneity and the size of the fiber in hardwoods [
56], the decomposition of hemicellulose and cellulose occurs between 573 and 673 K [
57], while the range for lignin is wider: 300–1200 K [
58]. In this study, this summarizes the thermal behavior of
C. aesculifolia wood. Another aspect observable in
Figure 3 is the left-to-right shift of the TGA-DTG curves as the heating rate (10–30 °C/min; 283–303 K/min) and temperature increase. The explanation of this phenomenon is based on heat transfer inside the biomass, which is inversely proportional to the heating rate [
59].
2.7. Kinetic Analysis of C. aesculifolia Wood: Ea and A
Based on the data obtained from the thermogravimetric analysis and the methods of chemical kinetics proposed by Friedman, Flynn–Wall–Ozawa, Kissinger–Akahira–Sunose, and Kissinger, it was possible to determine the principal kinetic parameters—that is, activation energy (Ea) and the pre-exponential factor (A). The reaction order, meanwhile, was determined using Avrami’s method, as explained in the materials and methods section. These methods are extremely important for improving our understanding of the pyrolytic process of
C. aesculifolia wood in an inert atmosphere. Since each one uses distinct aspects to determine the reaction mechanism, associated kinetics, and regression techniques according to the handling of the data obtained from the thermal process, distinct interpretations may be reached for each result, as occurs with Ea. Because each mathematical method applied in reaction kinetics studies represents diverse mechanisms, the values for the kinetic parameters (Ea and A) will show certain differences [
60].
Figure 4a–d and
Table 4 show the results for these parameters according to the four methods used to describe interrelations among the degree of advance, temperature, and velocity. They also show the calculation of the correlation coefficient, R
2 (
Table 4). This is particularly important because, according to the value obtained—preferably between 0.90 and 1.0—any graph of the kinetic methods proposed is considered acceptable. In this study, the graph of
C. aesculifolia presented low R
2 values (<0.80) when the degree of conversion (α) was ˃0.70. These values were not included in the analysis. It is important to mention that some research related to kinetic studies (Friedman, FWO, and KAS) of biomass have reported similar behavior, in fact with a lower degree of progress (α = 0.65). This situation can be attributed to the degree of complexity in the generation of carbon in the pyrolysis process [
61]. Recent research has had to adjust its results to values where the degree of advance is also 0.70 due to poor fitting of the data—for example, in Trapa natans peel biomass [
62]. However, this does not mean that the results are not valid for the study in question.
It is important to emphasize that the average R
2 value for the four methods was ˃0.95, a figure practically equal to those of the FWO (0.9938) and KAS (0.9865) methods. We should clarify that Kissinger’s method does not use the degree of conversion as its basis for analysis, although it must also fulfill the condition that R
2 is ˃0.90 in order to achieve more approximate values [
63]. According to
Table 4, the average Ea of the Friedman, FWO, and KAS methods was 132.03, 121.65, and 118.14 kJ/mol, respectively, while the figure for Kissinger’s method was 155.85 kJ/mol. The first three results are quite similar, with KAS having the lowest value and Kissinger’s method the highest. These differences are due to the mathematical bases employed in each method. Friedman’s method is differential, based on the rate of change in the degree of advance with respect to time (dα/dt), while FWO and KAS are based on the heating rate (β), and both are integral methods. Thus, we must keep in mind that these kinetic methods are complementary and should not be understood as being in competition [
64]. Also important, however, is that Kissinger’s method reports only a value for Ea. Since it is not possible to appreciate the degree of advance of the pyrolytic process, the complexity of that process is not observed clearly. From the phenomenological point of view, one of the variables used by these methods (Friedman, FWO, KAS, and Kissinger) is temperature. Such differences in activation energy can also be explained as a function of increasing temperature and advancing degree of advancement, causing a reduction in the molecular mobility of the biomass structure, resulting in an increase in activation energy. As the process continues, the activation energy is reduced as the kinetics progresses [
65]. It can be observed that the value of the activation energy in the FWO method is slightly higher compared to the KAS method; this may be related to the behavior of the bonds in the biomass structure, the weakest being the ones that are first randomly broken by scission effects [
66].
Table 4 also shows that the maximum values for Ea (Friedman, KAS, FWO) are within the range of the degree of advance (α)—that is, between 0.40 and 0.55. This means that, based on this value (0.55), the amount of mass available to participate in the thermal process will be lower. These maximum values can also be interpreted as a function of the greater decomposition of the biomass that, according to the TGA-DTG graph (
Figure 3), corresponds to the elimination of the HM and CE that make up the wood of
C. aesculifolia.
Another significant aspect is the variability in Ea over the range considered for the methods applied: 81.89–153.10, 74.79–141.78, and 69.89–138.95 kJ/mol for the Friedman, FWO, and KAS methods, respectively. These results mean that the pyrolytic process is complex, occurs in multiple stages, and presents diverse types of chemical reactions across the temperature range [
67]. In general, reports in the literature affirm that the range of Ea for lignocellulosic biomass made up of woody tissue lies between 60 and 170 kJ/mol [
68]. Our results are within this range. Other biomasses present similar values for Ea: nutshell (≈136 kJ/mol) [
69], Indian almond (≈133 kJ/mol) [
70], Brazil nut (≈137 kJ/mol) [
71], products of figs (≈160 kJ/mol) [
72], fistula cane (Cassia fistula L.) (≈137 kJ/mol), peach palm (≈112 kJ/mol) [
73], Manilkara zapota seeds (≈132 kJ/mol), Delonix regia (≈143 kJ/mol), and Cascabela thevetia (≈152 kJ/mol) [
74].
Another key parameter in kinetic analyses of thermogravimetric processes is the reaction order (
n). The variation in this parameter for
C. aesculifolia is shown in
Figure 5 and
Table 5. This was calculated using Avrami’s equation, as indicated in the Materials and Methods section. Calculating the reaction order requires considering the heating rate (β = 10–30 °C/min; 283–303 K/min) and the temperature range of 500–650 K for the most important event of pyrolysis (region 2)—that is, the TGA-DTG analysis. Clearly, according to the lines in
Figure 5, the value of the lineal correlation coefficient (R
2) lies between 0.95 and 0.99, with an average of 0.9903, which is acceptable for calculating the reaction order.
As
Table 5 shows, the reaction order initially increased from 0.3937 to 0.6141, before decreasing to 0.3895, with an average value of 0.4887. This result is similar to those obtained by authors who have analyzed pyrolysis using agricultural waste and applying Avrami’s equation, as in the cases of corn straw and rice husk, for which average reaction orders of 0.365 and 0.539 were reported, respectively. It is true that this value can change from one type of biomass to another. These differences may be related to the composition of ash content, since this exerts a significant effect on the alkaline metals present and, in turn, their impact on the thermal process [
75].
2.8. Thermodynamic Analysis of C. aesculifolia: A, ΔH, ΔG, and ΔS
The first parameter obtained in thermodynamic analyses is the pre-exponential factor, calculated using the Ea values throughout the range of the degree of advance (α = 0.10–0.70) and the temperature at the maximum peak of the DTG (
Figure 3).
Table 4 shows the pre-exponential factor (A) obtained at a low heating rate (β = 15 °C/min). This reflects the fact that the frequency of molecular collisions increases proportionally to the heating rate [
76]. In this regard, applying a low value of β reduces the degree of interaction among the constituents of the biomass in pyrolysis. Recent research (on cattle manure and pistachio shells) has also opted for the calculation of the pre-exponential factor at even lower heating rates (10 °C/min) [
77,
78]. It is worth mentioning that, with respect to different regions of mass degradation, more accurate results can be obtained by performing thermogravimetric analysis at low heating rates, even as low as 1 °C/min [
79]. In this regard, and based on the information reported in the aforementioned investigations, thermodynamic analysis of
C. aesculifolia is performed at a low heating rate. In this study,
A varied in the range of 10
4–10
11 s
−1, defining the occurrence of both reactions and complex structures in this process due to the thermal conversional of the biomass [
60]. This variability in
A may be conditioned by other factors, such as particle size, the presence of a catalyzer, and even the diverse heating rate applied. It is indicative of the collisions that occur among the particles during the thermal process. Here, the number of collisions increased with the increase in this value. The average values shown in
Table 4 for the KAS, FWO, Friedman, and Kissinger methods are 2.41E + 09, 4.30E + 09, 8.11E + 10, and 3.47E + 11 s
−1, respectively.
Similar results have been reported where A varied from 10
7 to 10
12 s
−1—for example, with red pepper, rice husk, and bran [
80]. Other reports affirm that, depending on whether the value of A is <10
9 s
−1 or ≥10
9 s
−1, either collisions or superficial reactions and the formation of a simple compound will occur, respectively. However, if this value lies between 10
10 and 10
12 s
−1, there is a possibility of activating rotations of certain compounds that were passive at first [
81]. This may cause the size of the complex formed to increase (unimolecular reaction) or remain unchanged (monomolecular reaction) in relation to the interaction with its neighbors [
82]. In this case, regardless of the method applied, the constant values of the pre-exponential factor (10
10–10
11 s
−1) represent the zone where the main loss of mass occurs—that is, degradation of HM and CE when the degree of advance (α) lies between 0.40 and 0.55 throughout the pyrolytic process. The highest value for
A (3.47E + 11 s
−1) was calculated by Kissinger’s method, perhaps reflecting the fact that the rotation of the active–reactive compound pair remained unchanged during the thermal process [
83]. However, in light of the disparity with the results of the Friedman, FWO, and KAS methods, it is likely that Kissinger’s method presented this result because of the method’s mathematical basis, as explained previously.
In addition to the pre-exponential factor, and as part of the thermodynamic study of the pyrolysis process of
C. aesculifolia,
Table 6 shows the variation in three other key parameters that make it possible to analyze energy behavior in terms of spontaneity, conservation, equilibrium, and quality [
84], namely, analyses of enthalpy (ΔH, kJ/mol), Gibbs free energy (ΔG, kJ/mol), and entropy (ΔS, J/mol. K). As mentioned above, these parameters were calculated at a low heating rate.
Table 6 displays the average values for these thermodynamic properties (ΔH, ΔG, and ΔS) according to the Friedman (127.12, 174.89, −78.83), FWO (116.75, 175.30, −96.61), KAS (113.23, 175.46, −102.68), and Kissinger (150.95, 173.97, −37.98) kinetic methods.
In the process of lignocellulosic biomass pyrolysis, i.e., the generation of compounds such as solids, liquids, and gases, the enthalpy can be defined as the energy that will be required by the biomass for the generation of these products [
85].
Table 6 shows that, regardless of the mathematical method applied, enthalpy begins with low values, increases proportionally with the reaction order up to high values, and then decreases. This behavior revealed that the biomass contained compounds that require greater energy to achieve their total transformation. The positive value of enthalpy is another important factor, for it indicates that the thermal degradation process of
C. aesculifolia involves an endothermic reaction [
86]. Recently, similar enthalpy values have been obtained from the biomass of Manilkara zapota seeds, with reports of an average value of 137 kJ/mol [
75]. Other studies, one with almond husk and another with acorn pericarp [
87], have reported values close to the ones found in our work. Here, we would emphasize that the difference between activation energy (
Table 4) and enthalpy (
Table 6) is <5 kJ/mol, so we can affirm that the results for the degradation of
C. aesculifolia, the formation of activated compounds, and the conversion into other compounds will all be favorable [
88]. This means that the energy barrier to carry out the pyrolysis process is low, and the formation of pyrolytic products will be favorable.
In lignocellulosic biomass pyrolysis processes, the possibility of generating activated complexes is defined according to the total energy available in the thermal system, i.e., Gibbs free energy (ΔG) [
82]. As
Table 6 shows, this parameter is positive throughout the range of the degree of advance (α = 0.10–0.70), indicating that the process does not develop directly or automatically but, rather, is oxidative in nature and requires external energy to achieve the pyrolysis reaction [
89]. This need for additional energy input may be a disadvantage; however, it is important to highlight that, in pyrolytic biomass processes, the results of entropy (negative) and ΔG > ΔH (
Table 6) imply consideration that a small amount of energy is surplus as input to the thermal system [
90]. It is also clear that an approximate value of 175 kJ/mol is maintained in all three iso-conversional methods (Friedman, FWO, KAS), as well as for Kissinger’s non-iso-conversional method. Recent studies of corncobs and pine wood report approximate ΔG values of 173 kJ/mol [
91] and 180 kJ/mol [
86], respectively. These values are relatively high for lignocellulosic biomass, which means that pyrolysis of
C. aesculifolia consumes a large amount of energy.
Theoretically, the entropy (ΔS) of a system represents its molecular disorder and randomness [
76]; however, in a pyrolytic system, this function may represent the level of order or disorder of the carbon layers formed in the thermal process [
82]. Results for this factor are shown in
Table 6, where negative values indicate that the system will undergo only relatively small physicochemical changes; that is, the level of disorder in the thermal process of
C. aesculifolia is relatively low compared to some of its degraded products [
92]. The highest ΔS value was found for the KAS method (−102.68 J/mol. K), while the lowest was found for Kissinger’s method (−37.98 J/mol. K). Studies of biomass waste (red pepper) have reported such negative values ranging from −100 J/mol. K to almost −250 J/mol. K over a wide conversion range. [
93]. This variability in entropy may be related to the mathematical method applied and to the precise components of the biomass, which can interfere with the thermal process, such as the presence of certain inorganic and alkaline earth compounds, or even the production method of the biomass [
94]. However, it has also been determined that such negative entropy values (i.e., reduction in the randomness of the system) as those reported in this study for
C. aesculifolia can be generated due to physicochemical aging processes, which can lead to a higher thermodynamic equilibrium [
95].
2.9. Fourier-Transform Infrared Analysis (FT-IR)
The FT-IR spectra of raw
C. aesculifolia biomass are portrayed in
Figure 6. There, we can identify the characteristic signals that correspond to the various functional groups that constitute the structure of wood, including polysaccharides, lignin, and other low-weight molecular compounds. The absorption band visible at 3354 cm
−1 corresponds to the stretching vibrations of hydroxyl groups (OH) that are interlaced intermolecularly. These vibrations are related to the main constituents of wood, including cellulose, hemicellulose, lignin, and some proteins [
75]. In another study, this functional group was reported for the bark of
C. pentandra at 3418 cm
−1 [
96]. There are reports in the literature that this functional group is associated with the constitutional water present in the cell walls of wood [
97]. Another functional group identified was CH, as has occurred with other species of this genus, like the bark of C. pentandra, where a similar signal to that of
C. aesculifolia was observed at 2935 cm
−1 [
96]. According to the literature, the signal at 2911 cm
−1 is associated with symmetric and non-symmetric vibrations, suggesting the presence of groups of aliphatic chains (CH
2, CH
3) derived from the elemental structure of the raw material, such as cellulose, hemicellulose, and lignin [
98]. The CH functional group has been found in other studies of hardwood species in the range of 2938–2933 cm
−1 [
99]. Likewise, for species of this genus, there are reports of a vibration in the region around 2918 cm
−1 in the cellulose of the fruit of
C. pentandra [
100]. Published reports also affirm that the bands centered around 2911 cm
−1 have been associated with extractable substances located on the surface of the primary cell wall. A study of
C. aesculifolia seeds found these kinds of substances at 2918 cm
−1, while another found them in C. speciosa fruit (also a member of this genus). In that case, the signals at 2900 cm
−1 were attributed to lipids [
101,
102].
The signal at 1734 cm
−1 is associated with stretching vibrations of the C=O type, such as the aldehyde, ketone, and ester functional groups present in cellulose, hemicellulose, and lignin [
103]. This finding is similar to those of a study on
C. aesculifolia seeds, which reported a signal at 1736 cm
−1 [
102], and to those of reports on hardwoods that found this functional group between 1740 and 1730 cm
−1 [
104]. In the absorption band centered at 1594 cm
−1, there were stretching vibrations characteristic of double bonds (C=C) in the aromatic skeleton of lignin [
105]. A study of
C. aesculifolia seeds reported this group (C=C) at 1594 cm
−1, a result consistent with the value reported here [
102]. The functional groups CH
2 and CH
3 have been found in studies of hardwood species at around 1464–1375 cm
−1 [
104]. The signal at 1238 cm
−1 is also linked to those groups, associated mainly with the CE and lignin in
C. aesculifolia [
106]. This signal has also been associated with the stretching of the C-O bonds in the xylene and syringyl ring of lignin and hemicellulose [
107]. A study of cellulose using
C. pentandra fruit identified C-O at 1249 cm
−1 [
99].
The peaks at 1150 cm
−1 are assigned to the stretching vibration of the C-O-C bridge of the characteristic esters of cellulose and hemicellulose. One study of
C. aesculifolia seeds reported this bond at 1157 cm
−1 [
100]. The functional groups C-O-C, C-OH, and C
4-OH were identified in the absorption band at 1025 cm
−1, where mainly β-glucopyranose was located in the cellulose [
101]. Similar values have been obtained for other hardwood species, including Prosopis laevigata [
101]. Some authors affirm that the signal centered at 1025 cm
−1 corresponds to the symmetric stretching of the C-OH groups in lignin, cellulose, and hemicellulose [
108].
There are reports that vibrations around 930 cm
−1 correspond to glycosidic bonds in cellulose and hemicellulose [
109]. A study of cellulose in
C. pentandra fruit identified this at 900 cm
−1 [
100]. Likewise, the vibration around 725 cm
−1 can be attributed to the balance of the CH
2 group in crystalline cellulose I [
110]. A study of
C. aesculifolia seeds located this at 710 cm
−1 [
102]. These peaks are characteristic of native cotton [
111]. Finally, the signal found with values near 600 cm
−1 was thought to be due to the deformation of OH groups [
112]. A study of cellulose in
C. pentandra fruit identified this at 614 cm
−1 [
100].