Kinetic triplet and parameters of activated complex for the process of thermal decomposition of calcite (Eq. 33 in Chapter 2).
This chapter provides the information about the concept of effective mass and effective velocity of the activated complex and its connection to the transition state theory. Therefore, these parameters are of essential importance for the field of homogenous as well as heterogeneous kinetics. They also prove to be useful for the calculation of many other properties of activated state, such as momentum, energetic density, mass flux, etc., as will be demonstrated on the example of thermal decomposition of calcite and aragonite. Since the activation energy and the momentum of activated state enable to complete the characterization of motion of this instanton (pseudoparticle) alongside the reaction coordinate, these parameters can be then considered as two quantum numbers of activated complex. The quantum numbers of activated state, that is, the activation energy and momentum, also explain the relation of activated complex to Planck energy, length and time, as well as to the Gravitational constant. This idea was also applied to derive the wave function of activated complex pseudoparticle, which is affected by the isotopic composition of the sample and polymorphism as well. Furthermore, the findings introduced in this chapter enable to derive and propose the modified Kissinger equation and experimental solution for the approximation parameter in the Doyle equation of temperature integral.
- activated complex
- activated state
- effective mass
- effective velocity
- mean lifetime
- group velocity
- phase velocity
- transition state theory
- activation energy
- quantum numbers
- Schrödinger equation
- wave function
- modified Kissinger equation
- temperature integral
- approximation parameter
- thermal decomposition of solids
1. Chapter introduction and basic assumptions
The definition of effective or reduced mass of activated complex (
where is the thermodynamic (absolute Kelvin scale) temperature of peak measured by thermal analysis (TA [3, 4, 5, 6, 7], such as DTA1, DTG2, DSC3…) with the heating rate (
where is the mass of activated molecule and the constant
More general expression for Eq. 1 is given by the following formula:
where is the constant exponent (power) of power function (please refer to Eq. 12) with the scaling factor :
For example, the DTG peak (numerical derivation of thermogravimetric experiment) for the thermal decomposition of calcite (Eq. 33 in Chapter 2) heated with the rate of 1 K∙min−1 under inert atmosphere of nitrogen is shown in Figure 1.
Since the activation energy of this reaction (Eq. 33 in Chapter 2) is already known (please refer to Figure 2(a) in Chapter 2), it is possible to calculate the mechanism of the process using the formula8 of Augis and Bennett :
Before we continue, it is also interesting to solve the nature of apparent change of kinetic coefficient with heating rate. This behavior results from the effect of heating rate on the peak temperature and full width at half maximum of peak (Eq. 7), while the value of activation energy stays constant (please check this assumption with regard to the discussion of Eqs. 19–24 in Chapter 1 and Footnote 44). Both dependencies can be accurately approximated by the power law:
where and are the scale factors (coefficients) and
where the kinetic coefficient, the fraction () as well as the correction term are all dimensionless and . It can also be derived, that the activation energy is:
Using the experimental data for the abovementioned example of the process of thermal decomposition of calcite, the dependence of , and on heating rate is shown in Figure 2(a).
It is obvious that the power dependence of the peak temperature (
where is the peak temperature calculated (Eq. 12) for the heating rate equal to the Euler number.
it can be treated with similar manner. Analogically to Eq. 14, it can be written:
Eqs. 8, 12, 14 and 17 then provide the set of important parameters, which can have interesting utilization in many branches of heterogeneous kinetics. For example, it enables to eliminate the constant term from Kissinger Equation [14, 15]:
This relation provides an easy experimental way for the determination of effective mass of activated complex from only two TA1 experiments. Furthermore, the modified Kissinger equation can be derived by similar way (please refer to Eq. (e) in Footnote 11).
in the following form12:
where the quantity of . The published value for the term
Since the average value of is then very close to the temperature-rate kinetic coefficient (
2. Effective mass and derived properties of activated complex
Since Eq. 23 contains three constants, it is also possible to write:
where (please refer to Eq. 39) and is the coefficient that describes how the partition function of the object changes with the change in temperature. Eq. 24 also means that the term is coded directly in the shape of peak of thermoanalytical curve. Since the unit of this term corresponds to the reciprocal temperature, for example, for applied example of the process of thermal decomposition of calcite, it is possible to calculate:
and to derive the temperature dependence of partition function on the temperature as follows:
|Parameter of activated complex|
The effective mass of activated state is then indirectly proportional to the reaction mechanism, the value of which is also affected by the values of
This relation then introduces the law, which enables to keep the change in mass when the reactants pass into an activated state. Small observed differences are most probably caused by the combination of different isotopic compositions (please refer to Section 6 in Chapter 2), different content of admixtures in the sample and an experimental error.
The value of
which means that the rate of the process is a linear function of the concentration of reactant [
However, the experimental value of
Eq. 1 can also be combined with Eq. 34 in Chapter 2 to provide the relation:
From that it can be derived that:
For the small difference between
This equation can next be transformed to the formula:
This equation has an important implication that:
These relations enable to describe
Since the dimension of term:
The combination with Eq. 36 in Chapter 2 then leads to the condition:
where , so the temperature of the reaction affects the thermodynamics (stability) of activated complex, that is, the reaction rate which is given by its decomposition, (Eq. 7 in Chapter 2), but the temperature does not alter the value of activation energy. It can be further derived that:
The density of activated state could be solved to:
where is the specific gas constant and (please refer to Eq. 80 and the discussion thereof). Since:
The effective diameter of activated complex can be calculated as follows:
If we consider that the velocity of activated state is equal to the term :
Assuming the parity that for the peak of energetic barrier it is 20, it is also possible to derive that:
Since the energetic density in activated volume is equal to the pressure:
more general form of Eq. 49 can be predicted:
which was first introduced in the previous work . This equation predicts the way, how the energetic density of activated state could be affected via the potential force field (
It is possible to derive the mass velocity (mass current density or mass flux) of activated complex () in solution from the dimension of the following term:
Since the activation energy of the process is considered to be constant, the velocity vector of activated state does not vary with time either:
That means that the divergence of velocity vector is zero and the net total flux through the surface of activated complex must be equal to zero.
Since the dimension of term:
has a physical meaning of length, it is possible to calculate the length of energetic barrier (
The mean lifetime of activated complex (), that is, the time, which is required to overcome the energetic barrier is then equal to:
where is the half-life of activated complex. That means that the concentration of activated complex can be written in terms of the exponential decay equation (please refer also to Eq. 7 in Chapter 2):
where term () and and are the concentrations of activated complexes at time
(please refer also to Eq. 20 in Chapter 1 and to discussion of Eq. 27 in Chapter 2). The effective rate of activated complex21 can be then formulated as follows:
for the peak of energetic barrier.
the bulk modulus of activated complex can also be calculated as follows:
The nature of velocities and can also be easily solved from the statistic probability function known as the Maxwell-Boltzmann distribution:
This relation can be further transformed to the following formula:
where the scale parameter of Maxwell–Boltzmann distribution (
and the mode of Maxwell-Boltzmann distribution (
that is, the diagonal of the square with the side length equal to (Figure 4).
Eq. 68 then also provides an important proof for the validity of Eq. 47. From this point of view, the velocities and are the scale factor and the mode, that is, the most probable speed, in Maxwell-Boltzmann distribution, respectively. Furthermore, it can also be written that:
Eq. 67 means that all activated complexes have constant values of the ratio:
The values of and are also related to the velocity of De Broglie’s “phase wave” (phase velocity ), which should occur in the phase with intrinsic particle periodic phenomenon . The relation between these velocities can then be expressed by the formula:
where is the invariant mass of activated complex (please refer to Eq. (a) in Footnote 21). Therefore, the ratio of (Eq. 70) has exactly the same value as the ratio of :
It is then obvious that is the diagonal of the square with the side length equal to . Since the ratio:
The values of , and for the example of thermal decomposition of calcite and aragonite applied in this book are:
529045.43: 1058090.86: 2116181.72;
539291.80: 1078583.61: 2157167.22;
respectively. That also enables to express the value of as the geometric ratio of and :
Therefore, the quantum numbers of activated complex, that is, the activation energy and the momentum of activated complex, enable to express the following three velocities, that is:
The group velocity of a supersonic activation wave, which is equal to the most probable velocity (mode) resulting from the Maxwell-Boltzmann distribution ().
Sonic dissipation (shock) wave, which corresponds to the velocity of sound. This velocity (scale parameter of Maxwell-Boltzmann distribution) is termed in this work as the effective speed of activated complex ().
Phase velocity, that is the subsonic De Broglie “phase wave” (), which is equal to the speed at which the phase of the wave propagates in given medium.
That behavior makes it possible to suggest the theory that the formation of activated complex (reversible process, please refer to Chapter 2) also generates the irreversible dissipative wave (shock wave). This discontinuity in the pressure, temperature and density preserves the energy, but increases the system entropy. That means that the formation of activated state (activated molecule and shock wave) causes increasing entropy of the system, even in the case, that activated complex decomposes back to the reactants.
Also, there is really interesting fact that which relates to the process of thermal decomposition of calcite (Table 1) and aragonite (Table 2), is of the magnitude, which corresponds to the rate of sound or ultrasound in liquids (please refer also to the discussion of Eq. 56), for example, 1497 m∙s−1 in distilled water at 25°C. Since only the longitudinal wave22 can be propagated in liquids and gases (fluids in general), the activated complex can oscillate (transfer energy) along the direction of the reaction coordinate only. It can be then derived:
|Parameter of Activated complex|
Then there is direct proportionality between and , where the value of proportionality constant is . From this, it can be calculated that the most probable rate of activated complex is equal to the speed of sound in air (346.3 m∙s−1, dry air at 25°C) and water (1497 m∙s−1) for equal to 230.1 and 994.7 K, respectively.
Since further increasing of temperature increases the rate of activated complex as well, so it becomes comparable with the speed of sound in the solids, it can then also be deduced that increasing the temperature of the process () increases the probability of oscillation in the direction perpendicular to the reaction coordinate, that is, the activated state can get additional degrees of freedom. The energy which is stored in those additional degrees of freedom then increases the energetic density of activated complex, which is on the contrary reduced by the value of
Therefore, there is a possibility to formally divide the reactions according to the behavior of activated state described above, as follows:
Reaction with most probable (group) velocity of activated state in fluid region, that is the value of is comparable to the speed of sound in gases and liquids. It can be potentially divided to the subregion of gases and liquids.
Reaction with most probable velocity of activated state in solid region, that is the value of is comparable to the speed of sound in solids.
3. Introducing the quantum numbers of activated complex
and to formulate the activation energy using the constitutive equation of the state of activated complex:
that is the effect of reaction temperature on enthalpy of activation (corresponds to the work of isobaric process) lays in the reduction of Δ
where . Eq. 84 could be further treated as follows:
Before we continue with quantum numbers of activated complex, it is interesting to mention that Eqs. (79) and (82) allow to solve the relation of these numbers (the activation energy and the momentum) to the Gravitational constant (
As results from footnote 29, there are other equations analogical to the Planck energy26 (
or to the Planck time26 (
and then the relation:
can be easily formulated. It can also be clearly seen that:
where is the Einstein’s constant (the coupling constant in the Einstein field equation), which is directly proportional to the gravitational constant, constant term
Using, for example Eq. 91, the relation can further be derived:
and Eq. 92:
It can also be verified that:
The motion of activated complex can then be fully characterized by two numbers, that is the activation energy and the momentum, which are directly proportional to physical action . Therefore, these two numbers can be considered as the quantum numbers of activated complex.
4. Formulation of wave function, reduced mass and relative velocity of activated complex
Using the momentum, the mass activation energy, etc., for one molecule of activated complex20,23, the Avogadro constant is left out from these relations. Furthermore, the quantum numbers of activated state and its mass can then be defined in terms of universal physical constants, for example, the Gravitational constant, reduced Planck constant and the speed of light in the vacuum.
The important consequences of relations 79 and 82 are the following:
The activated complex at the peak of energetic barrier, which represent the equilibrium, but not a stable state, has the kinetic energy (momentum) only.
The complete characterization of the motion of activated complex alongside the reaction coordinate requires the knowledge of two its quantum numbers:
Activation (kinetic) energy, it is obvious that activated processes have ;
Momentum of activated complex.
There is not any quantization either for the energy or for the momentum of activated state.
The shift from the equilibrium position (the peak of energetic barrier and the bottom of potential well at the same time) decreases the kinetic energy and increases the potential energy of activated complex.
The expression of
The substitution for
From the discussion above, it can then be concluded, that the definition of change in mass when the reactants pass into the activated state, enables to calculate many important parameters of activated complex. It is then feasible to use the Schrödinger equation29 for the description of its behavior, but the solution is in the same form as for classical physics:
is an angular frequency which is related to the “ordinary” universal frequency of activated complex (, please refer to Eq. 27 in Chapter 2 the reaction rate constant ratio ) via the relation:
It is obvious that the activated complex behaves as the pseudoparticle (instanton) that corresponds to one dimensional, that is the reaction coordinate (
where and the operator:
has the momentum:
The commutation of operators and then means that:
The solution of motion of activated complex (state) mentioned above enables to consider it and treat it as an instanton, that is a pseudoparticle. The pertinent wave function can be then used to describe its motion or to calculate the probability of tunneling through the energetic barrier, e.g. by means of the WKB (Wentzel-Kramers-Brillouin) approximation.
This behavior enables to use the following idea about the nature of activated complex (Figure 7(a)32). The activated complex is in equilibrium with products (
This behavior of activated state attached by spring to the product can be approximated by the Hooke’s law33 [34, 35], where the force
That also means that the universal frequency of activated state () is independent of
The acceleration of activated molecule during this oscillation is then given by the formula:
The potential energy, which is stored in the spring (activated complex) during the oscillation, is given by the relation:
This behavior is in agreement with previous conclusions derived from Eq. 82. Since the change in the potential energy is of constant rate, the relation can also be written:
where is the kinetic energy of activated complex during its oscillation around the energetic peak, so:
Therefore, the momentum of activated complex is equal to zero when (Figure 7(b)), that is for two spring limits, that is in amplitude positions:
where . Since there are no losses in the energy (the activated complex is an adiabatic system36), the activation energy has a constant value of:
Furthermore the value of is very small in comparison to the length of potential barrier (Eq. 54). In other words, the length of energetic barrier is much longer than the section which belongs to the oscillation of activated complex. The pertinent wave can then be described by the equation for the simple harmonic motion:
Since the magnitude of this function starts with the value , it is possible to shift the beginning of wave function (negative phase shift about radians38)39 to the peak of energetic barrier, where 0, then the wave function can be written:
The period of this equation is given by the relation:
This solution also provides an important insight to the nature of the mass and the velocity vector of activated state, which can be actually characterized as follows:
Reduced (effective32) mass of activated state ();
Relative velocity of the bodies before the collision ().
The activation can then be characterized as the change in kinetic energy during the perfectly elastic collision. Since there is no dissipation, if the kinetic energy is dissipated (the formation of activated state is a reversible process), the change of energy for this kind of collisions is:
where (elastic collision) is the coefficient of restitution (COR). Its value is defined as the ratio of the final to initial relative velocity or the square root ratio of the final () to initial energy () of two collided objects40.
The graphs of the wave function for the process of thermal decomposition of calcite (a) and aragonite (b) are shown in Figure 8. It can be seen that calcite with lower activation energy (Table 1 and Table 2) has longer period (2.15·10−15 s−1) than aragonite (2.0·10−15 s−1). Furthermore, the wavelength is then only a tiny quantity of length of the energetic barrier (Figure 7(a)).
The changes of the potential energy and of momentum of activated complex during the oscillation around the equilibrium position, that is the peak of energetic barrier (
For Eq. 126 (or Equation in Footnote 40) with the phase shift to the peak of energetic barrier, it is possible to write:
Since the activation energy per one activated complex is the sum of its kinetic () and potential energy ():
it can be written42:
Because the frequency of and is twice as high as the frequency of oscillation43, the middle (average) potential () and the kinetic energy () are of the same size during all periods. So, the relation of the middle activation energy () of the reaction can be written as follows:
Eq. 115 enables to use the principle of equivalency between undamped (the damping factor or ratio
In the series RCL circuit the inductance
If we conceive the idea that these oscillating systems are actually not undamped, but very slightly underdamped (factor
Since the rate of reaction is proportional to the decomposition of activated state, it should be proportional to the damping ratio as well. Therefore, the following approximation between the rate constant of reaction (
In other words, the exponential decay of underdamped oscillator is equal to
The value of so called
Using the example of the process of thermal decomposition of aragonite applied in this book, the lost energy can be calculated (Eq. 142) to 8.622 × 10−37 J per one cycle of activated molecule, which is only very tiny fraction of its value.
Since the loss of energy per each cycle increases with the energy of the activation impulse, the same idea can be used to explain the decomposition of activated system which collects the energy . In other words, if tunneling (please refer to Section 5 in Chapter 2) is not taken into the account, the energetic barrier behaves as an energy filter which allows only to the activated molecule with the energy of
5. Other parameters that affect the effective mass of activated state
As results from previous discussion, the effective mass as well as other derived parameters of activated complex are independent of the temperature at which the reaction proceeds at measurable rate44. Of course, that is possible only on the assumption that the change in temperature does not bring the change of reaction mechanism with different activation energy45 (Figure 10). In that case, the activation energy, the mass of activated state as well as other above mentioned parameters of activated complex must be changed too. The example of that behavior can be found in previous works, e.g. [36, 37].
Since the effective mass of activated state could also be affected by the isotopic composition of the specimen (please refer to Section 5 in Chapter 2), the activation energy for the most cases of the process is actually Relative Activation Energy. The value of which depends on the isotopic composition of the sample. Please do not be confused with the apparent activation energy. This term is usually applied to the kinetics of processes of unknown or uncertain mechanism.
The change in mass when reactants pass into the activated state is also affected by polymorphism46, that is by the crystal structure. For example, the aragonite has higher
The course of thermal decomposition of aragonite specimen (La Pesquera, Spain) is shown in Figure 11. In comparison with the process of thermal decomposition of calcite (Figure 1), the peak temperature of the process is higher. On the contrary, the full width at half maximum of peak (
The Arrhenius plot for the process of thermal decomposition of aragonite47 is shown in Figure 12. The process of thermal decomposition of aragonite requires slightly higher activation energy (Table 2) than that of Iceland spar (Table 1). The mechanism of the process of thermal decomposition, which includes zero or decreasing nucleation rate of new phase and the diffusion controlled growth of new phase, is also very similar for calcite and aragonite. Please compare the data in Table 1 (calcite) with those in Table 2 (aragonite) for the demonstration of effect of polymorphism on the kinetics of thermal decomposition and properties of activated complex. Further research has shown that there is actually not any significant effect of the polymorphism of calcium carbonate48 or its origin on the value of activation energy for the specimens of comparable purity .
The same cannot be recognized for polymorphs of SrCO3, which was investigated in previous work . The mass of activated state is higher for orthorhombic polymorph (0.165 kg·mol−1) than for hexagonal (0.148 kg·mol−1) strontium carbonate, but the effect of temperature change on the mechanism, that is the value of kinetic exponent
As was demonstrated in this chapter, the effective (reduced) mass of activated state is an important parameter of activated complex, which can be easily derived from the results of kinetic experiments. As was demonstrated on the example of calcite and aragonite, this change in mass when reactants pass into an activated state depends on real reaction mechanism, resp. on kinetic factor, which is often different from common transcript of the reaction. Therefore, the mass of activated state could be different from the value resulting from this equation, that is the sum of reactants multiplied by their stoichiometric coefficients.
The most important significance of this parameter lays also in the possibility of further definition of rate, density, energetic density, current mass density, momentum and many other properties of activated complex. Since the activation energy and the momentum of activated complex enable complete characterization of the motion of activated complex alongside the reaction coordinate, these parameters are its quantum numbers.
Furthermore, it is possible to introduce the idea to approximate the behavior of activated complex by the spring oscillation and to determine the nature of mass of activated complex as the reduced or effective mass of activated complex. This mass is also affected by the isotopic composition of the sample and by polymorphism.
Symbols and abbreviations
|j¯#||mass flux (mass current density) of activated complex|
|v¯p||the speed of the phase wave (phase velocity)|
|v¯x||the effective speed of activated complex, which corresponds to the scale parameter of Maxwell–Boltzmann distribution|
|ΔrG°||standard Gibbs free energy of the reaction|
|ΔrG||Gibbs free energy of the reaction|
|ΔG#||Gibbs energy of activation|
|ΔH#||enthalpy of activation|
|ΔS#||entropy of activation|
|ΔU#||internal energy of activation|
|ΔV#||volume of activation, the change of volume when reactants pass into the activated state|
|Bavx||the average value of approximation parameter Bx|
|CPT||Pythagorean triple constant for hypotenuses|
|Ei1||initial energy of two collided object|
|E0||energetic difference between energy of activated state and reactants|
|Ea||activation energy (Arrhenius activation energy)|
|Ep||Planck energy (1.956·109 J)|
|Eth||activation energy (theoretically calculated)|
|Ĥ||Hamiltonian (Evolutional) operator|
|K#||equilibrium constant (formation of activated complex)|
|K+||equilibrium constant of activation (K#=K+)|
|Kb||the bulk modulus|
|Lp||Planck length (1.616∙10−35 m)|
|M#||the mass of activated state|
|NA||Avogadro constant (6.022140857∙10−23 mol−1)|
|Q0#||partition function of activated complex|
|Qi||partition function of reactant molecules|
|T#||period of activated complex wave function|
|T̂||kinetic energy operator|
|T´||activated complex temperature term (Ea/R)|
|V#||molar volume of activated state (V#=RT´/p)|
|V̂||potential energy operator|
|Vp,1#||potential energy of activated complex|
|Zp||the size of particle|
|a#||acceleration of activated molecule|
|gj||degeneracy factor, that is the number of allowed equimolar quantum microstates|
|k#||activated state rate constant (coefficient)|
|kB||Boltzmann constant (1.38064852∙10−23 J∙K−1)|
|kexp||rate constant (experimentally determined)|
|kth||rate constant (theoretically calculated)|
|m1#,0||the invariant mass of activated complex|
|r#||specific gas constant|
|rAC#||the effective diameter of activated complex|
|t1/2||the half-life of activated complex|
|tp||Planck time (5.391∙10−44 s)|
|tδ=1/A||the mean lifetime of activated complex, that is the time required to overcome energetic barrier|
|v¯||the most probable speed of activated complex (the mode of Maxwell-Boltzmann distribution), which corresponds to the group velocity of the activation wave|
|w1/2||full width at half maximum of peak (FWHM)|
|αQ#||coefficient of response of partition function to a change in temperature|
|βv||dimensionless rate ratio of activated complex to the speed of light|
|εj||energy level of j-microstate|
|λ#||de Broglie wavelength of activated complex|
|ν#||universal frequency (ν#=kBT/h), that is the frequency of decomposition of activated complex|
|νAC#||universal frequency of activated complex (νAC#=kBT´/h)|
|ρ#||density of activated state|
|ω#||angular frequency of activated complex|
|∆=∇2||Laplace operator (Laplacian)|
|℘||coefficient of restitution (COR)|
|ΔV||molar change in volume during reaction|
|A||frequency or pre-exponential Factor, sometimes prefactor|
|A´||temperature independent constant|
|ART||Absolute reaction rates theory|
|Bx||approximation parameter in the Doyle equation for px|
|E||energy, usually reaction energy|
|G||gravitational constant (Universal constant or Newton’s constant, 6.67408 × 10−11 m3·kg−1·s−2)|
|K||equilibrium constant of reaction|
|P||Ad hoc quantity, that is “fudge factor”|
|Q||partition function of the molecule and Q factor|
|R||Gas constant, also Molar, Universal or Ideal gas constant (8.3144598 J∙K−1∙mol−1)|
|STR||Special theory of relativity|
|T||thermodynamic or absolute temperature|
|TST||Transition state theory|
|W||energetic difference between reactants and activated complex in basic state|
|Y||the constant term (Y≈7.0383×10−45m5/s3)|
|c||speed of light in vacuum (299,792,458 m·s−1)|
|const.´´´||temperature-rate kinetic coefficient|
|e||the base of natural logarithms (Euler’s number, Napier’s constant)|
|h||Planck constant (6.626070040 × 10−34 J·s)|
|k||reaction rate constant (coefficient)|
|px||temperature integral, where the quantity of x=Ea/RT|
|q||the heat of reaction|
|r||reaction rate (rate or speed of reaction) or shift from equilibrium position (according to the context)|
|∇||Nabla (del) operator|
|β||thermodynamic beta (occasionally perk)|
|γ||Lorentz factor (term)|
|δ||the length of energetic barrier|
|κ||proportionality constant (transmission coefficient)|
|ν||frequency of harmonic oscillator|
|ψ||wave function (psi)|
Truhlar DG, Garrett BC, Klippenstein SJ. Current status of transition-state theory. The Journal of Physical Chemistry. 1996; 100(31):12771-12800. DOI: 10.1021/jp953748q
Ptáček P, Opravil T, Šoukal F. Introduction of novel kinetic approach to calculation of activation energy and its application to the sinter-crystallization of strontian feldspar. Ceramics International. 2016; 42(15):16969-16980. DOI: 10.1016/j.ceramint.2016.07.203
Lever T, Haines P, Rouquerol J, Charsley EL, Van Eckeren P, Burlett DJ. ICTAC nomenclature of thermal analysis (IUPAC recommendations 2014). Pure and Applied Chemistry. 2014; 86(4):545-553. DOI: 10.1515/pac-2012-0609
Kostyrko K, Skoczylas M, Klee A. Certified reference materials for thermal analysis. Journal of Thermal Analysis. 1988; 33(1):351-357. DOI: 10.1007/BF01914623
Andon RJL, Connett JE. Calibrants for thermal analysis. Measurement of their enthalpies of fusion by adiabatic calorimetry. Thermochimica Acta. 1980; 42(2):241-247. DOI: 10.1016/0040-6031(80)87106-1
Robbins SA, Rupard RG, Weddle BJ, Maull TR, Gallagher PK. Some observations on the use of strontium carbonate as a temperature standard for DTA. Thermochimica Acta. 1995; 269-270:43-49. DOI: 10.1016/0040-6031(95)02602-9
Šesták J, Hubík P, Mareš JJ. Thermal physics and thermal analysis: From macro to micro, highlighting thermodynamics, kinetics and nanomaterials. Hot topics in thermal analysis and calorimetry. Springer. 2017; 11:567 p. DOI: 10.1007/978-3-319-45899-1
Augis JA, Bennett JE. Calculation of the Avrami parameters for heterogeneous solid state reaction using a modification of the Kissinger method. Journal of Thermal Analysis. 1978; 13(2):283-292. DOI: 10.1007/BF01912301
Avrami M. Kinetics of phase change. I. General theory. Journal of Chemical Physics. 1939; 7(12):1103-1112. DOI: 10.1063/1.1750380
Avrami M. Kinetics of phase change. II. Transformation-time relations for random distribution of nuclei. Journal of Chemical Physics. 1940; 8(2):212-224. DOI: 10.1063/1.1750631
Kolmogorov AN. A statistical theory for the recrystallization of metals. Izvestiya Akademii Nauk SSSR, Ser. Math. 1937; 1(3):355-359
Johnson WA, Mehl RF. Reaction kinetics in processes of nucleation and growth. Transactions of the American Institute of Mining & Metallurgical Engineers. 1939; 135:416-458. DOI: 10.1007/s11661-010-0421-1
Christian JW. The Theory of Transformations in Metals and Alloys. 2nd ed. New York: Pergamon Press; 1975. 586 p
Kissinger HE. Reaction kinetics in differential thermal analysis. Analytical Chemistry. 1957; 29(11):1702-1706. DOI: 10.1021/ac60131a045
Ptáček P, Frajkorová F, Šoukal F, Opravil T. Kinetics and mechanism of three stages of thermal transformation of kaolinite to Metakaolinite. Powder Technology. 2014; 264:439-445. DOI: 10.1016/j.powtec.2014.05.047
Doyle CD. Estimating isothermal life from thermogravimetric data. Journal of Applied Polymer Science. 1962; 6(24):639-642. DOI: 10.1002/app.1962.070062406
Lyon RE. An integral method of nonisothermal kinetic analysis. Thermochimica Acta. 1997; 297(1–2):117-124. DOI: 10.1016/S0040-6031(97)00158-5
Mallakpour S, Taghavi M. The accuracy of approximation equations in the study of thermal decomposition behaviour of some synthesized optically active polyamides. Iranian Polymer Journal. 2009; 18(11):857-872
Flynn JH, Wall LA. General treatment of the Thermogravimetry of polymers. Journal of Research of the National Bureau of Standards-A, Physics and Chemistry. 1966; 70A(6):487-523. DOI: 10.6028/jres.070A.043
Coats AW, Redfern JP. Kinetic parameters from thermogravimetric data. Nature. 1964; 201(4914):68-69. DOI: 10.1038/201068a0
Balarin M. Improved approximations of the exponential integral in tempering kinetics. Journal of Thermal Analysis and Calorimetry. 1977; 12(2):169-177. DOI: 10.1007/BF01909473
Senum GI, Yang RT. Rational approximations of the integral of the Arrhenius function. Journal of Thermal Analysis and Calorimetry. 1977; 11(3):445-447. DOI: 10.1007/BF01903696
Doyle CD. Kinetic analysis of thermogravimetric data. Journal of Applied Polymer Science. 1961; 5(15):285-292. DOI: 10.1002/app.1961.070051506
Málek J. The shape of Thermoanalytical curves as a function of the reaction kinetics. Thermochimica Acta. 1993; 222(1):105-113. DOI: 10.1016/0040-6031(93)80544-K
Gorbatchev VM. Remarks on the application of the combined Kolmogorov - Erofeev - Kazeev - Avrami - Mampel equation in the kinetics of non-isothermal transformations. Journal of Thermal Analysis and Calorimetry. 1978; 13(3). DOI: 509-514. DOI: 10.1007/BF01912390
Connors KA. Chemical kinetics: The study of reaction rates in solution. 2nd ed. John Wiley & Sons. 1990; 480 p
Beruto D, Searcy AW. Use of the Langmuir method for kinetic studies of decomposition reactions: Calcite (CaCO3). Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases. 1974; 70:2145-2153. DOI: 10.1039/f19747002145
Maciejewski M. 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. 2000; 355(1–2):145-154. DOI: 10.1016/S0040-6031(00)00444-5
L’vov BV, Polzik LK, Ugolkov VL. Decomposition kinetics of calcite: A new approach to the old problem. Thermochimica Acta. 2002; 390(1–2):5-19. DOI: 10.1016/S0040-6031(02)00080-1
Navarro CR, Agudo ER, Luque A. Thermal decomposition of calcite: Mechanisms of formation and textural evolution of CaO nanocrystals. American Mineralogist. 2009; 94(4):578-593. DOI: 10.2138/am.2009.3021
Halikia I, Zoumpoulakis L, Christodoulou E, Prattis D. Kinetic study of the thermal decomposition of calcium carbonate by isothermal methods of analysis. Technical Note. The European Journal of Mineral Processing and Environmental Protection. 2001; 1(2):89-102. DOI: 10.1179/174328507X198591
Matvienko AA, Chizhik SA, Sidel’nikov AA. A new kinetic model of calcite thermal decomposition. Doklady Physical Chemistry. 2013; 451(2):184-186 (Original Russian text published in Doklady Akademii Nauk 2013; 45(6):647–649)
De Broglie L. Waves and quanta. Nature. 1923; 112:540. DOI: 10.1038/112540a0
Hooke R. Lectures de Potentia Restitutiva, Or of spring explaining the power of springing bodies. London: John Martyn; 1678. 56 p
Roesler J, Harders H, Baeker M. Mechanical Behaviour of Engineering Materials: Metals, Ceramics, Polymers, and Composites. Springer Science & Business Media; 2007. 534 p
Ptáček P, Křečková M, Šoukal F, Opravil T, Havlica J, Brandštetr J. The kinetics and mechanism of kaolin powder sintering I. The Dilatometric CRH study of sinter-crystallization of Mullite and Cristobalite. Powder Technology. 2012; 232:24-30. DOI: 10.1016/j.powtec.2012.07.060
Ptáček P, Šoukal F, Opravil T, Havlica J, Brandštetr J. Crystallization of spinel phase from Metakaoline: The nonisothermal Thermodilatometric CRH study. Powder Technology. 2013; 243:40-45. DOI: 10.1016/j.powtec.2013.03.031
Ptáček P, Šoukal F, Opravil T. Kinetics and mechanism of thermal decomposition of calcite and aragonite. Journal of Mineral, Metal and Material Engineering. 2017; 3:71-79
Ptáček P, Bartoníčková E, Švec J, Opravil T, Šoukal F, Frajkorová F. The kinetics and mechanism of thermal decomposition of SrCO3 polymorphs. Ceramics International. 2015; 41(1):115-126. DOI: 10.1016/j.ceramint.2014.08.043
- Differential thermal analysis (DTA) is the TA method where the difference between heat flow rates into a sample and inert reference material (usually alumina) is measured . The calibration of DTA and DSC3 instruments uses melting and phase transition of temperature standard reference materials such as pure metals (In, Sn, Zn, Ag, Au…) or salts KNO3, KClO4, Ag2SO4, K2CrO4, quartz, K2SO4, BaCO3 and SrCO3. A comprehensive effort related to standardization and nomenclature of TA methods was launched in 1965 by the International Confederation for Thermal Analysis (ICTA) [4, 6]. In 1992, the name was changed to International Confederation for Thermal Analysis and Calorimetry (ICTAC) in order to reflect close relationship between TA and calorimetry.
- DTG, where adjective derivative “D” is pertaining 1st derivative (mathematical) of TG (thermogravimetric) or TGA (Thermogravimetric Analysis) curve (thermoanalytical or ICTAC discouraged collocation thermal curve) .
- Differential scanning calorimetry (DSC) is the TA method where the heat flow rate difference between sample and reference material is measured .
- The numerical value of const.″ to 20 decimal places is 1.13255326439210664113… The decimal approximation for const.‴ (Eq.2) is then 1.06421485818988011399…
- In order to avoid the confusion with the rate constant of reaction, this work uses the abbreviation “const.”, instead of letter “k”.
- Swiss mathematician, physicist, astronomer and engineer Leonhard Euler (1707–1783).
- The constant (irrational transcendental number) was actually discovered by the Swiss mathematician Jacob Bernoulli (1655–1654), who solved (1683) the value of the formula:limn→∞1+1nn=e.(a)
- With the exception of possible change in the reaction mechanism of the process, which may take place with increasing temperature (please refer to discussion of Figure 10), the value of kinetic coefficient should stay constant over the given interval of temperature (please refer to Footnote 9). Within the same interval it should be also insensitive to the shift of peak temperature with heating rate (Eqs.8 and 12). In the case that the average value of FWHM over the investigated interval of heating rate is applied for the calculation of kinetic exponent, the obtained dependence can be accurately approximated by the power law. Despite of the fact, that coefficient of variation is usually lower than 10%, the non-zero value of skewness and kurtosis of data means that results are not normally distributed. Therefore, the calculation of kinetic exponent with the average value of FWHM over the investigated interval of heating rate cannot be recommended.
- Together with the activation energy and the frequency factor, the kinetic coefficient is the part of so-called kinetic triplet. There is no clear physical interpretation for this constant (n, reaction “order”, nucleation rate, growth morphology, etc.), as well as for so-called effective overall reaction rate constant or temperature-dependent factor (k, which depends on the nucleation as well as on the growth rate) in the Avrami equation:α=1−exp−ktn;
- The ratio of derivation (Eq. 13) to its function (Eq. 12) is:cdΘd−1cΘd=dΘ(a)
- It is then possible to write the relation:
- In further work, we dealt with the most important formulation of modified equation, where the approach of compensation of the mutual influence between the transport of heat and the transport of mass was applied.
- The processes of thermal decomposition of calcite (Figure 1) and aragonite (Figure 11) show B(x) = 1.043 and 1.040 (please refer also to the footnote 51). When using the experimental data from previous work focused on the kinetics of thermal decomposition of strontium carbonate :
- Well known approximations of p(x) are provided by, e.g. Doyle :
- The concertation of one reactant in the second-order reaction of the type:
- The substitution of data from Figure 1(a) and Table 1 to Eq. 1 enables to calculate the activation energy:
- For example, using calcite with higher sizes of crystallites leads to the increase of Tm and w1/2 as well. Therefore, the nature of applied sample, its treatment (purification, intensive milling process, etc.) and applied conditions of analysis may also affect the mechanism of the investigated process, that is, the value of parameter M# as well.
- The solution of the limit (refer to the right side of Eq.33):
- An interesting consequence of these relations is to be solved in the next article.
- The details can be found in the previous work .
- It is also interesting to apply the Albert Einstein’s (1879–1955) special theory of relativity (STR), where it is possible to calculate the invariant mass of activated state as follows:
- On the contrary, the longitudinal as well as transversal waves (oscillate perpendicular to the direction of energy transfer) can propagate in the solids.
- The calculation of momentum per particle then requires to divide Eq. 79 by NA (Footnote 11 in Chapter 2).
- Eq. 79 then corresponds to the de Broglie matter waves (Louis Victor Pierre Raymond de Broglie, 1892–1987, Nobel Prize in Physics in 1929).
- Named after English physicist Paul Dirac (1902–1984). Dirac was awarded the Nobel Prize in Physics in 1933.
- Reduced Planck constant is the ratio:
- Please refer also to Eq. 11 in Chapter 2.
- For example, the molar volume of ideal gas at the temperature of 966.55 K (according to Figure 1 and Figure 2(b) this temperature corresponds to thermal decomposition of calcite at the heating rate of 1 K·min−1) is Vm=RT/p=7.931∙10−2m3∙mol−1, so Eq. 83 provides the value Ea = −101,325·(7.931·10−2 - 1.833) + 8.314·966.55 = 185.73·103 kJ·mol−1. Please compare to the value in Table 1.
- Erwin Rudolf Josef Alexander Schrödinger (1887–1961). The Schrödinger’s equation describes the behavior (evolution with time) of a physical system in which the quantum effects, such as the particle-wave duality take place:
- The vector differential operator “nabla” in three dimensional coordinate system with the basis vectors i¯,j¯ and k¯ and is written as:∇=i¯∂∂x+j¯∂∂y+k¯∂∂z.
- The differential operator given by the divergence of gradient of a function in Euclidian space:∆=∇∙∇=∇2=∂2∂x2+∂2∂y2+∂2∂z2;
- For the peak of energetic barrier, it can then be written:
- English natural philosopher, architect and polymath Robert Hooke (1635–1703).
- The other end of spring, that is, the reactants, is considered for the fixed position.
- The net energy of the activated complex is constant during the time.
- The “walls” of this system do not allow the transport of matter and heat.
- The formula is the solution of an ordinary differential equation (ODE):
- Like e and π, the π/2 is the transcendental number. This value corresponds to 90°, i.e. to 1 quadrant. The positive phase shift about 3π/2 = 4.712388… radians, i.e. 270° (3 quadrants), which leads to the equation:
- There is a very important fact that this phase shift corresponds to complete elliptic integral of the second kind:
- From this point of view, please see also the text (Van’t Hoff statement) related to the discussion of Eq. 9 in Chapter 1.
- In the case that there is not any phase shift.
- Please refer also to Eq. 115 and to Figure 9.
- On assumption that the activation energy is independent of the temperature. As was firstly assumed by van’t Hoff (please refer to discussion of Eq. 18 in Chapter 1), the activation energy could be a function of temperature. From this point of view, please refer also to the work of S. Vyazovkin and B.V. L’vov. Also according to our investigation, the activation energy is most probably the function of temperature, but within the temperature interval, where the reaction takes place in measurable rate, this change is usually smaller than the uncertainty of experimental results. We plan to publish our solutions on this topic in the future work.
- In this case, we can observe two or more linear parts in Arrhenius plot which are pertinent to the change of reaction mechanism. The mechanism crossing is usually gradual without an abrupt change of the slope. The effect of temperature on the reaction mechanism can be, for example, illustrated by the effect of temperature on the process of nucleation and diffusion. The increasing temperature makes the transport phenomenon’s, for example, the diffusion toward growing nucleus of new phase, easier, but it also reduces overcooling, for example, driving force of nucleation, as well.
- Polymorphism is the ability of solids to crystallize in various structures in different intervals of temperature and pressure. In the case of chemical elements, the same ability is termed as allotropism (allotropy). For example, the calcite and aragonite are two polymorphs of CaCO3 and graphite and diamond are two of allotropes of carbon.
- The chemical equation is the same as for calcite (Eq. 33 in Chapter 2). The Arrhenius and Eyring plot for the process of thermal decomposition of calcite can be found in Figure 2 in Chapter 2.
- This behavior can be explained by very small value of enthalpy of calcite (trigonal) ↔ aragonite (orthorhombic) phase transition (please refer also to the footnotes 18 and 19 in Chapter 2). The same cannot be said for vaterite, but this phase is the synthetic hexagonal calcium carbonate polymorph, that is, vaterite cannot be classified as a mineral.