Next Article in Journal
Molecular Informatics Studies of the Iron-Dependent Regulator (ideR) Reveal Potential Novel Anti-Mycobacterium ulcerans Natural Product-Derived Compounds
Next Article in Special Issue
Thermal Decomposition of Maya Blue: Extraction of Indigo Thermal Decomposition Steps from a Multistep Heterogeneous Reaction Using a Kinetic Deconvolution Analysis
Previous Article in Journal
Bridged Nucleic Acids Reloaded
Previous Article in Special Issue
Continuous Monitoring of Shelf Lives of Materials by Application of Data Loggers with Implemented Kinetic Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Critical Appraisal of Kinetic Calculation Methods Applied to Overlapping Multistep Reactions

1
Energetic Materials Laboratory, Semenov Institute of Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Russia
2
Chemistry Laboratory, Department of Science Education, Graduate School of Education, Hiroshima University, 1-1-1 Kagamiyama, Higashi-Hiroshima 739-8524, Japan
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(12), 2298; https://doi.org/10.3390/molecules24122298
Submission received: 23 May 2019 / Revised: 18 June 2019 / Accepted: 19 June 2019 / Published: 21 June 2019
(This article belongs to the Special Issue Thermal Analysis Kinetics for Understanding Materials Behavior)

Abstract

:
Thermal decomposition of solids often includes simultaneous occurrence of the overlapping processes with unequal contributions in a specific physical quantity variation during the overall reaction (e.g., the opposite effects of decomposition and evaporation on the caloric signal). Kinetic analysis for such reactions is not a straightforward, while the applicability of common kinetic calculation methods to the particular complex processes has to be justified. This study focused on the critical analysis of the available kinetic calculation methods applied to the mathematically simulated thermogravimetry (TG) and differential scanning calorimetry (DSC) data. Comparing the calculated kinetic parameters with true kinetic parameters (used to simulate the thermoanalytical curves), some caveats in the application of the Kissinger, isoconversional Friedman, Vyazovkin and Flynn–Wall–Ozawa methods, mathematical and kinetic deconvolution approaches and formal kinetic description were highlighted. The model-fitting approach using simultaneously TG and DSC data was found to be the most useful for the complex processes assumed in the study.

1. Introduction

Thermally-induced transformations in heterogeneous system usually do not obey the idealized single-step kinetic pattern but are comprised of consecutive or concurrent reaction steps. The reaction behavior is further complicated by the contribution of some physical phenomena (melting, diffusion, etc.) [1,2,3]. Thermal decomposition of solids represents a typical example of such complex process [4,5,6,7,8]. Determination of the kinetic parameters for the thermal decomposition is not a straightforward, but large experimental efforts should be paid for finding a possible way to attain the rigorous solution [9,10,11]. Alongside, examination of the applicability of different kinetic calculation methods through analyzing the mathematically simulated kinetic curves assuming the specific complex process is the possible method to obtain a guideline for the successful kinetic calculations [12,13,14]. For example, Svoboda et al. [15] analyzed the simulated process with independent reactions and revealed that the commonly used Kissinger method [16,17] provides the good estimate for the activation energy of the dominant reaction. Vyazovkin et al. applied isoconversional methods to the simulated parallel independent [18] and consecutive [19] reactions, and summarized the typical shapes of the activation energy dependency on the conversion degree [20]. Some examples of the complex reactions where the isoconversional hypothesis is not fulfilled were also reported [21,22,23]. Burnham [24] proposed the formal kinetic methods to describe the complex processes, but usually some preliminary insights are necessary to give the close initial guesses for many fitted parameters. The attempt to consider simultaneously various kinetic calculation methods to underline their advantages and shortcomings was performed by the Kinetic committee of ICTAC (International Confederation for Thermal Analysis and Calorimetry) [25]. Since then several thermokinetic techniques have been proposed [26,27], however, the critical assessment of the applicability of the modern thermokinetic methods for analyzing the complex processes is still missing.
If a partially overlapping multistep reaction is considered, the contributions of each component reaction step to the overall process determined by the different thermoanalytical (TA) techniques are not equal in general, because of the different physical quantities subjected to individual TA measurements. In some cases, the respective reaction steps can bring the opposite contributions to the overall process as measured using a TA technique. Experimental studies concerning the overall reaction composed of more than one reaction step with oppositely signed TA signals have been reported, i.e., partially overlapping mass loss and mass gain [3,28,29], or exothermic and endothermic events [30,31,32]. In practice, such complex processes are observed during thermal decomposition of energetic materials [5] and high-temperature operations with ionic liquids [33]. This particular situation, e.g., two overlapping mass-loss processes corresponding to endothermic and exothermic effects, respectively, educe that the apparent kinetic curves recorded using differential scanning calorimetry (DSC) can be different from that obtained using thermogravimetry (TG). The phenomenon is caused by different nature of these TA signals and was observed experimentally for isopropylammonium nitrate decomposition [30] and kerogen pyrolysis [34]. A sophisticated kinetic deconvolution technique have been applied recently to this type of solid-state reactions, where the kinetic parameters obtained for the constituted single processes have been linked with the physicochemical and physicogeometrical features of the transformations [31,32].
The aim of the present study is to evaluate the applicability of the available kinetic calculation methods to the simulated complex process characterized as the partially overlapping multistep reaction and to provide a guideline of successful kinetic analysis to the processes of this type. The kinetic curves were simulated mathematically (thus, correct kinetic parameters are available for verification), and processed using different kinetic calculation methods including Kissinger method [16], isoconversional methods [12,35,36,37], mathematical [27] and kinetic [31,38] deconvolution methods, and formal kinetic analysis [39,40]. Through the kinetic calculation using different methods, the resulting apparent kinetic parameters determined for simulated TG and DSC curves are mutually compared for finding the possible way to obtain relevant kinetic information of the real-life complex reaction process, for which one has limited information a priori.

2. Theoretical

2.1. Simulation of Successive Reactions

The thermoanalytical data were simulated for the following theoretical process comprising two successive reaction steps:
A ( s ) k 1 B ( s ) + C ( g ) ,   B ( s ) k 2 D ( s ) + E ( g )
To model the existence ratio of solid components – A(s), B(s), and D(s) at a time, a set of differential kinetic equations has to be solved. When TA signals under linear nonisothermal conditions at a heating rate β are simulated, a set of kinetic equations is expressed using conversion degrees, α1 and α2, for the respective component reaction steps [41]:
d α 1 d T = A 1 β e E a 1 / ( R T ) f 1 ( α 1 ) ,   d α 2 d T = A 2 β e E a 2 / ( R T ) f 2 ( α 1 , α 2 ) ,  
where A and Ea are the apparent Arrhenius parameters and f(α) is the kinetic model function. The subscripts 1 or 2 indicate the first and second component reaction steps. Considering the derivation of the kinetic rate data from derivative TG (DTG), the normalized overall reaction rate, (dα/dt)DTG, is expressed:
( d α d t ) DTG = η d α 1 d t + ( 1 η ) d α 2 d t ,  
where α is the conversion degree for the overall reaction and η is the ratio of the total mass-loss value during the first reaction process, Δm1, to the total mass-loss for the overall reaction, ΔmΣ: η = Δm1mΣ. On the other hand, the normalized overall reaction rate derived from DSC, (dα/dt)DSC, is represented as follows [19]:
( d α d t ) DSC = γ d α 1 d t + ( 1 γ ) d α 2 d t ,  
where γ is the ratio of the total thermal effect for the first reaction, Q1, to the total thermal effect for the overall reaction, QΣ: γ = Q1/QΣ. It must be noted from Equations (3) and (4) that DTG and DSC curves for the overall process are generally different because of different definitions of η and γ.
For simplicity, a series of kinetic curves were simulated mathematically by assuming a successive two first-order reactions; therefore, f1(α1) = 1 − α1 and f2(α1, α2) = α1α2. In addition, equivalent total mass loss for the first and second reaction processes was assumed, thus η = 0.5. The thermal effects of the component reaction processes were assumed to have opposite signs, and the absolute value of the total thermal effect of the first reaction process (exothermic) was set to be twice as large as that for the second reaction process (endothermic), i.e., γ = 2. Based on the aforementioned assumptions, four different cases of the temperature dependences of the rate constants of each component reaction step were simulated:
  • Case 1: the rate constant for the first reaction, k1, is larger than that of the second reaction, k2, and this difference decreases with increasing temperature, T;
  • Case 2: the value of k1 is larger than k2, and its ratio, k1/k2, is constant independent of T;
  • Case 3: the value of k1 is smaller than k2, and this difference increases with T;
  • Case 4: the value of k1 is smaller than k2, and k1/k2 is constant independent of T.
Variation of k1/k2 with temperature is plotted for cases 1–4 in Figure S1, and the apparent Arrhenius parameters assumed for simulating the kinetic curves are listed in Table 1. The simulation of the kinetic curves at different β values were performed using Mathcad (ver. 15.0), in which the fourth-order Senum and Yang approximation [42] was employed for the approximation of exponential temperature integral.

2.2. Kinetic Calculation Methods

Kinetic calculations for the mathematically simulated kinetic curves were performed using self-developed thermokinetic code THINKS [43], which comprises various calculation methods including Kissinger (ASTM E698 [17]), Friedman [35], Flynn–Wall–Ozawa (FWO) [36], Starink [44], and Vyazovkin [45] methods, combined kinetic analysis (CKA) [26], and model-fitting approaches, kinetic deconvolution [31,38] and formal kinetic analysis [39,40]. For analyzing the partially overlapping thermal events, the statistical peak separation technique, i.e., mathematical deconvolution [27], was applied. The fundamentals of these kinetic calculation methods can be found in original papers, the basic ideas relevant to the present study are implemented in the text.

3. Results and Discussion

3.1. Apparent Features of the Simulated Kinetic Datas

The simulated kinetic curves at β = 2 K min−1 for the Cases 1–4 are shown in Figure 1. All the simulated kinetic curves at different β values for each case are presented in Figures S2–S5. Although various patterns of overlapping two mass-loss processes with exothermic and endothermic thermal effects are expected, the kinetic curves simulated by assuming four different cases in this study are believed to cover typical apparent features of the kinetic curves possibly obtained by TG and DSC measurements.
In Case 1 (Figure 1a), the kinetic curve of DTG indicates a single peak, whereas overlapping two positive peaks with distinguishable peak tops are observed. At first glance, the process is misinterpreted as a single-step reaction from TG/DTG and as overlapping two exothermic reactions from DSC. However, the kinetic features assumed to simulate the kinetic rate can be deduced by comparing the kinetic curves derived from DTG and DSC: the peak top of DTG and the minimum between the two apparent positive peak tops in DSC appear at the same temperature. The assumed overlapping feature of the complex process is easily deduced in Case 2 (Figure 1b): DTG indicates overlapping of two peaks with well separated peak tops, while DSC indicates one positive and one negative peak. In Case 3 (Figure 1c), the smooth single peak is present on DTG curve. However, because of clearly separated positive and negative peaks on the kinetic curve of DSC, the assumed overlapping feature of the process is easily recognized. Case 4 is the most difficult one for deducing the complex features of the process from DTG and DSC data (Figure 1d). Because the kinetic curves derived from DTG and DSC are practically identical, it is the case usually interpreted from TG and DSC measurements as the single reaction process. In the present case, the situation appears as the result of the practically identical reaction temperature region and the same rate behavior assumed; however, similar situation is also expected when the contribution of a component reaction to the overall mass loss and thermal effect is very limited.

3.2. Kissinger Method

For the kinetic analysis of the single-step reaction, the peak maximum of DTG and DSC curves, i.e., point where (d2α/dt2) = 0, contain important kinetic information [25] as formalized in the Kissinger method [16]: Ea value is evaluated from the shift of the peak temperature with varying β value. Because of its simplicity and robustness, the Kissinger method is one of the most widely used kinetic calculation methods for the determination of Ea value and also is implemented in ASTM E698-05 [17]. In general, application of the Kissinger method to the complex process as assumed in this study lacks theoretical validity. However, as demonstrated by Svoboda and Málek [15] for the kinetic analysis of the independent overlapping reactions, the results of Kissinger technique often reveal the relevant Ea value for the rate-limiting step.
Table 1 lists the apparent Arrhenius parameters for simulated kinetic data, calculated from the Kissinger plots and assuming the first-order reaction behavior. Irrespective of DTG and DSC data, superficial Arrhenius parameters were obtained for Case 1. Herein, the peak maximums cannot be correlated directly to the constituting single reactions as assumed by Kissinger method, but are determined by the superposition of two processes, which relative contributions are varied with β (Figure 1a, Figure S2). Table 1 also lists the variation range of the conversion degree at peak maximum, Δαp, among the curves with different β. This variation should be negligible for rigorous application of Kissinger method [25], so its magnitude can serve as an indicator of the reaction complexity. In Case 2, both DTG and DSC data indicate well separated two peak maxima (Figure 1b) because of limited overlapping of two constituent reactions. The limited influence on each peak maximum from another stage results in nearly constant αp with β, i.e., negligible Δαp (Table 1). Therefore, the Arrhenius parameters determined by Kissinger method for the first and second peaks of DTG and DSC data fully agree with the correct kinetic parameters. Application of Kissinger method in Case 3 brings the incorrect A, Ea values due to high overlapping of stages, that is recognized by a considerable Δαp values (Table 1). As was seen in Figure 1d, Case 4 is very specific with the same Ea values for stages, but lower rate constant k1 throughout the process. The Δαp values for DSC and DTG data are negligible supporting the value of Kissinger analysis. As a result, the Ea values obtained from DTG and DSC data closely correspond to that assumed to simulate the kinetic data. The preexponent value calculated from the overall DTG and DSC data correspond to the assumed A value for the first reaction step, because the first reaction step is the rate-limiting step in Case 4.
Summarizing the above findings, the Kissinger method becomes less applicable with increasing the interrelationship between the constituent processes. One should look on the magnitude of peak shift Δαp as a first indicator of the reaction complexity.

3.3. Isoconversional Methods

The kinetic analysis of the TA data using the isoconversional methods is recommended by the kinetic committee of ICTAC [25]. For the ideal single-step reaction recorded under ideal measurement conditions, a constant Ea value should be obtained at different α during the course of reaction, irrespective of the kind of isoconversional methods, i.e., integral method of Flynn-Wall-Ozawa [36], differential method of Friedman [35], and advanced integral method of Vyazovkin [45]. However, in practice, the constant Ea is often not attained even for the relatively simple processes, due to the gradual changes in reaction conditions, reaction geometry, and rate-limiting step as the reaction advances. Therefore, along with the examples of successful applications to complex processes [37], some superficial Ea values and observations of the kinetic parameters variation in an exotic manner during reaction course have been reported [22]. Moreover, when the isoconversional methods are applied to the multistep processes, the physical meaning of conversion degree α should be carefully reconsidered: while, in the fundamental kinetic equation, α is defined as the fractional reaction for the single reaction step [25], for the multistep process the α value is determined experimentally as the fraction of the total changes in the physical quantities during the overall process. Keeping in mind these theoretical limitations, several isoconversional methods were applied to the simulated kinetic data of Cases 1–4.
Application of differential Friedman [35] and integral FWO [36] methods reveals large (up to 35%) discrepancies of calculated Ea values for the same input data (Figure S6). A possible reason of inaccuracy of the Doyle approximation of temperature integral used in FWO method [25] was disregarded after coincidence of the FWO result with the output of Starink isoconversional method [44] that uses more precise representation of the temperature integral (Figures S6 and S7). Thus, the difference is inherent for the integral isoconversional methods and caused by considerable dependency of Ea on α, while in calculation the apparent Ea value is smoothed over 0–α region resulting in systematic error [45]. To prove this, we used the advanced integral method developed by Vyazovkin [45] that is based on the optimization of incremental change in Ea by iterative calculation over small α − (α + Δα) and therefore eliminates the “smoothing” effect. Indeed, it reveals the kinetic parameters equal to the output of differential Friedman method (Figure S6).
The results of the isoconversional Friedman analysis applied to the simulated data of Cases 1–4 are shown in Figure 2 along with the assumed values for the constituting steps (dashed lines). For Case 1, the Ea value calculated from DTG indicates the correct value of Ea1 at the beginning and increases toward Ea2 in α range characterized by overlapping of two reaction processes but goes back to Ea1 at the final stage (Figure 2a). As for isoconversional dependency built on DSC data, its variation is superficial with no meaning for the lowering of Ea, although the values at the beginning and end of the reaction correspond to Ea1 as in the DTG.
For Case 2 that describe two partially overlapping processes with similar Ea values, the Ea values determined from DTG are close to the correct Ea and approximately constant during the course of reaction. Due to distinction between the exothermic and endothermic effects on DSC curve of Case 2, the conversion, derived normally as a partial area under curve, will not only increase but decrease as the reaction advances. This fact results in α value larger than unity or smaller than zero depending on the magnitude relation between the exothermic and endothermic effects (Points designated as circles in Figure 1b,c correspond to the moment when the total conversion, calculated from DSC, is equal to unity). Thus, “normal” isoconversional procedure will treat the truncated to 1 part of DSC data offering the partial and distorted kinetic information. One of possible empirical procedures to deal with the “α > 1” issue is the separate consideration of the exothermic and endothermic parts [46]. The breakpoint between two parts is where DSC signal changes its sign (denoted by square in Figure 1b). To put the isoconversional kinetic parameters on a single plot we adjust the first “exo“ part as 0..0.7 (= Qexo/(Qexo + |Qendo|)) range of overall conversion, while the second “endo” part – to 0.7..1 region (Figure 2b). The results of the “separation” procedure correspond well to the isoconversional data calculated on DTG and the correct Ea for both stages. Another possible approach to the “α > 1” issue is representing of the overall α as the index of the reaction advancement calculated as the partial area under the absolute magnitude of DSC data. It should be noted that the procedure lacks theoretical validity, because the sum of the absolute values of exothermic and endothermic peak areas possibly changes accompanied by the changes of overlapping degree of the exothermic and endothermic processes with β. Even in such theoretically invalid situation, the apparent Ea values calculated from DSC of Case 2 are close to the correct Ea value (Figure 2b, |DSC| curve).
The results for DTG of Case 3 seem to be more relevant as the barriers for stages are different. Although DTG curve looks like a smooth single peak (Figure 1c), the isoconversional analysis catches the variation of Ea values during the process (Figure 2c). The calculated Ea values at the beginning and the end of the overall reaction indicate the correct Ea values for the first and second mass-loss processes, respectively. For DSC data the “α > 1” issue arises due to clear exothermic and endothermic parts of the heat flow. The separate isoconversional analysis of exo- and endothermic parts reveals the activation energy for the first stage about 175 kJ mol−1 (exact value is 185 kJ mol−1), for the second endothermic stage Ea value approaches to 120 kJ mol−1, which is the correct value. Treatment of the absolute magnitude of DSC data as in Case 2, seems doubtful due to significant change in the sum of absolute peak areas, i.e., at β = 10 K min−1 it is larger in 1.23 times than at β = 1 K min−1. However, the Ea values calculated from |DSC| data indicate the similar variation as the reaction advances with that from DTG (Figure 2c).
Case 4 represents two stages that are overlapping in a smooth single peak either in DSC or DTG (Figure 1d). Because the Ea values assumed for the first and second reaction steps are identical, the results of Friedman method represent the close correspondence to the correct Ea value during the course of reaction both for DTG and DSC data (Figure 2d). The multistep reaction behavior cannot be deduced from the TA curves and the results of isoconversional analysis in this instance.
These results show the strength of isoconversional technique in catching the change of Ea in course of reaction and its weakness in case of significant overlapping. To evaluate the isoconversional results further, we compare the prediction of thermal behavior based on isoconversional kinetic parameters with the output of the exact model. Indeed, the results of the isoconversional analysis appear as the mathematical representation of the process. They can be easily extended outside of the temperature region where the kinetic parameters have been determined, i.e., for prediction purposes (e.g., [47]). What are the limits of this procedure, and does the isoconversional parameters derived for considered complex process allow successful prediction? Answering these questions, we compare the model output (exact solution) for heating rates from 10−9 to 109 K min−1 with the prediction based on isoconversional results (obtained with data at 1–10 K min−1 rates). In Table S1 the results of the isoconversional-based simulation and exact solution are compared in terms of the peak temperatures of the conversion rate, Figure S8 depicts the reaction profiles. For Cases 2 and 4, i.e., where the Ea values for both stages are the same, the isoconversional prediction based on TG data is fully consistent with the true outputs from the assumed model. Isoconversional results for Case 2, where the overall α calculated from DSC goes above unity, reveal that for 0–1 region the prediction is correct (e.g., Figure S8c). To account the rest of the process, we apply two modified approaches, i.e., analysis of separate exo- and endothermic parts, and absolute magnitude of DSC, both results in inaccurate prediction. For Cases 1 and 3 with more complex kinetic pattern the isoconversional technique gives the results differing from the exact output of the assumed model in terms of the peak temperatures up to 260 K using DSC data and 75 K using TG data. Therefore, in case of high reaction complexity, the predictions based on the isoconversional kinetic parameters have to be performed with great caution.

3.4. Mathematical Deconvolution and Subsequent Combined Kinetic Analysis

An alternative approach is to extract (deconvolute) the constituent peaks from the complex profile for their subsequent kinetic analysis. This approach is based on the mathematical deconvolution—the empirical peak separation by summing up mathematically fitted component peaks using a fitting function, F(t) [48]:
d α d t = i = 1 N F i ( t ) ,  
where N is the total number of component reaction steps. The fitting function with the asymmetric peak shape is generally recommended, e.g., Fraser–Suzuki function [27]:
F ( t ) = a 0 exp ( ln 2 [ ln ( 1 + 2 a 3 t a 1 a 2 ) a 3 ] 2 ) ,  
where a0, a1, a2, a3 are the amplitude, position, half-width and asymmetry of the curve, respectively. This approach has been successfully applied for real-life complex processes [49,50].
The deconvolution procedure for DSC data of Case 1, where the endothermic event is not apparent, yields the superficial results as shown in Figure 3a. The kinetic data derived from DSC (thick line) is equally well, with the correlation coefficient higher than 0.997, described by the superposition of two exothermic processes (dotted lines) or superposition of exothermic and endothermic events (dashed lines). Moreover, for the latter (correct choice) the deconvoluted exothermic and endothermic peaks are largely different from the assumed contribution of stages, i.e., γ (dα1/dt) and (1 − γ) (dα2/dt) (solid lines, Figure 3a). The problem of selecting the peak combination is easily eliminated for Cases 2 and 3, where the overall process clearly represents the combination of the exothermic and endothermic steps. Truly, in Case 2, where both stages have the same Ea value, but different A values, the deconvoluted peaks closely resemble the true reactions conversions (Figure S9a). However, in more realistic Case 3, when the Ea values for component reactions differ, the deconvoluted peaks although being exothermic and endothermic, not follow true reaction rates (Figure S9b). For DTG data, one is caught between deconvolution with two peaks (dashed lines, Figure 3b) or single peak (dotted line). However, the deconvoluted curves again are different from the true ones, i.e., η(dα1/dt) and (1 − η)(dα2/dt) (solid lines, Figure 3b).
In the practical kinetic analysis, usually it is not possible to recognize the superficial mathematical deconvolution as was seen in Figure 3. For critical evaluation of the results we investigate how these superficial deconvoluted reactions for Cases 1–4 (Figure 3, Figure S9) behave being further subjected to the kinetic analysis. The kinetic analysis assumed a single-step reaction, so called combined kinetic analysis (CKA) [26], is based on the optimization of the linearized equation:
ln d α d t ln [ α m ( 1 α ) n ] = ln ( c A ) E a R T ,  
During the optimization we use for m and n the ranges of –1 ≤ m ≤ 2 and 0 ≤ n ≤ 2 as the reasonable limits [10]. The calculated linear dependency of Equation (7) is shown in Figure 4a for the first reaction in Case 1, deconvoluted as the combination of two exothermic processes on DSC (incorrect scheme!). Note, that the nice linearity not allows one to suspect the superficial nature of analyzed data. After the determination of the optimum m and n values, the reaction type have been discriminated by comparing the reduced f(α)/f(0.5) function with that for ideal reaction types. The analyzed reaction follows three-dimensional diffusion of Ginstling–Brounshtein model (designated as D4 in Figure 4b)!
Table S3 summarizes the apparent kinetic results obtained via the mathematical peak deconvolution and subsequent CKA treatment. Although the kinetic calculation using CKA method is superficially successful as was seen in Figure 4, the correct kinetic parameters were obtained only for Cases 2 and 4 with the same Ea for both component reactions. Friedman isoconversional analysis of the deconvoluted peaks (Figure S10) supports the conclusion: improper deconvolution disturbs the results, although some kinetic features of the initial kinetic model can still be noticed. Overall, the advantage of the mathematical deconvolution—high flexibility and excellent fit quality—appears as its drawback when the superficial peaks are reconstructed from the reaction profiles of the complex process.

3.5. Kinetic Deconvolution Analysis

Another deconvolution approach – kinetic deconvolution analysis (KDA) [10,48] –represents the input data as a superposition of the independent stages, i, with contributions ci:
d α d t = i n c i A i exp ( E a , i R T ) f i ( α i ) .  
An empirical reaction model f(α) is usually taken in a flexible Šesták–Berggren form [51]:
f ( α ) = α m ( 1 α ) n [ ln ( 1 α ) ] p .  
Nonlinear regression is carried out for fitting the calculated data based on Equation (8) to the experimental data giving the optimized kinetic parameters ci, Ea, A, m, n, p for each reaction step.
Whereas the mathematical deconvolution-based analysis and the KDA are based on the similar idea of superposition with several peaks, these two approaches are largely different in the light of methodological procedures. KDA accepts the “kinetic” shape of peaks according to Equations (8), (9), and this restriction results in more sound kinetic results (vide infra). However, the considerable number of optimized parameters imposes heavy demands on the reliability of the initial values for kinetic parameters. Therefore, the detailed kinetic information collected using a range of physicochemical and microscopic techniques is necessary to the successful KDA, in addition to preliminary kinetic calculations for kinetic curves recorded systematically [3,7,32].
To assess the value of the KDA for analyzed data and allow the nonlinear regression to converge, we take the Šesták-Berggren model in truncated form, i.e., with p = 0, and use the recently proposed approach for scheme of two partially overlapping reactions [31]. Therein, combining Equations (3) and (4), one obtains:
P ( T ) = ( d α d t ) DSC ( d α d t ) DTG = γ d α 1 d t + ( 1 γ ) d α 2 d t η d α 1 d t + ( 1 η ) d α 2 d t .  
Analyzing the P(T) dependency, we look at the beginning of the process, when dα2/dt ≈ 0 and Pγ/η. In fact, in all Cases 1–4, the P(T) curve starts from value of 4 (Figure 5a). In Cases 2 and 3, where at the end of the process k1 exceeds k2, we obtain P ≈ (1 − γ)/(1 − η) = −2 since dα1/dt ≈ 0. Therefore, in specific instances when the contribution of either reaction step is negligible at the beginning and the end of the process, values of γ and η can be calculated from starting and final P values. These γ and η values were used as initial guesses in regression according to Equation (8).
Results of KDA show its advantage and robustness in revealing the correct reaction model, i.e., first-order reaction, using either DSC or DTG data (Table S3). As for kinetic pairs (Ea, lgA) the correct results were obtained for simple Cases 2 and 4 where the ratio between the rate constants are independent on temperature. Considering the Cases 1 and 3, the statistical and visual (Figure 5b) quality of the fit is still high, but the Ea values differ from that in exact model. However, the A values compensate this difference in accordance with the kinetic compensation effect (KCE). It seems that this KCE is caused by the conceptual feature of the above deconvolution techniques: the conversion rate in form of Equation (8) in general is not equal to that parameter for the successive reactions due to unaccounted dependency of the component kinetic processes as considered in Equations (2) and (3) by the dependence of dα2/dt on α1 or in case of parallel processes because of the same interdependence of the conversion degrees for the components.

3.6. Formal Kinetic Analysis

The last considered method is formal kinetic analysis. It combines the strength of the model-fitting techniques with allowing for various connections between the process stages (e.g., successive, parallel and independent steps). This approach seems to be the most advantageous for the complex multistep process especially with exothermic and endothermic reactions in line with observations by Burnham [40]. The kinetic calculations have been performed with THINKS software by assuming single and various combinations of two-step reactions.
While performing the formal kinetic analysis, one faces with even higher (than for KDA) demands of the appropriate initial guesses and a need of the selection of the reaction scheme prior to calculation. Till now the selection of the reaction pattern depends on the experience of the researcher and the supplementary physicochemical insights to the process that are available. Recently, Muravyev and Pivkina propose to use the artificial neural networks to make a “recognition” of the kinetic scheme, however, the results only for the single-step reactions have been obtained [52]. Herein, we first tried the correct scheme with two subsequent first-order reactions to obtain the kinetic parameters and compare them with exact values. Overall, the kinetic parameters that were optimized throughout nonlinear regression based on DSC or DTG are close to that in assumed model (Table S4). But can we show that this kinetic scheme is indeed better than the above KDA model with two independent reactions, or model with the single-step reaction (e.g., for single DTG peak, Case 1)? From our calculations, all these schemes offer equally high correlation coefficients in certain cases and to select statistically best model we use the recently proposed advanced tool, the Bayes information criteria (BIC) [53,54]. Table 2 and Table S5 list the BIC values with the minimal one corresponding to the statistically best model. Evidently, true kinetic scheme can be identified by this kind of analysis. Also, we perform the calculations on joint DSC and DTG input data as was proposed recently by authors [31], i.e., the same reaction scheme but various contributions of stages for gravimetric and caloric data. Results show that using both datasets simultaneously increases the accuracy in kinetic parameters determination and allows easier distinguishing of the appropriate kinetic scheme (Figures S11 and S12).
High performance of the model-fitting analysis is evidenced by these results; however, the analyzed kinetic curves were drawn by assuming overlapping two-step reactions that are characterized by homogeneous-like first-order reactions with a simple interaction between the reaction steps as expressed by f1(α1) = 1 − α1 and f2(α1, α2) = α1α2. Further careful examinations with more complex input data are required to further assess the strength of the formal kinetic approach and its applicability to the real heterogeneous reactions.

4. Summary and Outlook

The kinetic parameters calculated by above methods for both reaction steps of Case 1 are brought together in Figure 6 as lgA dependency on Ea. The exact values (used to simulate the data) are marked as red circles on the plot. Results of Kissinger method, kinetic parameters for deconvoluted peaks and the output of the formal kinetic analysis fall onto the straight line drawn through the exact values of (lgA, Ea). Same trend is observed for Case 3 (Figure S13). This phenomenon is known as kinetic compensation effect and it has been extensively discussed in literature [55,56,57]. Application of the incorrect reaction model f(α) or inappropriate computational technique to fit the data will manifest in results as KCE [57,58,59,60,61]. Apparently, our results shown in Figure 6 and Figure S13 reveal this, mathematical, type of KCE.
Summarizing the findings of present study, we can highlight the guidelines for designing the kinetic approach to the considered type of complex processes. First of all, it should be noted that the kinetic data resolved using TG and DSC are generally different for the multistep reactions and the observation of both signals is advantageous. Application of the isoconversional method can be useful for preliminary kinetic assessment as it often reflects some kinetic features of the process, e.g., the starting/final Ea values. However, the resulting changes of Ea during reaction course can easily become superficial when significant overlapping of stages takes place. The situation will propagate in isoconversional prediction; thus, it has to be performed with great precautions. The mathematical deconvolution technique shows high flexibility in catching of the reaction profile, while the flip side is extremely unstable kinetic results. Kinetic deconvolution analysis uses the same idea of superposition, but its “kinetic shape” of peaks increases robustness of the approach drastically. KDA method is recommended for the preliminary evaluation of the kinetic parameters, number of stages, and the reaction types. The formal kinetic analysis is shown to be the most useful in present evaluation since it provides the correct parameters for each reaction step in considered process. Its main problem is selection of an appropriate reaction scheme, must be solved in each particular case with the help of the preliminary kinetic analysis and additional physicochemical insights to the process.
There appears to be no single approach leading to the correct kinetic triplet in all situations but having the clear idea about the considered process and having analyzed the whole available experimental data (DSC and TG, structure changes, evolved gases etc.) will help one to select an appropriate way to perform the kinetic analysis and obtain correct results. As a general framework that eases the selection of particular kinetic calculation method, we use the following logical approach (Figure 7). The synthetic kinetic data that reveal the specific features of the real data was first generated. Then possible kinetic calculation methods are applied and justified by comparison of their output with exact kinetic parameters used in model. Last step is performing the kinetic analysis using selected kinetic calculation method over the real experimental data (i.e., with unknown answer). Present study and previous experimental work [31] taken altogether show the value of recommended approach.

Supplementary Materials

The following are available online, Figures S1–S13. Tables S1–S5.

Author Contributions

Conceptualization, N.V.M. and N.K.; methodology, N.V.M. and N.K.; software, N.V.M.; validation, N.V.M., A.N.P., and N.K.; formal analysis, N.V.M. and A.N.P.; investigation, N.V.M., resources, A.N.P.; data curation, N.V.M.; writing—original draft, N.V.M.; writing—review & editing, N.V.M. and N.K., visualization, N.V.M. and N.K., project administration, N.V.M. and A.N.P.; supervision, A.N.P.

Funding

This research was funded by the State research project 0082-2018-0002, AAAA-A18-118031490034-6 to Semenov Institute of Chemical Physics.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kaisersberger, E.; Opfermann, J. Kinetic evaluation of exothermal reactions measured by DSC. Thermochim. Acta 1991, 187, 151–158. [Google Scholar] [CrossRef]
  2. Burnham, A.K.; Weese, R.K.; Weeks, B.L. A Distributed Activation Energy Model of Thermodynamically Inhibited Nucleation and Growth Reactions and Its Application to the β−δ Phase Transition of HMX. J. Phys. Chem. B 2004, 108, 19432–19441. [Google Scholar] [CrossRef]
  3. Kitabayashi, S.; Koga, N. Thermal Decomposition of Tin(II) Oxyhydroxide and Subsequent Oxidation in Air: Kinetic Deconvolution of Overlapping Heterogeneous Processes. J. Phys. Chem. C 2015, 119, 16188–16199. [Google Scholar] [CrossRef] [Green Version]
  4. Galwey, A.K.; Brown, M.E. Thermal Decomposition of Ionic Solids: Chemical Properties and Reactivities of Ionic Crystalline Phases; Elsevier: Burlington, NJ, USA, 1999; ISBN 978-0-08-054279-9. [Google Scholar]
  5. Long, G.T.; Vyazovkin, S.; Brems, B.A.; Wight, C.A. Competitive Vaporization and Decomposition of Liquid RDX. J. Phys. Chem. B 2000, 104, 2570–2574. [Google Scholar] [CrossRef] [Green Version]
  6. Muravyev, N.V.; Monogarov, K.A.; Asachenko, A.F.; Nechaev, M.S.; Ananyev, I.V.; Fomenkov, I.V.; Kiselev, V.G.; Pivkina, A.N. Pursuing reliable thermal analysis techniques for energetic materials: Decomposition kinetics and thermal stability of dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate (TKX-50). Phys. Chem. Chem. Phys. 2017, 19, 436–449. [Google Scholar] [CrossRef] [PubMed]
  7. Koga, N.; Kameno, N.; Tsuboi, Y.; Fujiwara, T.; Nakano, M.; Nishikawa, K.; Iwasaki Murata, A. Multistep thermal decomposition of granular sodium perborate tetrahydrate: A kinetic approach to complex reactions in solid–gas systems. Phys. Chem. Chem. Phys. 2018, 20, 12557–12573. [Google Scholar] [CrossRef] [PubMed]
  8. Kameno, N.; Koga, N. Heterogeneous Kinetic Features of the Overlapping Thermal Dehydration and Melting of Thermal Energy Storage Material: Sodium Thiosulfate Pentahydrate. J. Phys. Chem. C 2018, 122, 8480–8490. [Google Scholar] [CrossRef]
  9. Opfermann, J.; Hädrich, W. Prediction of the thermal response of hazardous materials during storage using an improved technique. Thermochim. Acta 1995, 263, 29–50. [Google Scholar] [CrossRef]
  10. Sánchez-Jiménez, P.E.; Perejón, A.; Criado, J.M.; Diánez, M.J.; Pérez-Maqueda, L.A. Kinetic model for thermal dehydrochlorination of poly(vinyl chloride). Polymer 2010, 51, 3998–4007. [Google Scholar] [CrossRef] [Green Version]
  11. Koga, N.; Yamada, S.; Kimura, T. Thermal Decomposition of Silver Carbonate: Phenomenology and Physicogeometrical Kinetics. J. Phys. Chem. C 2013, 117, 326–336. [Google Scholar] [CrossRef]
  12. Flynn, J.H.; Wall, L.A. General treatment of the thermogravimetry of polymers. J. Res. Natl. Bur. Stand. Phys. Chem. 1966, 70A, 487. [Google Scholar] [CrossRef]
  13. Vyazovkin, S.; Linert, W. The Application of Isoconversional Methods for Analyzing Isokinetic Relationships Occurring at Thermal Decomposition of Solids. J. Solid State Chem. 1995, 114, 392–398. [Google Scholar] [CrossRef]
  14. Agrawal, R.K. Analysis of non-isothermal reaction kinetics: Part 2. Complex reactions. Thermochim. Acta 1992, 203, 111–125. [Google Scholar] [CrossRef]
  15. Svoboda, R.; Málek, J. Is the original Kissinger equation obsolete today? J. Therm. Anal. Calorim. 2014, 115, 1961–1967. [Google Scholar] [CrossRef]
  16. Kissinger, H.E. Reaction Kinetics in Differential Thermal Analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  17. ASTM. ASTM E698-05 Standard Test Method for Arrhenius Kinetic Constants for Thermally Unstable Materials; ASTM International: West Conshohocken, PA, USA, 2005. [Google Scholar]
  18. Vyazovkin, S.V.; Goryachko, V.I.; Lesnikovich, A.I. An approach to the solution of the inverse kinetic problem in the case of complex processes. Part III. Parallel independent reactions. Thermochim. Acta 1992, 197, 41–51. [Google Scholar] [CrossRef]
  19. Vyazovkin, S. Conversion dependence of activation energy for model DSC curves of consecutive reactions. Thermochim. Acta 1994, 236, 1–13. [Google Scholar] [CrossRef]
  20. Vyazovkin, S.V.; Lesnikovich, A.I. An approach to the solution of the inverse kinetic problem in the case of complex processes. Thermochim. Acta 1990, 165, 273–280. [Google Scholar] [CrossRef]
  21. Criado, J.M.; Sánchez-Jiménez, P.E.; Pérez-Maqueda, L.A. Critical study of the isoconversional methods of kinetic analysis. J. Therm. Anal. Calorim. 2008, 92, 199–203. [Google Scholar] [CrossRef] [Green Version]
  22. Cai, J.; Wu, W.; Liu, R. Isoconversional Kinetic Analysis of Complex Solid-State Processes: Parallel and Successive Reactions. Ind. Eng. Chem. Res. 2012, 51, 16157–16161. [Google Scholar] [CrossRef]
  23. Moukhina, E. Determination of kinetic mechanisms for reactions measured with thermoanalytical instruments. J. Therm. Anal. Calorim. 2012, 109, 1203–1214. [Google Scholar] [CrossRef] [Green Version]
  24. Burnham, A.K. Application of the Šesták-Berggren Equation to Organic and Inorganic Materials of Practical Interest. J. Therm. Anal. Calorim. 2000, 60, 895–908. [Google Scholar] [CrossRef]
  25. Vyazovkin, S.; Burnham, A.K.; Criado, J.M.; Pérez-Maqueda, L.A.; Popescu, C.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim. Acta 2011, 520, 1–19. [Google Scholar] [CrossRef]
  26. Pérez-Maqueda, L.A.; Criado, J.M.; Sánchez-Jiménez, P.E. Combined Kinetic Analysis of Solid-State Reactions: A Powerful Tool for the Simultaneous Determination of Kinetic Parameters and the Kinetic Model without Previous Assumptions on the Reaction Mechanism. J. Phys. Chem. A 2006, 110, 12456–12462. [Google Scholar] [CrossRef]
  27. Perejón, A.; Sánchez-Jiménez, P.E.; Criado, J.M.; Pérez-Maqueda, L.A. Kinetic Analysis of Complex Solid-State Reactions. A New Deconvolution Procedure. J. Phys. Chem. B 2011, 115, 1780–1791. [Google Scholar] [CrossRef] [Green Version]
  28. Noda, Y.; Koga, N. Phenomenological Kinetics of the Carbonation Reaction of Lithium Hydroxide Monohydrate: Role of Surface Product Layer and Possible Existence of a Liquid Phase. J. Phys. Chem. C 2014, 118, 5424–5436. [Google Scholar] [CrossRef]
  29. Koga, N.; Kodani, S. Thermally induced carbonation of Ca(OH)2 in a CO2 atmosphere: Kinetic simulation of overlapping mass-loss and mass-gain processes in a solid–gas system. Phys. Chem. Chem. Phys. 2018, 20, 26173–26189. [Google Scholar] [CrossRef]
  30. Schmid, H.; Eisenreich, N.; Krause, C.; Pfeil, A. Analysis of complex thermoanalytical curves: The thermodynamic and kinetic parameters of isopropylammonium nitrate. J. Therm. Anal. Calorim. 1989, 35, 569–576. [Google Scholar] [CrossRef]
  31. Muravyev, N.V.; Koga, N.; Meerov, D.B.; Pivkina, A.N. Kinetic analysis of overlapping multistep thermal decomposition comprising exothermic and endothermic processes: Thermolysis of ammonium dinitramide. Phys. Chem. Chem. Phys. 2017, 19, 3254–3264. [Google Scholar] [CrossRef]
  32. Nakano, M.; Wada, T.; Koga, N. Exothermic Behavior of Thermal Decomposition of Sodium Percarbonate: Kinetic Deconvolution of Successive Endothermic and Exothermic Processes. J. Phys. Chem. A 2015, 119, 9761–9769. [Google Scholar] [CrossRef] [Green Version]
  33. Heym, F.; Etzold, B.J.M.; Kern, C.; Jess, A. Analysis of evaporation and thermal decomposition of ionic liquids by thermogravimetrical analysis at ambient pressure and high vacuum. Green Chem. 2011, 13, 1453. [Google Scholar] [CrossRef]
  34. Skala, D.; Sokić, M.; Tomić, J.; Kopsch, H. Kinetic analysis of consecutive reactions using TG and DSC techniques: Theory and application. J. Therm. Anal. Calorim. 1989, 35, 1441–1458. [Google Scholar] [CrossRef]
  35. Friedman, H.L. Kinetics of thermal degradation of char-forming plastics from thermogravimetry. Application to a phenolic plastic. J. Polym. Sci. Part C Polym. Symp. 1964, 6, 183–195. [Google Scholar] [CrossRef]
  36. Ozawa, T. A New Method of Analyzing Thermogravimetric Data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef] [Green Version]
  37. Vyazovkin, S. Isoconversional Kinetics of Thermally Stimulated Processes; Springer International Publishing: Cham, Switzerland, 2015; ISBN 978-3-319-14174-9. [Google Scholar]
  38. Sánchez-Jiménez, P.E.; Pérez-Maqueda, L.A.; Perejón, A.; Criado, J.M. A new model for the kinetic analysis of thermal degradation of polymers driven by random scission. Polym. Degrad. Stab. 2010, 95, 733–739. [Google Scholar] [CrossRef]
  39. Opfermann, J. Kinetic Analysis Using Multivariate Non-linear Regression. I. Basic concepts. J. Therm. Anal. Calorim. 2000, 60, 641–658. [Google Scholar] [CrossRef]
  40. Burnham, A.K.; Dinh, L.N. A comparison of isoconversional and model-fitting approaches to kinetic parameter estimation and application predictions. J. Therm. Anal. Calorim. 2007, 89, 479–490. [Google Scholar] [CrossRef] [Green Version]
  41. Budrugeac, P.; Homentcovschi, D.; Segal, E. Critical Considerations on the Isoconversional Methods. III. On the evaluation of the activation energy from non-isothermal data. J. Therm. Anal. Calorim. 2001, 66, 557–565. [Google Scholar] [CrossRef]
  42. Pérez-Maqueda, L.A.; Criado, J.M. The Accuracy of Senum and Yang’s Approximations to the Arrhenius Integral. J. Therm. Anal. Calorim. 2000, 60, 909–915. [Google Scholar] [CrossRef]
  43. Muravyev, N.V. THINKS—Thermokinetic Software, Moscow, 2016.
  44. Starink, M.J. The determination of activation energy from linear heating rate experiments: A comparison of the accuracy of isoconversion methods. Thermochim. Acta 2003, 404, 163–176. [Google Scholar] [CrossRef]
  45. Vyazovkin, S. Modification of the integral isoconversional method to account for variation in the activation energy. J. Comput. Chem. 2001, 22, 178–183. [Google Scholar] [CrossRef]
  46. Muravyev, N.V.; Monogarov, K.A.; Bragin, A.A.; Fomenkov, I.V.; Pivkina, A.N. HP-DSC study of energetic materials. Part I. Overview of pressure influence on thermal behavior. Thermochim. Acta 2016, 631, 1–7. [Google Scholar] [CrossRef]
  47. Roduit, B.; Borgeat, C.; Berger, B.; Folly, P.; Alonso, B.; Aebischer, J.N. The prediction of thermal stability of self-reactive chemicals: From milligrams to tons. J. Therm. Anal. Calorim. 2005, 80, 91–102. [Google Scholar] [CrossRef]
  48. Koga, N.; Goshi, Y.; Yamada, S.; Pérez-Maqueda, L.A. Kinetic approach to partially overlapped thermal decomposition processes: Co-precipitated zinc carbonates. J. Therm. Anal. Calorim. 2013, 111, 1463–1474. [Google Scholar] [CrossRef]
  49. Yan, Q.-L.; Zeman, S.; Zhang, J.-G.; Qi, X.-F.; Li, T.; Musil, T. Multistep Thermolysis Mechanisms of Azido-s-triazine Derivatives and Kinetic Compensation Effects for the Rate-Limiting Processes. J. Phys. Chem. C 2015, 119, 14861–14872. [Google Scholar] [CrossRef]
  50. Svoboda, R.; Málek, J. Applicability of Fraser–Suzuki function in kinetic analysis of complex crystallization processes. J. Therm. Anal. Calorim. 2013, 111, 1045–1056. [Google Scholar] [CrossRef]
  51. Šesták, J.; Berggren, G. Study of the kinetics of the mechanism of solid-state reactions at increasing temperatures. Thermochim. Acta 1971, 3, 1–12. [Google Scholar] [CrossRef]
  52. Muravyev, N.V.; Pivkina, A.N. New concept of thermokinetic analysis with artificial neural networks. Thermochim. Acta 2016, 637, 69–73. [Google Scholar] [CrossRef]
  53. Bohn, M.A. Assessment of description quality of models by information theoretical criteria based on Akaike and Schwarz-Bayes applied with stability data of energetic materials. In Proceedings of the 46th International Annual Conference of the Fraunhofer ICT, Karlsruhe, Germany, 23–26 June 2015; pp. 6.1–6.23. [Google Scholar]
  54. Roduit, B.; Hartmann, M.; Folly, P.; Sarbach, A.; Baltensperger, R. Prediction of thermal stability of materials by modified kinetic and model selection approaches based on limited amount of experimental points. Thermochim. Acta 2014, 579, 31–39. [Google Scholar] [CrossRef]
  55. Garn, P.D. An examination of the kinetic compensation effect. J. Therm. Anal. Calorim. 1975, 7, 475–478. [Google Scholar] [CrossRef]
  56. Agrawal, R.K. On the compensation effect. J. Therm. Anal. Calorim. 1986, 31, 73–86. [Google Scholar] [CrossRef]
  57. Koga, N. A review of the mutual dependence of Arrhenius parameters evaluated by the thermoanalytical study of solid-state reactions: The kinetic compensation effect. Thermochim. Acta 1994, 244, 1–20. [Google Scholar] [CrossRef]
  58. Koga, N.; Šesták, J. Kinetic compensation effect as a mathematical consequence of the exponential rate constant. Thermochim. Acta 1991, 182, 201–208. [Google Scholar] [CrossRef]
  59. Brown, M.E.; Galwey, A.K. The significance of “compensation effects” appearing in data published in “computational aspects of kinetic analysis”: ICTAC project, 2000. Thermochim. Acta 2002, 387, 173–183. [Google Scholar] [CrossRef]
  60. Koga, N.; Šesták, J.; Málek, J. Distortion of the Arrhenius parameters by the inappropriate kinetic model function. Thermochim. Acta 1991, 188, 333–336. [Google Scholar] [CrossRef]
  61. Koga, N.; Šesták, J. Further aspects of the kinetic compensation effect. J. Therm. Anal. Calorim. 1991, 37, 1103–1108. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds are available from the authors.
Figure 1. Comparison of the simulated kinetic curves at β = 2 K min−1 for each component reaction step calculated according to Equation (2), overall kinetic curve calculated according to Equation (3) (from DTG), and overall kinetic curve calculated according to Equation (4) (from DSC): (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Figure 1. Comparison of the simulated kinetic curves at β = 2 K min−1 for each component reaction step calculated according to Equation (2), overall kinetic curve calculated according to Equation (3) (from DTG), and overall kinetic curve calculated according to Equation (4) (from DSC): (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Molecules 24 02298 g001
Figure 2. Isoconversional Friedman analysis of the simulated DTG and DSC data at β of 1–10 K min−1: (a) Case 1, (b) Case 2 (for clarity the uncertainties are not shown), (c) Case 3, and (d) Case 4.
Figure 2. Isoconversional Friedman analysis of the simulated DTG and DSC data at β of 1–10 K min−1: (a) Case 1, (b) Case 2 (for clarity the uncertainties are not shown), (c) Case 3, and (d) Case 4.
Molecules 24 02298 g002
Figure 3. Selection of the peak combinations for Case 1, β = 1 K min−1: (a) superposition of two exothermic processes (dotted lines) or exothermic and endothermic ones (dashed lines) results in perfect correlation with DSC signal (black solid) but completely different from true heat evolution data, (b) DTG peaks are described well both by two peaks (dashed lines) and by single peak (dotted line), but both are different from true signals.
Figure 3. Selection of the peak combinations for Case 1, β = 1 K min−1: (a) superposition of two exothermic processes (dotted lines) or exothermic and endothermic ones (dashed lines) results in perfect correlation with DSC signal (black solid) but completely different from true heat evolution data, (b) DTG peaks are described well both by two peaks (dashed lines) and by single peak (dotted line), but both are different from true signals.
Molecules 24 02298 g003
Figure 4. CKA analysis of the first reaction step of Case 1 deconvoluted from DSC curve by assuming two successive exothermic processes. (a) CKA plot, (b) Comparison of theoretical master plots of different kinetic models with calculated data (circles).
Figure 4. CKA analysis of the first reaction step of Case 1 deconvoluted from DSC curve by assuming two successive exothermic processes. (a) CKA plot, (b) Comparison of theoretical master plots of different kinetic models with calculated data (circles).
Molecules 24 02298 g004
Figure 5. Kinetic deconvolution analysis: (a) Ratio of the conversion rates from DSC to DTG dependences on temperature for cases 1–4 at β = 2 K min−1, (b) Fit of DSC data for Case 3 under various heating rates (points) with Equation (8) with optimized parameters (thick red line).
Figure 5. Kinetic deconvolution analysis: (a) Ratio of the conversion rates from DSC to DTG dependences on temperature for cases 1–4 at β = 2 K min−1, (b) Fit of DSC data for Case 3 under various heating rates (points) with Equation (8) with optimized parameters (thick red line).
Molecules 24 02298 g005
Figure 6. Kinetic compensation effect on the kinetic parameters calculated for Case 1 data.
Figure 6. Kinetic compensation effect on the kinetic parameters calculated for Case 1 data.
Molecules 24 02298 g006
Figure 7. Proposed framework for the kinetic analysis of complex multistep processes.
Figure 7. Proposed framework for the kinetic analysis of complex multistep processes.
Molecules 24 02298 g007
Table 1. Kinetic parameters calculated by the Kissinger plot for the simulated successive reactions of Cases 1–4.
Table 1. Kinetic parameters calculated by the Kissinger plot for the simulated successive reactions of Cases 1–4.
CaseStepAssumed Kinetic Parameters 1Ea/kJ mol−12lg(A/s−1) 2Δαp 3
Ea/kJ mol−1lg(A/s−1)DTGDSCDTGDSCDTGDSC
1112010138 ± 4148 ± 112.9 ± 0.413.7 ± 0.10.06 0.19
218517.5160 ± 1014.0 ± 1.10.19
2111510.4115 ± 1115 ± 110.3 ± 0.110.4 ± 0.10.0020.007
21159.4115 ± 1115 ± 19.4 ± 0.019.4 ± 0.10.0020.003
3118517.5203 ± 4177 ± 219.3 ± 0.416.7 ± 0.30.110.11
212010159 ± 314.0 ± 0.30.16
411159.5116 ± 3116 ± 29.6 ± 0.39.6 ± 0.20.030.02
211511
1 Stands for the kinetic parameters used in Equations (2)–(4) to calculate the input data. 2 Correlation coefficients of linear fit in Kissinger plot is higher than 0.999 for all cases. 3 Variation of αp with β (1–10 K min−1).
Table 2. Bayes information criteria for the formal kinetic analysis of Case 1 data with several reaction schemes.
Table 2. Bayes information criteria for the formal kinetic analysis of Case 1 data with several reaction schemes.
Kinetic SchemeDSC DataDTG DataDSC+DTG Data
Single-step reaction−4671−6339−23010
Two parallel reactions−4888−5355−18186
Two consecutive reactions−9873−10733−32343
Two independent reactions (KDA)−6791−7611−24676

Share and Cite

MDPI and ACS Style

Muravyev, N.V.; Pivkina, A.N.; Koga, N. Critical Appraisal of Kinetic Calculation Methods Applied to Overlapping Multistep Reactions. Molecules 2019, 24, 2298. https://doi.org/10.3390/molecules24122298

AMA Style

Muravyev NV, Pivkina AN, Koga N. Critical Appraisal of Kinetic Calculation Methods Applied to Overlapping Multistep Reactions. Molecules. 2019; 24(12):2298. https://doi.org/10.3390/molecules24122298

Chicago/Turabian Style

Muravyev, Nikita V., Alla N. Pivkina, and Nobuyoshi Koga. 2019. "Critical Appraisal of Kinetic Calculation Methods Applied to Overlapping Multistep Reactions" Molecules 24, no. 12: 2298. https://doi.org/10.3390/molecules24122298

APA Style

Muravyev, N. V., Pivkina, A. N., & Koga, N. (2019). Critical Appraisal of Kinetic Calculation Methods Applied to Overlapping Multistep Reactions. Molecules, 24(12), 2298. https://doi.org/10.3390/molecules24122298

Article Metrics

Back to TopTop