Research Article  Open Access
Velyana Georgieva, Lyubomir Vlaev, Kalinka Gyurova, "NonIsothermal Degradation Kinetics of CaCO_{3} from Different Origin", Journal of Chemistry, vol. 2013, Article ID 872981, 12 pages, 2013. https://doi.org/10.1155/2013/872981
NonIsothermal Degradation Kinetics of CaCO_{3} from Different Origin
Abstract
Thermogravimetric studies on two varieties of calcium carbonate, namely, analytical reagentgrade and in situ product from thermal degradation of calcium oxalate monohydrate, were carried out at four rates of linear increase of the temperature. The kinetics and mechanism of their solidstate thermal decomposition reaction were evaluated from the TG data using four calculation procedures and isoconversion method, as well as 27 mechanism functions. The comparison of the results obtained with these calculation procedures showed that they strongly depend on the selection of proper mechanism function for the process. Therefore, it is very important to determine the most probable mechanism function. In this respect the isoconversion calculation procedure turned out to be more appropriate. In the present work, the values of apparent activation energy E, preexponential factor A in Arrhenius equation, as well as the changes of entropy , enthalpy , and free Gibbs energy for the formation of the activated complex from the reagent are calculated. All calculations were performed using programs compiled by ourselves.
1. Introduction
The thermal decomposition of calcium carbonate has been intensively studied over the years [1–12]. It is explained by great technological importance due to his various industrial applications such as manufacturing of lime and additives or fillers in medicine, cosmetics, foods, plastics, printing inc and an apparent simplicity of the reaction. It is a high temperature reversible decomposition reaction involving a relatively large massloss associated with the evolution of carbon dioxide according to the reaction: Admittedly the studying and measurements are complicated by the reversibility of this reaction and its strongly endothermic nature; however, these complications should only add interest to the resulting comparison.
Calcium carbonate mineral filler has been one of the most popular fillers used in the thermoplastic industry [13, 16]. can be generally found in three distinct crystalline structures, that is, calcite, aragonite, and vaterite [17, 18]. Calcite is the most stable and most commonly found in the nature. Aragonite, however, can only be found in precipitate or seashells whereas vaterite is found from synthesized and does not occur naturally [17, 18]. Generally, aragonite has higher density and hardness than calcite and vaterite which makes it a valuable inorganic material that can be used as a filler for armature plastic, rubber, paper, glass fiber, print ink, paint pigment, and composite [13, 16, 17]. From another side, very interesting is the comparative kinetics studying of the thermal degradation of in situ generated calcium carbonate, obtained from the calcium oxalate monohydrate thermal degradation [14, 15].
Apart from mechanical properties, thermal stability is also one of the most important factors for the processing of polymeric materials. Thermal decomposition kinetics of the polymeric materials depend from the thermal stability of the used fillers, which may be characterized with its kinetic parameters such as reaction order , activation energy , and preexponential or frequency factor in the Arrhenius equation [19–21]. According to the results obtained from the different authors the apparent magnitudes of the Arrhenius parameters range from 110 to 3800 kJ for the activation energy, and from 10^{2} to for the frequency factor [1–22]. The reason for widely varying values are different—the shape of the crystalline phase of used calcium carbonate, the mean size of the particles, sample mass, isothermal or nonisothermal heating, heating rate, static or dynamic atmosphere around the studied sample, furnace atmosphere, the partial pressure of carbon dioxide, and the used calculation procedure [2, 14, 21, 23, 24].
The aim of the present paper is to compare the parameters characterizing the kinetics of nonisothermal degradation of analytical reagentgrade commercial sample of and insitu produced calcium carbonate from , obtained on the basis of thermogravimetrical data using different kinetic equations and calculation procedures.
2. Experimental
2.1. Material and Measurement
The thermal decomposition on two calcium carbonate, namely, analytical reagentgrade commercial sample and in situ produced calcium carbonate from calcium oxalate monohydrate, were carried out. The samples of (puriss, Fluka) and (puriss, Aldrich) were used after vigorously grounding in agate vibration mortar. The thermogravimetrical measurements were carried out in a flow of nitrogen (99.999%) at a rate of 25 min^{−1} under nonisothermal conditions on an instrument STA 449 F3 Jupiter (Nietzsch, Germany) with its high temperature furnace. Samples of about mg mass were used for the experiments varied out at hearing rates of 3, 6, 9, and 12 K min^{−1} up to 1000°C. The samples were loaded without pressing into an open 6 mm diameter and 3 mm high platinum crucible, without using of a standard reference material. The TG, DSC, and DTG curves were recorded simultaneously with 0.1 mg sensitivity.
The crystal structures of studied samples were determined by Xray diffraction (XRD). The XRD intensity of the samples were collected from a JEOL/JDX 3530 diffractometer equipped with generator ( nm) at voltage of 32 kV, Ni filter, and a filament current of 30 mA. The scans were taken in the range of the diffraction angle = 5–60°.
2.2. Mathematical Background
The kinetics of such solidstate reactions is described by various equations taking into account the special features of their mechanisms. The reaction rate can be expressed through the degree of conversion according to the formula: where , , and are the initial, final, and current sample mass at the moment , respectively. Generally, the kinetic equation of the process can be written as follows [1, 3, 13–15, 19–21, 23]
The temperature dependence of the rate constant for the process is described by the Arrhenius equation: where is the preexponential factor, is the apparent activation energy, is the absolute temperature, and is the gas constant. Substitution of (3) in (2) gives
When the temperature increases at a constant rate, therefore:
The conversion function for a solidstate reaction depends on the reaction mechanism and can generally be considered to be as follows: where , , and are empirically obtained exponent factors, one of them always being zero [25].
The combinations of different values of , , and make it possible to describe various probable mechanisms.
After substitution of (7) in (6), separation of variables and integration, the following general equation was obtained
The solutions of the left hand side integral depend on the explicit expression of the function and are denoted as . The formal expressions of the functions depend on the conversion mechanism and its mathematical model [13–15]. The latter usually represents the limiting stage of the reaction—the chemical reactions; random nucleation and nuclei growth; phase boundary reaction or diffusion. Algebraic expressions of functions of the most common reaction mechanisms operating in solidstate reactions are presented in Table 1 [13–15].

Several authors [1, 3, 14, 26, 27] suggested different ways to solve the right hand side integral. For the present study, one calculation procedure was based on Coats and Redfern equation [1]. Data from TG and DTG curves in the decomposition range were used to determine the kinetic parameters of the process in all used calculation procedures. The integral method of Coats and Redfern has been mostly and successfully used for studying of the kinetics of dehydration and decomposition of different solid substances [1, 14]. The kinetic parameters can be derived using a linear form of modified Coats and Redfern equation: where is an integral form of the conversion function, the expression of which depends on the kinetic model of the occurring reaction. If the correct function is used, a plot of against should give a straight line from which the values of the activation energy and the preexponential factor in Arrhenius equation can be calculated.
Later, several authors [14, 28–31] suggested different solutions of the temperature integral in (8), sustaining the opinion that this increases the precision of the kinetic parameters being calculated. For instance, MadhysudananKrishnanNinan [28, 29] suggested the following equation:
Tang et al. [30] suggested another kinetic equation: and Wanjun et al. [31] suggested the equation:
Equations (9)–(12) imply that there would be differences in the calculated values of the activation energy and preexponential factor even when the same function is used. To find which calculation procedure would turn out to be the most suitable for the calculations, they were estimated by the criterion “the best” correlation coefficient of the linear regression for (9)–(12). The advantage of these equations is that the values of and can be calculated on the basis of single rate thermogravimetric curves and the type of the most probable mechanism function of the studied reaction can be determined. For the calculations of the kinetic parameters a computer program was developed for all the data manipulations.
The second approach used for the calculations was based on multiple rates thermogravimetric curves and socalled isoconversion calculation procedures [13, 15, 21, 23, 24, 32–35].
2.3. Isoconversional Methods
Isoconversion methods employ multiple temperature programs (e.g., different heating rates) in order to obtain data on varying rates at a constant extent of conversion. Thus, isoconversion methods allow complex (i.e., multistep) processes to be detected via a variation of with . As has results from some critical analyses [15, 23, 24, 32–35], the correct determination of nonisothermal kinetic parameters involves the use of experimental data recorded at several heating rates. These data have allowed applying the isoconversion (modelfree) methods in accessing the activation energy on the conversion degree that can be correlated with the investigated process mechanism. These methods are recommended from ICTAC kinetics committee for performing kinetic computations on thermal analytical data [24].
2.3.1. Calculation of Activation Energy by Iterative Procedure
Thermal decomposition kinetic of OzawaFlynnWall (OFW) is an isoconversion or modelfree method, which can calculate the values of the activation energy E and frequency factor A in Arrhenius equation through following equation [13, 27, 32, 33, 35]: and KissingerAkahiraSunose (KAS) equation [32–35]:
These two methods of plotting a linear regressive curve were used at fraction conversion , different heating rates , and sample masses. The plots of versus (Equation (13)) and versus (see (14)) have been proved to give the values of the apparent activation energies for the different stages of decomposition of and at different values of . According these equations, the reaction mechanism and the shape of function cannot affect the calculation of the activation energies of the different stages. Iterative procedure was used to calculate the values of approximating the exact value, according to the next equations [32, 35]: Because and where is expressed by the fourth Senum and Yang approximation formulae [32, 34, 35]: where and is equal to [32, 35]
The iterative procedure performed involved the following steps. (i) Assume or to estimate the initial value of the activation energy . The conventional isoconversional methods stop the calculation at this step. (ii) Using , calculate a new value of for the activation energy from the plot of versus or versus . (iii) Repeat step (ii), replacing with . When kJ mol^{−1}, the last value of was considered to be the exact value of the activation energy of the reaction. These plots are model independent since the estimation of the apparent activation energy does not require the selection of particular kinetic model (type of function). Therefore, the activation energy values obtained by this method are usually regarded as more reliable than these obtained by a single TG curve.
2.3.2. Determination of the Most Probably Mechanism Function
The following equation was used to estimate the most correct reaction mechanism, that is, function, [32, 34]:
Plotting versus and using a linear regressive of least square method, if the mechanism studied conforms to certain function, the slope of the straight line should be equal to −1.0000 and the linear correlation coefficient should be equal to unity. The values of and do not influence the shape of the most correct reaction mechanism function determined. For determination of the most probable mechanism function, the values of the conversion corresponding to multiple rates taken at the same temperature were put into the left side of (19) and all the thirtyfive types of mechanism functions presented in Table 1 were tested. The slope and correlation coefficient were obtained from the plot of verus . The most probable mechanism function was assumed to be the one for which the value of the slope of the straight line was closest to −1.0000 and the correlation coefficient was highest. If several functions comply with this requirement, the values of conversion corresponding to multiple rates at the same temperature were applied to calculate the probable mechanism by the same method. Thus, the most probable mechanism function was considered to be the one for which the slope of the straight line was closest to −1.0000 and the linear correlation coefficient was near to unity.
2.3.3. Calculation of Preexponential Factor in Arrhenius Equation
The preexponential factor A can be estimated from the intercept of the plots of (15), inserting the most probable function determined. All calculations were performed using a programs compiled by ourselves.
The values of the preexponential factor in Arrhenius equation for solid phase reactions are expected to be in a wide range (six or seven orders of magnitude), even after the effect of surface area is taken into account [15, 36–40]. For first order reactions, the preexponential factor may vary from 10^{5} to 10^{18} s^{−1}. The low factors will often indicate a surface reaction, but if the reactions are not dependent on surface area, the low factor may indicate a “tight” complex. The high factors will usually indicate a “loose” complex [15, 36]. Even higher factors (after correction for surface area) can be obtained for complexes having free translation on the surface. Since the concentrations in solids are not controllable in many cases, it would have been convenient if the magnitude of the preexponential factor indicated for reaction molecularity. However, this appears to be true only for nonsurfacecontrolled reactions having low (<) preexponential factors. Such reactions (if elementary) can only be bimolecular [36, 37].
From the theory of the activated complex (transition state) of Eyring [15, 36–40], the following general equation may be written where is the Neper number; :transition factor, which is unity for monomolecular reactions; :Boltzmann constant; :Plank constant, and is the peak temperature of DTG curve. The change of the entropy may be calculated according to the formula Since the changes of the enthalpy and Gibbs free energy for the activated complex formation from the reagent can be calculated using the wellknown thermodynamical equation: The values of , , and were calculated at ( is the DTG peak temperature at the corresponding stage), since this temperature characterizes the highest rate of the process, and therefore, is its important parameter.
3. Results and Discussion
Figure 1 shows the Xray diffraction patterns of reagentgrade and producedin situ as a result from the thermal degradation of calcium oxalate monohydrate.
The Xray diffraction pattern of obtained in situ from calcium oxalate monohydrate corresponds to calcite which has hexagonal crystal lattice with the following parameters of the unit cell: a—4.99013 Å, b—4.99013 Å, c—17.0690 Å, α—90°, β—90°γ—120°, and volume—368.09 Å^{3}. Beside the reflexes for calcite observed in the Xray diffraction pattern of the initial reagentgrade (Fluka) used, two other reflexes characteristic for aragonite were registered which have orthorhombic crystal lattice with the following parameters of the unit cell: a—4.9632 Å, b—7.9673 Å, c—5.7442 Å, α—90°, β—90°, γ—90°, and volume227.15 Å^{3}. The mean crystalline size of the samples studied was calculated according to the DebyeScherrer formula: where is a dimensionless constant or shape factor that may range from 0.89 to 1.39 depending on the specific geometry of the scattering objects (in our case equal to 0.9), is the Xray wavelength (0.154 nm) for Cu anode, is the peak width at half the maximum intensity reflex in radians, and is the Bragg angle. For the first sample, the mean size of the crystallites was 920 Å and for the second one −2310 Å, respectively.
The TG and DTG curves of thermal degradation of reagentgrade and generated in situ from calcium oxalate monohydrate are presented in Figure 2.
In the DTG curve of CaCO_{3}, there is only one peak at 742.8°C corresponding to the process of substance decarbonation. The thermal decomposition of proceeded in three reaction steps observed at increased temperature, corresponding to the splitting of , , and CO_{2}, respectively [15]:
The stages of dehydration, decarbonylation, and decarbonation of are all well separated, though the kinetic characteristics are affected by the presence of , , and , as well as by the reaction conditions, including heating rate, sample mass, sample size, and type of sample container [13, 15]. In the DTG curve of decomposition, there are three peaks at 180.3, 488.7, and 749.5°C corresponding to the dehydration, decarbonylation and decarbonation processes, respectively. In both cases, the massloss values measured were all in good agreement with the value of 44% calculated according Stage III.
The TG curves of thermal decomposition of the initial reagentgrade CaCO_{3} obtained at different heating rates are presented in Figure 3
Using CoatsRedfern calculation procedure and functions with different values of ,the values of the activation energy , and frequency factor were calculated at four heating rates of thermal decomposition of the reagentgrade . The dependencies of the coefficient of linear regression on the values of n at the used heating rates are presented in Figure 4.
As can be seen from Figure 4, each curve had a local maximum. The curves presented can be successfully described by empirical polynomials of second power and different coefficients. Differentiating and nullifying these polynomials, the value of n at which has maximum value can be calculated. Despite the heating rate, it turned out that the maximal value of was obtained at . It means that, in all cases, the most appropriate mechanism function describing the kinetics of decomposition of in integral form is and in differential form which belongs to the mechanism of phase boundary reaction (contracting sphere). The same tendency was established for the other calculation procedures. The values of , , and for the thermal decomposition of reagentgrade obtained at different heating rates by different calculation procedures are presented in Table 2 for comparison.

As can be seen from Table 2, regardless of the calculation procedure used, the values of E and A regularly increased with the heating rate while these of decreased. In our opinion, the tendency of increase of the activation energy values with the increase of the heating rate observed was due to the process sensibility to the heat transfer to the reaction interface. The average values of the activation energy obtained from the calculation procedures used was 326.4 2.5 kJ/mol.
The same approach was used for the studies on the kinetics of the thermal decomposition ofin situ produced . It was established that the most appropriate mechanism function is the same— and in differential form which belongs to the mechanism of phase boundary reaction or contracting sphere. The values of , , and for the thermal decomposition of in situproduced obtained at different heating rates and by different calculation procedures are shown for comparison in Table 3.

At the thermal decomposition of in situ produced is observed the same tendency of increasing of the values of the activation energy and frequency factor with the increasing of the heating rate, as well as at the thermal decomposition of reagentgrade . In this case, however, the values of and are considerably lower at the same heating rate. According to us, the reason for this is the fact that the reagent had higher reactivity at the in situ generated CaCO_{3} because it was in “status nascendi” and the mean size of the crystallites according to the Xray diffraction data was smaller.
Figure 5 presents the linear relationships between and for the thermal degradation of reagentgrade according to OFW calculation procedure.
At constant heating rate, the kinetic parameters can be calculated by using the slope and cutoff from the ordinate axis. The results obtained are presented in Table 4.

The linear relationships between and for the thermal degradation of reagentgrade drawn according to KAS calculation procedure are presented in Figure 6.
The activation energy was directly evaluated from the slopes of these plots and frequency factor —from the cutoff from the ordinate axis. The results obtained are summarized in Table 4.
As can be seen from Table 4 for all the cases shown, the values of the calculated activation energy regularly increase with the increase of the degree of conversion . According to some authors [20], three rate controlling processes are possible: heat transfer to the reaction interface, chemical reaction, and diffusion through the product layer. In our case, diffusion control of the reaction seems to be plausible; therefore, as the thickness of the products layer increases, the resistance against CO_{2} diffusion increases, and, as a consequence, the activation energy also increases. In corroboration of the suggestion that diffusion through the product layer is the rate controlling process comes the fact that the value of the activation energy of the thermal decomposition of CaCO_{3} in vacuum is smaller than that obtained in dry argon atmosphere [22]. The values for the activation energy of the thermal decomposition of the reagentgrade CaCO_{3} obtained were close to these reported in [20, 22].
Comparing the values of the kinetic parameters calculated on the base of single curve methods and the isoconversion procedures, it can be seen that the latter were twice smaller. Obviously, the type of the calculation procedure used had significant effect on the values of the kinetic parameters of certain reaction which can explain the great difference between the values of the calculated activation energy E (from 110 to 3800 kJ mol^{−1}) and from 10^{2} to 10^{157} s^{−1} for the frequency factor [1–22]. According to ICTAC kinetics committee recommendations for performing kinetic computations on thermal analysis data, more appropriate calculation procedures are these based on the isoconversion procedure [24].
Figure 7 shows the straight lines drawn according to (19) using different functions.
As can be seen from Figure 7, the lines have different slopes. Line 1 corresponds to D_{4}, line 2 to (), and line 3 to D_{3} mechanism function, respectively. The values of the slope and obtained for different mechanism functions are presented in Table 5 for comparison.

As can be seen from Table 5 for all the cases studied, the values of were very high but the slope of the line was closest to −1.0000 only for mechanism (−1.0086 and −1.0102 ). For the D_{3} mechanism, it was higher and for D_{4} smaller than −1.0000. This was the reason to be concluded that the most appropriate mechanism describing the thermal degradation of is the phase boundary reaction —contracting sphere or spherical symmetry. The values of the parameters characterizing the kinetic of thermal decomposition of CaCO_{3} from different origins are presented in Table 6.

The change of entropy for the formation of the activated complex from the reagent reflects how near the system is to its thermodynamic equilibrium. Lower activation entropy means that the material has just passed through some kind of physical or chemical rearrangement of the initial structure, bringing it to a state near its own thermodynamic equilibrium. In this situation, the material shows little reactivity, increasing the time necessary to form the activated complex. On the other hand, when highactivation entropy values are observed, the material is far from its own thermodynamic equilibrium. In this case, the reactivity is higher and the system can react faster to produce the activated complex and, consequently, short reaction times are observed. In particular, for example, the negative values of would indicate that the formation of activated complex is connected with the decrease of entropy, that is, the activated complex is “more organized” structure compared to the initial substance and such reactions are classified as “slow” [41]. The more negative values of for the in situ produced CaCO_{3} showed that its structure is far from its own thermodynamic equilibrium compared to reagentgrade .
Using (22) and (23), the values of the change of enthalpy and the Gibbs free energy for the formation of the activated complex from the reagent can be calculated. The change of the activation enthalpy shows the energy difference between the reagent and activated complex. If this difference is small, the formation of the activated complex is favored because the potential energy barrier is low. The change of the Gibbs free energy reveals the total energy increase of the system at the approach of the reagents and the formation of the activated complex. As can be seen from Table 6, they are practically equal for both studied compounds, but there are significant differences between the values of and .
The linear relationship between and for all values of and corresponding values of obtained for both samples with different calculation procedures is presented in Figure 8.
The kinetic compensation effect was developed in order to determine the effect from the different specimens or experimental conditions when calculating the change of the activation energy. Under these circumstances, the kinetic parameter of preexponential factor A would vary with the activation energy . A good linear dependence between and E was observed for these parameters, obtained from all the calculation procedures used. The straight line presented in Figure 8 is described with the following empirical equation: According to the some authors [39, 41–44] the socalled kinetic compensation effect, isokinetic effect, or θrule may be described with the equation: On the basis of (25) and (26) the values of and may be calculated. In our case they are 0.059 min^{−1} and 983 K, respectively. This temperature belongs in the temperature interval of the beginning and the end of the thermal decomposition of the both studied samples of .
4. Conclusion
Thermogravimetric studies on two varieties of calcium carbonate namely, analytical reagentgrade andin situproduced from thermal degradation of calcium oxalate monohydrate, were carried out at four rates of linear increase of the temperature. The kinetics and mechanism of their solidstate thermal decomposition reaction were evaluated from the TG data using four calculation procedures and isoconversion method, as well as 27 mechanism functions. For all the sets of results, the and values ofin situ generated were marginally lower. On this basis, we can infer that the sample history has a definite effect on the computed values of and , even though the magnitude of the effect is very small. The observation can be explained as follows. In the phase boundary reaction the surface nucleation is instantaneous and hence the rate controlling process is the movement of the interface towards the interior. Thus, even though the numbers of nucleation sites for the both types of samples are likely to be largely different, it will have no effect on the reaction rate. However, marginal decrease of the values of the kinetics parameters for thein situproduced can be attributed to the smaller average crystallite size and because this sample is in “status nascendi” and for this reason more reactive.
Independently from the high values of the correlation coefficient of linear regression obtained on the basis of single kinetic TGcurves, the values of , calculated according to isoconversion methods are two times less and strongly depends on the heating rate and degree of conversion. For this reason it may be concluded that will be more correct if the values of obtained on isoconversion calculation procedures are more representative.
References
 A. W. Coats and J. P. Redfern, “Kinetic parameters from thermogravimetric data,” Nature, vol. 201, no. 4914, pp. 68–69, 1964. View at: Publisher Site  Google Scholar
 Z. Adonyi, “Correlation between kinetic constants and parameters of differential thermogravimetry, in the decomposition of calcium carbonate,” Periodica Polytechnica, vol. 11, pp. 325–336, 1967. View at: Google Scholar
 T. V. Mesteryakova and N. D. Topor, “Determination of thermal stability solid inorganic compounds and minerals using thermoravimetric method,” Vestnik Moskovskogo Universiteta, vol. 3, pp. 73–78, 1967 (Russian). View at: Google Scholar
 J. Zsakó and H. E. Arz, “Kinetic analysis of thermogravimetric data—VII. Thermal decomposition of calcium carbonate,” Journal of Thermal Analysis, vol. 6, no. 6, pp. 651–656, 1974. View at: Publisher Site  Google Scholar
 P. K. Gallagher and D. W. Johnson, “Kinetics of the thermal decomposition of CaCO_{3} in Co_{2} and some observations on the kinetic compensation effect,” Thermochimica Acta, vol. 14, no. 3, pp. 255–261, 1976. View at: Google Scholar
 J. P. Elder and V. B. Ready, “The kinetics of the thermal degradation of calcium carbonate,” Journal of Thermal Analysis, vol. 31, no. 2, pp. 395–405, 1986. View at: Publisher Site  Google Scholar
 Y. Wang and W. J. Thomson, “The effect of sample preparation on the thermal decomposition of CaCO_{3},” Thermochimica Acta, vol. 255, no. C, pp. 383–390, 1995. View at: Google Scholar
 Y. Wang and W. J. Thomson, “The effects of steam and carbon dioxide on calcite decomposition using dynamic Xray diffraction,” Chemical Engineering Science, vol. 50, no. 9, pp. 1373–1382, 1995. View at: Google Scholar
 B. V. L'vov, “Mechanism of thermal decomposition of alkalineearth carbonates,” Thermochimica Acta, vol. 303, no. 2, pp. 161–170, 1997. View at: Google Scholar
 M. Maciejewski, “Computational aspects of kinetic analysis. Part B: the ICTAC Kinetics Project—the decomposition kinetics of calcium carbonate revisited, or some tips on survival in the kinetic minefield,” Thermochimica Acta, vol. 355, no. 12, pp. 145–154, 2000. View at: Publisher Site  Google Scholar
 Y. T. Zhao, T. S. Sun, and B. Sun, “The thermal decomposition of calcium carbonate,” Chinese Chemical Letters, vol. 12, no. 8, pp. 745–746, 2001. View at: Google Scholar
 B. V. L'Vov, “Mechanism and kinetics of thermal decomposition of carbonates,” Thermochimica Acta, vol. 386, no. 1, pp. 1–16, 2002. View at: Publisher Site  Google Scholar
 S. Thumsorn, K. Yamada, Y. W. Leong, and H. Hamada, “Development of cockleshedderiwed CaCO_{3} for flame retardancy of recycled PET/recycled PP blend,” Materials Sciences and Applications, vol. 2, pp. 59–69, 2011. View at: Google Scholar
 K. N. Ninan, K. Krishnan, and V. N. Krishnamurthy, “Kinetics and mechanism of thermal decomposition of insitu generated calcium carbonate,” Journal of Thermal Analysis, vol. 37, no. 7, pp. 1533–1543, 1991. View at: Publisher Site  Google Scholar
 L. Vlaev, N. Nedelchev, K. Gyurova, and M. Zagorcheva, “A comparative study of nonisothermal kinetics of decomposition of calcium oxalate monohydrate,” Journal of Analytical and Applied Pyrolysis, vol. 81, no. 2, pp. 253–262, 2008. View at: Publisher Site  Google Scholar
 R. Gachter and H. Muller, Plastics Aditives Handbook, Hanser, Munich, Germany, 1990.
 N. H. de Leeuw and S. C. Parker, “Surface structure and morphology of calcium carbonate polymorphs calcite, aragonite, and vaterite: an atomistic approach,” Journal of Physical Chemistry B, vol. 102, no. 16, pp. 2914–2922, 1998. View at: Google Scholar
 Y. Mori, T. Enomae, and A. Isogai, “Preparation of pure vaterite by simple mechanical mixing of two aqueous salt solutions,” Materials Science and Engineering C, vol. 29, no. 4, pp. 1409–1414, 2009. View at: Publisher Site  Google Scholar
 B. V. L'Vov, L. K. Polzik, and V. L. Ugolkov, “Decomposition kinetics of calcite: a new approach to the old problem,” Thermochimica Acta, vol. 390, no. 12, pp. 5–19, 2002. View at: Publisher Site  Google Scholar
 C. RodriguezNavarro, E. RuizAgudo, A. Luque, A. B. RodriguezNavarro, and M. OrtegaHuertas, “Thermal decomposition of calcite: mechanisms of formation and textural evolution of CaO nanocrystals,” American Mineralogist, vol. 94, no. 4, pp. 578–593, 2009. View at: Publisher Site  Google Scholar
 N. Koga, Y. Yamane, and T. Kimura, “Thermally induced transformations of calcium carbonate polymorphs precipitated selectively in ethanol/water solutions,” Thermochimica Acta, vol. 512, no. 12, pp. 13–21, 2011. View at: Publisher Site  Google Scholar
 J. P. Sanders and P. K. Gallagher, “Kinetic analyses using simultaneous TG/DSC measurements. Part I: decomposition of calcium carbonate in argon,” Thermochimica Acta, vol. 388, no. 12, pp. 115–128, 2002. View at: Publisher Site  Google Scholar
 P. Budrugeac, “An iterative modelfree method to determine the activation energy of nonisothermal heterogeneous processes,” Thermochimica Acta, vol. 511, no. 12, pp. 8–16, 2010. View at: Publisher Site  Google Scholar
 S. Vyazovkin, A. K. Burnham, J. M. Criado, L. A. PérezMaqueda, C. Popescu, and N. Sbirrazzuoli, “ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data,” Thermochimica Acta, vol. 520, no. 12, pp. 1–19, 2011. View at: Publisher Site  Google Scholar
 J. Šesták and G. Berggren, “Study of the kinetics of the mechanism of solidstate reactions at increasing temperatures,” Thermochimica Acta, vol. 3, no. 1, pp. 1–12, 1971. View at: Google Scholar
 H. H. Horowitz and G. Metzger, “A new analysis of thermogravimetric traces,” Analytical Chemistry, vol. 35, no. 10, pp. 1464–1468, 1963. View at: Google Scholar
 T. Ozawa, “A new method of analyzing thermogravimetric data,” Bulletin of the Chemical Society of Japan, vol. 38, pp. 1881–1886, 1965. View at: Google Scholar
 P. M. Madhusudanan, K. Krishnan, and K. N. Ninan, “New approximation for the p(x) function in the evaluation of nonisothermal kinetic data,” Thermochimica Acta, vol. 97, pp. 189–201, 1986. View at: Publisher Site  Google Scholar
 P. M. Madhusudanan, K. Krishnan, and K. N. Ninan, “New equations for kinetic analysis of nonisothermal reactions,” Thermochimica Acta, vol. 221, no. 1, pp. 13–21, 1993. View at: Google Scholar
 W. Tang, Y. Liu, H. Zhang, and C. Wang, “New approximate formula for Arrhenius temperature integral,” Thermochimica Acta, vol. 408, no. 12, pp. 39–43, 2003. View at: Publisher Site  Google Scholar
 T. Wanjun, L. Yuwen, Z. Hen, W. Zhiyong, and W. Cunxin, “New temperature integral approximate formula for nonisothermal kinetic analysis,” Journal of Thermal Analysis and Calorimetry, vol. 74, no. 1, pp. 309–315, 2003. View at: Publisher Site  Google Scholar
 L. Liqing and C. Donghua, “Application of ISOtemperature method of multiple rate to kinetic analysis: dehydration for calcium oxalate monohydrate,” Journal of Thermal Analysis and Calorimetry, vol. 78, no. 1, pp. 283–293, 2004. View at: Publisher Site  Google Scholar
 P. Budrugeac and E. Segal, “Some methodological problems concerning nonisothermal kinetic analysis of heterogeneous solidgas reactions,” International Journal of Chemical Kinetics, vol. 33, no. 10, pp. 564–573, 2001. View at: Publisher Site  Google Scholar
 Z. Gao, I. Amasaki, and M. Nakada, “A description of kinetics of thermal decomposition of calcium oxalate monohydrate by means of the accommodated Rn model,” Thermochimica Acta, vol. 385, no. 12, pp. 95–103, 2002. View at: Publisher Site  Google Scholar
 T. T. Su, H. Jiang, and H. Gong, “Thermal stabilities and the thermal degradation kinetics of poly(εcaprolactone),” Polymer, vol. 47, no. 4, pp. 398–403, 2008. View at: Publisher Site  Google Scholar
 D. V. Sokolskii and V. A. Druz, Introduction in Theory Heterogeneous Catalysis, Vyshaya Shkola, Moscow, Russia, 1981.
 J. Šesták, Thermodynamical Poperties of Slids, Academia, Prague, Czech Republic, 1984.
 F. Tudos and P. K. David, “Comments on the eringinterpretation of the compensation effect,” Polymer Degradation and Stability, vol. 50, no. 2, pp. 159–161, 1995. View at: Google Scholar
 R. Bigda and A. Mianowski, “Influence of heating rate on kinetic quantities of solid phase thermal decomposition,” Journal of Thermal Analysis and Calorimetry, vol. 84, no. 2, pp. 453–465, 2006. View at: Publisher Site  Google Scholar
 B. Boonchom and S. Puttawong, “Thermodynamics and kinetics of the dehydration reaction of ${\text{FePO}}_{4}\xb72{\text{H}}_{2}\text{O}$,” Physica B, vol. 405, no. 9, pp. 2350–2355, 2010. View at: Publisher Site  Google Scholar
 A. Atanassov, S. Genieva, and L. Vlaev, “Study on the thermooxidative degradation kinetics of tetrafluoroethyleneethylene copolymer filled with rice husks ash,” Polymer, vol. 49, no. 6, pp. 541–554, 2010. View at: Publisher Site  Google Scholar
 V. M. Gorbachev, “The compensation effect in the kinetics of the thermal decomposition of calcium carbonate,” Journal of Thermal Analysis, vol. 9, no. 1, pp. 121–123, 1976. View at: Publisher Site  Google Scholar
 H. S. Ray, “The kinetic compensation effect in the decomposition of calcium carbonate,” Journal of Thermal Analysis, vol. 24, no. 1, pp. 35–41, 1982. View at: Publisher Site  Google Scholar
 S. C. Turmanova, S. D. Genieva, A. S. Dimitrova, and L. T. Vlaev, “Nonisothermal degradation kinetics of filled with rise husk ash polypropene composites,” Express Polymer Letters, vol. 2, no. 2, pp. 133–146, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Velyana Georgieva et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.