- Open Access
Detailed kinetic and chemometric study of the cellulose thermal breakdown in artificially aged and non aged commercial paper. Different methods for computing activation energy as an assessment model in archaeometric applications
Chemistry Central Journalvolume 6, Article number: S7 (2012)
The thermal oxidative degradation of aged and non aged cellulose samples of commercial paper was studied using thermogravimetry and derivative thermogravimetry under a forced air flow up to 800°C.
TG and DTG data were processed using two non-isothermal-based model-fitting methods and one based on linear least squares to calculate Ea trend values, measured as a function of artificially induced sample age. The Ea trends thus obtained were compared in order to assess their potential for yielding archaeometric curves. As the trends of first two methods show an inversion of the direction between non aged cellulose samples and artificially aged samples, while the third method does not, an in-depth study was carried out using a multilinearity assumption.
The results are discussed and the outcomes indicate that the above cited inversion is real and not linked to the method. Additionally, it was evidenced that the number of points used for the estimation of linear least squares model parameters is of capital importance.
It was postulated that in the paper or wood samples respectively aged naturally or artificially, the state of the cellulose polymer chains is damaged by well-known oxidative, photo-oxidative and hydrolytic processes, which occur during ageing . It was therefore attempted to correlate ageing, or the age, of the cellulosic samples with age-linked parameters that could be measured experimentally, if possible using relatively simple procedures and thus avoiding expensive apparatus [5–7].
In essence, it was endeavoured to detect any relationship between the breakdown of the cellulose and the physico–chemical parameters by which this could be evidenced. Consequently, the thermal oxidative degradation of non aged, or artificially aged, cellulose samples of commercial paper as such, or UV irradiated in a weatherometer, were studied using thermogravimetry (TG) and first-order derivative thermogravimetry (DTG) . In this framework the activation energy (Ea) values, which are related to cellulose fiber packing, seem to be more promising for the purpose of constructing archaeometric curves . However, there appears to be no recommended literature method for calculating the Ea of the cellulose thermal degradation process. Thus, in a previous study , in order to obtain reliable activation energy values and to identify the method yielding the most reliable data, a comparison was made between the results obtained by applying three different methods for calculating the Ea of the oxidative thermal breakdown of the cellulose contained in several different samples of aged and non aged paper. In particular, two non isothermal model-fitting methods – the differential Arrhenius [9, 11] and the integral Satava [9, 12] methods - were applied and the results obtained compared with those determined using the Wyden-Widmann method . All three methods are based on the linearization of suitable equations using the least squares method, from the slope of which the Ea values are determined. Nevertheless, a good correlation was found between Ea trends obtained using Arrhenius’ and Satava’s calculation methods and only partially with those obtained using Wyden-Widmann’s method.
The Ea trends obtained so far from TG/DTG data as a function of sample age for artificially aged extra-strong paper samples using the first two methods have not always been found to be completely linear. However, an inversion of the direction of these trends was found above all on going from fresh samples to aged samples, while the trends seem to be sufficiently linear whenever only samples of different artificial ageing times are considered.
For this purpose, from a previous study , it was found possible only to conclude that the cause of the inversion in the direction of the trend of Ea values observed in non aged cellulose samples, compared with aged samples, may be different and is often multiple in nature. For the time being it is not possible to conclude what the most plausible, or at least the most prevalent, reason might be that a degree of synergism possibly occurs also between several causes of this phenomenon. Furthermore, it must be determined to what extent the kinetic method used to calculate the Ea values actually affects the trends obtained. This point obviously calls for a more thorough investigation in order to to achieve full elucidation, which is precisely what has been attempted in the present note.
Results and discussion
As stated above, a reasonable correlation was found in a previous paper  between the trends of Ea values obtained using Arrhenius’ and Satava’s methods and partly also from the trend found using the Wyden-Widmann least squares method (MLR) applied to the data obtained from the TG and DTG curves of artificially aged extra-strong paper shown in Fig 1, while some difficulty was found when the Ea value of cellulose thermal degradation of non aged paper was compared with that referring to the same process in non aged paper found using Arrhenius’ or Satava’s methods (the TG-DTG curve of which is also reported in Fig. 1).
In the latter case, the representative Ea point was quite anomalous and differed as to the trends obtained (see Fig. 2), so that in the latter case non negligible differences over the entire trend of Ea values were found between the first two model-fitting methods and that obtained using Wyden-Widmann’s method.
In the interpretation of these outcomes, it was not clear to what extent they depended on the real differences in the physico-chemical characteristics of the cellulose in the non aged sample versus the aged samples and, conversely, to what extent the kinetic method used to calculate the Ea values actually affected the Ea trends obtained using the equations (6, 7 and 9) for the Arrhenius, Satava and Wyden-Widmann methods, respectively. In the present research Wyden-Widmann’s method, which is very often applied automatically using computer implemented software coupled with the thermogravimetric apparatus, was thoroughly tested. Initially, computation of Wyden-Widmann model parameters was also tried using a non-least squares iterative method based on continuous-valued genetic algorithms  but results statistically comparable to those ones deriving from the application of the MLR method were obtained, both when applied to the linearized version and when nonlinear solution was evaluated directly. However, as also discussed in the Experimental section, we observed how the estimated value of Ea can change significantly depending on the temperature range (and, as a consequence, on the number of data points) selected for use in the fitting procedure. Accordingly, the computation of model parameters in the Wyden and Widmann method was repeated by including a different number of data points in the fitting procedure in order to account for this observation.
When examining the results obtained we realized that a reasonable correlation existed (see Fig. 3) between the trends obtained using the Arrhenius and Satava methods with that found by the Wyden-Widmann method when a sufficiently large number of experimental points was used (at least > 60).
Indeed, a thorough investigation, as shown in Fig. 4, revealed that the data corresponding to the TG/DTG curves of extra-strong paper, after a logarithmic transformation (see equation (9)), fit the multilinearity assumption only over a limited temperature range. Hence, the quality of the fitting process is highly dependent on the number of points used to compute the regression model, as shown in Fig. 5.
In particular, when Wyden - Widmann’s method is used, a trend similar to those found using Satava’s and Arrhenius’ methods is obtained only when a number of points lying between 62 and 100, before and after the DTG peak are processed in the fitting. In practice it was thus possible to conclude that: (i) the most correct Ea trend, as a function of the irradiation time in the weatherometer, seems to be the one obtained using Arrhenius’ or Satava’s methods, based on the identification of the model function that best fits the experimental trend of the DTG peak value - in practice the type of reaction mechanism [11, 12]. However, a sufficiently correct Ea trend (i.e. in sufficient agreement with the Satava and Wyden – Widmann trend) may be obtained also by applying the Wyden-Widmann method but using only an appropriate number of points right and left of the DTG peak, (ii) It therefore seems that the inversion of trend found on going from the non aged extra-strong paper sample to the artificially aged samples is real and linked to physical – chemical differences between the non aged and aged paper matrixes (at least for the Ea values referring to the thermal degradation of the commercial extra-strong paper studied in the present research), (iii) When the results of Ea values obtained by applying all three methods are compared in this way (see Fig. 3), it seems that the most probable Ea values lie in the range of 150-180 kJ mol–1 and that the fitting procedure indicates that the most suitable mechanism of all the functions reported in Table 1 is the F1 function (first-order mechanism).
Furthermore, if the same procedures to calculate Ea, as those used in the three methods examined are repeated for values referring to all the extra-strong paper samples tested at different irradiation times, excluding only the non irradiated sample and using all three methods considered alternately, it may be concluded (see Fig. 6) that, in this case, the function F1 produces the expected quasi monotonically decreasing trend and at the same time Ea values having an order of magnitude equal to those previously found and that also Wyden-Widmann’s method, again using a number of experimental points approximately lying between 40 and 60, in this case produces a trend in Ea values very similar to those of the two Model Fitting methods.
Once again this seems to indicate that the anomalous value obtained for the non irradiated sample does not depend on the method used to calculate Ea, but on physico-chemical differences inherent in the sample itself vis-à-vis the irradiated samples. Further confirmation of this conclusion was obtained by observing the Ea trends obtained using each of the three calculation methods described herein but referring to the TG-DTG curves shown in Fig. 7 recorded for another type of commercial Italian paper - “Fabriano”, both artificially aged and non aged.
As shown in Fig. 8 good agreement in this case is found among the Ea value trends obtained using the customary three calculation methods only when, in the case of the Wyden-Widmann method, the number of values processed around the DTG peak is ≥ about 50. Here the anomaly observed for the Ea value of the non irradiated sample of extra-strong paper does not occur. This amounts to a definitive confirmation that this anomaly is to be attributed to the non irradiated sample and not to the computation method used.
Lastly, the decrease in the total water content in aged, in comparison with non aged, paper samples, detected by us using TG measurements (see Fig. 9) appears to be related to some extent to the decrease in microcrystalline domains in the paper samples and to the increase in the amorphous domains which seems to occur with increasing paper sample irradiation time. Our thermogravimetric method findings actually seem to be in good agreement with the results of NMR studies on several paper samples [15–17]. These studies demonstrated that water is wholly bound to the cellulose lattice, but only semi-crystalline fibrous cellulose is rich in water, while the amorphous variety has a low water content. On the other hand, comparative NMR measurement also of the total water amount in paper samples revealed a different degree of degradation linked to total water variation. It was thus possible to conclude that the loss of water is associated with increasing cellulose degradation of polymorphs, which, in turn, may be associated with the increase in the amorphous fraction. Moreover, the results of NMR measurements agree with those obtained using the SANS method , that is, in both cases clusters (or microdomains) in which the highest concentration of water contained in the aged paper samples were detected.
Samples and apparatus
Paper ageing was achieved by photoirradiation at λ = 310 nm in a weatherometer (Model QUV/spray Q-Panel LAB-Products) for up to 750 hours. The apparatus was equipped with 8 mercury lamps (40 W each) able to simulate sunlight perfectly up to 370 nm.
Samples of extra-strong paper for photocopiers and laser printers (white paper for office photocopiers, Copy blu, Code 13092 - 6IN - format 21x29.7 cm, A4, basic weight 80 g m–2) were studied either as such or artificially aged in a weatherometer for different ageing periods up to a maximum of 750 h.
The TG/DTG measurements were carried out on a Mettler 50 thermobalance coupled with a Mettler TG 10-TA processor system and to a Swiss dot matrix printer.
For the purpose of performing thermogravimetric experiments, from each extra-strong paper sample 5 small circular disks (5.0±0.1 mm in diameter) were cut out, each weighing 1.7±0.1 mg. The diskettes were then placed in a special cylindrical alumina crucible, 0.7 cm in diameter, and carefully compressed using a small piston before being subjected to TG analysis under a forced air flow (0.1 dm3 min–1) at a heating rate of 10 °C min–1.
Another typical and well known commercially available Italian paper (“Fabriano” paper), both non aged or artificially aged in a weatherometer, was also analysed. In this case all the paper samples were minced finely before analysis. The size of individual paper fragments was of the order of 1.0 ± 0.2 mm2 on average.
In the course of artificial ageing in the weatherometer, Fabriano paper samples were taken at several irradiation times up to a maximum of 200 hours.
Determination of activation energy using the Arrhenius and Satava model-fitting procedure
The thermal breakdown of the cellulose contained in the paper and wood samples was studied by treating the TG/DTG data referring to this process. The kinetic study of the reactions undergone by a compound in the solid state is generally based on the following relation:
where α is the degree of conversion, calculated according to the following equation:
where mi, and mf are the initial and final sample mass, while mT is the sample mass at a given temperature T; k is the process kinetic constant and f(α) is a function defined by different mathematical expressions depending on the type of reaction mechanism [19, 20]. The explicit dependence of the reaction constant on temperature is obtained by expressing the kinetic constant k using the Arrhenius equation:
which, when included in Eq. (1), gives:
in which A (pre-exponential factor) and Ea (activation energy), are the Arrhenius parameters and R is the universal gas constant. When the sample is subjected to controlled heating at a constant rate β, where β = dT/dt, the explicit dependence of reaction rate dα/dt on time in equation (4) is eliminated, as:
Multiplying both sides of equation (5) by β · f(α)–1 and calculating the natural logarithms, the following equation is obtained:
In the two methods used, the left-hand side of equations (6) and (7) is assumed to be linear with respect to the reciprocal of the temperature (1/T) only for the functions f(α) and g(α) which represent the best mechanism for the reaction considered. From the slope and intercept of the regression straight lines described in equations (6) and (7) the values of the activation energy and the pre-exponential factor may be calculated. The differential and integral f(α) and g(α) functions that best fit the reaction mechanism were chosen from among those reported in literature  on the basis of experimental parameters referring to the shapes of DTG peaks related to the thermal oxidative process of cellulose contained in the wood or paper samples, as described in detail in previous papers [9, 22]. The model functions reported in literature  and selected in previous papers for similar cellulosic samples  are denoted by the symbols A (Avrami-Erofeev models) and F (nth-order decay models).
Determination of activation energy using the Wyden-Widmann method
Using this method the activation energy values were always obtained starting from the TG and DTG data, which were however processed by means of the multiple linear regression method proposed by Wyden-Widmann . Also using this method the activation energy Ea of the cellulose decomposition step was calculated on the basis of the values derived from a single thermogravimetric experiment carried out at a constant heating rate, using the nth order multiple linear regression method and adopting the following Arrhenius type equation:
where n is the order of reaction.
By applying the least squares method, the sum of the squares of the differences between the values of dα/dt calculated by equation (8) and those derived using TG/DTG measurements (taking into account that α is calculated using equation (2) and dα/dt = (dm/dt) /Δmtot where dm/dt is the rate of mass loss and Δmtot= mf–mi) attains its minimum value for given values of A, n and Ea, which are thereby identified. In particular, the value of the activation energy Ea is thus determined.
Further investigation: equation for multilinearity assumption
The left-hand side of equation (9) is linear with respect to both the reciprocal temperature and log(1-α). Therefore, it is customary for the values of the kinetic parameters to be computed following a multiple linear regression (MLR) approach. However, when the values of log(dα/dt), T and log(1-α) are obtained from the TG/DTG curves according to equation (9) and log(dα/dt) is plotted against (1/T) and log(1-α), not all the fitted data points are found to lie on the same plane, as is assumed in the same equation. Therefore, in order to obtain a reliable estimation of kinetic parameters, a suitable number of data points, sufficient to give an accurate fit but not so many as to include deviations from the linearity assumption, must be considered. Accordingly, in this investigation the effect of the number of data points to be used for computing model parameters in Wyden-Widmann’s approach on the quality of the final solution was also studied. This goal was obtained for each TG/DTG curve by computing MLR models based on an increasing number of data points from 21 to 201, symmetrically distributed around the temperature, corresponding to the minimum in DTG cellulose peak.
It may be concluded that: (i) the most correct Ea trend, as a function of the irradiation time in the weatherometer, seems to be the one obtained using the Satava and Arrhenius methods to identify the model function that best fits the experimental trend of the DTG peak values, (ii) A sufficiently correct trend may be obtained also by applying the Wyden–Widmann method but only using a suitable number of points right and left of the DTG peak, (iii) It is possible to choose the appropriate number of points in Wyden-Widmann’s approach by taking into account the multilinearity assumption and observing the temperature interval over which it is followed, (iv) It therefore seems that the inversion of trend found on going from the non aged extra-strong paper sample to the artificially aged samples is real and linked to the non aged paper sample and not to the method used to compute the Ea values, (iv) the choice of the best functional form (in the case of model fitting methods) and of the number of points to use for the computation of the regression parameters in the case of Wyden-Widmann’s method, is of capital importance for a correct final result.
Tomassetti M, Campanella L, Tomellini R: Thermogravimetric analysis of ancient and fresh woods. Thermochim Acta. 1990, 170: 51-65.
Koutznetsov DA, Ivanov AA, Veletsky PR: Analysis of cellulose chemical modification: a potentially promising technique for characterizing cellulose archaeological textiles. J Archaeol Sci. 1996, 23 (1): 23-34. 10.1006/jasc.1996.0003.
Davidson RS: The photodegradation of some naturally occurring polymers. J. Photochem. Photobiol. B:Biol. 1996, 33: 3-25. 10.1016/1011-1344(95)07262-4.
Agarwal UP, Mcsweeny JD: Photoyellowing of thermo mechanical pulps: looking beyond a-carbonyl and ethylenic groups as the initiating structures. J. Wood Chem. Technol. 1997, 17 (1&2): 1-26.
Campanella L, Costanza C, Tomassetti M: Determination of environmental persistence and thermogravimetric analysis of paper artificially aged by photoirradiation. Ann Chim. 2006, 96 (9-10): 575-585. 10.1002/adic.200690059.
Campanella L, Antonelli A, Favero G, Tomassetti M: New archeometric method for wood based on an enzymatic biosensor. L'Actualite' Chimique. 2001, 10: 14-20.
Campanella L, Battilotti M, Costanza C: Studio della degradazione fotochimica della cellulosa della carta. Ambiente, Risorse e Salute. 2004, 99:
Campanella L, Chicco F, Favero G, Gatta T, Tomassetti M: Further applications of a new biosensor method for dating cellulosic finds. Ann Chim. 2005, 95 (3-4): 133-141. 10.1002/adic.200590015.
Campanella L, Favero G, Rodante F, Tomassetti M, Vecchio S: Thermogravimetric and kinetic methods to date wood finds. First results. Ann Chim. 2003, 93 (11): 897-907.
Campanella L, Costanza C, Tomassetti M, Vecchio S: Kinetic processing of thermoanalytical curves of cellulose samples and finds, assessment of possible archaeometric applications. Curr Anal Chem. 2010, 6 (4): 260-268.
Brown ME, Dollimore D, Galwey AK: Reactions in the Solid State. Comprehensive Chemical Kinetics. 1980, Elsevier: Amsterdam, 22: 9786611777760
Satava V: Mechanism and kinetics from non-isothermal thermogravimetric traces. Thermochim Acta. 1971, 2 (5): 423-428. 10.1016/0040-6031(71)85018-9.
Wyden H, Widmann G: Kinetic evaluation of TG measurements with multiple regression. Experentia. Supplementum. Angew. Chem. Thermodyn. Thermoanal. 1979, 37: 284-290.
Wright A: Genetic algorithms for real parameter optimization. Foundation of Genetic Algorithms. Edited by: G. J. Rawlins. 1991, San Mateo, CA: Morgan Kaufmann, 205-218. 1558601708
Alesiani M, Proietti F, Capuani S, Maraviglia B: 13C CPMAS NMR spectroscopic analysis applied to wood characterisation. Appl Magn Reson. 2005, 29 (2): 177-184. 10.1007/BF03167005.
Capitani D, Emanuele MC, Bella J, Segre AL, Attanasio D, Focher B, Capretti G: 1H NMR relaxation study of cellulose and water interaction in paper. TAPPI J. 1999, 82 (9): 117-124.
Paci M, Federici C, Capitani D, Perenze N, Segre AL: NMR study of paper. Carbohydr Polym. 1995, 26 (4): 289-297. 10.1016/0144-8617(95)00011-U.
Missori M, Mondelli C, De Spirito M, Castellano C, Bicchieri M, Schweins R, Arcovito G, Papi M, Congiu Castellano A: Modifications of the mesoscopic structure of cellulose in paper degradation. Phys Rev Lett. 2006, 97 (23): 238001-1-238001-4.
Rodante F, Vecchio S, Catalani G, Guidotti M: Thermal analysis and non-isothermal kinetic study of some pesticides. Part II. Chlorinate derivatives. J Therm Anal Calorim. 2000, 60 (2): 605-622. 10.1023/A:1010159424733.
Garn PD: Thermoanalitical Methods of Investigation. 1965, Academic Press, London & New York, 9780122760501
Doyle CD: Estimating isothermal life from thermogravimetric data. J Appl Polym Sci. 1962, 6 (24): 639-642. 10.1002/app.1962.070062406.
Sestak J: Thermophysical Properties of Solids, Comprehensive Analytical Chemistry. 1984, Elsevier, Amsterdam, 12 D: 9780444996534
Dollimore D, Evans TA, Lee YF, Wilburn FW: Correlation between the shape of a TG/DTG curve and the form of the kinetic mechanism which is applying. Thermochim Acta. 1992, 198 (2): 249-257. 10.1016/0040-6031(92)85081-6.
Authors are indebted with Dr. C. Costanza and Prof. L. Campanella for making available the paper samples analysed which were artificially aged in a weatherometer. This work was funded by “Sapienza University of Rome”, Ateneo and University Projects.
This article has been published as part of Chemistry Central Journal Volume 6 Supplement 2, 2012: Proceedings of CMA4CH 2010: Application of Multivariate Analysis and Chemometry to Cultural Heritage and Environment. The full contents of the supplement are available online at http://journal.chemistrycentral.com/supplements/6/S2.
The authors declare that they have no competing interests.
FM performed chemometric data analysis using kinetic data processed by TG. SV processed TG/DTG data according to model-fitting methods and provided the related kinetic parameters. MT carried out TG/DTG measurements, coordinated the study. MT and SV edited the text and prepared the final draft of the paper. All the authors have read and approved the final manuscript.