American Journal of Analytical Chemistry
Vol.4 No.12(2013), Article ID:40790,7 pages DOI:10.4236/ajac.2013.412091

Thermal Stability and Decomposition Kinetics of Polysuccinimide

Li Zhang, Mingxing Huang, Cairong Zhou

1School of Chemical Engineering and Energy, Zhengzhou University, Zhengzhou, China

Email: zhanglizibo@163.com, zhoucairong@zzu.edu.cn

Copyright © 2013 Li Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received November 5, 2013; revised December 1, 2013; accepted December 9, 2013

Keywords: Polysuccinimide; Thermal Gravimetric Analysis; Thermal Stability; Decomposition Kinetics

ABSTRACT

The thermal stability and decomposition kinetics of polysuccinimide (PSI) were investigated using analyzer DTG-60 under high purity nitrogen atmosphere at different heating rates (3, 6, 9, 12 K/min). The thermal decomposition mechanism of PSI was determined by Coats-Redfern method. The kinetic parameters such as activation energy (E), pre-exponential factor (A) and reaction order (n) were calculated by Flynn-Wall-Ozawa and Kissinger methods. The results show that the thermal decomposition of PSI under nitrogen atmosphere mainly occurs in the temperature range of 619.15 - 693.15 K, the reaction order (n) was, the activation energy (E) and pre-exponential factor (A) were obtained to be 106.585 kJ/mol and 4.644 × 109 min−1, the integral and differential forms of the thermal decomposition mechanism of PSI were found to be and, respectively. The results play an important role in understanding the thermodynamic properties of polysuccinimide.

1. Introduction

Polyaspartic acid (PASP, CAS181828-06-8) has amino and carboxyl groups, which belongs to the biological macromolecule material and is a kind of polymer of amino acids. Since it has the characteristics of good biocompatibility and biodegradability, PASP has been widely used in industrial and medical fields as a new type of green chemicals [1]. Polysuccinimide (PSI, CAS6899- 03-2) is the intermediate of PASP. The experimental formula and relative molecular mass of its monomer are C4H3O2N and 97.074 g/mol, respectively. The chemical structure of the monomer can be written as:

In view of this, PSI is a kind of linear polyimide with high activity and it is easy to open ring changing into poly asparagine with side chains. With the special property, many kinds of derivatives that were used as drug carriers have been prepared [2].

The thermal stability and decomposition kinetics of Polysuccinimide (PSI) were investigated by TG-DTA method. Kinetic parameters such as activation energy (E) and pre-exponential factor (A) were calculated by FlynnWall-Ozawa (F-W-O) and Kissinger methods. The kinetic mechanism function of thermal decomposition of PSI was established by Coats-Redfern method. Using TG-DTA method to study the stability of drugs has the advantages of less sample dose, short experimental period and reliable results. These data not only play an active role in understanding the thermodynamic properties of PSI but also provide a theoretical basis for practical application.

2. Experimental

2.1. Materials and Instruments

Polysuccinimide (the mass fraction was higher than 99.5%) was purchased from Henan Xinlianxin Fertilizer Limited Co. α-Al2O3 (standard material, Shimadu Company in Japan) was used as standard material in the process of thermal analysis.

The TG-DTA analyzer (type DTG-60, Shimadzu Corporation, Japan) was used to determine the TG-DTA curves of the sample. SPN-500-type nitrogen generator (Hewleet-Packard, Beijing Institute of Technology, China) was used to provide a high purity nitrogen atmosphere for the experimental system of thermal analysis. Fourier Transform Infrared Spectrometer (type WQF, Beijing Beifen-Ruili Analytical Instrument (Group) Co, Ltd) was used to analyze PSI. Gel permeation chromatography (type Agilent1100, Agilent Corporation, America) was used to analyze PSI’s purity and the number-average molecular weight (Mn) and polydispesity index (Mw/Mn).

2.2. Experiment Methods

The PSI sample was dried in the vacuum oven at 378.15 K before analysis. The thermogravimetric measurements were carried out at different heating rates (3, 6, 9, 12 K/min) from room temperature to 873.15 K under high purity nitrogen atmosphere (20 ml/min). Mass of each powdered sample was about 4 - 5 mg.

2.3. Characterization of PSI

2.3.1. FTIR Analysis of PSI

FTIR spectra of PSI were recorded in the wave number range of 4000 - 500 cm−1. The result is shown in Figure 1. From Figure 1, we can know that the carboxyl group of branched chain in the opened ring appears at 1400 cm−1, the carbon-carbon bond appears at 1165 cm−1, the coupled carbonyl group appears at 1720 cm−1 and the carbonyl group in the ring obviously appears at 1797 cm−1.

2.3.2. Gel Permeation Chromatography Analysis of PSI

The analysis diagram of PSI by gel permeation chromatography is shown in Figure 2. The analysis results show

Figure 1. FTIR spectra of PSI.

Figure 2. The analysis diagram of PSI by gel permeation chromatography.

that the number-average molecular weight (Mn) of PSI is 2.7941 × 104 g/mol, the polydispersity index of relative molecular mass (Mw/Mn) of PSI is 21.644 and Polymerization degree of PSI is obtained to be 287.5.

2.4. Theoretical Analysis

There are two major categories to study the thermal decomposition kinetics of polymer, which are differential method and integral method. The differential method contains Kissinger [3-5], Caroll-Freeman [6] and Friedman [7] methods while Coats-Redfern [8-12], Doyle [13] and Flynn-Wall-Ozawa (F-W-O) [14-16] methods belong to the integral method. Coats-Redfern, Flynn-Wall-Ozawa (F-W-O) and Doyle methods are usually used for investigating thermal decomposition kinetics of polymers.

In general, the thermal decomposition of a solid polymer under inert atmosphere can be summarized as: Bsolid → Csolid + Dgas. Polymer is finally decomposed into solid residue (C) and gaseous matter (D).

The kinetic analysis of solid-state sample is usually given by Equation (1) [17]

(1)

where is the rates of conversion; is the conversion degree that can be defined as in which m0 and mf are the initial and final masses of the sample, respectively; mt is the mass of the sample at time t (or temperature T) of the decomposition process, mg. In Equation (1), is the temperature dependent rates constant and is normally assumed to obey the Arrhenius equation:

(2)

where A is the pre-exponential factor (min−1), E is the activation energy (kJ/mol) of the kinetic process, R is the gas constant (8.314 J/(mol∙K)) and T is the absolute temperature (K).

Moreover, taking into account the heating rates β =

dT/dt under non-isothermal condition, can be described by Equation (3):

(3)

Equation (4) can be obtained by combing Equations (1)-(3), which describes the thermal decomposition kinetics. Based on TGA data, the kinetic parameters can be calculated from Equation (4). [18]:

(4)

2.4.1. Kissinger Method

The formula of Kissinger method is given by Equation (5) [3-5]:

(5)

where i = 1, 2, 3, 4 (or even more); Tpi is the peak temperatures of different DTA curve at different heating rates. By plotting versus 1/Tpi, the activation energy Ek and pre-exponential factor Ak can be calculated based on its slope (−Ek/R) and intercept ln(AkR/Ek), respectively.

2.4.2. Flynn-Wall-Ozawa (F-W-O) Method

F-W-O method can be used directly to calculate the activation energy E. The integral formula of F-W-O method [14-16] is showed in Equation (6):

(6)

In Equation (6), since the value of

is approximately constant when the values of α are the same at the different heating rates βi, so it is easy to obtain values of E by plotting lgβ against 1/T at the certain conversion degree α.

2.4.3. Coats-Redfern Method

Coats-Redfern method can be described as Equation (7). [8-10,12]:

(7)

where g(α) comes from one of 34 forms of integral formula in the literature which are shown in Table 1 [11]. From Equation (7), the values of both Ec and Ac can be obtained for any selected g(α) and fixed βi (i = 1, 2, 3, 4, even more). The calculating steps in detail are: 1) choose the same αj for each heating rates βi and calculate the corresponding g(αij); 2) give the corresponding temperature Tij according to αij; 3) describe the sketch of

versus 1/Ti for each fixed βi, and calculate the values of both Ec and ln(Ac) according to the slope and intercept of the line, respectively; 4) determine the mechanism function of the thermal decomposition process. Generally, the selected mechanism functions should meet all the conditions, which are: 1) 0 < Ec < 400 kJ/mol; 2)

where E0 come from F-W-O method 3) in which ln(Ak) are obtained from Kissinger method. If g(α) meet the requirements mentioned above, g(α) can be regarded as the probable mechanism function for the thermal decomposition process.

2.5. Determination of High Temperature Heat-Resistance of PSI

The activation energy data obtained by above-mentioned methods can be used to evaluate the high temperature heat-resistance of PSI. The relationship between activa-

Table 1. 34 types of thermal decomposition mechanism functions.

tion energy and lifetime of polymer can be expressed by Equation (8) [19]:

(8)

where tf is the lifetime of PSI at the temperature of Tf, Tα is the temperature at the conversion degree α. Combined Equations (4) and (8), Tf can be shown as:

(9)

3. Results and Discussion

3.1. Thermal Decomposition of PSI

TG and DTA curves of PSI sample are shown in Figures 3 and 4. The thermal analysis data are summarized in Table 2. In Figure 3, the thermal decomposition temperature at different heating rates increases with the increasing of. That indicates the decomposition temperature is affected by the heating rate. Besides, each TG curve of PSI under nitrogen atmosphere has an obvious mass loss process where percentage of mass loss

Figure 3. The TG curves of PSI in nitrogen atmosphere.

Figure 4. The DTA curves of PSI in nitrogen atmosphere.

Table 2. Basic data of kinetics of PSI from TG.

increases gradually with the increasing of heating rates and TG curve moves toward to the right.

From Figures 3 and 4, the thermal decomposition of PSI under nitrogen atmosphere mainly occurs in the temperature range of 619.15 to 693.15 K. Mass loss is accompanied by heat absorption, so thermal decomposition process and solid rearrangement reaction may simultaneously happen in this process.

3.2. Non-Isothermal Kinetic of PSI

Plots of against 1/Tp of PSI by Kissinger method is shown in Figure 5. The values of activation energy Ek and exponential factor ln(Ak) calculated by Kissinger method are 143.874 kJ/mol and 23.897 min-1, respectively. The linear correlation coefficient (R2) is 0.9852.

The activation energies calculated at different conversion degree using Flynn-Wall-Ozawa method are shown in Table 3 and the relationship between E0 and conversion degree α are showed in Figure 6.

The results show that the correlation coefficients (R2) of FWO method is better, the values of activation energy is between 81.509 and 122.205 kJ/mol and increase with the increasing of conversion degree. From Table 3 and Figure 6, the average value of activation energy Eo is 104.202 kJ/mol and the pre-exponential factor lgAo is 9.967.

The kinetic parameters calculated by Coats-Redfern method are listed in Table 4. The values of activation energy and exponential factor calculated by Coats-Redfern method compare respectively with the average value of activation energy calculated using FWO method and the value of pre-exponential factor obtained by Kissinger method, the results are listed in Table 5.

The thermal decomposition process of PSI in the temperature stage of 619.15 to 693.15 K is consistent with the sequence number 14 in Table 1, this is because Ec and ln(Ac) meet better with the conditions of both

and

than others; besides, the relative coefficients (R2) are much better. The integral and differential forms of the mechanism function are

and, respectively. The values of Ec and ln(Ac) are obtained to be 108.967 kJ/mol and 16.802 min−1. So the activation energy and pre-exponential factor of thermal decomposition process of PSI can be regarded as 106.585 kJ/mol and 4.644 × 109 min−1 (i.e. average values obtained by the Coats-Redfern and FWO methods). The thermal decomposition kinetic equation of PSI can be described as:

Figure 5. Plots of and 1/Tp of PSI of Kissinger method.

Table 3. The activation energy Eo and exponential factor lgAo of PSI.

Figure 6. FWO curves of PSI at different conversion.

Table 4. Results of 34 types of kinetic equations of PSI calculated with Coats-Redfern method.

Table 5. Calculation of 6 types of kinetics equations for PSI with Coats-Redfern method.

(10)

3.3. High Temperature Heat-Resistance of PSI

Selecting tf = 60 s and α = 15% as the evaluation indexes of the high temperature heat-resistance of the polymer. According to Equation (12), value of Tf of PSI in nitrogen atmosphere is 259.27˚C. The result shows that the high temperature heat-resistance of PSI is not very well.

4. Conclusion

The thermal behavior of polysuccinimide under nonisothermal condition was investigated using TG-DTA method at different heating rates in nitrogen atmosphere. The results show that the thermal decomposition of PSI under nitrogen atmosphere mainly occurs in the temperature range of 619.15 - 693.15 K, the reaction order (n) was 3/4, the activation energy (E) and pre-exponential factor (A) were obtained to be 106.585 kJ/mol and 4.644 × 109 min−1, the integral and differential forms of the thermal decomposition mechanism of PSI were found to be and, respectively.

REFERENCES

  1. A. W. Yang, G. P. Cao and M. H. Zhang, “Synthesis of Polysuccinimide and Determination of the Intrinsic Viscosity,” Polymer Materials Science & Engineering, Vol. 26, 2010, pp. 4-7.
  2. J. H. Jeong, H. S. Kang, S. R. Yang and J. D. Kim, “Polymer Micelle-Like Aggregates of Novel Amphiphilic Biodegradable Poly(Asparagine) Grafted with Poly(Caprolactone),” Polymer, Vol. 44, No. 3, 2003, pp. 583-591. http://dx.doi.org/10.1016/S0032-3861(02)00816-9
  3. A. Rotaru, M. Anca, G. Popa, P. Rotaru and E. Segal, “Non-Isothermal Kinetics of 2-Allyl-4-((4-(4-Methylbenzyloxy)Phenyl) Diazenyl) Phenol in Air Flow,” Journal of Thermal Analysis and Calorimetry, Vol. 97, No. 2, 2009, pp. 485-491. http://dx.doi.org/10.1007/s10973-009-0214-1
  4. H. E. Kissinger, “Variation of Peak Temperature with Heating Rate in Different Rate in Differential Thermal Analysis,” Journal of Research of the National Bureau of Standards, Vol. 57, No. 4, 1956, pp. 217-221. http://dx.doi.org/10.6028/jres.057.026
  5. H. E. Kissinger, “Reaction Kinetic in Differential Thermal Analysis,” Analytical Chemistry, Vol. 29, No. 11, 1957, pp. 1702-1706. http://dx.doi.org/10.1021/ac60131a045
  6. E. S. Freeman and B. Carroll, “The Application of Thermoanalytical Technique to Reaction Kinetics,” Journal of Physical Chemistry, Vol. 3, 1958, pp. 394-397.
  7. H. L. Friedman, “Kinetics and Gaseous Products of Thermal Decomposition of Polymers,” Journal of Macromolecular Science: Part A—Chemistry, Vol. 1, No. 1, 1967, pp. 57-59. http://dx.doi.org/10.1080/10601326708053917
  8. A. W. Coats and J. P. Redfern, “Kinetic Parameters from Thermogravimetric Data,” Nature, Vol. 201, 1964, pp. 68-69. http://dx.doi.org/10.1038/201068a0
  9. C. R. Zhou, Q. H. Li and H. F. Wang, “Thermal Analysis for the Thermal Decomposition of Methylsulfonate Tin,” Journal of Chemical Engineering of Chinese Universities, Vol. 20, 2006, pp. 669-672.
  10. C. R. Zhou, X. H. Shi, H. F. Wang and D. G. Jiang, “Thermal Decomposition and the Non-Isothermal Decomposition Kinetics of DL-2-Naproxen,” Journal of Chemical Engineering of Chinese Universities, Vol. 25, 2011, pp. 442-446.
  11. L. G. Lu, Q. Zhang, X. N. Xu, X. L. Dong and D. W. Wang, “Thermal Degradation Kinetics of Novel Intumescent Flame Retardant Polypropylene,” China Plastics, Vol. 23, 2009, pp. 53-60.
  12. C. Y. Ou, C. H. Zhang, S. D. Li, L. Yang and J. J. Dong, “Thermal Degradation Kinetics of Chitosan-Cobalt Complex as Studied by Thermogravimetric Analysis,” Carbohydrate Polymers, Vol. 82, No. 4, 2010, pp. 1284-1289. http://dx.doi.org/10.1016/j.carbpol.2010.07.010
  13. C. D. Doyle, “Kinetic Analysis of Thermogravimetric Data,” Journal of Applied Polymer Science, Vol. 5, No. 15, 1961, pp. 285-292. http://dx.doi.org/10.1002/app.1961.070051506
  14. T. Ozawa, “Kinetic Analysis of Derivative Curves in Thermal Analysis,” Journal of Thermal Analysis, Vol. 2, No. 3, 1970, pp. 301-310. http://dx.doi.org/10.1007/BF01911411
  15. F. X. Chen, C. R. Zhou and G. P. Li, “Study on Thermal Decomposition and the Non-Isothermal Decomposition Kinetics of Glyphosate,” Journal of Thermal Analysis and Calorimetry, Vol. 109, No. 3, 2012, pp. 1457-1462.
  16. Q. F. Wang, L. Wang, X. W. Zhang and Z. T. Mi, “Thermal Stability and Kinetic of Decomposition of Nitrated HTPB,” Journal of Hazardous Materials, Vol. 172, No. 2-3, 2009, pp. 1659-1664. http://dx.doi.org/10.1016/j.jhazmat.2009.08.040
  17. C. D. Gamlin, N. K. Dutta, N. R. Choudhury, D. Kehoe and J. Matisons, “Evaluation of Kinetic Parameters of Thermal and Oxidative Decomposition of Base Oils by Conventional, Isothermal and Modulated TGA, and Pressure DSC,” Thermochimica Acta, Vol. 392-393, 2002, pp. 357-369. http://dx.doi.org/10.1016/S0040-6031(02)00121-1
  18. X. Y. Li, Y. Q. Wu, D. H. Gu and F. X. Gan, “Thermal Decomposition Kinetics of Nickel(II) and Cobalt(II) Azo Barbituric Acid Complexes,” Thermochimica Acta, Vol. 493, No. 1-2, 2009, pp. 85-89. http://dx.doi.org/10.1016/j.tca.2009.04.010
  19. Z. W. Zhou and Q. X. Wu, “Studies on Thermal Properties of Poly(Phenylene Sulfide Amide),” Journal of Applied Polymer Science, Vol. 66, No. 7, 1997, pp. 1227- 1230. http://dx.doi.org/10.1002/(SICI)1097-4628(19971114)66:7<1227::AID-APP2>3.0.CO;2-I