- Research article
- Open Access
- Published:

# Non-isothermal kinetics of thermal degradation of chitosan

*Chemistry Central Journal***volume 6**, Article number: 81 (2012)

## Abstract

### Background

Chitosan is the second most abundant nitrogen containing biopolymer in nature, obtained from the shells of crustaceans, particularly crabs, shrimp and lobsters, which are waste products of seafood processing industries. It has great potential application in the areas of biotechnology, biomedicine, food industries, and cosmetics. Chitosan is also capable of adsorbing a number of metal ions as its amino groups can serve as chelation sites. Grafted functional groups such as hydroxyl, carboxyl, sulfate, phosphate, and amino groups on the chitosan have been reported to be responsible for metal binding and sorption of dyes and pigments. The knowledge of their thermal stability and pyrolysis may help to better understand and plan their industrial processing.

### Results

Thermogravimetric studies of chitosan in air atmosphere were carried out at six rates of linear increasing of the temperature. The kinetics and mechanism of the thermal decomposition reaction were evaluated from the TG data using recommended from ICTAC kinetics committee iso-conversional calculation procedure of Kissinger-Akahira-Sunose, as well as 27 mechanism functions. The comparison of the obtained results showed that they strongly depend on the selection of proper mechanism function for the process. Therefore, it is very important to determine the most probable mechanism function. In this respect the iso-conversional calculation procedure turned out to be the most appropriate.

### Conclusion

Chitosan have excellent properties such as hydrophilicity, biocompatibility, biodegradability, antibacterial, non-toxicity, adsorption application. The thermal degradation of chitosan occurs in two stages. The most probable mechanism function for both stages is determined and it was best described by kinetic equations of *n*^{-th} order (*F*_{n} mechanism). For the first stage, it was established that *n* is equal to 3.0 and for the second stage – to 1.0 respectively. The values of the apparent activation energy *E*, pre-exponential factor *A* in Arrhenius equation, as well as the changes of entropy Δ*S*^{≠}, enthalpy Δ*H*^{≠} and free Gibbs energy Δ*G*^{≠} for the formation of the activated complex from the reagent are calculated.

## Background

A number of reviews [1–9] and articles [10–15] have been dedicated to chitin and the products obtained from its chemical treatment at different conditions. Chitosan is the second most abundant nitrogen containing biopolymer in nature, obtained from the shells of crustaceans, particularly crabs, shrimp and lobsters, which are waste products of seafood processing industries [13, 15].

Chitosan, namely poly-β(1–4)-2-amino-2-deoxy-*D*-glucopyranose is a natural high molecular mass biopolymer which is generally obtained by extensive deacetilation of chitin [poly-β(1–4)-2-acetamide-2-deoxy-*D*-glucopyranose] [10]. It occurs as a component of the cell wall of some fungi but it is generally produced by carrying out the deacetilation of chitin. A sharp nomenclature with respect to the degree of N-deacetilation has not been defined between chitin and chitosan. The major procedure for obtaining chitosan is based on the alkaline deacetylation of chitin with strong alkaline solution. When deacetilation is incomplete a mixture of acetamide-polymer and amino-polymer are obtained producing materials with different properties called chitosans. Degree of deacetilation of 70-80% is very common. The degree of deacetylation and the crystallinity of chitosans are the principal characteristics affecting both solubility in aqueous medium and capacity for forming complexes [14]. The structure of cellulose, chitin and chitosan are shown in Figure 1.

Chitin and chitosan are of commercial interest due to their high percentage of nitrogen (6.89%) compare to synthetically substituted cellulose (1.25%). In this respect, chitin and chitosan are recommended as suitable functional materials, because these natural polymers have excellent properties such as hydrophilicity, biocompatibility, biodegradability, antibacterial, non-toxicity, adsorption properties, and remarkable affinity for many biomacromolecules. Chitosan has been widely studied for biosensors, tissue engineering, separation membrane, waste water treatment, because of its multiple functional groups [15].

Recently, much attention has been paid to chitin and chitosan as a potential polysaccharide resource. Due to their special chemical and biological properties and widespread availability chitins and their derivatives have extensive applications in many industrial, and agricultural fields [15, 16]. Chitosan have great potential application in the areas of biotechnology, biomedicine, food industries, and cosmetics. Chitosan is also capable of adsorbing a number of metal ions as its amino groups can serve as chelation sites. Due to their high nitrogen content and porosity, chitosan-based sorbents have exhibited relatively high sorption capacities and kinetics for most heavy metal ions. This could be easily explained by the presence of amino groups in the polymer matrix, which can interact with metal ions in the solution by ion exchange and complexation reactions [8, 15]. The high content of amino groups also makes possible many chemical modifications in the polymer with the purpose of improving its sorptive features, such as selectivity and adsorption capacity. Grafted functional groups such as hydroxyl, carboxyl, sulfate, phosphate, and amino groups on the chitosan have been reported to be responsible for metal binding and sorption of dyes and pigments [17–22]. The knowledge of their thermal stability and pyrolysis may help to better understand and plan their industrial processing.

Studies on mechanism and kinetics of reactions involving solid compounds is challenging and difficult task with complexity results from the great variety of factors with diverse effects, e.g. reconstruction of solid state crystal lattice, formation and growth of new crystallization nuclei, diffusion of gaseous reagents or reaction products materials heat conductance, static or dynamic character of the environment, physical state of the reagents – dispersity, layer thickness, specific area and porosity, type, amount and distribution of the active centers on solid state surface, etc. [23]. Recently, the methods of thermogravimetry (TG), differential thermal analysis (DTA) and differential scanning calorimetry (DSC) are quite useful, since they provide reliable information on the physico-chemical parameters, characterizing the processes of transformation of solids or participation of solids in processes of isothermal or non-isothermal heating [23–28]. The kinetic triplet (i.e., *E* *A* and *f*(*α*)) is typical outcome of the regular kinetic analysis. Practically, the kinetic triplet is needed to provide mathematical description of the process. If the kinetic triplet is determined correctly, it can be used reproduce the original kinetic data as well as to predict the process kinetics outside the experimental temperature region. For example, may be predicting at any desired temperature, the time to reach any extent of conversion [29]. The results, obtained on this basis can be directly applied in materials science for the preparation of various metals and alloys, cements, ceramics glasses, enamels, glazes, synthetic or natural polymers and composite materials [23]. In the case of chitosan, several studies concerning the thermal degradation by means of TG, DTA and DSC techniques were reported and different kinetic values were obtained depending on the experimental conditions in which essay were performed [12–14, 30, 31]. The knowledge of thermal degradation kinetic may help a better understanding and planning of processes to recover metals or metal oxides previously adsorption on chitosan.

The aim of this paper is to assess the kinetic models of chitosan pyrolysis and to estimate the kinetics triplet (*E*, *A* and the shape of the most appropriate *f*(*α*)-function) of this process by means of thermogravimetric data, using different calculation procedures.

## Methods

### Materials

The materials used for thermal degradation chitosan obtained from crab shells were commercially obtained from SIGMA-ALDRICH (Cat. No. C3646), deacetylated with ≥75% and were used without further purification. Before using, chitosan was vigorously grounded in agate mortar and dried in air at 60°C for 4 hours.

### Thermal analysis

The thermogravimetrical measurements (TG-DTG-DTA) were carried out in a flow of synthetic air at a rate of 25 cm^{3} min^{–1} under non-isothermal conditions on an instrument STA 449 F3 Jupiter (Nietzsch, Germany) with its high temperature furnace. Samples of about 7 ±0.1 mg mass were used for the experiments varied out at hearing rates of 3, 6, 9, 12, 15 and 18°C min^{-1} up to 800°C. The samples were loaded without pressing into an open 6 mm diameter and 3 mm high platinum crucible, without using of a standard reference material. The TG, DTG and DSC curves were recorded simultaneously with 0.1 mg sensitivity.

### Theoretical approach and calculation procedures

TG/DTG technique is very useful for the determination of decomposition temperature, decomposition steps and kinetics triplet for solid substances. The kinetics of thermal degradation reactions is described with various equations taking into account the special features of their mechanisms. The reaction rate can be expressed through the degree of conversion *α*, according to the Equation (1):

where: *m*_{i}*m*_{f} and *m*_{t} are the initial, final and current sample mass at the moment *t*, respectively. The rate of many condensed phase reactions can be conveniently parameterized as a function of the temperature *T*, and the extent of the reaction conversion *α*. In non-isothermal kinetics of heterogeneous condensed phase reactions, it is usually accepted that the reaction rate is given by [23–25]:

where: *α* is the degree of conversion, *T* – absolute temperature, *t* – time, *f*(*α*) – the differential conversion function, *R* – the gas constant (8.314 Jmol^{-1} K^{-1}), *q* is the linear constant heating rate *q* = d*T*/d*t* and *A* and *E* are the pre-exponential factor or frequency factor and the activation energy respectively, given by the Arrhenius equation. The conversion function f(*α*) for a solid-state reactions depend on the reaction mechanism and can generally be considered to be as follows [32]:

where *m* *n* and *p* are empirically obtained exponent factors, one of them always being zero. The combination of different values of *m* *n* and *p* make it possible to describe various probable mechanisms [23, 25, 32].

After substitution of Eq. (3) in Eq. (2), separation of variables and integration, the following general equation was obtained:

The solution of left-hand side integral depend on the explicit expression of the function *f*(*α*) and are denoted as *g*(*α*). The formal expression of the function *g*(*α*) depend on the conversion mechanism and its algebraic expression. The latter usually represents the limiting stage of the reaction – the chemical reactions; random nucleation and nuclei growth; phase boundary reaction or diffusion. The algebraic expression of functions of the most common reaction mechanisms operating in solid-state reactions are presented in some papers [23–25, 27, 29].

The right-hand side integral or so called “temperature integral” in Eq. (4) has no exact analytical solution, and can be numerically calculated, making the substitution *u* = *E*/(*RT*) and using the relation:

this difficulty can be overcome [28, 33, 34] may be rewritten as follows:

where *p*(*u*) is the exponential integral. Several author [35–38] suggested different ways to solve this exponential integral.

### Coats-Redfern calculation procedure

Coats and Redfern [35] making some approximations, proposed the next linear equation:

For most values of *E* and for the temperature range over which reactions generally occurs, and because 2*RT*/*E* is much lower than unity, we have reason to write the right side of Eq. (7). If the correct *g*(*α*) function is used, a plot of ln*g*(*α*)/*T*^{2} against 1/*T* should give a straight line from which the values of the activation energy *E* and the pre-exponential factor *A* in Arrhenius equation can be calculated. They can be calculated from the slope and intercept respectively. The model that gives the best linear fit is selected as the chosen model. The integral method of Coats-Redfern has been mostly and successfully used for studying of the kinetics of dehydration and decomposition of different solid substances when 20 < *E*/*RT* < 60 [23, 24, 26–28]. This approach is applied and for single TG curves.

### Iso-conversional method

Iso-conversional methodology in non-isothermal experiments is recommended from ICTAC kinetics committee [39]. Iso-conversional methodology in non-isothermal experiments assumes that for a given degree of conversion α, the reaction mechanism does not depend on the heating rate. In this respect Kissinger-Akahira-Sunose [28, 40] proposed own model-free method and the next linear equation:

For each conversion degree *α*, the linear plot of ln(*q*/*T*^{2}) versus 1/*T* enables *E* and ln*AR*/*g*(*α*)*E*to be determined from the slope and the intercept respectively. If the reaction model *g*(*α*) is a priory known, the corresponding pre-exponential factor can be calculated for each conversion degree. If the reaction model *g*(*α*) is unknown, the following equation may be used to estimate the most probably mechanism function [23, 25]:

where *h*(*u*) is expressed by the fourth Senum and Yang [41] approximation formulae:

Plotting ln*g*(*α*) versus ln*q* and using linear regressive of least square method, if the mechanism studied conforms to certain *g*(*α*) function, the slope of the straight line should be equal to −1.0000 and the linear correlation coefficient *R*^{2} should be equal to unity. Obviously, the values of *E* and *A* do not influenced on the shape of this obtained straight lines.

The Ozawa-Flynn-Wall and Popescu integral methods [42–44], based on Doyle’s approach [28] was not used because it gives similar results to the KAS calculation procedure.

From the theory of activated complex (transition state) of Eyring [23, 24, 26, 27] the following general equation may be written:

where *e* = 2.7183 is the Neper number, *χ* is the transition factor, which is unity for monomolecular reactions, *k*_{B} is the Boltzmann constant, *h* is the Plank constant and *T*_{p} is the peak temperature of DTG curve. The change of the entropy at the formation of the activated complex from the reagent may be calculated according to the formula [23, 24, 26, 27]:

Since

the changes of the enthalpy Δ*H*^{≠} and Gibbs free energy Δ*G*^{≠} for the activated complex formation can be calculated, using well known thermodynamic equation:

The values of Δ*S*^{≠}, Δ*H*^{≠} and Δ*G*^{≠} were calculated at *T* = *T*_{p} (*T*_{p} is the peak temperature at the corresponding stage), since this temperature characterizes the highest rate of the process, and therefore, is its important parameter [23, 24, 26, 27].

### Estimation of lifetime

Lifetime estimations are very useful in the development or selection of polymers for different applications. The lifetime is usually determined by accelerated aging, like air oven aging studies, which require long time periods. Chitosan, being natural biopolymer have own lifetime, which value is very important parameter. The apparent kinetic parameters calculated according above described manner has been used to calculate the value of lifetime for chitosan. The estimated lifetime or the time accelerating ageing *t*_{f} of a chitosan to failure has been defined as the time when the mass loss reaches 5 mass%, i.e. *α* = 0.05 [28, 45–47] and can be estimated by the following equation:

or

where the value of the reaction order *n* is obtained previously. With these equations, the time to equivalent damage at different temperatures can be calculated.

## Results and discussion

The TG, DTG and DTA curves, obtained by the thermal degradation of chitosan at heating rate of 15°C min^{-1} are presented on Figure 2.

Two steps can be observed in the TG-curve: the first one at about 90°C connected with 5.8% mass loss was accompanied by endothermic effect and was attributed to the evaporation of water absorbed in the inner polymer. The second one, beginning at about 245°C and ending at 580°C was connected with 78.5% mass loss and was indicated for vaporization and burning of volatile compounds produced from the thermal degradation of polymeric chain. The pyrolysis of polysaccharides structure starts by a random split of the glycosidic bonds followed by a further decomposition to form acetic and butyric acids and a series of lower fatty acids where C2, C3 and C6 predominate [14]. At temperatures above 380°C, significant change was observed in the course of the TG-curve. It is may be due to the change of the structure of the material and the change of the mechanism of its thermal degradation process. In the DTG and DTA-curves, two peaks were observed – at 294.5 and 544.5°C, respectively. These two stages are strongly exothermic. The first stage ending at 400°C is connected with 43.5% mass loss while the second one – with 35% mass loss. According to some authors [12–14] the first stage is connected with deacetylation and depolymeryzation of chitosan. The second one corresponds to the residual cross-liked degradation chitosan [14].

Using Coats-Redfern calculation procedure and on the basis of different f(*α*) functions, which algebraic expression is presented in our previous work [48] were plotted the dependences of ln(g(*α*)/*T*^{2}) on 1/*T* according to Eq. (7). The values of the coefficient of linear regression *R*^{2} strongly depends form the shape of used f(*α*) function. On Figure 3 are presented the dependences of *R*^{2} versus *n* for *F*_{n} mechanism functions with different values of *n* at a heating rate of 15°C min^{-1}.

As can be seen from Figure 3, the both stages of thermal degradation strongly depends from the values of *n* and are characterized with different most appropriate mechanism functions. The curves presented can be successfully described by empirical polynomials of second power and different coefficients. Differentiating and nullifying these polynomials, the value of *n* at which *R*^{2} has maximum value can be calculated. For the first stage of thermal decomposition was established that the maximal value of *R*^{2} was obtained at *n* = 3, and for the second stage – at *n* = 1 respectively. It means that, the kinetics of the both stages is different and the shapes of the most appropriate mechanism function are different. The same tendency was established for the all used heating rates. The values of *E* and *A* for the both stages of the thermal degradation of chitosan obtained at different heating rates are presented in Table 1 for comparison.

As can be seen from Table 1 the average value of *E* and *A* for the first stage of the thermal degradation of chitosan is higher than these for the second one. Because the obtained value of the apparent activation energy *E* for the second stage are near to 100 kJ mol^{-1}, may be say that the kinetic of this stage is diffusion controlled process.

The TG and DTA-curves obtained at different heating rates are presented in Figure 4 for comparison.

As can be seen from Figure 4, the TG and DTA-curves are shifted to higher temperatures with the increase of the heating rates. Furthermore, two exo-effects were observed, the second of them stronger. These effects correspond to a step which abruptly changes its course at temperatures above 400°C. It means that the kinetic model describing the thermal decomposition (the algebraic expression of *f*(*α*)-function) will also change.

To perform the kinetic computations on the thermal analytical data, the Kissinger-Akahira-Sunose calculation procedure was used, recommended from ICTAC kinetics committee [39]. Taking subsets of the TG-curves at certain values of *α* and taking the corresponding temperatures, the dependence of ln(*q*/*T*^{2}) on 1/*T* was drawn according to Eq.(8) for both stages (Figure 5).

The apparent activation energy *E* was directly evaluated from the slope of these plots and frequency factor *A* – from the cut-off from the ordinate axis respectively.

For determination of the most probable mechanism function *g*(*α*), Eq.(9) was used. Plotting ln[*g*(*α*)] against ln*q* with different algebraic expressions of *g*(*α*)-function and using a linear regression of the least squares method, we were looking for a straight line with slope equal to −1.0000 for which the linear correlation coefficient *R*^{2} is close to unity. The graphical dependence of ln[*g*(*α*)] against ln*q* for both stages is presented in Figure 6.

Lines with slope equal to −1.0000 and linear correlation coefficient *R*^{2} close to unity were obtained with F_{n} mechanism functions at different values of *n*. For the first stage, it was established that *n* is equal to 3.0 and for the second stage - to 1.0 respectively. A reaction order *n* > 2 is mathematically equivalent to Gamma distribution of frequency factor. The interpretation of the high values of reaction order has been discussed in the literature [49]. However, the reason why this parameter has such a high value in relation to the described phenomena has not been completely explained.

The values of the apparent activation energy *E* and pre-exponential factor *A* calculated for both stages of the thermal degradation of chitosan are presented in Table 2. The values of the apparent activation energy are in good agreement with earlier reported ones [13, 14, 49]. Using Eqs. (12), (13) and (14), the values of the change of entropy Δ*S*^{≠}, enthalpy Δ*H*^{≠} and the Gibbs free energy Δ*G*^{≠} for the formation of the activated complex from the reagent were calculated. The change of the activation enthalpy *H*^{≠} shows the energy difference between the reagent and activated complex. If this difference is small, the formation of the activated complex is favored because the potential energy barrier is low. The change of the Gibbs free energy Δ*G*^{≠} reveals the total energy increase of the system at the approach of the reagents and the formation of the activated complex.

As can be see from Table 2, the average value of the apparent activation energy *E* and pre-exponential factor *A* in the Arrhenius equation is higher for the first stage of the thermal degradation of chitosan. The values of the pre-exponential factor for a solid phase reactions are expected to be in a wide range (six or seven orders of magnitude), even after the effect of surface area is taken into account [24]. For first order reactions, the pre-exponential factor may vary from 10^{5} to 10^{16} min^{-1}. The low factors will often indicate a surface reaction, but if the reactions are not dependent on surface area, the low factor may indicate a “tight” complex. The high factors will usually indicate a “loose” complex [23]. Even higher factors (after correction for surface area) can be obtained for complexes having free translation on the surface. Since the concentrations in solids are not controllable in many cases, it would have been convenient if the magnitude of the preexponential factor indicated for reaction molecularity. However, this appears to be true only for non-surface-controlled reactions having low (<10^{8} min^{-1}) pre-exponential factors. Such reactions (if elementary) can only be bimolecular.

The change of entropy for the formation of the activated complex from the reagent Δ*S*^{≠} reflects how close the system is to its own thermodynamic equilibrium. Lower activation entropy means that the material has just passed through some kind of physical or chemical rearrangement of the initial structure, bringing it to a state close to its own thermodynamic equilibrium. In this situation, the material shows little reactivity, increasing the time necessary to form the activated complex. On the other hand, when high activation entropy values are observed, the material is far from its own thermodynamic equilibrium. In this case, the reactivity is higher and the system can react faster to produce the activated complex, and consequently, short reaction times are observed. In particular, for example, the negative values of Δ*S*^{≠} would indicate that the formation of activated complex is connected with decrease of entropy, i.e. the activated complex is “more organized” structure compared to the initial substance and such reactions are classified as “slow” [23]. The negative values of Δ*S*^{≠} obtained for the second stage of the thermal degradation of chitosan showed that its structure is far from its own thermodynamic equilibrium, comparing with the second one.

On Table 3 are given for comparison the kinetic parameters, characterizing thermal degradation of chitin and chitosan.

As can be seen from Table 3 chitin is more stable than chitosan. The values of the apparent activation energy and pre-exponential factor for both stages of thermal degradation of chitin are higher than these of chitosan. The same tendency was established from other authors [12, 49]. The negative values of Δ*S*^{≠} for both compounds showed that the formation of the activated complex from the reagents is connected with a decreasing of entropy, i.e. the activated complex is “more organized” structure and the formation process may be classify as “slow” [23].

The curves representing the dependence of lifetime on the temperature at 5% conversion for both stages of thermal degradation of chitosan are shown in Figure 7.

It is obvious from Figure 7 that the lifetime is a parameter strongly depending on the temperature and decreases exponentially with the increase of the temperature. The lifetime is more sensitive concerning the temperature for the first stage of the thermal decomposition of chitin. The same tendency was established for the lifetime of chitin [48].

## Conclusion

Chitosan have excellent properties such as hydrophilicity, biocompatibility, biodegradability, antibacterial, non-toxicity, adsorption application. The thermal degradation of chitosan occurs in two stages. The kinetics and mechanism of the thermal decomposition reaction were established using iso-conversional calculation procedure. The most probable mechanism function for both stages is determined and it was best described by kinetic equations of *n*^{-th} order (*F*_{n} mechanism). For the first stage, it was established that *n* is equal to 3.0 and for the second stage – to 1.0 respectively. The values of the apparent activation energy *E*, pre-exponential factor *A* in Arrhenius equation, as well as the changes of entropy Δ*S*^{≠}, enthalpy Δ*H*^{≠} and free Gibbs energy Δ*G*^{≠} for the formation of the activated complex from the reagent are calculated.

## References

- 1.
Ravi Kumar MNV: A review of chitin and chitosan applications. React Funct Polym. 2000, 46: 1-27.

- 2.
Goycoolea FM, Argüelles W, Peniche C, Higuera-Ciapara I: Chitin and chitosan. Developments in Food Science. 2000, 41: 265-308.

- 3.
Shahidi F, Abuzaytoun R: Chitin, chitosan and co-products: Chemistry, production, applications and health effects. Adv Food Nutr Res. 2005, 49: 93-135.

- 4.
Rinaudo M: Chitin and chitosan: Properties and applications. Prog Polym Sci. 2006, 31: 603-632.

- 5.
Guibal E: Heterogeneous catalysts on chitosan-based materials: a review. Prog Polym Sci. 2007, 30: 71-109.

- 6.
Prashanth KVH, Tharanathan RN: Chiti/chitosan: modifications and their unlimited application potential – an overview. Trends Food Sci Technol. 2007, 18: 117-131.

- 7.
Tokura S, Tamura H: Analysis of Glycans/Polysaccharide Functional Properties. Comprehensive Glycoscience from Chemistry to Systems Biology. 2007, 2: 449-475.

- 8.
Miretzky P, Fernandez Cirelli A: Hg(II) removal from water by chitosan and chitosan derivatives: A review. J. Hazard. Mater. 2009, 167 (1–3): 10-23.

- 9.
Pillai CKS, Paul W, Sharma CP: Chitin and chitosan polymers: Chemistry, solubility and fiber formation. Prog Polym Sci. 2009, 34: 641-678.

- 10.
Peniche-Covas C, Argüelles-Monal W, Roman JS: A kinetic study of the thermal degradation of chitosan and a mercaptan derivative of chitosan. Polym Degrad Stab. 1993, 39: 21-28.

- 11.
Pawlak A, Mucha M: Thermogravimetric and FTIR studies of chitosan blends. Thermochm. Acta. 2003, 396: 153-166.

- 12.
Wanjun T, Cunxin W, Donghua C: Kinetic studies on the pyrolysis of chitin and chitosan. Polymer Degrad. Stabil. 2005, 87: 389-394.

- 13.
De Britto D, Campana-Filho SP: Kinetics of thermal degradation of chitosan. Thermochim Acta. 2007, 465: 73-82.

- 14.
Lopez FA, Merce ALR, Alguacil FJ, Lopez-Delgado A: A kinetic study on the thermal behaviour of chitosan. J. Therm. Anal. Cal. 2008, 91 (2): 633-639.

- 15.
Zhou L, Wang Y, Liu Z, Huang Q: Characteristics of equilibrium, kinetics studies for adsorption of Hg(II), Cu(II), and Ni(II) ions by thiourea-modified magnetic chitosan microspheres. J Hazard Mater. 2009, 161: 995-1002.

- 16.
Stawski D, Rabiej S, Herczynski L, Draczynski Z: Thermogravimetric analysis of chitin of different origin. J. Thermal. Anal. Calorim. 2008, 93 (2): 489-494.

- 17.
Varma AJ, Deshpande SV, Kennedy JF: Metal complexation by chitosan and its derivatives: a review. Carbohydr Polym. 2004, 55: 77-93.

- 18.
Guibal E, Toutaud E, Roussy J: Chitosan interactions with metal ions and dyes: dissolved-state vs. solid-state application. World Journal of Microbiolgy & Biotechnology. 2005, 21: 913-920.

- 19.
Crini G: Non-conventional low-cost adsorbents for dye removal: A review. Bioresource Technol. 2006, 97|: 1061-1085.

- 20.
Annadurai G, Ling LY, Lee J-F: Adsorption of reactive dye from aqueous solution by chitosan: isotherm, kinetic and thermodynamic analysis. J Hazard Mater. 2008, 152: 337-346.

- 21.
Gupta VK: Suhas. Application of low-cost adsorbents for dye removal – A review. J. Environ. Manag. 2009, 90: 2313-2342.

- 22.
Kyzas GZ, Kostoglou M, Lazaridis NK: Copper and chromium(VI) removal by chitosan derivatives – Equilibrium and kinetic studies. Chem Eng J. 2009, 152 (2–3): 440-448.

- 23.
Vlaev L, Nedelchev N, Gyurova K, Zagorcheva M: A comparative study of non-isothermal kinetics of decomposition of calcium oxalate monohydrate. J Anal Appl Pyrolysis. 2008, 81: 253-262.

- 24.
Turmanova SC, Genieva SD, Dimitrova AS, Vlaev LT: Non-isothernal degradation kinetics of filled with rice husk ash polypropene composites. eXPRESS Polym. Letters. 2008, 2 (2): 133-146.

- 25.
Atanassov A, Genieva S, Vlaev L: Study of the thermooxidative degradation kinetics of tetrafluoroethylene-ethylene copolymer filled with rice husks ash. Polym.-Plast. Technol. Eng. 2010, 49: 541-554.

- 26.
Boonchom B, Puttawong S: Thermodynamics and kinetics of the dehydration reaction of FePO4·2H2O. Physica B. 2010, 405: 2350-2355.

- 27.
Boonchom B, Thongkam M: Kinetics and thermodynamics of the formation of MnFeP4O12. J. Chem. Eng. Data. 2010, 55: 211-216.

- 28.
He W, Deng F, Liao G-X, Lin W, Jiang Y-Y, Jian X-G: Kinetics of thermal degradation of poly(aryl ether) containing phthalazinone and life estimation. J. Thermal. Anal. Calorim. 2010, 100 (3): 1055-1062.

- 29.
Vyazovkin S: Model-free kinetics. Staying free of multiplying entities without necessity. J. Therm. Anal. Calorim. 2006, 83 (1): 45-51.

- 30.
Nietto JM, Peniche-Covas C, Padron G: Characterization of chitosan by pyrolysis-mass spectrometry, thermal analysis and differential scanning calorimetry. Thermochim. Acta. 1991, 176: 63-68.

- 31.
Cardenas G, Bernal L, Tagle LH: Thermogravimetric studies of chitosan derivatives. Thermochim. Acta. 1992, 195: 33-38.

- 32.
Šesták J, Berggren G: Study of the kinetics of the mechanism of solid-state reactions at increasing temperatures. Thermochim. Acta. 1971, 3 (1): 1-12.

- 33.
Elbeyli IJ, Piskin S: Combustion and pyrolysis characteristics of tuncbilek lignite. J. Therm. Anal. Calorim. 2006, 83 (3): 721-726.

- 34.
Cadenato A, Morancho JM, Fernandez-Francos X, Salla JM, Ramis X: Comparative kinetic study of the non-isothermal curing of bis-GMA/TEGDMA systems. J. Therm. Anal. Calorim. 2007, 89 (1): 233-244.

- 35.
Coats AW, Redfern JP: Kinetic parameters from thermogravimetric data. Nature (London). 1964, 201: 68-69.

- 36.
Madhusudanan PM, Krishnan K, Ninan KN: New approximation for thr p(x) function in the evaluation of non-isothermal kinetic data. Thermochim. Acta. 1986, 97: 189-201.

- 37.
Wanjun T, Juwen L, Hen Z, Zhiyong W, Cunxin W: New temperature integral approximate formula for non-isothermal kinetic analysis. J. Therm. Anal. Calorim. 2003, 74: 309-315.

- 38.
Tang W, Liu Y, Zhang H, Wang C: New approximate formula for Arrhenius temperature integral. Thermochim. Acta. 2003, 408: 39-43.

- 39.
Vyazovkin S, Burnham AK, Criado JM, Perez-Maqueda LA, Popescu C, Sbirrazzuoli N: ICTAC kinetics committee recommendations for performing kinetic computations on thermal analysis data. Thermochim. Acta. 2011, 520 (1–2): 1-19.

- 40.
Kissinger HE: Reaction kinetics in different thermal analysis. J. Anal. Chem. 1957, 29: 1702-1706.

- 41.
Senum GI, Yang RT: Rational approximations of the integral of the Arrhenius function. J. Therm. Anal. 1977, 11: 445-447.

- 42.
Ozawa TA: A new method of analyzing thermogravimetric data. Bul. Chem. Soc. Japan. 1965, 38: 1881-1886.

- 43.
Flynn JH: The “temperature integral”- its use and abuse. Thermochim. Acta. 1997, 300: 83-92.

- 44.
Popescu C: Integral method to analyze the kinetics of heterogeneous reactions under nonisothermal conditions. A variant on the Ozawa-Flynn-wall method. Thermochim. Acta. 1996, 285: 309-323.

- 45.
Roy PK, Surekha P, Rajagopal C, Choudhary V: Thermal degradation of LDPE containing cobalt stearate as pro-oxidant. eXPRESS Polym. Letters. 2007, 1 (4): 208-216.

- 46.
Paik P, Kar KK: Kinetics of thermal degradation and estimation of lifetime for polypropylene particles: Effects of particle size. Polymer Degrad. Stab. 2008, 93: 24-35.

- 47.
Wang H, Tao X, Newton E: Thermal degradation kinetics and lifetime prediction of a luminescent conducting polymer. Polym. Int. 2004, 53 (1): 20-26.

- 48.
Georgieva V, Zvezdova D, Vlaev L: Non-isothermal kinetics of thermal degradation of chitin. J Therm Anal Calorim. 2012, DOI 10.1007/s10973-012-2359-6.

- 49.
Stolarek P, Ledakowicz S: Pyrolysis kinetics of chitin by non-isothermal thermogravimetry. Thermochim. Acta. 2005, 433: 200-208.

## Acknowledgements

The authors would like to express their gratitude to Operating Programme Human Resources Development; Grant agreement: BG051PO001/3.3-05-001 “Science and Business” Funded by*:* EC, European Social Fund for their financial support.

## Author information

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

LV has formulated the research idea and planed the experiment. DZ has carried out the collection of data. VG has processed results and has prepared the figure and tables and has finalized the manuscript. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

## About this article

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Chitosan
- Thermal degradation
- Non-isothermal kinetics
- Kinetics triplet