4 Kinetic Analysis Based on Thermogravimetry Data

DOI

10.34663/9783945561171-06

Citation

Arhangel‘skii, Igor V., Dunaev, Alexander V., Makarenko, Irina V., Tikhonov, Nikolay A. and Tarasov, Andrey V. (2013). Kinetic Analysis Based on Thermogravimetry Data. In: Non-Isothermal Kinetic Methods: Workbook and Laboratory Manual. Berlin: Max-Planck-Gesellschaft zur Förderung der Wissenschaften.

4.1 Dehydration of Calcium Oxalate Monohydrate

Let us consider the dehydration as the reaction that occurs by the scheme , that is,

4.1

Figure 4.1 shows the TG curves of the dehydration of calcium oxalate monohydrate.

Fig. 4.1: TG and DTG curves of the dehydration of calcium oxalate monohydrate at different heating rates.

Fig. 4.1: TG and DTG curves of the dehydration of calcium oxalate monohydrate at different heating rates.

The CaC2O4 · H2O dehydration was studied on a Netzsch TG 209 F 3 Tarsus thermo-microbalance. The experiment was run at three heating rates: 5, 7.5, and 10 K/min. Three measurements were taken at each heating rate, other conditions being identical. Standard aluminum crucibles without lids were used as holders. The process was carried out in a dry air flow at a rate of 200 mL/min. The initial reagent was freshly precipitated calcium oxalate with a particle size of 15–20 μm. Weighed portions of the reagent were 5–6 mg for each heating rate.

4.2 Computational Procedure. Solution of the Inverse and Direct Kinetic Problems. Quasi-One-Stage Process

As follows from Figure 4.1, the dehydration in the given temperature range can be considered a quasi-one-stage reaction at all heating rates used. The experimental data obtained on NETZSCH equipment are processed with the NETZSCH Proteus software.

For further work with the experimental data using the NETZSCH Thermokinetics software, it is necessary to export data from the Proteus program in the tabulated form (measured signal as a function of temperature or time) as an ASCII file. To do this, the user must select the desired curve and click the Extras Export data button in the Proteus toolbar. The user must enter the lower and upper limits of the data range to be exported. To correctly specify the limits, the derivative of the selected curve is used. The left- and right-hand limits are chosen in the ranges where the derivative becomes zero (Figure 4.2).

Fig. 4.2: Selection of the number of points while exporting data.

Fig. 4.2: Selection of the number of points while exporting data.

Remember that the derivative of the selected curve can be obtained by clicking the corresponding icon in the NETZSCH Proteus program window.

An important step for further calculations is to select the number of points on the TG curve used for subsequent processing. This can be done in two ways. While exporting the data from the NETZSCH Proteus software, the user enters either the number of points of the curve to be exported (input window No.Pts., Figure 4.2) or the temperature increment between neighboring points (window Steps, Figure 4.2) from which the program calculates the number of points. In our case, from the parallel measurements performed under identical conditions, the root-mean-square deviation of the temperature for the rate curve minimum (DTG) was calculated to be S . For the selected temperature range, 100 points for each curve are used. It is worth noting that the use of a smaller increment is physically unreasonable since the S value reflects the properties of the system in question, but is not the temperature measurement accuracy provided by an instrument, which is an order of lower magnitude. In this case, specific features of the process should be taken into account.

If necessary, the experimental curves can be smoothed before exporting data, using the NETZSCH Proteus Analysis program (Figures 4.3, 4.4). The smooth level (Figure 4.4) is considered sufficient when the DTG curve has no spikes.

Fig. 4.3: Curve smoothing.

Fig. 4.3: Curve smoothing.

Fig. 4.4: Selection of the curve smooth level.

Fig. 4.4: Selection of the curve smooth level.

After all data has been imported as ASCII files, they can be used for kinetic analysis with the NETZSCH Thermokinetics software. To do this, a new project is created (Figure 4.5). In the initialization dialog box, the project name, the sample, the number of measured curves (scans), and the type of measurement (TG, DSC, etc.) are entered.

Fig. 4.5: Creation of a new project in the NETZSCH Thermokinetics software.

Fig. 4.5: Creation of a new project in the NETZSCH Thermokinetics software.

Fig. 4.6: Initiation of a new project in the NETZSCH Thermokinetics software.

Fig. 4.6: Initiation of a new project in the NETZSCH Thermokinetics software.

The user imports then the data into the program. To do this, the type of loaded data should be specified (Figure 4.7). In the dialog box (Figure 4.8), the type of measurement is specified. In the case of thermogravimetry, GTA TG is selected.

Fig. 4.7: Selection of the type of measurement for loading in the NETZSCH Thermokinetics software.

Fig. 4.7: Selection of the type of measurement for loading in the NETZSCH Thermokinetics software.

Fig. 4.8: Selection of the measurement type for loading in the NETZSCH Thermokinetics software.

Fig. 4.8: Selection of the measurement type for loading in the NETZSCH Thermokinetics software.

The data is loaded by clicking the Load ASCII file icon (Figure 4.9). In the popup window, the user specifies the desired data file.

Fig. 4.9: Loading data into the NETZSCH Thermokinetics software.

Fig. 4.9: Loading data into the NETZSCH Thermokinetics software.

If necessary, the user corrects the evaluation range limits (Figure 4.10).

Fig. 4.10: Selection of the evaluation range limits in the NETZSCH Thermokinetics software.

Fig. 4.10: Selection of the evaluation range limits in the NETZSCH Thermokinetics software.

The loaded files are checked by clicking the Check scans icon (Figure 4.11).

Fig. 4.11: Check of the loaded data in the NETZSCH Thermokinetics software.

Fig. 4.11: Check of the loaded data in the NETZSCH Thermokinetics software.

The next step is to perform the model-free analysis. To do this, the user should enter the Model free menu (Figure 4.12).

Fig. 4.12: Model-free analysis.

Fig. 4.12: Model-free analysis.

Fig. 4.13: Logarithm of heating rate versus inverse temperature.

Fig. 4.13: Logarithm of heating rate versus inverse temperature.

The Ozawa–Flynn–Wall analysis is started by pressing the corresponding button to show the graph of the logarithm of heating rate versus inverse temperature (Figure 4.13).

Fig. 4.14: Activation energy and preexponential factor versus conversion.

Fig. 4.14: Activation energy and preexponential factor versus conversion.

The activation energy and preexponential factor are then calculated as a function of conversion (Figure 4.14). Analogous calculations can be performed by the Friedman method (Figure 4.15).

Fig. 4.15: Model-free Friedman analysis.

Fig. 4.15: Model-free Friedman analysis.

The resulting activation energies and preexponential factors are used as a zero approximation to solve the direct kinetic problem. To do this, the Linear regr.-toolbar is used (Figure 4.16). This procedure is applicable to one-step processes.

Fig. 4.16: Linear regression.

Fig. 4.16: Linear regression.

Fig. 4.17: Model selection in case of linear regression.

Fig. 4.17: Model selection in case of linear regression.

Fig. 4.18: Graphical comparison of the computation results (solid line) and experimental data (symbols) for TG curves of the CaC2O4∙H2O dehydration. Three heating rates; one-step reaction.

Fig. 4.18: Graphical comparison of the computation results (solid line) and experimental data (symbols) for TG curves of the CaC2O4∙H2O dehydration. Three heating rates; one-step reaction.

Fig. 4.19: Kinetic parameters of the CaC2O4 · H2O dehydration.

Fig. 4.19: Kinetic parameters of the CaC2O4 · H2O dehydration.

If the kinetic model of a process is a priori unknown, it is advisable to perform calculations using all models (Figure 4.17). To do this, in the Selection box, the user should click on ‘Yes’ in the corresponding cell of the column ‘Selected’. The program outputs a graph of measured and calculated curves, as well as tables.

4.3 Analysis of Computation Results

Let us consider the computation results obtained by the linear regression method for the CaC2O4 · H2O dehydration (Figure 4.19).

Figure 4.19 presents the Arrhenius parameters and the form and characteristics of the function best fitting the experimental results (from the statistical viewpoint). For the reaction under consideration, the best fitting function is the Prout–Tompkins equation with autocatalysis (the Bna code), which is indicated at the top left of the table. However, before discussing the meaning of the results obtained, let us consider the F-test: Fit Quality window (Figure 4.20).

Fig. 4.20: Fit quality for different models.

Fig. 4.20: Fit quality for different models.

4.3.1 The F-test: Fit quality and F-test

Step significance windows present the statistical analysis of the fit quality for different models. This allows us to determine using the statistical methods which of the models provides the best fit for the experimental data.

To perform such an analysis, Fisher’s exact test is used. In general, Fisher’s test is a variance ratio which makes it possible to verify whether the difference between two independent estimates of the variance of some data samples is significant. To do this, the ratio of these two variances is compared with the corresponding tabulated value of the Fisher distribution for a given number of degrees of freedom and significance level. If the ratio of two variances exceeds the corresponding theoretical Fisher test value, the difference between the variances is significant.

In the Thermokinetics software, Fisher’s test is used for comparing the fit qualities ensured by different models. The best-fit model, that is, the model with the minimal sum of squared deviations, is taken as a reference (conventionally denoted as model 1). Then, each model is compared to the reference model. If the Fisher test value does not exceed the critical value, the difference between current model 2 and reference model 1 is insignificant. There is no reason to then believe that model 1 provides a more adequate description of the experiment in comparison to model 2.

The Fexp value is estimated by means of Fisher’s test:

4.2

The Fexp value is compared with the Fisher distribution Fcrit(0.95) for the significance level of 0.95 and the corresponding number of degrees of freedom.

Figure 4.20 shows that three models meet this requirement: the Prout–Tompkins nth-order equation with autocatalysis (Bna), the nth-order equation with autocatalysis (Cnb), and the nth-order reaction (Fn) (see Table 3.1). To discriminate between these functions, a nonlinear regression calculation is performed for each function separately under the assumption of a one-step process. To do this, the Nonlin regr. toolbar is opened (Figure 4.21).

Fig. 4.21: Nonlinear regression.

Fig. 4.21: Nonlinear regression.

Fig. 4.22: Menu for model selection.

Fig. 4.22: Menu for model selection.

Fig. 4.23: Selection of the model and function type in the nonlinear regression menu.

Fig. 4.23: Selection of the model and function type in the nonlinear regression menu.

Fig. 4.24: Table of initial parameters.

Fig. 4.24: Table of initial parameters.

In this toolbar, the Model Selection button is pressed (Figure 4.22) to select the process model. In the opened window (Figure 4.23), we select the one-step process and the Cn model (nth order reaction with autocatalysis). The program is started by clicking on the green triangular button (Figure 4.22). The table of initial parameters that should be verified then appears (Figure 4.24).

In the Const. column, the option ‘false’ is set for the parameters that should be varied and the option ‘true’ is chosen for the parameters that remain constant. The three columns to the right of this column are intended for imposing constraints on the selected values. The computation results are presented in Figure 4.25.

Fig. 4.25: Computation results for the Cnb function.

Fig. 4.25: Computation results for the Cnb function.

The calculation for the Fn function is performed analogously. The results are shown in Figure 4.26. To compare the parameter values obtained for three statistically equivalent functions, these parameters should be evaluated with allowance for correlation of the initial data. In the Statistics menu (Figure 4.25), the Durbin–Watson test value is presented (4.3) which should be added to the calculated tS error, where t is the Student’s test for a given number of degrees of freedom, and S is the mean-square deviation. Thus, we obtain the following set of parameters and their errors (Table 4.1) for all statistically equivalent functions.

Fig. 4.26: Computation results for the Fn function.

Fig. 4.26: Computation results for the Fn function.

The following conclusions can be drawn from Table 4.1: first, the autocatalysis parameters for the Bna and CnB functions are almost zero, that is, all is reduced to the Fn function. Second, the error of the Arrhenius parameters for Fn are minimal. Hence, the calcium oxalate dehydration is better described by the nth-order function. The reaction order can be considered to be 1/3, that is, the process is described by the “contracting sphere equation.” This means that the sample consists of spherical particles of the same size and that the dehydration is a homothetic process, that is, the particles during decomposition undergo a self-similar decrease in size. Such a mechanism is inherent in thermolysis of inorganic crystal hydrates. Thus, the problem of CaC2O4·H2O dehydration macrokinetics can be thought of as solved.

Function code log A E, kJ/mol Reaction order n log K cat 1 Exp a1
Bna 6.8±0.6 73±6 0.34±0.25 0.06±0.14
CnB 6.8±1 74±8 0.41±0.48 -0.45±1.9
Fn 7.0±0.1 75.2±0.8 0.34±0.25

Tab. 4.1: Kinetic parameters of the CaC2O4∙H2O dehydration process

Tab. 4.1: Kinetic parameters of the CaC2O4∙H2O dehydration process

The project created in the NETZSCH Thermokinetics software is saved by clicking on the common button (Figure 4.27).

Fig. 4.27: Saving the project in the NETZSCH Thermokinetics software.

Fig. 4.27: Saving the project in the NETZSCH Thermokinetics software.

The output data of the NETZSCH Thermokinetics software can be exported as an ASCII file by clicking on a button in the Nonlin regr. toolbar (Figure 4.28). In the opened dialog box, the user can select data that must be exported: kinetic parameters, statistics, and so on. (Figure 4.29).

Fig. 4.28: Export of the results from NETZSCH Thermokinetics as an ASCII file.

Fig. 4.28: Export of the results from NETZSCH Thermokinetics as an ASCII file.

Fig. 4.29: Selection of data for export as an ASCII file.

Fig. 4.29: Selection of data for export as an ASCII file.