Reaction Rate Models for the Thermal Decomposition of Ibuprofen Crystals

Kinetics of the decomposition of racemic ibuprofen crystals were studied by non-isothermal analysis. Thermogravimetric analysis revealed that ibuprofen is thermally stable up to 152.6 ̊C and the initial loss of mass was due to evaporation only. Activation energy, pre-exponential factor, activation entropy and Gibbs free energy for the decomposition of ibuprofen were determined using the integral method of Coats-Redfern (CR). Geometrical contraction models were found to be the best fits. The Arrheinus equation for the thermal decomposition of ibuprofen is k = (1.1 × 107) e–79125/RT sec–1.


Introduction
Thermal analysis is one of the most frequently used instrumental techniques in the pharmaceutical research.TGA/ DTA and differential scanning calorimetry help to solve pharmaceutical problems like the purity level, qualitative and quantitative analysis of drug formulations and polymorphism [1] [2].Also, determination of the kinetic parameters helps to estimate the thermal stability [3] and hence the shelf-life of drugs.
Kinetic analysis has been the subject of interest for many investigators involved in the field of thermal decomposition.Solid-state kinetics studies have increasing importance in thermal analysis and being helpful in the calculation of the parameters of Arrhenius equation and to determine the mechanism of decomposition reaction.Various methods were proposed for the kinetic study of the thermal decomposition of materials and are generally classified under model-fitting and model-free methods.In each case, isothermal and non-isothermal data can be used.These analyses can provide valuable information about the shelf-life period and conditions of storage.
Ibuprofen or 2-(4-isobutylphenyl)-propionic acid belongs to a group of drugs called nonsteroidal anti-in-flammatory drugs (NSAIDs).For many years, it has been used in the treatment of rheumatoid arthritis and other inflammatory diseases.Recently, renewed interest in this compound has arisen as new therapeutic properties have been discovered.Ibuprofen was shown to protect neurons from glutamate toxicity in-vitro [4].The kinetic parameters of ibuprofen were determined by Bogdan Tita et al. [5] earlier.A renewed attempt is made by the authors.Since single crystal is the purest form of a compound, ibuprofen was synthesized from crystalline powder.Presently ibuprofen single crystals were grown by gel growth method.The grown crystal was confirmed using single crystal X-ray diffraction analysis.The kinetic parameters were determined for various solid state reaction models using the CR method under non-isothermal condition.The best fitting and suited model was determined using these parameters.Gibbs free energy and the entropy of activation were also calculated.

Experimental Details
Ibuprofen (99.9%) was purchased from Gopish Pharma Ltd., India.Silica gel prepared from an aqueous solution of sodium metasilicate (Na 2 SiO 3 .5H 2 O) (s.d.fine-Chem Ltd, Mumbai) was used as crystal growth medium.Organic solvent assays of 99% (Spectrum Reagent & Chemicals Pvt, Ltd., Cochin) were used for the preparation of binary mixtures.An aqueous solution of sodium metasilicate of specific gravity 1.06g/cm 3 was prepared and the pH of the solution was adjusted using 5M acetic acid (Nice Chemicals Pvt.Ltd, Cochin).2% w/v of ibuprofen in methanol-diethyl ether mixture (0.5:1) solution was placed over the set gel taken in a test tube.The pH of the gel solution was varied to get the best conditions for crystal growth and the crystal growth experiments were carried out at ambient temperature (≈ 27˚C).Transparent, platy crystals crystallized (Figure 1).The crystals were removed after a growth period of four weeks.Enraf-Nonius CAD-4 Diffractometer, with MoKα (λ = 0.71073 Å) radiation was used to obtain the accurate cell parameters of the grown ibuprofen single crystal at room temperature.Cell parameters were obtained from least-squares refinement of the setting angles of 87 reflections.
TGA/DTA of ibuprofen (Figure 2) was carried out using a Perkin Elmer Pyris Diamond thermogravimetric/differential thermal analyzer instrument with an alumina crucible, applying a heating rate 10˚C min -1 in the temperature range of 40˚C -610˚C, under nitrogen atmosphere with a flow rate 40 mL min -1 .Crystalline powder of mass 10.776 mg of the grown crystals was used for the thermal analysis.
The TGA/DTA curve is shown in Figure 2. Two non-symmetrical endothermic peaks were observed in the DTA curve at 76.2˚C and 255˚C.The endothermic peak at 76.2˚C, and the absence of loss of mass in the TGA indicates that it is the melting point of ibuprofen.Also the melting point of racemic ibuprofen reported earlier [7] agrees with the experimental data.The other peak at 255˚C indicates the boiling point.The evaporation of ibuprofen is intiated at 152.6˚C, and is followed by the exothermic behaviour between 189.7˚C and 255˚C.
Thermodynamical parameters, enthalpy of fusion (∆H f ) at the melting point and enthalpy of vaporization(∆H vap ) were determined using the Gibbs function and by Trouton's rule as 28.10 kJ/mol and 46.47 kJ/mol  respectively.It is found that the activation energy of ibuprofen is twice that for enthalpy of vaporization (Table 1).This suggests that ibuprofen continue the dimer structure in liquid state and becomes a monomer in vapour state.The dimer structure is also confirmed by single crystal X-ray diffraction data.Hence the loss of mass of 50% may be attributed to evaporation of ibuprofen molecules in liquid phase, which is complete at a temperature nearly 245˚C.This behavior is also observered in the work of Lerdkanchanaporn et al [8].The thermal decomposition of ibuprofen starts at 255˚C, transforming into a monomer and decompose with the release of CO 2 molecule and hydrogen from the carboxylic acid chain accounting for a loss of mass of 21.8% that happens in the neighbourhood of the boiling point.Further decomposition results in a loss of mass of 7.76% and is due to the release of CH 4 molecule at 274˚C, leaving carbonaceous residue behind.The material is found to be thermally stable up to 152.6˚C.

Kinetic Parameters
Kinetic parameters primarily depend on the rate of reaction.The TGA data is useful in the determination of the extent of reaction [9].Kinetic analysis of decomposition process involves the reaction model and the Arrhenius equation.The reaction models are classified as nucleation model, diffusion model and reaction order model.Each model depends on the function extent of reaction (α).The extent of reaction was determined from the TGA/DTA (Figure 2) using the relation, where m t is the mass of the sample at particular temperature.m o and m f are the initial and final masses of the reaction respectively.Pre-exponential factor and activation energy were calculated for every model and the goodness of fit is analyzed.Kinetic parameters were determined from solid state rate reaction and integral rate law.The general form of solid state rate reaction Equation (2) and integral rate law Equation (3) are: ( ) ( ) ( ) where t is time, T-temperature, f(α)-reaction model and g(α) is the integral reaction model.The rate constant k is determined using Arrhenius equation.Equations ( 2) & (3) are modified using the Arrhenius equation as, where A and E are the pre-exponential factor and activation energy; R is the gas constant.
The Arrhenius parameters (A and E) together with the reaction model are called as the kinetics triplet.Under non-isothermal condition equations incorporating β the heating rate ( 4) & ( 5) can be written as, x" and transforming the limits in the Equation (7), ( ) The exponential integral is replaced with p(x) which is a standard one, has no analytical solution but has approximations.

( ) ( )
Kinetics parameters may be obtained from non-isothermal rate laws by both model-fitting and model-free methods.Model-fitting methods involve fitting of different models to the extent of reaction-temperature curves and simultaneous determination of the activation energy and frequency factor.Though there are several non-isothermal model-fitting methods, one of the most popular being the CR method [10] was used.It is an integral method that assumes various order of reaction and compares the linearity in each case to select the correct order.This method uses the asymptotic series expansion for approximating the exponential integral in Equation (9) and hence becomes, ( ) T is the mean experimental temperature.Plotting ( ) (10), gives E and A from the slope and intercept respectively.The model that gives the best linear fit is selected as the chosen model for the decomposition.Different reaction models available in the literature are listed in Table 2 [11] [12].Arrhenius Parameters for non-isothermal decomposition of Ibuprofen were also determined for various models.Kinetic analysis of ibuprofen was made using the TGA/DTA (Figure 1) and the thermal degradation ibuprofen in the range 152.6˚C to 255˚C was chosen as it is linear.The graphs between ln{g(α)/T 2 } vs. 1/T of Equation (10) for the models 1 -3 and 4 -12 are plotted and depicted in Figures 3(a) & 3(b) respectively.Two Table 2. Solid state reaction models.

M. No. Reaction Model f(α) g(α)
Nucleation models Diffusion models Reaction order and geometrical contraction models separate plots were presented to view the response distinctly.The activation energy and pre-exponential factor were determined from the graph and are listed in Table 1.The activation energy should be positive for the absence of phase transition [13].The calculated activation energies are positive for all the models, indicating that no phase transition took place in the selected temperature range.Moreover the best fitting model is decided from the correlation coefficient (r 2 ) and it should approach 1.Though all the models seem to be better fits for α < 0.80, the reaction order models 8, 11 and 12 seem to be the best fitted models for the decomposition of ibuprofen in the chosen temperature range.But the Arrhenius parameters obtained from the models 11 and 12 almost matched with that of the reported values [5] [14].Hence the mean value of activation energy and preexponential factors these models are presented and are: 79.1 kJ/mol and 1.162 × 10 7 sec -1 respectively.The Arrhenius equation for the thermal decomposition of ibuprofen is k = (1.1 × 10 7 ) e -79125/RT sec -1 .This equation agrees with that of the values reported earlier [15].Further a graph (Figure 4) is drawn between the activation energy (E) and pre exponential factor (A) obtained in Table 1.This indicates that there is a linear relation between E and A.

Thermodynamic Parameters
Entropy of activation is defined as the difference between the entropy of the transition state and the sum of entropies of the reactants.Free energy of change of decomposition is the driving force of a chemical reaction and hence is the difference between the enthalpy of the transition state and the sum of enthalpies of the reactants in the ground state.The Thermodynamic parameters of ibuprofen: entropy of activation (∆S), enthalpy of activation (∆H) and Gibbs free energy (∆G) were calculated using Equations ( 11)-( 13).
where h and k are Planck's constant and Boltzmann constant respectively.

H E RT
The calculated rate constant and thermodynamic parameters or different models are listed in Table 3.It is observed that the entropy of activation increases with the increase of activation energy.The positive value of ∆G shows that the reaction involved in the decomposition of ibuprofen is not spontaneous.It is noted that the activation energy is nearly equal to enthalpy of activation (Table 3) suggesting that ibuprofen is in condensed phase between 152.6˚C to 255˚C.The mean values of ∆S, ∆H and ∆G for the models 11 & 12 are:39.60J/K mol, 75.17 kJ/mol and 56.32 kJ/mol respectively.

Conclusion
Single crystals of ibuprofen were grown by gel growth method.Monoclinic form of the crystal was confirmed  using single crystal X-ray diffraction analysis.Thermal decomposition kinetics of ibuprofen was analyzed.Ibuprofen is thermally stable up to 152.6˚C.It is concluded that ibuprofen exists in dimer before the boiling point and in monomer after the boiling point.The thermal decomposition of ibuprofen occurred in a single step after the melting point.This happened in the form of evaporation only up to the boiling point.The kinetic parameters were determined using CR method.Thermodynamic parameters also confirmed that the decomposition of ibuprofen is due to non-spontaneous reaction.Geometrical contraction models were found to be the best fits.

Figure 4 .
Figure 4. Plot of E vs. ln A.

Table 1 .
Arrhenius parameters for non-isothermal decomposition of ibuprofen.

Table 3 .
Thermodynamic parameters of ibuprofen