Development of software for the calculation of kinetic parameters from solid phase reaction data, application to synthesized N’-(4-Chlorobenzoyl)-N,N-di-N-butylthiourea ligand and its nickel (II) complex

In gas and liquid phase reactions, the conversion of starting materials to products is directly controlled by the rate of the reaction, whereas in solid state reactions, it is determined by the rate of diffusion. When working with thermal analysis methods, tens of thousands of data to be obtained from measured differential thermal analysis (DTA), thermal gravimetry (TG), and differential thermal gravimetry (DTG) curves must first be digitized, and then these digitized data must be converted to concentrations and thermal values using appropriate software programs. Companies that produce devices using thermal analysis techniques generally adapt the encrypted software programs they use according to the needs of the user companies. Since the devices work like a black box, it is impossible for users to make any changes to the packaged programs. In addition, since the program cannot be seen, its logic cannot be understood most of the time. The program developed in this study will be put at the service of every researcher, and the system of the program can be used on any Windows Software compatible computer. While preparing the program, an easy-to-understand and flexible MATLAB environment was used. The thermal analysis data of synthesized N’-(4-chlorobenzoyl)-N,N-di-n-butylthiourea (HL) and its Ni(II) complex have been digitalized by measuring with thermal analyses such as DTA, TG with hardware technical equipment. By entering the raw experimental data taken into the prepared software program, thermodynamic data such as enthalpy and entropy can be calculated as well as kinetic parameters such as activation energy, reaction order, and rate constant. In this study, a software program was developed that can be used in the calculation of the desired thermodynamic parameters by using the experimentally measured ΔT, Δm, and dΔm values. The accuracy of the results was compared with the thermal results of the NaN3 reference material.

programming languages. MATLAB also includes a set of ready-made solution packages and functions written for specific applications, called 'toolboxes' . Toolboxes consist of M-files needed to solve a problem. The algorithms used in MATLAB are based on the "least squares method" used in numerical integral and derivative calculations with the help of computers. It is possible to use MATLAB on PC, Macintosh, SUN, UNIX, VAX, Apollo, HP, DECStation, SGI, RS/6000, and Convex computer systems.
The prepared program was applied for calculating kinetic parameters of thermal decomposition of newly synthesized thiourea derivative and its complex. N'-(4-chlorobenzoyl)-N,N-di-n-butylthiourea and the Ni(II) complex of these ligands were synthesized and characterized. Decomposition knetic studies were carried out on the program prepared using thermal analysis data [8].

Calculation of thermal data of the ligand and complex
The ligand and its metal complex were synthesized and characterized as given in the literature [9]. HL and its nickel complex were studied by thermogravimetric analysis from room temperature to 1300 K in nitrogen atmosphere. DTA/ TG/DTG diagrams of HL and NiL 2 complex are shown in Figures 1a and 1b, respectively. From the TG curve of HL, it appeared that the sample decomposes in two stages in the temperature range of 403 to 802 K. The first decomposition occurs between 403 and 459 K, with a mass loss of 10.4%; the second decomposition starts at 459 K, and ends at 802 K, with an 88.9% mass loss. From the corresponding DTA profile, three endothermic peaks are noted, the first between 355 and 390 K, with a maximum at 362 K; the second between 390 and 443 K, with a maximum at 425 K; and the third between 443 and 600 K, with a maximum at 508 K. The first endothermic effect is related to the melting of HL (362 K). The other effects are due to the decomposition of the related compound. The TG curve of the NiL 2 complex shows an initial mass loss in the temperature range of 501-612 K, corresponding to the decomposition of the complex to Ni(SCN) 2 . The mass loss at this stage is attributed to the evolved moieties di-n-butylbenzamide (theoretical mass loss: 75.4%, experimental mass loss: 75.7%). These agree with the literature data [9]. The last decomposition step occurs in the temperature range of 612-1273 K and it corresponds to the formation of Ni 3 S 2 (theoretical mass loss: 88.7%, experimental mass loss: 88.4%) [9][10][11][12][13][14][15][16].

PC program for thermal decomposition kinetics 2.1.1. Content of the program
Based on the Broido kinetic method in the program software, the graph of 1/T versus lnln(1/y) was drawn and a linear graph was obtained [17]. The desired kinetic parameters can be calculated from the obtained graph. The part up to the drawing of this graph is explained step by step.
Pure solids decompose when heated under vacuum, and at least some of the resulting decomposition products are volatile. The reaction can be monitored by continuously measuring the sample mass. Mass changes are evaluated with the following equality: y=N/N 0 =(W t -W s ) /( W 0 -W s ) (2.21) In this equation, W 0 is the initial mass, W t is the mass at any time t, W s is the final mass, and y is the mass fraction of the remaining material. If the pyrolysis is carried out in an isothermal environment, the rate of reaction is given by the following relation: d y / d t = -ky n (2.22) In this equation, n is the order of the reaction and k is the rate constant of the reaction. Rate constant according to Arhenius is given by the following equation: y ∫ 1 dy/y n = -(A/u ) y ∫ 1 e -E/RT dT (2.26) has the expression. Pyrolysis reactions are generally first order. Equation (2.26) can be regulated as follows: y ∫ 1 dy/y n = y ∫ 1 dy/y = -ln y = ln(1/y ) (2.27) The integration of the right-hand side of the equation was explored by Vallet and published as a monograph in 1961 [25]. According to this, the following :

Flow chart of PC program
The PC program was written and compiled using C++ programming language in the MATLAB 5.2 environment. Preparation of the program was in easy and flexible MATLAB concept. By using the TG data of the synthesized ligands and chelates, the program can perform the kinetic analysis of the pyrolysis reactions with the help of the Broido kinetic method [17].
In the program, only the percent amount of mass change and the relevant temperature from the TG analysis results are used as data. By using the prepared software program, raw experimental data and thermodynamic data such as entropy can be calculated as well as kinetic parameters such as activation energy, reaction order, and rate constant. The program accepts the reaction degree as '1' for the Broido kinetic method and the flow diagram is given in Figure 1c.

Application of PC program to NaNO3
In this study, the kinetics of the pyrolysis reaction of NaNO 3 was examined as a reference and the program written in the 'MATLAB 5.2' environment was tested. TG/DTG/DTA analyses were performed in a differential thermal analyzer with analytical grade NaNO 3 and α-Al 2 O 3 reference material obtained from Merck.
Seventeen milligrams of sample was used for the measurement made at 10 K min −1 N 2 gas flow rate and dynamic nitrogen atmosphere using a Pt reaction vessel. The TG/DTG/DTA diagrams of the NaNO 3 compound are given in Figure  1d and the graph of the pyrolysis reaction according to the Broido method is given in Figure 1e. It was determined that the energy value calculated for NaNO 3 (51.2 kJ/mol) was consistent with that in the literature (45.3 kJ/mol).

Implementation of the software program for the thermal analysis data of N,N-di-n-butyl-N'-(4-chloro-benzoyl) thiourea ligand (HL) according to the Broido method
TG/DTA/DTG and GC/MS analyses were interpreted and it was determined that the ligand was decomposed in accordance with the reaction (3.1, 3.1a, 3.1b The decomposition kinetic data of the HL ligand calculated in the software program using the Broido method are given in Table 1 and the Broido graph is given in Figure 2a.

3.2.b
The decomposition kinetic data of the [NiL 2 ] complex, calculated in the software program using the Broido method, are given in Table 2 and the Broido graph is given in Figure 2b.

Discussion
The calculation of thermodynamic parameters using the Broido method is frequently encountered in the literature. For example, in the studies by Meena and Sharma, copper(II) soap complex was synthesized and kinetic parameters were calculated by using this equation. The present study revealed that the values of energy of activation for all the equations applied follow the order: Step III > Step II > Step I [18]. Hai et al. synthesized N-acryloyl-N-phenylthiourea. Thermal degradation kinetic parameters are determined for polymer samples From TGA curves using Broido's, Coats-Redfern and Horowitz-Mitzger methods which provide overall kinetic data. The present study revealed that the values of energy of activation for all the equations applied follow the order: Step II > Step I [19]. A series of functionalized phenolformaldehyde polymer resins have been synthesized by the reactionof 2,4-dihydroxyacetophenone-formaldehyde resin with theamines, such as ethanolamine, aminophenol, ethylenediamine, and propylenediamine in dichloromethane. The calculated values for the activation energy of decomposition are 13.86, 9.98, 14.55, and 9.98 kJ/mol for DAPF-ea , DAPF-ap, DAPF-en, and DAPF-pn, respectively [20]. In this study, the decomposition mechanism and kinetics of the previously synthesized and characterized N,N-di-nbutyl-N'-(4-chloro-benzoyl)thiourea (HL) ligand and the bis(N,N-di-n-butyl-N'-(4-chloro-benzoyl)thioureato)nickel(II) [NiL 2 ] complex were investigated using TG/DTG/DTA and GC-MS combined systems.
TG/DTG/DTA analyses were performed in a differential thermal analyzer. Analytical grade NaNO 3 and α-Al 2 O 3 were used as reference materials. Seventeen milligrams of sample was used and the measurements were made under 10 mL/min N 2 gas flow rate using a Pt reaction pan. It is seen that the calculated energy value (51.2 kJ/mol) for NaNO 3 is compatible with that in the literature (45.3 kJ/mol). It was determined that there is an absolute error of 1.1% between the values in the literature and the calculated values of the NaNO 3 compound. This error is thought to be caused by the calibration of the device and the purity level of the substance. Thus, the accuracy of the kinetic parameters calculated in the program written in MATLAB 5.2 environment according to the Broido method of the decomposition reactions with the NaNO 3 compound was tested.
Calculated kinetic parameters (E*, S, A) of the HL ligand and its Ni(II) complex obtained with the prepared PC software are given in Tables 1 and 2. The calculation method of the prepared PC software can be used with other integral equations besides the Broido equation with some modifications to the software. In addition, the results obtained for the compounds examined in this study can be compared with the results obtained with other methods.