Relevance of particle size distribution to kinetic analysis: the case of thermal dehydroxylation of kaolinite

processes Article Relevance of Particle Size Distribution to Kinetic Analysis: The Case of Thermal Dehydroxylation of Kaolinite Juan Arcenegui-Troya 1,*, Pedro E. Sánchez-Jiménez 1,2,*, Antonio Perejón1,2,* and Luis A. Pérez-Maqueda 1,*   Citation: Arcenegui-Troya, J.; Sánchez-Jiménez, P.E.; Perejón, A.; Pérez-Maqueda, L.A. Relevance of Particle Size Distribution to Kinetic Analysis: The Case of Thermal Dehydroxylation of Kaolinite. Processes 2021,9, 1852. https:// doi.org/10.3390/pr9101852 Academic Editor: Gurutze Arzamendi Received: 30 September 2021 Accepted: 14 October 2021 Published: 19 October 2021 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1Instituto de Ciencia de Materiales de Sevilla, CSIC-Universidad de Sevilla, C. Américo Vespucio No. 49, 41092 Sevilla, Spain 2Departamento de Química Inorgánica, Facultad de Química, Universidad de Sevilla, 41012 Sevilla, Spain *Correspondence: [email protected] (J.A.-T.); pedr[email protected] (P.E.S.-J.); [email protected] (A.P.); [email protected] (L.A.P.-M.) Abstract: Kinetic models used for the kinetic analysis of solid-state reactions assume ideal conditions that are very rarely fulfilled by real processes. One of the assumptions of these ideal models is that all sample particles have an identical size, while most real samples have an inherent particle size distribution (PSD). In this study, the influence of particle size distribution, including bimodal PSD, in kinetic analysis is investigated. Thus, it is observed that PSD can mislead the identification of the kinetic model followed by the reaction and even induce complex thermoanalytical curves that could be misinterpreted in terms of complex kinetics or intermediate species. For instance, in the case of a bimodal PSD, kinetics is affected up to the point that the process resembles a reaction driven by a multi-step mechanism. A procedure for considering the PSD in the kinetic analysis is presented and evaluated experimentally by studying the thermal dehydroxylation of kaolinite. This process, which does not fit any of the common ideal kinetic models proposed in the literature, was analyzed considering PSD influence. However, when PSD is taken into account, the process can be successfully described by a 3-D diffusion model (Jander’s equation). Therefore, it is concluded that the deviations from ideal models for this dehydroxylation process could be explained in terms of PSD. Keywords: kinetics; particle size distribution; kaolinite 1. Introduction Solid-state kinetic models give a mathematical description of a process. These algebraic expressions of simplified physical models are based on both mechanistic assumptions and geometrical considerations [ 1 – 5 ]. These ideal conditions are rarely met in real experiments, which has motivated criticism towards the use of kinetic models [ 6 ]. Thus, for example, real samples are rarely constituted by particles with identical shape and size. In fact, usually, they have a certain particle size distribution (PSD), which is a sample’s inherent feature that leads to different values of the fractional reaction α during the process depending on particle size. Moreover, it has been reported that PSD plays a significant role in the thermal behavior and kinetics of solid-state processes [ 7 – 9 ]. The modification of ideal kinetic models f(α) when considering PSD has been previously studied for some diffusion and interface reaction models [ 10 – 12 ]. Therefore, the modification of the kinetic models caused by the presence in the sample of particles with different sizes needs to be considered, as a kinetic analysis based on the shape of f(α) without considering the PSD might lead to incorrect conclusions. However, PSD is seldom considered in kinetics analyses. The first part of this paper emphasizes the importance of considering PSD when performing kinetic analysis. Consequently, we study how PSD affects the kinetic analysis in terms of the modification of the shape of interface reaction and diffusion models. Moreover, the specific case of bimodal particle size distributions is also considered. In the second part, as an example, we study the kinetics of dehydroxylation of kaolinite (Al 2 Si 2 O 5 (OH) 4 ). This is a very important process with applications in a number of Processes 2021,9, 1852. https://doi.org/10.3390/pr9101852 https://www.mdpi.com/journal/processes Processes 2021,9, 1852 2 of 16 industries (e.g., ceramics manufacturing, preparation of sialon, zeolite synthesis, geopolymer obtention, pozzolanic product manufacturing, paper production, pharmaceuticals, and environmental decontamination) [ 13 – 20 ], and it has largely been investigated without attention paid to the effect of PSD [ 21 – 28 ]. Moreover, the conclusions of these studies do not lack controversy. The values of the reported activation energies lie within a wide range, from 140 to 250 kJ/mol [ 29 – 32 ]. Besides, there is no consensus on the kinetic model that best describes this reaction. Some authors concluded that the reaction can be interpreted in terms of a first- or n-order kinetic model [ 26 , 33 – 38 ]. Other researchers have described it using Johnson–Mehl–Avrami or even diffusion models [ 21 , 23 , 39 – 41 ]. However, in [ 39 ], the authors observed a deviation from the diffusion models for values of α> 0.6. The complexity of the dehydroxylation has led some authors to suggest a change in mechanism during the process [ 42 , 43 ]. Other authors have pointed out that the process involves many individual steps with variations in the activation energy throughout it [ 25 ]. Here, we revisit the dehydroxylation of kaolinite considering PSD in this kinetic analysis, concluding that the apparent complexity of the process could be explained in terms of its PSD. 2. Theoretical Approach The general equation that describes a solid-state reaction is: dα dt =Aexp−E RT f(α)(1) where α stands for the extent of the reaction, t is the time, A is a constant termed preexponential factor, Eis the apparent activation energy, R= 8.31 J · mol −1· K −1 is the gas constant, and T is the temperature of the process. The function f(α) represents the kinetic model relating the rate of the reaction to α . Under isothermal conditions, this equation can be integrated to obtain [44]: α Z 0 dα0 f(α0)=Aexp−E RT t Z 0 dt0(2) Using the notation g(α)=Rα 0dα0 f(α0) , k=Aexp−E RT  and integrating the right side of Equation (2), we can write: g(α)=kt (3) The dependence of kinetics on the particle size r lies on k (Equation (3)). In general, we can write: k=k0S(r)(4) where k0 is a constant and S(r) is a function of the particle size. Table 1shows the expressions for S(r) for the different ideal models studied in this paper. Substituting Equation (4) in (3) and ordering terms, we get: g(α)−k0S(r)t=0 (5) Table 1. Kinetic models of diffusion and interface reaction studied in this work. Symbol Particle Shape Meaning of r S(r)g(α) 2D diffusion D2 Cylinder Base diameter 1/r2α+ (1−α)·ln(1−α) 3-D diffusion (Jander) D3 Sphere Diameter 1/r21−(1−α)1/32 3D diffusion (Ginstling–Brounshtein) D4 Sphere Diameter 1/r2 1 −2α 3−(1−α)2/3 2D interface reaction R2 Cylinder Base diameter 1/r1−(1−α)1/2 3D interface reaction R3 Sphere Diameter 1/r1−(1−α)1/3 Processes 2021,9, 1852 3 of 16 Expressions for g(α)are given in the right column in Table 1[1]. In general, Equation (5) can be numerically solved for any kinetic model to obtain the extent of the reaction α as a function of time for a given value of r . In the case of an R3 model, Equation (5) takes the form (Table 1): 1−(1−αr)1/3 −k0 rt=0 (6) whose solution is: αr=1−1−k0 rt3 (7) This latter function is plotted in Figure 1a, with k0= 2.8 × 10 −12 µ m · s −1 , for different particle sizes. As expected, the time required to complete the reaction increases with the size of the particle. In fact, larger particles start to react at temperatures when the smallest ones are almost completely converted. This result has been substantiated by experimental investigations on the dehydroxylation of fractions of pyrophyllite with different particle sizes, which showed that the smaller the particles, the lower its average dehydroxylation temperature [45]. Figure 1. ( a ) Fractional reaction as a function of normalized time for different particle sizes. The overall values for the sample are plotted as a pink solid line. (b) Lognormal PSD with σ=1 and µ=ln 10−5. The overall values of the extent of the reaction, shown as a pink solid line in Figure 1a, were calculated according to: α=∑ r αrV(r)∆r(8) where V(r)∆r represents the volume fraction occupied by the particles whose size is r , with ∆r being the interval of sizes in which the volume fraction is considered to be constant. In this study, we use a lognormal-type PSD: V(r)=1 rσ√2πexp −(ln r−µ)2 2σ2!(9) Specifically, the results of the simulation plotted in Figure 1a were obtained using the PSD shown in Figure 1b, with σ= 1 and µ=ln 10 −5 , and the particle size ranging from 0 to 100 µ m. The whole range was discretized into intervals of ∆r = 1 µ m. As can be observed, the shape of the curve that represents the temporal evolution of the overall Processes 2021,9, 1852 4 of 16 fractional reaction, considering the PSD, differs from the shape of the curve corresponding to a single particle with a specific size. 3. Experimental Section A low-defect kaolinite sample from Washington County, Georgia (KGa-1 from the Source Clay Mineral Repository, University of Missouri, Columbia, MO, USA), was used for the present study. Dehydroxylation experiments were conducted in a thermogravimetric analyzer (TGA). The experiments were conducted in small samples (approx. 10 mg) at high vacuum (residual pressure: 3 · 10 −5 millibar) to minimize mass transfer phenomena. The series of experiments were conducted under conventional linear heating conditions at 1, 5, and 10 K · min −1 and non-conventional sample-controlled thermal analysis (SCTA) at a constant reaction rate of 4.6 · 10 −3 min −1 . In the latter case, feedback from the thermogravimetric signal is used as an input in the algorithm commanding the furnace control in such a way that the total reaction rate remains constant over the entire process [46–49]. Particle size distribution of the kaolinite sample used here was measured using a low-angle laser light scattering instrument (Mastersizer Malvern Instruments). 4. Results and Discussion 4.1. Effect of PSD in Simulated Linear Heating Experiments Data plotted in Figure 1a can be used to derive the kinetic model that describes a 3D interface reaction occurring in a sample with the PSD shown in Figure 1b. Indeed, according to Equation (1), this can be achieved by differentiating the curve plotted as the pink solid line as follows: f(α) f(0.5)= dα dt dα dt 0.5 (10) For the sake of clarity and ease of comparison with other models in the literature, the kinetic model was normalized to its value for α= 0.5. The normalized kinetic model is represented as a function of the extent of the reaction in Figure 2. The ideal model R3 is also plotted in Figure 2. Consistently with the results shown in Figure 1a, the kinetic model is significantly modified when we take PSD into account. Figure 2. Normalized kinetic models. The dashed green line represents the ideal model R3, whilst the continuous red line corresponds to the kinetic model obtained when PSD is taken into account. Processes 2021,9, 1852 5 of 16 Using the kinetic model plotted in Figure 2, we simulated linear heating experiments intended to study the kinetics of a thermally induced reaction. The results of this simulation are shown in Figure 3a. To simulate the experiments, we solved the following system of equations using the Runge–Kutta method with the initial conditions T(t=0)= 275 K and α(t=0)=10−4: dα dt =Aexp−E RT f(α)dT dt =β(11) where β represents the heating rates. Four different heating rates were considered: 1, 2, 5, and 10 K · min −1 . The pre-exponential factor used was A= 10 10 s −1 , and the activation energy was set to E= 100 kJ·mol−1. Figure 3. ( a ) Curves simulated under linear heating conditions using the kinetic model R3 with the PSD shown in Figure 1b. ( b ) Values of activation energy as a function of the fractional reaction obtained by the Friedman isoconversional method. (c) Combined kinetic analysis. Processes 2021,9, 1852 6 of 16 Results of the Friedman isoconversional method applied to data in Figure 3a are depicted in Figure 3b. As expected, the values of activation energy remain constant for all the values of conversion. Thus, if this were an analysis of experimental data collected in the laboratory, the conclusions would be that this process can be described with a sole value of activation energy, and there is only one reaction kinetic mechanism [50,51] To discriminate the kinetic model followed by the process, the combined kinetic analysis, which simultaneously analyzes all experimental data obtained under any heating conditions, was used. This analysis is based on the general kinetic Equation (11) that after rearranging terms can be written in logarithmic form as follows: ln dα dt f(α)!=ln A−E RT (12) Thus, only the correct kinetic model, f( α ), would fit all of the experimental data yielding a straight line from whose slope and intercept the activation energy and preexponential factor can be determined, respectively [ 52 ]. Moreover, a modified Sestak–Berggren equation is often used in the simultaneous combined kinetic analysis of experiments conducted under different experimental conditions: f(α)=C(1−α)nαm(13) where C , n , and m are the fitting parameters. The original Sestak–Berggren equation includes a term of the form (−ln(1−α))p [ 53 ]; however, it has been observed that this modified, simplified version given by Equation (13) can fit every f( α ) of the ideal kinetic models most extensively used in the literature [ 54 ]. Using this expression, Equation (12) can be rewritten as follows: lndα/dt C(1−α)nαm=ln A−E RT (14) Therefore, the model that best describes the reaction linearizes lndα/dt/C(1−α)nαm as a function of 1 /T . Thus, the parameters of Equation (13) are obtained by an optimization procedure that yields the parameters that maximize the linear correlation coefficient [ 54 ]. Then, the activation energy and the pre-exponential factor are both calculated from the slope and intercept of the plot of lndα/dt/C(1−α)nαm as a function of 1 /T , respectively. Figure 3c shows the results of combined kinetic analysis applied to data presented in Figure 3a. The values of the fitting parameters and correlation coefficient are given in the same graph. The values of E and A obtained are close to those used in the simulation, which validates the use of the combined kinetic analysis in this case. Nonetheless, although Equation (13) can be used to reconstruct the ideal kinetic models accurately with the appropriate choice of values for C , n , and m , in general, Equation (13) has no physical meaning for an arbitrary selection of the fitting parameters. A practice commonly employed to provide the obtained equation with a physical significance consists of a graphical comparison with the ideal kinetic models from the literature. This comparison is made in Figure 4. The ideal kinetic models depicted are A0.5, F1, A2, A3 (nucleation), R2, R3 (interface reaction), D2, D3, and D4 (diffusion). Among these models, the one that most resembles the model obtained from the combined kinetic analysis is A0.5, which corresponds to a particular type of Avrami–Erofeev model used to describe a process of nucleation and growth. Therefore, these results lead to the conclusion that the reaction studied follows a nucleation and growth model instead of an interface reaction mechanism modified by the particle size distribution, which is the actual case. Processes 2021,9, 1852 7 of 16 Figure 4. Comparison between the modified Sestak–Berggren equation obtained from combined analysis (open black squares) and some ideal kinetic models from the literature (lines) [1]. 4.2. Effect of PSD on Diffusion and Interface Reaction Models The study on how PSD modifies the shape of the R3 kinetic model was extended to other ideal models, such as interface growth and diffusion-controlled models (Figure 5 shows the results). In this study, three possible PSDs with σ = 0.1, 0.5, and 1, while maintaining µ=ln 10 −5 constant, were considered in the calculus, (see the resulting PSD curves in Figure 5a). As observed in Figure 5b–d, in the cases of diffusion models, it is hard to discern if the reaction is described by one model or another when PSD comes into play. For instance, a sample with a PSD characterized by µ=ln 10 −5 and σ= 0.5 that reacts according to D2 could be wrongly assigned to a D3 model (Figure 5c), while D3 can be confused with a nucleation A0.5 model (Figure 5b) when considering PSD. As observed in Figure 5e,f, the impact seems to be more noticeable in the case of interface reaction models, which can be easily mistaken with a diffusion-type model or a nucleation A0.5 model. 4.3. Effect of a Bimodal PSD in Kinetic Curves When measuring particle sizes, a PSD that exhibits two separated peaks or humps is commonly found, as observed in the graphs shown in Figure 6. The PSDs shown in this figure was calculated as the sum of two lognormal functions, with µ=ln10−5 and σ= 0.5, whose peaks are separated at a variable distance of a. For the sake of clarity, the peaks are plotted separately in each graph as lines, while the sum of both is represented by open blue circles. The calculation was performed assuming that every individual particle in the sample reacts according to the ideal model R3. The extent of the reaction and its derivative as a function of temperature are plotted in Figure 7. These data were simulated assuming linear heating conditions with a heating rate of β= 1 K · min −1 . As in the case of Figure 6, the contribution of each peak is depicted separately as a line, while the total is plotted as open blue circles. As can be observed, the separation between peaks has a clear impact on the shape of the curves α -T (Figure 7a–c) and derivatives (Figure 7d–f). For the largest value of aused, the derivative clearly exhibits two humps. This can be explained by the fact that the smaller particles in the left peak of the PSD require much less time to attain a complete reaction than those particles in the right peak of the PSD. Such behavior was observed in Figure 1. Therefore, when the smaller particles complete the reaction, at low values of temperature, a considerable portion of the volume occupied by the larger ones remains unreacted. This result could lead to the misleading conclusion that the reaction is driven by two processes or mechanisms Processes 2021,9, 1852 8 of 16 separated by temperature. Therefore, taking PSD into account when performing kinetic analysis is crucial to avoid such mistakes. Figure 5. Modification of ideal kinetic models when PSD is taken into consideration. ( a ) Different PSD used. ( b ) D3, ( c ) D2, (d) D4, (e) R2 and (f) R3. Figure 8shows the kinetic models modified by the consideration of the bimodal PSDs. Again, for the sake of clarity, the contributions of each peak to the kinetic model are presented separately. When ais very large, the kinetic model obtained deviates largely from the ideal one for α< 0.1, whilst both coincide for higher values of the fractional reaction. As can be observed in Figure 6, the relative relevance of σ decreases as aincreases, and for a = 150 µ m, the second peak can be practically considered as single particle size distribution. This reflects in the graphs of Figure 8, where the contribution of the second peak to f(α) tends to be the ideal kinetic model used in the simulation (R3) as aincreases. For a = 150 µm, the contribution of the second peak and the ideal kinetic model overlap. Processes 2021,9, 1852 9 of 16 Figure 6. Bimodal particle size distributions using the simulations. The total PSD was calculated as the sum of two lognormal functions, with µ=ln10−5 and σ= 0.5 being the amplitude of the second peak ten times the amplitude of the first one. The separation between peaks was set at (a) 2 µm, (b) 20 µm, and (c) 150 µm. Figure 7. ( a – c ) Curves of α as a function of temperature, simulated using the R3 kinetic models modified by bimodal PSD, with β =1K · min −1 and a= 2, 20 and 150 µ m, respectively. ( d – f ) Derivative of α with time as a function of temperature for a= 2, 20 and 150 µm, respectively. Figure 8. Kinetic models modified by the consideration of the bimodal PSDs. The curves that represent the contribution of the peaks are plotted as lines, while the total kinetic model is depicted as a set of blue open circles. The separation between the two peaks of the bimodal PSDs was set at (a) 2 µm, (b) 20 µm and (c) 150 µm. Processes 2021,9, 1852 16 of 16 55. Franco, F.; Cecila, J.; Pérez-Maqueda, L.; Pérez-Rodríguez, J.; Gomes, C. Particle-size reduction of dickite by ultrasound treatments: Effect on the structure, shape and particle-size distribution. Appl. Clay Sci. 2007,35, 119–127. [CrossRef] 56. Wiewióra, A.; Pérez-Rodrıguez, J.L.; Perez-Maqueda, L.A.; Drapała, J. Particle size distribution in sonicated high- and low-charge vermiculites. Appl. Clay Sci. 2003,24, 51–58. [CrossRef] 57. Perez-Maqueda, L.A.; Blanes, J.M.M.; Pascual, J.; Pérez-Rodríguez, J.L. The influence of sonication on the thermal behavior of muscovite and biotite. J. Eur. Ceram. Soc. 2004,24, 2793–2801. [CrossRef] 58. Pérez-Maqueda, L.A.; Montes, O.M.; González-Macias, E.M.; Franco, F.; Poyato, J.; Pérez-Rodrıguez, J.L. Thermal transformations of sonicated pyrophyllite. Appl. Clay Sci. 2004,24, 201–207. [CrossRef] 59. Vyazovkin, S.; Burnham, A.; Criado, J.M.; Perez-Maqueda, L.A.; Popescu, C.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim. Acta 2011,520, 1–19. [CrossRef] 60. Otero-Arean, C.; Letellier, M.; Gerstein, B.C.; Fripiat, J.J. Protonic structure on kaolinite during dehydroxylation studied by proton nuclear magnetic resonance. Summary 1982,105, 499. [CrossRef]