- Short Report
- Open Access
Clarifications regarding the use of model-fitting methods of kinetic analysis for determining the activation energy from a single non-isothermal curve
Chemistry Central Journal volume 7, Article number: 25 (2013)
This paper provides some clarifications regarding the use of model-fitting methods of kinetic analysis for estimating the activation energy of a process, in response to some results recently published in Chemistry Central journal.
The model fitting methods of Arrhenius and Savata are used to determine the activation energy of a single simulated curve. It is shown that most kinetic models correctly fit the data, each providing a different value for the activation energy. Therefore it is not really possible to determine the correct activation energy from a single non-isothermal curve. On the other hand, when a set of curves are recorded under different heating schedules are used, the correct kinetic parameters can be clearly discerned.
Here, it is shown that the activation energy and the kinetic model cannot be unambiguously determined from a single experimental curve recorded under non isothermal conditions. Thus, the use of a set of curves recorded under different heating schedules is mandatory if model-fitting methods are employed.
In a paper recently published in the Chemistry Central Journal, a kinetic study of the thermal decomposition of both aged and non-aged commercial cellulosic paper was presented, and the apparent activation energy (Ea) of the degradation reaction was determined for each case . According to the authors, the Ea of the process is related to the breakdown of cellulose chains and, since the apparent activation energy of the process is found to decrease with the aging time of the cellulose paper, it is proposed that such evolution could be used to construct archaeometric curves. Three different model-fitting methods were used to determine the activation energy: the differential Arrhenius method, the integral Savata method and the Wyden-Widmann method. Also, the authors speculate with the possibility that the kinetic method selected influences the obtained Ea. Actually, when using the Wyden-Widmann method, it is observed that only a limited number of data points around the DTG peak should be employed or else the Ea obtained would not fit that obtained by the other kinetic methods. Finally, it is concluded that a first order model is the most suitable for describing the cellulose decomposition reaction. However, recent works have reached to different conclusions, finding a chain scission model to be far more appropriate . Such discrepancy is due to some fundamental misconceptions in the manner the kinetic methods are employed in Marini’s work. Firstly, the activation energy for every sample studied was obtained by means of applying a model-fitting method to experimental data proceeding from a single non-isothermal run. Secondly, only the fit to two kinetic models, F1 and A2, were explored in the analysis. Basically, model-fitting methods of kinetic analysis consist of fitting the experimental data to a series of theoretical kinetic models, which are algebraic functions that reflect the relationship between reaction rate and degree of conversion and can be related to the reaction mechanism. The model providing the best linear fit is usually regarded as the correct one, and the activation energy is deduced from the slope of the fit. Unfortunately, it has been long established that the activation energy cannot be reliably determined from a single non isothermal curve because the experimental data almost always provides a reasonable fit regardless the kinetic model selected [3–5]. Despite that significant flaw, such inappropriate practice is still nevertheless widely used. As a result, it is common that nth order models are incorrectly selected because they are often tested as the first option for simplicity and a good fit is usually obtained. Here, we attempt to throw some light on the use of model-fitting methods and clarify such still widespread misuses.
Results and Discussion
Figure 1 includes a simulated α-T curve, constructed assuming a heating rate of 10 K min-1 and the following kinetic parameters: Ea=150 kJ/mol, A=1010 s-1 and a F1 (first order) kinetic model. The model-fitting methods of Arrhenius and Savata, those used in Marini’s work , were selected to determine the activation energy of the simulated curve. Thus, from the simulated data, the left hand of Eqs (1) and (2), corresponding to the Arrhenius and Savata methods respectively and shown in the Methods section, are plotted vs. the inverse of the temperature considering several of the most usual models in the literature. The f(α) and g(α) functions are listed in Table 1. The resulting activation energies, as obtained from the slope of the plots, and the regression coefficients showing the quality of the fits are included for the Arrhenius and Savata methods in Tables 2 and 3, respectively. Additionally, Figures 2 and 3 shows a selection of the plots resulting from the fitting of the data to the different kinetic models, so that the quality of the fit is clearly illustrated. Thus, when the simulated curve is analyzed by the Arrhenius method, six out of nine models deliver excellent fits to the data, with regression coefficients over 0.99. Therefore, it is not really possible to effectively discern the correct model with this procedure. Moreover, as it can be observed by the values in Table 2, the activation energy obtained from the analysis is highly dependent on the kinetic model assumed, with only the fit to F1 yielding the correct value. Consequently, without further evidence regarding the correct kinetic model, the activation energy cannot be established. The results are even more problematic when the Savata model is employed (Table 3 and Figure 3). Using such method the fit to all models tested are excellent and no clear candidate can be appropriately selected.
Consequently, neither the model nor the activation energy can be deduced from the use of a model-fitting method of kinetic analysis and a single non-isothermal curve. Note that this conclusion is reached after analyzing simulated, error-free data. When using experimental data it is probable that even a higher percentage of the tested models will adequately fit such data. Thus, in order to unambiguously determine the kinetic parameters by a model-fitting procedure, a set of experimental curves, each recorded under different heating schedules, must be employed [3, 6]. Then, a set of four curves simulated assuming the aforementioned kinetic parameters and heating rates of 1, 2, 10 and 20 K min-1 were analyzed simultaneously using the Arrhenius method. The resulting plots are shown in Figure 4. As it can be clearly noticed, only when the tested model is the right one, F1 in this case, are the plots positioned along a straight line. Therefore, using a set of curves recorded under different heating schedules, it is possible to unambiguously determine the kinetic model and, consequently, the activation energy.
As a final note, it should be considered that these methods are proposed for single step reactions, which can be described by a single kinetic triplet. When more than one process is taking place, each of them is expected to be defined by a different triplet. The thermogravimetric curves included in Marini’s work explicitly show a complex, multistep process and, therefore, such model-fitting methods cannot be employed. It is then recommended to resort to isoconversional methods or attempt the deconvolution of the contributing steps in order to study them independently [3, 7–9].
The uncertain results provided by the Wyden-Widmann method in Marini’s work can be likewise explained. It is still a model-fitting method since it assumes the process is driven by an nth order kinetic model. As it has been reported elsewhere, an nth order model is described by a mathematical function that cannot replicate the initial induction period typical of the function describing a chain scission model . Thus, being the wrong model to describe the reaction, it is understandable that the required linearity is not achieved along the entire temperature range, as described in the paper. Had the experimental data been fitted to the right model, such dependence of the activation energy on the number of data points considered would most probably have not been found.
It has been shown that the activation energy cannot be reliably determined by applying model-fitting methods of kinetic analysis to data obtained under non isothermal experimental conditions. Thus, the use of a set of curves recorded under different heating schedules is necessary, as recently recommended by the ICTAC Kinetics Committee . It is also important to consider the nature of the reaction under study since only one step processes can be analyzed by this methodology. More complex or multiple step reactions require the use of isoconversional methods or the deconvolution of the individual steps.
The simulated curves were constructed using a Runge–Kutta 4th order numerical integration method by means of the Mathcad software (Mathsoft, Needham, MA, USA).
The Arrhenius method is based in the following equation:
where dα/dt is the reaction rate, f(α) is the kinetic model, A the preexponential factor, E the activation energy and T the temperature in Kelvin. The Savata method relies in the following equation, which is obtained after integrating the expression above and reordenating terms:
where g(α) is the integral form of the kinetic model. For determining the activation energy, the left-hand side of Eqs (1) and (2) are plotted against the inverse of the temperature. The value of the activation energy is then deducted from the slope of such plot.
Marini F, Tomassetti M, Vecchio S: Detailed kinetic and chemometric study of the cellulose thermal breakdown in artificially aged and non aged commercial paper. Different methods for computing activation energy as an assessment model in archaeometric applications. Chem Cent J. 2012, 6 (Suppl 2): S7-10.1186/1752-153X-6-S2-S7.
Sanchez-Jimenez PE, Perez-Maqueda LA, Perejon A, Pascual-Cosp J, Benitez-Guerrero M, Criado JM: An improved model for the kinetic description of the thermal degradation of cellulose. Cellulose. 2011, 18 (6): 1487-1498. 10.1007/s10570-011-9602-3.
Vyazovkin S, Burnham AK, Criado JM, Perez-Maqueda LA, Popescu C, Sbirrazzuoli N: ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim Acta. 2011, 520 (1–2): 1-19.
Brown ME: Steps in a minefield - Some kinetic aspects of thermal analysis. J Therm Anal. 1997, 49 (1): 17-32. 10.1007/BF01987418.
Criado JM, Ortega A: Remarks on the discrimination of the kinetics of solid-state reactions from a single non-isothermal trace. J Therm Anal. 1984, 29 (6): 1225-1236.
Sanchez-Jimenez PE, Perez-Maqueda LA, Perejon A, Criado JM: Combined kinetic analysis of thermal degradation of polymeric materials under any thermal pathway. Polym Degrad Stab. 2009, 94 (11): 2079-2085. 10.1016/j.polymdegradstab.2009.07.006.
Perejon A, Sanchez-Jimenez PE, Criado JM, PerezMaqueda LA: Kinetic Analysis of Complex Solid State Reactions. A New Deconvolution Procedure. J Phys Chem B. 2011, 115: 1780-1791. 10.1021/jp110895z.
Vyazovkin S: Model-free kinetics - Staying free of multiplying entities without necessity. J Therm Anal Calorim. 2006, 83 (1): 45-51. 10.1007/s10973-005-7044-6.
Sanchez-Jimenez PE, Perejon A, Criado JM, Dianez MJ, Perez-Maqueda LA: Kinetic model for thermal dehydrochlorination of poly(vinyl chloride). Polymer. 2010, 51 (17): 3998-4007. 10.1016/j.polymer.2010.06.020.
Sanchez-Jimenez PE, Perez-Maqueda LA, Perejon A, Criado JM: A new model for the kinetic analysis of thermal degradation of polymers driven by random scission. Polym Degrad Stab. 2010, 95 (5): 733-739. 10.1016/j.polymdegradstab.2010.02.017.
Financial support from project CTQ2011-27626 (Spanish Ministerio de Economía y Competitividad) and FEDER funds is acknowledged. Additionally, one of the authors (PESJ) is supported by a JAE-Doc grant (CSIC-FSE). We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
The authors declare that they have no competing interests.
PSJ and AP performed the kinetic analysis of the simulated curves presented in this paper. The draft was prepared by PSJ and edited into final form by LPM. JC coordinated the study. All authors read and approved the final manuscript.
About this article
Cite this article
Sánchez-Jiménez, P.E., Pérez-Maqueda, L.A., Perejón, A. et al. Clarifications regarding the use of model-fitting methods of kinetic analysis for determining the activation energy from a single non-isothermal curve. Chemistry Central Journal 7, 25 (2013) doi:10.1186/1752-153X-7-25
- Activation Energy
- Kinetic Model
- Isoconversional Method
- Simulated Curve
- Heating Schedule