DOI:
10.1039/D3RA07045J
(Paper)
RSC Adv., 2023,
13, 34972-34986
Insights into a co-precursor driven solid-state thermal reaction of ferrocene carboxaldehyde leading to hematite nanomaterial: a reaction kinetic study†
Received
16th October 2023
, Accepted 17th November 2023
First published on 1st December 2023
Abstract
Thermal decomposition of a mixture of ferrocene carboxaldehyde and oxalic acid dihydrate in O2 atmosphere produced rod-like hematite nanomaterial. The decomposition reaction was complex as evident from the overlapped multistep reaction steps in the non-isothermal thermogravimetry (TG) profiles obtained in the 300–700 K range. A peak deconvolution method was applied to separate the overlapped reaction steps. The multistep TG profiles were successfully deconvoluted, which showed that the decomposition occurs in six individual steps. However, it was found that only the last three reaction steps were responsible for the production of hematite. To estimate the activation energy values for these thermal reactions, six model-free integral isoconversional methods were used. The activation energy value significantly depends on the extent of conversion in each step; however, the nature of its dependence significantly different for each step. The most probable stepwise reaction mechanism functions for the solid-state reactions were obtained using the master plot method. The reaction mechanism was found to be different for different steps. Utilizing the activation energy and reaction mechanism function, the reaction rates of decomposition for each step were determined. To substantiate the validity of the assumed kinetic models, the experimental conversion curves were compared with the constructed ones, and the agreement was quite reasonable. The conversion-dependent thermodynamic parameters were obtained utilising the estimated kinetic parameters. Role of the co-precursor in the thermal reaction of the precursor was plausibly revealed. The present study describes how the use of a co-precursor significantly enhances the thermal decomposition of the precursor, how hematite nanomaterials can be synthesized from a co-precursor driven solid state reaction at low temperatures, and how the kinetic calculations facilitate the understanding of the solid-state reaction process. This study proposes the use of a suitable combination of precursor and co-precursor for solid-state thermal synthesis of iron-based nanoparticles using organo-iron compounds as precursor and also illustrates the effective application of the thermal analysis technique to understand the decomposition reaction.
1. Introduction
Among several naturally occurring iron oxides, hematite (α-Fe2O3) possesses the highest thermodynamic stability. Depending on the size and shape, hematite at nanoscale shows a large number of applications.1–4 Recent studies demonstrated that fundamental properties of bulk materials can be effectively improved and customized by tuning the size, shape, composition, and crystallinity in nanomaterials.3–9 Hematite is an environmentally compatible, cost-effective, and good corrosion-resistive narrow band gap semiconductor. These striking properties led to the fabrication of various nanostructures of hematite (e.g., nanorods, nanotubes, nanocubes, nanospheres, and nanoflakes) by employing a variety of synthesis methods to allow for major applications in magnetic storage devices,10 biomedicine,11 engineering and industrial fields,3 lithium-ion batteries,12 gas sensors,13 water treatment,14 and magnetic resonance imaging (MRI).15 However, the physical property of these nanomaterials strongly depend on the size, shape, and crystallinity, which are dependent on the synthesis method.16 There are numerous well-known chemical, physical, and biological methods17 that have been used to synthesize nanomaterials, among which ∼90% of the chemical methods (e.g., hydrothermal, microemulsion, sonochemical, chemical precipitation, electrochemical, and thermal decomposition)4,18,19 are deployed for the preparation of hematite nanoparticles in order to control their size, shape, and crystallinity. Undoubtedly, every method has its merits and demerits. However, among the methods reported so far, the thermal decomposition of organoiron solid precursors is a popular, less expensive technique which is easy to handle.20–23 The thermal decomposition of a solid (upon heating a sample) is a process of redistribution of bonds and formation of new products different from the reactants.24 This technique has a number of merits such as synthesis at comparatively low temperature, control by reaction time, temperature and environment, faster reaction, and use of a variety of iron-bearing compounds as precursors and organic compounds as co-precursors.20–23
Application of kinetic theory of thermal analysis is important to know the underlying process in a solid-state thermal decomposition reaction leading to the formation of new material(s). Among several thermo-analytical methods, the non-isothermal thermogravimetry (TG) is a popular method to explore solid-state thermal reactions,21,25,26 where analysis of the TG profiles (i.e., mass vs. temperature plots obtained at constant heating rates) based on iso-conversional procedure allows the evaluation of the kinetic parameters (activation energy, reaction mechanism function and reaction rate) and predicts the mechanism of the rate-controlled conversion reaction. The iso-conversional model free integral methods21,23,25,27–32 are mostly employed to evaluate the activation energy quite accurately over a range of temperatures utilizing the TG data obtained at different heating rates, and these integral methods are dependent on certain approximations used in solving specific temperature integrals.
Among various organo-iron compounds, ferrocene, (C5H5)2Fe and its derivatives are good sources of iron for the formation of iron oxide nanoparticles through thermal decomposition.33–37 Ferrocene decomposes completely at ∼450 K, but produces a hematite nanomaterial when decomposed in presence of oxalic acid.38 Acetyl ferrocene,33 on decomposition in a nitrogen atmosphere, produced hematite nanoparticles too. Similarly, decomposition of 1-ferrocenyl ethanol in air led to hematite nanoparticles,21 while ferrocene carboxaldehyde on thermal decomposition in oxygen atmosphere yielded a much smaller size (av. diameter ∼5 nm) of dot-shaped hematite nanoparticles.37 Apart from the reaction environment, additional presence of a co-precursor influences the decomposition process of ferrocene materials, affecting the reaction kinetics.23,39,40 For a co-precursor driven thermal decomposition, the co-precursor is chosen in such a way that it decomposes prior to the decomposition of the precursor, so that the decomposed products may react with the intermediate decomposed products from the precursor. In a previous article,25 the kinetics of thermal decomposition of ferrocene carboxaldehyde in an oxidative atmosphere was reported.
Presently, in order to find the effect of a co-precursor on the nature of solid-state thermal decomposition of ferrocene carboxaldehyde (C5H4CHO)Fe(C5H5) (Scheme 1), as well as on the reaction product, we studied the thermal decomposition of this precursor in the presence of oxalic acid dihydrate as a co-precursor in O2 atmosphere. Remarkably, oxalic acid dihydrate decomposes to CO2, CO, and H2O with complete mass loss when heated beyond ∼373 K,35 whereas ferrocene carboxaldehyde starts to decompose at ∼380 K.25 Here, the TG profiles with multiple heating rates of a 1
:
1 mixture of ferrocene carboxaldehyde and oxalic acid dihydrate were obtained in O2 atmosphere, and the decomposed product was identified as a hematite nanomaterial. The multi-step differential thermogravimetry (DTG) profiles were resolved into individual reaction steps by peak-deconvolution technique, and the reaction steps responsible for producing hematite were identified. Model-free non-isothermal methodology was utilized to evaluate the reaction kinetic parameters for each step of the reaction, following the recommendation of the International Confederation of Thermal Analysis and Calorimetry Kinetics Committee (ICTAC).41,42 Relevant thermodynamic parameters (ΔS, ΔH and ΔG) were found, and their temperature dependence was studied. A substantial resemblance between the experimentally obtained thermal decomposition profiles and those constructed using the estimated kinetic parameters established the accuracy of the modelling approach applied herein and the methodology adopted. The present work highlights how a co-precursor can affect the decomposition reaction of a precursor to lead to materials of nanometric size. Such results might be extremely useful for researching on the thermal synthesis methodology of nanomaterials of different size, starting with a precursor but by varying the co-precursor.
![image file: d3ra07045j-s1.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-s1.gif) |
| Scheme 1 A schematic diagram of ferrocene carboxaldehyde (C5H4CHO)Fe(C5H5). | |
2. Experimental and analytical methods
2.1 Material and methods
Ferrocene carboxaldehyde (hereafter, FcCHO) and oxalic acid dihydrate, (COOH)2·2H2O (hereafter, OxA) were procured from Sigma, and were used as precursor and co-precursor, respectively. Both compounds were mixed in a 1
:
1 weight ratio, and then finely powdered using a mortar pastle. Thermal decomposition of the mixture (hereafter, FO11) was performed in a thermogravimetric analyser, TGA (STA 449 F3 Jupiter, Netzsch, Germany) in UHP (99.999%) O2 gas environment (flow rate: 50 ml min−1) with UHP (99.999%) N2 gas as the protective gas (flow rate: 20 ml min−1). Heating rate-dependent thermogravimetry (TG) profiles were obtained when the sample of ∼5 mg mass was scanned with an alumina-made sample and reference crucibles.
A powder X-ray diffraction (XRD) study of the decomposed material was performed with a diffractometer D8 Advance of Brucker for the Cu-Kα radiation source within the 2θ values ranging from 20° to 80°. MATCH (free) software was used for the powder XRD data analysis. Field Effect Scanning Electron Microscope (FE-SEM) images were collected using Gemini SEM 450 of Zeiss, where the Energy Dispersive X-ray Analysis (EDX) study was made with the Model-Z2 of AMETEK. ImageJ software was used for the analysis of the SEM images. For the analysis of the observed thermogravimetric data, the ORIGIN software was used, whereas all calculations for estimation of the reaction kinetic and thermodynamic parameters were done with the help of a program compiled in MATLAB.
2.2.1 Theoretical background of kinetic analysis. In a kinetic study of the thermal solid-state reaction, the rate equation of a single-step thermal decomposition is assumed as,25 |
![image file: d3ra07045j-t1.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t1.gif) | (1) |
where α is the extent of conversion (
, mi is the initial mass, mf is the final mass, mT is the mass at temperature T), k(T) = Arrhenius rate function, f(α) = conversion function, Aα and Eα are the conversion-dependent reaction rate and activation energy, respectively, while R is the universal gas constant. For a non-isothermal reaction process with constant linear heating rate β
, eqn (1) can be written as, |
![image file: d3ra07045j-t4.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t4.gif) | (2) |
which on rearrangement, reduces to:
|
![image file: d3ra07045j-t5.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t5.gif) | (3) |
Integrating eqn (3), one obtains the expression for the mechanism function:
|
![image file: d3ra07045j-t6.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t6.gif) | (4) |
where
![image file: d3ra07045j-t7.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t7.gif)
,
![image file: d3ra07045j-t8.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t8.gif)
and integral function
![image file: d3ra07045j-t9.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t9.gif)
. From
eqn (4), it is clearly visible that determination of
g(
α) is dependent on the approximation of
p(
x). After using the best approximation of
p(
x) for a suitable range of
x, it is possible to calculate the kinetic parameters –
Eα,
g(
α) and
Aα.
2.2.2 Kinetic parameters. The minimum energy required for the formation of the new bonds in a new complex is known as the activation energy. It is possible to solve the above stated integral p(x) using modern integral methods utilizing numerical integration and best approximation of p(x) suitable for a range of x,28–30 and hence to evaluate the activation energy Eα. There are several well-known integral methods, like Flynn–Wall–Ozawa (FWO),43,44 Kissinger–Akahira–Sunose (KAS),45,46 Tang,47 Starink,29 Akbi–Mekki–Rafai (AMR),48 and Vyazovkin,28 defined as:FWO:43,44
|
![image file: d3ra07045j-t10.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t10.gif) | (5a) |
KAS:45,46
|
![image file: d3ra07045j-t11.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t11.gif) | (5b) |
Tang:47
|
![image file: d3ra07045j-t12.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t12.gif) | (5c) |
Starink:29
|
![image file: d3ra07045j-t13.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t13.gif) | (5d) |
AMR:48
|
![image file: d3ra07045j-t14.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t14.gif) | (5e) |
Vyazovkin:28
|
![image file: d3ra07045j-t15.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t15.gif) | (5f) |
The slope of the linear fit obtained from ln(β/Tn) vs. 1/T plots following eqn (5a–d) leads to the estimation of the Eα value. The Eα value can also be derived from eqn (5e) utilizing the method of iteration.48 Use of eqn (5f) leads to the estimation of the relatively more accurate value of the activation energy from the minimum value of the objective function Ω vs. Eα plots.28 Herein, the equations eqn (5a–f) are applied to evaluate and compare the activation energy values of thermal decomposition.
The mechanism function describes the reaction path followed by the reactants during the product formation, while the reaction rate indicates the frequency of collisions, and hence gives an idea about the reactivity. Different reaction mechanism functions49 g(α) are presented in Table S1,† where models like the chemical reaction, random nucleation, phase boundary reaction and diffusion, and subsequent growth are some of the frequently utilized reaction mechanisms. Usually, the master plot method50 is used to determine the most probable mechanism function. A master plot is the reference theoretical curve based on the kinetic model function. This method is based on the comparison of theoretical master plots with the experimental master plots. According to eqn (4) with α = 0.5 as a reference point, the following equation is obtained:
|
![image file: d3ra07045j-t16.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t16.gif) | (6) |
with
Aαa = the average value of reaction rate
Aα values,
Eαa = the average value of activation energy
Eα values,
p(0.5) =
Eαa/
RT0.5, and
T0.5 is the temperature corresponding to 50% conversion. Dividing
eqn (4) by
eqn (6), one obtains
|
![image file: d3ra07045j-t17.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t17.gif) | (7) |
The plot of [g(α)/g(0.5)] as a function of α originates the theoretical master plot for a range of g(α) functions. The experimental master plots are obtained using a suitable approximation of p(x). When the most probable model is used, eqn (7) shows that for a particular α, the experimentally obtained value of p(x)/p(0.5) and the theoretically calculated value of g(α)/g(0.5) are equal. There are many approximations of p(x). Among those, the Senum–Young approximation51 is a good approximation of p(x). Using this approximation, p(x) is written as,
|
![image file: d3ra07045j-t18.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t18.gif) | (8) |
After determining Eα using eqn (5a–f) and g(α) using eqn (7), the reaction rate Aα is estimated as a function of α for different heating rates using the equation:48
|
![image file: d3ra07045j-t19.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t19.gif) | (9) |
2.2.3 Thermodynamic parameters. The thermodynamic parameters – changes in entropy (ΔS), enthalpy (ΔH) and Gibbs free energy (ΔG) for the activated complex formation during the thermal decomposition can be estimated using the following relations:21 |
ΔS = R ln(Aαh/eχkBTp)
| (10a) |
where e = 2.7183 (base of natural logarithm); χ is the transition factor; kB is the Boltzmann constant; h is Planck's constant, and Tp is the peak temperature obtainable from the
vs. T curves.
2.3 Peak deconvolution
Eqn (1) represents a single-step decomposition reaction. However, most of the time, thermal decomposition of the iron-based organic solid precursors is a multistep complex reaction process, i.e., they follow multiple consecutive overlapping reaction steps.21,23,25,52 For a multistep complex process, the overall rate equation is the combination of multiple single-step reactions.53 To separate the individual steps, the peak deconvolution method is a suitable process. For multiple overlapping peaks, the overall
can be written as |
![image file: d3ra07045j-t22.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t22.gif) | (11) |
where
and
. Here, n represents the total number of steps, and ci denotes the contribution of the ith process. Ai, Eα,i and fi(αi) are the Arrhenius parameters of process i, and α represents the extent of reaction of the overall process. For a multistep process, the peak deconvolution method is a good approach to separate the overlapping peaks. Perejón et al.54 used Gaussian, Lorentzian, Weibull and Fraser Suzuki functions to fit
. Considering the asymmetry, the Weibull and Fraser Suzuki functions21 were fitted to the experimental
vs. T curves. However, in the present study, the fitted curves for the Fraser Suzuki function produced better accuracy than the Weibull function (Fig. S1†). Thus, using the Fraser Suzuki function,21,23,25 eqn (11) is written as |
![image file: d3ra07045j-t27.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t27.gif) | (12) |
where
is the conversion rate of the overall process, and hi, pi, si and wi are the parameters corresponding to the peak amplitude, position, asymmetry and half-width of the ith peak, respectively. After setting the initial values of the fit parameters, an optimization operation is done to minimize F defined as |
![image file: d3ra07045j-t29.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t29.gif) | (13) |
where M is the total number of data points under a constant heating rate. The lower value of F provides the better accuracy in fitting. On successful deconvolution, the rate equation for each individual and independent steps of thermal decomposition can be represented by eqn (12) with the derived fit parameters.
3. Results and discussion
3.1 Analysis of thermogravimetry results
Fig. 1a compares the observed TG profile of FO11 in an O2 atmosphere with those obtained for its components, OxA and FcCHO. The TG profile of OxA obtained under 5 K min−1 heating rate shows that it starts to decompose at ∼320 K and decomposed completely at ∼474 K following a three-step process, while the TG profile of FcCHO obtained under 3 K min−1 heating rate shows that decomposition started at ∼380 K and was completed at ∼678 K with 35% of residual mass, seemingly following a three-step process. On the other hand, thermal decomposition of FO11 starts at ∼305 K and is completed at ∼625 K with 20% residual mass, exhibiting an apparent four-step process. The mass loss associated with each decomposition step in the TG profiles is due to the release of volatile material(s), owing to the progress of the solid-state decomposition reaction. With the progress of reaction, the reaction steps occur at different temperatures, and the reaction is completed when a thermally stable compound is obtained. Comparing the TG profiles, as FcCHO is seen to be stable up to ∼380 K from the TG profile, the initial mass loss of FO11 observed in the 305–380 K range should be due to OxA decomposition. After OxA decomposes completely (at ∼474 K), any mass loss in FO11 should be due to decomposition of FcCHO. Hence, as FO11 is a physical mixture of FcCHO and OxA, one can expect that the TG profile of FO11 should be similar to that of FcCHO above ∼474 K. However, this is not the result observed, as reflected in Fig. 1a. The nature of TG profile of FO11 is quite different from that of FcCHO (shifted towards lower temperature compared to FcCHO) and the decomposition of FO11 is completed at ∼625 K, which is about 50 K lower than that of FcCHO.
![image file: d3ra07045j-f1.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f1.gif) |
| Fig. 1 A comparison of the TG (a) and DTG (b) profiles observed for FO11, FcCHO, and OxA. | |
Fig. 1b shows the differential thermogravimetry (DTG) profiles obtained from the TG data for OxA, FcCHO and FO11 shown in Fig. 1a, where the existence of multiple peaks/humps in a single profile indicates the occurrence of multi-step thermal decomposition reactions and the positions of the peaks/humps are the temperatures where maximum decomposition takes place in a particular reaction step. Here, OxA shows two distinct peaks at 375 K and 467 K, along with a hump at ∼404 K. Thermal decomposition of FcCHO shows three peaks, along with small humps at 433 K, 465 K (hump), 482 K (hump), 493 K (hump), 601 K and 644 K, in agreement with a previous report.25 Apparently, the DTG profile of FO11 shows three large peaks at 364 K, 433 K and 529 K, along with shoulder peaks/humps at 408 K, 465 K and 564 K. Thus, a comparison of the DTG profiles of FO11 and its constituents shows that there is a clear shift in the observed DTG peak positions when FcCHO is decomposed in the presence of OxA. Hence, the observations made from the analysis of the TG and DTG profiles of OxA, FcCHO and FO11 (shown in Fig. 1) undoubtedly establish that the presence of OxA in the solid reaction atmosphere significantly influences the thermal decomposition of FcCHO.
Fig. 1 showed that thermal decomposition of FO11 in an O2 atmosphere is a complex multistep process. Hence, the peak deconvolution method following the Fraser Suzuki function, as discussed earlier, was adopted to separate the overlapping DTG peaks of FO11, FcCHO and OxA. The fit parameters are available as supplementary material (Table S2†). Fig. 2 compares the deconvoluted
vs. T plots for FO11, FcCHO and OxA with the experimental ones for a single heating rate, where the agreement is very good. The peak deconvolution method shows that thermal decomposition of FO11 takes place in six different steps, where the steps are quite overlapped, while OxA and FcCHO decompose in three and six different steps, respectively, in agreement with earlier reports.23,25 Table 1 compares the temperature corresponding to the maximum decomposition (Tp) for each reaction step and the temperature range of steps (ΔT) for the decomposition of FO11, FcCHO and OxA in an O2 atmosphere. From Fig. 2 and Table 1, when the reaction steps of individual FcCHO and OxA are compared with those of FO11, it is imminent that peak-1 and peak-2 observed for FO11 are due to the evaporation of adsorbed water molecules, followed by partial decomposition of OxA, and peak-3 is due to the initial decomposition of FcCHO. OxA starts to decompose much before FcCHO, and Step 1 of OxA completes at ∼387 K, while FcCHO starts its Step 1 of decomposition at around the same temperature. Peak-4 observed for FO11 is low in intensity. Notably, peak-3 of OxA and peak-2 of FcCHO appear at about the same temperature at lower heating rates. Eventually, peak-2, 3, 4 of lone FcCHO may be present in the same temperature range. However, considering the level of error, any further attempt to deconvolute the low intensity curve (i.e., peak-4 for FO11) was avoided. Heating rate-dependent
vs. T plots of OxA (see Fig. S2†) show that the intensity of peak-2 increases with increasing heating rate, and peak-3 is sustained at higher heating rates when lone OxA is heated in an O2 atmosphere. This discards the possible merger of peak-3 of OxA with peak-2 of FcCHO as peak-4 in FO11. The disappearance of peak-3 of OxA in FO11 is inquisitive. On the other hand, with increasing heating rate, peak-2 of lone FcCHO decreases faster in intensity and shifts towards higher temperatures.25 It will be shown in the later discussions that with increasing heating rate, peak-4 observed for FO11 gradually decreased in intensity and shifted towards higher temperature. The positions of Tp for Step 5 and -6 for FO11 are significantly different than those of lone FcCHO. Thus, it may be presumed that in FO11 under O2 atmosphere, after the Step 2 decomposition of OxA, the gaseous products that are formed react with those formed by Step 1 decomposition of FcCHO with increasing temperature, leading to Step 4, -5 and -6 of reactions in FO11, as well as the residual product. This causes the different TG and DTG profiles for thermal decomposition of FO11 in comparison to that of FcCHO.
![image file: d3ra07045j-f2.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f2.gif) |
| Fig. 2 The deconvoluted vs. T profiles of FO11, FcCHO, and OxA. | |
Table 1 Values of the temperature corresponding to maximum decomposition (Tp) for each reaction step and the temperature range of steps (ΔT) for decomposition of FO11, ferrocene carboxaldehyde, and oxalic acid in O2 atmosphere
Material |
Parameter |
Step 1 |
Step 2 |
Step 3 |
Step 4 |
Step 5 |
Step 6 |
Heating rate. |
FO11 @ 3a K min−1 |
Tp (K) |
364 K |
408 K |
433 K |
465 K |
530 K |
564 K |
ΔT (K) |
300–390 K |
332–432 K |
382–468 K |
439–515 K |
440–568 K |
515–625 K |
FcCHO @ 3a K min−1 |
Tp (K) |
432 K |
465 K |
482 K |
493 K |
601 K |
643 K |
ΔT (K) |
380–440 K |
446–479 K |
464–494 K |
472–512 K |
452–641 K |
585–678 K |
OxA @ 5a K min−1 |
Tp (K) |
375 K |
402 K |
469 K |
|
|
|
ΔT (K) |
304–387 K |
378–414 K |
380–474 K |
|
|
|
3.2 Identification of the decomposed material
In order to identify the residual product obtained on thermal decomposition of FO11 in O2 atmosphere, a room temperature powder XRD pattern was obtained in the 2θ range of 20°–80°. Fig. 3 shows the observed powder XRD pattern along with the Rietveld refinement results, where the dot represents the observed data, the solid line corresponds to the calculated pattern, and the vertical lines indicate the Bragg diffraction positions. Clear distinct Bragg peaks are observed at 24.20°, 33.28°, 35.69°, 40.89°, 49.61°, 54.17°, 62.52° and 63.99°. Close agreement between the observed and calculated powder XRD patterns recognized that the thermally decomposed material obtained from FO11 in the O2 atmosphere is pure hematite (α-Fe2O3) with a hexagonal crystal structure with a = b = 5.0428 Å, c = 13.7326 Å and α = γ = 90° and β = 120° (COD No: 96-901-5066).55 The sharp and intense peaks of α-Fe2O3 imply the absence of impurities and good crystalline nature of the sample.
![image file: d3ra07045j-f3.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f3.gif) |
| Fig. 3 The powder XRD profile of the decomposed material obtained on thermal decomposition of FO11 in O2 atmosphere. The lower panel shows the difference between the observed and fitted XRD profiles. | |
The average particle size of the hematite material derived using the Scherrer method56 and Wagner–Halder method57 is 12.51 ± 0.02 nm and 10.72 ± 0.03 nm, respectively. So, the average particle size values estimated using both methods are consistent with each other and lie in the nanometric range. Thus, it is recognized that on thermal decomposition of ferrocene carboxaldehyde in the presence of oxalic acid dihydrate in an oxidative atmosphere completed at 625 K, pure hematite nanoparticles can be obtained. At this point, it may be recalled that the thermal decomposition of lone ferrocene carboxaldehyde completed at ∼678 K in oxygen atmosphere produced dot-shaped hematite nanomaterials of ∼5 nm av. diameter.37
FE-SEM image was used to analyse the morphology of the prepared hematite sample (Fig. 4). The images revealed that the hematite material has rod-like forms, interestingly different from the spherical nanoparticles obtained when bare FcCHO was decomposed under the same oxidative atmosphere. Some nanorods were randomly selected and then processed with ImageJ software to obtain the size distribution of the resulting particles (shown as insert of Fig. 4). The average length and diameter of the rods obtained from the histogram are 148 ± 6 nm and 83 ± 3 nm, respectively. Next, EDX analysis was used to measure the proportion of elemental species contributing to the synthesized hematite nanorods. The observed EDX spectrum is shown in Fig. S3,† which verified the presence of 57.0% of iron and 43.0% of oxygen (atomic weight percentage) species, confirming that the thermally synthesized hematite is made up of only iron and oxygen with no impurities.
![image file: d3ra07045j-f4.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f4.gif) |
| Fig. 4 FE-SEM image of the decomposed material of FO11. The inserts show the magnified image and histograms. | |
3.3 Estimation of kinetic parameters
Here, we try to understand the reaction processes involved by estimating the reaction kinetic parameters. The TG profiles of FO11 under O2 atmosphere at different heating rates were obtained, where the DTG (see Fig. S4†) showed multiple overlapping peaks/humps. The peak-deconvolution method, as discussed earlier, is applied to separate the individual peaks/humps, and the results obtained are shown in Fig. 5. The fit parameters are available as supplementary material in Table S3.† In this table, the different values of w for different heating rates indicate the different reactivities for the separated reactions. The different values of h with each heating rate signify the different intensities. A positive value of s signifies the reaction is right symmetric, while a negative value of s signifies the reaction is left symmetric.58 Fig. 5 compares the deconvoluted
vs. T plots for FO11 with the experimental ones obtained from the TG data at different heating rates, where the agreement between the calculated (cumulative) and the experimentally observed
vs. T plots is very good. The observed shift in Tp towards the higher temperature with increasing heating rate may be due to the lowered heat transfer rate between the furnace of TGA and the sample. From the estimated ΔT and Tp values (Table S3†), thermal decomposition of FO11 in O2 occurs in six overlapped steps in the 300–419 K, 332–465 K, 382–488 K, 439–545 K, 440–586 K and 515–635 K temperature ranges having Tp values (averaged over heating rates) equal to ∼371 K, ∼420 K, ∼447 K, ∼488 K, ∼545 K and ∼573 K, respectively. Thus, as stated above, the first three peaks (i.e., first three steps) that appeared for FO11 are due to decomposition of FcCHO and OxA in the O2 atmosphere, whereas the last three peaks (i.e., last three steps) are for the thermal reactions leading to hematite. Considering the aim of the present study, we will focus only on the last three reaction steps.
![image file: d3ra07045j-f5.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f5.gif) |
| Fig. 5 The results of the peak deconvolution analysis ( vs. T plots) for each separate step of the thermal decomposition of FO11 at different heating rates under O2 atmosphere. | |
For kinetics analysis, the extent of conversion α for each step of decomposition of FO11 was calculated as a function of temperature using the relation
,25 where paTemp,T and paTotal are the area of the dα/dT vs. T curve up to a temperature T and the area under the individual dα/dT vs. T curve, respectively for a particular heating rate. The α vs. T plots thus obtained for Step 4, -5 and -6 of FO11 are shown in Fig. 6.
![image file: d3ra07045j-f6.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f6.gif) |
| Fig. 6 The α–T curves for Steps 4, 5, and 6 for the thermal decomposition of FO11. | |
The Eα values for 0.1 ≤ α ≤ 0.9 for the thermal decomposition of FO11 for Step 4, -5 and -6 are estimated by the FWO, KAS, Tang, Starink, AMR and Vyazovkin methods using eqn (5a–f). The ln(β/Tn) vs. 1/T curves obtained for Step 4, -5 and -6 are compared in Fig. S5,† while Fig. S6† shows the Ω(Eα) curves according to eqn (5f) for the same steps. Table S4† presents the results of linear fitting of the experimental data with eqn (5a–f). The Eα values thus estimated using eqn (5a–f) are presented in Table S5.† Fig. 7 shows the variations of these Eα values for thermal decomposition of FO11 for Step 4, Step 5 and Step 6 with the extent of conversion α in an O2 atmosphere. For step 4, the Eα values estimated using the FWO method lie between 57 and 71 kJ mol−1, whereas for the rest of the methods, the Eα values lie between 52 and 67 kJ mol−1. Meanwhile, for Step 5 and Step 6, the Eα values estimated following different methods are quite close and lie within 82–137 kJ mol−1 and 128–207 kJ mol−1 for Step 5 and -6, respectively. As the Eα values for each reaction step using eqn (5a–f) are quite close to each other and particularly, those estimated using eqn (5b–f) are the same, the Eα values estimated using eqn (5f) will be preferentially used for further calculations of the kinetic parameters.
![image file: d3ra07045j-f7.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f7.gif) |
| Fig. 7 Eα vs. α curves for Step 4, 5, and 6 of FO11 using different methods eqn (5a–f). Solid lines are the fits to the estimated values. | |
The nature of dependence of Eα values on α obtained following different equations are identical for each individual step. The Eα values significantly depend on the extent of conversion in each step. However, the nature of dependence of Eα on α is significantly different from step to step. The Eα values are estimated for some distinct α. Since Eα is continuous, we interpolate Eα values for all α by fitting with suitable functions. In Step 4, the Eα value initially decreases up to α = 0.5, and then increases following a polynomial equation of second order. In Step 5, the Eα value increases with an increase in α, following an exponential growth equation. Meanwhile, in Step 6, the Eα value increases with an increase in α, following a polynomial equation of second order. The functions are available in Table S6.† Since the activation energy is the energy barrier to overcome for a reaction process to occur, the higher values of Eα with the increase of reaction step (i.e., temperature) indicate that more energy is required for the formation of the corresponding products with increasing temperature. In other words, it can be stated that the higher activation energy indicates the high thermal stability of the activated complex. An increasing trend of the activation energy with increasing temperature indicates competing, independent or consecutive reactions.59 Generally, this type of reaction occurs during a solid to (solid + gas) transformation. Usually, the gas–solid reaction is a heterogenous reaction existing with multiple parallel reactions. If the dominant reaction path or the reaction control steps change, the activation energy may vary significantly.60 The decreasing trend of the activation energy with increasing temperature indicates an endothermic reversible reaction, followed by an irreversible reaction.61
The most probable reaction mechanism functions g(α) for the present solid-state reactions are obtained following the master plot method described earlier using all 35 mechanisms (shown in Table S1†). Fig. 8 shows the comparison of the theoretical and experimental master plots for FO11 in Step 4, 5 and 6. From Fig. 8, it is clear that none of the experimental plots exactly coincide with the theoretical master plot. Hence, to determine the most probable reaction mechanism function, the standard deviation between the theoretical master data and experimental data50,52 was calculated, and then the function for which the Σδ became minimum was chosen as the most probable reaction function. The expression for Σδ52 is stated as
|
![image file: d3ra07045j-t37.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t37.gif) | (14) |
where
m = number of data points and
n = number of heating rates, and
δ = averaged square value of the deviation between
p(
x)/
p(
x0.5) calculated from the experimental data and
gk(
αi)/
gk(0.5) calculated for any theoretical model.
![image file: d3ra07045j-f8.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f8.gif) |
| Fig. 8 The master plot curves of FO11 for Step 4, 5, and 6. | |
The most probable mechanism functions for thermal decomposition of FO11 (described in Table S7†) are the sigmoidal reaction (A2), sigmoidal reaction (A3/2) and chemical reaction (F3/2) for Step 4, -5 and -6, respectively. So, the thermal decomposition of FO11 in an O2 atmosphere follows different reaction mechanisms for different steps. The chemical reaction represents the bond redistribution between the reactants and products. The acceleratory rate equations indicate nucleation without considering any restriction of nuclear growth. For these kinds of reactions, the rate of the reaction increases with the conversion. Sigmoidal equations represent the random nucleation and growth of new nuclei. In the nucleation process, two types of chemical change may be involved – (i) the chemical transformation of one or more constituents of the reactant into the constituents of the product, and (ii) recrystallization of the reacted material to the lattice structure of the product phase. The effect of nucleation causes the active reactant/product interface generation. The reaction occurs in these interfaces during its subsequent advance into an unchanged reactant as the nucleus grows. Growth is thus the maintained reaction within the interface between the reactant and product phases. For such reactions, the nucleation first occurs at local nucleation sites, and then the nucleation sites are eliminated through the growth of new nuclei.62 Hence, the rate of the reaction first increases and then decreases with the conversion of the reaction.
The reaction rate values of Step 4, -5 and -6 for thermal decomposition of FO11 as a function of α are calculated from the estimated Eα and g(α) using eqn (9), and are presented in Table S8.† Fig. 9 shows the variations of lnAα with α for Step 4, -5 and -6 of FO11. The α-dependent Aα values of for Step 4, -5 and -6 lie within the range 8.4 × 104–2.7 × 106 min−1, 5.6 × 106–2.7 × 1012 min−1 and 9.9 × 109–1.8 × 1018 min−1, respectively. Reaction rate Aα represents the frequency of the collision in the activated complex. For surface-independent reactions, the higher values of Aα indicate a ‘loose’ complex (>109 s−1) indicate a highly reactive system, while a low value of Aα indicates a ‘tight’ complex and obviously a less reactive system.63 The high values of the reaction rate for Step 6 indicates the high reactivity of the activated complex. The increasing trend of Aα with increasing temperature indicates that the reactivity of the activated complex increases with the extent of conversion. The nature of variation of lnAα with α for each step is identical with the corresponding Eα vs. α plots.
![image file: d3ra07045j-f9.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f9.gif) |
| Fig. 9 The variations of the reaction rate with the extent of conversion for Step 4, 5, and 6 of FO11. Solid lines are the fits to the estimated values. | |
3.4 Reliability of the kinetic models used
When the estimated kinetic parameters are combined, the observed step-wise solid-state thermal decomposition of FO11 in an O2 atmosphere can be summarised as follows:
Step 4:
|
![image file: d3ra07045j-t38.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t38.gif) | (15a) |
Step 5:
|
![image file: d3ra07045j-t39.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t39.gif) | (15b) |
Step 6:
|
![image file: d3ra07045j-t40.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t40.gif) | (15c) |
To verify the reliability of the used kinetic models as well as the validity of the estimated reaction kinetic parameters for the thermal decomposition of FO11 leading to hematite nanoparticles, the conversion curves (α vs. T plots) under different heating rates are made following eqn (14), and are then compared with those obtained from experimental data. Fig. 10 shows such constructed α vs. T plots for different heating rates with the corresponding experimentally obtained ones. There is a reasonable agreement between these α vs. T plots, which specifies the precision and utility of the methodology implemented for the kinetic analysis and the choice of appropriate models for different steps of reaction in the present study.
![image file: d3ra07045j-f10.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f10.gif) |
| Fig. 10 A comparison of α vs. T plots for Step 4, 5, and 6 obtained using the kinetic parameters and experimentally observed for thermal decomposition of FO11 at different heating rates. | |
3.5 Estimation of thermodynamic parameters
The values of the thermodynamic parameters ΔS, ΔH, and ΔG for Step 4, 5, and 6 of thermal decomposition of FO11 in O2 atmosphere are calculated using eqn (10a–c) using the estimated kinetic parameters. The estimated values are found to be heating rate-independent. Thus, their values averaged for different heating rates were considered for discussion. The values of ΔS, ΔH and ΔG are presented in Table S9,† and are plotted in Fig. 11 as a function of α.
![image file: d3ra07045j-f11.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-f11.gif) |
| Fig. 11 Variations of ΔS (a), ΔH (b) and ΔG (c) with the extent of conversion for Step 4, 5 and 6 of FO11. | |
From Fig. 11a, for Step 4, the ΔS values decrease with increasing α up to α = 0.5, followed by an increase up to α = 0.9 and lie within −138.5 to −134 kJ mol−1 K−1. For Step 5, the ΔS values linearly increase with increasing α and lie within −129.1 to −20.2 kJ mol−1 K−1. For Step 6, the ΔS values increase with increasing α up to α = 0.8, followed by a decrease up to α = 0.9 and lie within the range of −67.3 to –90.7 kJ mol−1 K−1. ΔS implies the change in the entropy values while going from the intermediate reactant state to the transition state, and provides a measure of the degree of randomness in the system. A negative value of ΔS corresponds to a low degree of freedom in the activated complex and vice versa. Consequently, these can be referred to as the “slow” and “fast” stages of the reaction, respectively. A lower entropy value implies that the substance is closer to its equilibrium state. Consequently, it will be less reactive and a longer time is required to form the activated complex. Equally, a higher entropy value implies that the system is far-off from its equilibrium. Thus, it is highly reactive and will not take a longer time to form an activated complex.
The energy required by a system for the reaction to occur is the enthalpy. The difference in the potential energy between the activated complex and the precursor is denoted as the change in enthalpy ΔH. Fig. 11b reveals that the estimated ΔH values are α-dependent, and that the dependence of ΔH on α follows analogous form to that of ΔS and activation energy. The average values of ΔH (52.86, 101.34, 175.0 kJ mol−1) for Step 4, -5 and -6 are lower than the average activation energy values (56.92, 105.87 and 179.76 kJ mol−1) of the corresponding steps, which indicate the formation of the activated complexes. The small difference between the activation energy and change in enthalpy suggests a lower energy barrier between the products produced and the activated complex.63 Hence, there is a smaller required amount of thermal energy for the occurrence of the reactions. The variation of ΔG with α for the above-stated reactions is shown in Fig. 11c. ΔG is the difference in the free energies between the activated complex and the precursor. The sign of ΔG points to the spontaneity of the reactions. Thus, ΔG is the key factor that directs a chemical reaction. Presently, for Step 4, -5 and -6, the ΔG values lie within 128.5 to 128.7 kJ mol−1, 144.5 to 140.9 kJ mol−1 and 159.6 to 152.6 kJ mol−1, respectively. As the ΔG values are higher than the ΔH values for each step, the energy supplied to the system is partially used for the creation of products. For all of the reaction steps, the positive ΔG values specify that the present thermal reactions are non-spontaneous and require energy to continue.
3.6 Role of the co-precursor
The thermal decomposition of bare ferrocene carboxaldehyde in oxygen atmosphere is a multi-step process completed at ∼678 K, leaving pure hematite as the residual mass. Conversely, the thermal decomposition of the same compound in nitrogen atmosphere completed at ∼560 K yields a mixture of hematite, cementite and magnetite.37 Thus, the significant effect of the oxygen atmosphere on the nature of the thermal decomposition of ferrocene carboxaldehyde and also on the resultant products are established. Herein, a mixture of ferrocene carboxaldehyde and oxalic acid is thermally decomposed in an oxygen atmosphere, and produced hematite at 625 K, as confirmed from the XRD study, which also recognized the effect of co-precursor on the reaction process of ferrocene carboxaldehyde in oxygen atmosphere. Ferrocene carboxaldehyde is a substituted derivative of ferrocene (C5H5)2Fe. It has two parallel cyclopentadienyl rings with an iron ion in the centre between the rings. It is reported that ferrocene sublimates at ∼448 K. Above ∼773 K, gaseous ferrocene decomposes spontaneously, giving rise to metallic iron along with some volatile gases [(C5H5)2Fe → Fe + H2 + CH4 + C5H6 + …].64 At such high temperatures, solid or liquid-like iron particles exist in the reaction medium. However, when ferrocene was heated at ∼453 K together with oxalic acid, hematite was formed35 [2(C5H5)2Fe + (COOH)2·2H2O → Fe2O3 (s) + CO2(g) + CO(g) + 13H2(g) +22C], where oxalic acid first decomposed to CO2 and CO and H2O. The CO2 and CO are deoxidized to oxygen and carbon by metallic iron, and it then reacted with the available oxygen, giving rise to hematite.35 Presently, considering the strength of the bonds present in ferrocene carboxaldehyde, it is obvious that the –CH
O bond will break prior to the breakage of the Fe–C bonds, and the resulting metallic iron will react with the oxygen available to give rise to hematite.37
There are many predictions on the thermal decomposition of oxalic acid. According to Lapidus et al.,65 the thermal decomposition of oxalic acid, at temperatures up to 430 K, yielded only CO2 and HCOOH. Multiphoton infrared decomposition study of oxalic acid by Yamamoto and Back66 pointed out that though two primary processes yielding CO2 + HCOOH and CO2 + CO + H2O may occur, but CO and H2O are more vital products than HCOOH considering the energy required. From the ab initio calculations over the temperature range of 300–1300 K, Kakumoto et al.67 considered the decomposition of oxalic acid to occur through several reaction channels as follows:
|
(COOH)2 → HCOOH + CO2
| (a) |
|
(COOH)2 → CO2 + CO + H2O
| (b) |
|
![image file: d3ra07045j-t41.tif](/image/article/2023/RA/d3ra07045j/d3ra07045j-t41.gif) | (e) |
where it was predicted that the lowest energy path was (COOH)
2 → CO
2 + CO + H
2O (in bold), as the dihydroxycarbene [C(OH)
2] produced on decomposition of oxalic acid in path (e) mostly decomposes to CO + H
2O. Kakumoto
et al.67 performed thermal decomposition of oxalic acid vapor in an inert atmosphere in 850–1300 K, which led to CO
2, CO, and H
2O as the major products. This is consistent with their proposed reaction path obtained from the
ab initio calculations, and also the infrared multiphoton study by Yamamoto and Back.
66 According to Higgins
et al.,
68 although rapid formation of CO
2 and HCOOH in gas phase oxalic acid thermal decomposition is possible, with increasing temperature, the production of CO
2, CO, and H
2O becomes more significant.
Presently, the 1st peak that appeared for the decomposition of the bare oxalic acid dihydrate arises due to evaporation of absorbed water present, and the dehydrated oxalic acid further decomposes in two-steps. It should be noted that during the thermal decomposition of oxalic acid dihydrate in the presence of ferrocene carboxaldehyde, the third step of oxalic acid thermal decomposition did not appear. Thus, Step 2 of the oxalic acid decomposition process may be recognized as the formation of C(OH)2 and CO2, following Kakumoto et al.67 The [C(OH)2] produced on further heating decomposes to CO + H2O. Notably, thermal decomposition of ferrocene carboxaldehyde is also initiated at around the same temperature. Looking at the complexity of the reaction, it may plausibly be inferred that the metallic iron produced (as discussed above) will react with the oxygen available in the reaction environment to form hematite. Additionally, the metallic iron may deoxidize CO2 into carbon and oxygen,34 and react with this additional oxygen as an additional reaction, forming hematite. As no traceable carbon or iron carbide is detected in the residual product, it is likely that the carbon produced as an intermediate might have been oxidized to CO2/CO.69
Table 2 compares the temperature range, activation energy, reaction mechanism functions and reaction rate observed for thermal decomposition of FO11 and FcCHO in an O2 atmosphere. It is clear from this table that the use of oxalic acid as a co-precursor has significantly reduced the temperature at which the step-wise thermal decomposition reaction occurred in FcCHO, modified the activation energy of the thermal decomposition reaction of FcCHO particularly in Step 5, and enhanced the reaction rates by several orders in Step 5 and -6, but did not meaningfully change the reaction mechanisms in Step 4 and -5. When the changes involved in the thermodynamic parameters are compared, no significant change in them are observed for Step 4, whereas a huge change in those parameters are observed for Step 5 in comparison to those in Step 6. Perhaps in FO11, the reactivity of the reactants after Step 3 is higher, which led to the early completion of the reaction for FO11 with respect to bare FcCHO on the temperature scale. Thus, the oxalic acid co-precursor acts as a reaction enhancer70 in the thermal decomposition of ferrocene carboxaldehyde. It may be recalled that ferrocene sublimates at ∼773 K. When ferrocene is heated together with oxalic acid, it produces hematite at much lower temperature (∼453 K)35 as is seen in the present case.
Table 2 A comparison of the kinetic parameters (Eα, g(α), and Aα) and thermodynamic parameters (ΔS, ΔH, and ΔG) for Step 4, 5, and 6 observed for the thermal decomposition of FO11 and ferrocene carboxaldehyde (FcCHO)
Parameters |
Material |
Step 4 |
Step 5 |
Step 6 |
25. For symbols used, refer to Table S1. |
Steps (K) |
FO11 |
439–515 |
440–568 |
515–625 |
FcCHOa |
472–512 |
452–641 |
585–678 |
Eα (kJ mol−1) |
FO11 |
52–67 |
82–137 |
128–209 |
FcCHOa |
40–81 |
45 – 105 |
165 – 235 |
g(α)b |
FO11 |
[−ln (1 − α)]1/2, A2 |
[−ln (1 − α)]2/3, A3/2 |
[−ln (1 − α)]2/3, F3/2 |
FcCHOa |
[−ln (1 − α)]1/3, A3 |
[−ln (1 − α)]2/3, A3/2 |
1 − (1 − α)1/2, R2/F1/2 |
Aα at α = 0.5 (min−1) |
FO11 |
8.4 × 104 |
3.8 × 109 |
2.1 × 1016 |
FcCHOa |
2.3 × 105 |
9.2 × 105 |
4.5 × 1013 |
Av. ΔS (J mol−1 K−1) |
FO11 |
−154.8 |
−74.4 |
36.6 |
FcCHOa |
−150.6 |
−126.8 |
−15.3 |
Av. ΔH (kJ mol−1) |
FO11 |
52.9 |
101.3 |
175.0 |
FcCHOa |
58.7 |
83.3 |
181.6 |
Av. ΔG (kJ mol−1) |
FO11 |
128.4 |
141.9 |
154.1 |
FcCHOa |
137.6 |
162.6 |
191.8 |
4. Conclusion
Thermal decomposition of ferrocene carboxaldehyde in the presence of oxalic acid under an oxygen atmosphere was systematically studied with the help of non-isothermal thermogravimetry. The decomposed product was identified to be pure hematite of nano-rod shape. The multistep thermogravimetry profiles were successfully deconvoluted, which showed that the decomposition occurs in six individual steps. It was analysed that the last three reaction steps were responsible for the production of hematite. The reaction kinetics for these three steps of decomposition were studied. Step-wise systematic dependence of the activation energy on the extent of conversion was observed. Here, the thermal decomposition progressed through different reaction mechanisms and with different range of reaction rates. Reactions between the gaseous products formed due to the decomposition of oxalic acid, and the intermediate product obtained on decomposition of ferrocene carboxaldehyde were plausibly held accountable for the production of hematite. Use of a co-precursor significantly enhanced the thermal reaction, lowered the reaction completion temperature, and modified the activation energy and reaction rate. Probably, during such reactions, the reactivity of the reactant particles is modified by factors like defect formation, particle disintegration, development of crystal strain, etc. Presently, by increasing the temperature during the decomposition, these factors might have been severely affected, which resulted in the different reaction mechanism functions for each step. This affected the dependence of the activation energy and reaction rate on the extent of conversion for each step. The present study describes the reaction kinetics of a co-precursor-driven solid state reaction, leading to hematite nanomaterial, and proposes that a suitable blend of precursor and co-precursor may produce important iron oxide nanomaterials of various sizes through solid-state thermal decomposition.
Author contributions
Material preparation, data collection and analysis were performed by M. C., whereas S. K. was involved in some data collection and analysis. A. B. conceptualized the problem and designed the study. The first draft of the manuscript was written by M. C. and finalized by A. B. All authors reviewed the final manuscript.
Conflicts of interest
The authors declare that they have no conflict interest.
Acknowledgements
S. K. is thankful to DST-INSPIRE, Government of India for providing a fellowship. Financial support for the thermogravimetry analyser (STA 449 F3 Jupiter) and powder X-ray diffractometer (Ultima IV of Rigaku) from the Department of Science and Technology (DST), Govt. of India through a grant (Ref. No. SR/FIST/PSI-157/2010) to the Department of Physics, Visva-Bharati University is gratefully acknowledged. The authors thank DST-PURSE Facility, Visva-Bharati University for providing the FE–SEM and EDX facilities.
References
- M. A. Shah and K. A. Shah, Nanotechnology: The Science of Small, Wiley Hoboken, NJ, 2019, vol 200 Search PubMed.
- K. A. Shah and F. A. Khanday, Nanoscale Electronic Devices and Their Applications, CRC Press, 2020 Search PubMed.
- Y. Li, J. Dang, Y. Ma and H. Ma, Hematite: A Good Catalyst for the Thermal Decomposition of Energetic Materials and the Application in Nano-Thermite, Molecules, 2023, 28(5), 2035 CrossRef CAS PubMed.
- A. Ali, H. Zafar, M. Zia, I. ul Haq, A. R. Phull, J. S. Ali and A. Hussain, Synthesis, Characterization, Applications, and Challenges of Iron Oxide Nanoparticles, Nanotechnol., Sci. Appl., 2016, 9, 49 CrossRef CAS PubMed.
- Z. L. Wang, Structural Analysis of Self-assembling Nanocrystal Superlattices, Adv. Mater., 1998, 10(1), 13–30 CrossRef CAS.
- S. J. Oldenburg, R. D. Averitt, S. L. Westcott and N. J. Halas, Nanoengineering of Optical Resonances, Chem. Phys. Lett., 1998, 288(2–4), 243–247 CrossRef CAS.
- C. J. Orendorff, T. K. Sau and C. J. Murphy, Shape-Dependent Plasmon-Resonant Gold Nanoparticles, Small, 2006, 2(5), 636–639 CrossRef CAS PubMed.
- M. A. Shah, Hydrothermal Synthesis of a–Fe2O3 Nanorods Prepared by a New and a Fast Route, Orient. J. Chem., 2009, 25(3), 489 CAS.
- B. Xu, B. Huang, H. Cheng, Z. Wang, X. Qin, X. Zhang and Y. Dai, α-Fe2O3 Hollow Structures: Formation of Single Crystalline Thin Shells, Chem. Commun., 2012, 48(52), 6529–6531 RSC.
- P. Tartaj, M. P. Morales, T. Gonzalez-Carreño, S. Veintemillas-Verdaguer and C. J. Serna, The Iron Oxides Strike Back: From Biomedical Applications to Energy Storage Devices and Photoelectrochemical Water Splitting, Adv. Mater., 2011, 23, 5243–5249 CrossRef CAS PubMed.
- P. Tartaj, M. del Puerto Morales, S. Veintemillas-Verdaguer, T. González-Carreño and C. J. Serna, The Preparation of Magnetic Nanoparticles for Applications in Biomedicine, J. Phys. D: Appl. Phys., 2003, 36(13), R182 CrossRef CAS.
- C. Wu, P. Yin, X. Zhu, C. OuYang and Y. Xie, Synthesis of Hematite (α-Fe2O3) Nanorods: Diameter-Size and Shape Effects on Their Applications in Magnetism, Lithium Ion Battery, and Gas Sensors, J. Phys. Chem. B, 2006, 110(36), 17806–17812, DOI:10.1021/jp0633906.
- Z. Wu, K. Yu, S. Zhang and Y. Xie, Hematite Hollow Spheres with a Mesoporous Shell: Controlled Synthesis and Applications in Gas Sensor and Lithium Ion Batteries, J. Phys. Chem. C, 2008, 112(30), 11307–11313 CrossRef CAS.
- P. N. Dave and L. V. Chopda, Application of Iron Oxide Nanomaterials for the Removal of Heavy Metals, J. Nanotechnol., 2014, 1–14, DOI:10.1155/2014/398569.
- H. Wan, P. Rong, X. Liu, L. Yang, Y. Jiang, N. Zhang, R. Ma, S. Liang, H. Wang and G. Qiu, Morphological Evolution and Magnetic Property of Rare-earth-doped Hematite Nanoparticles: Promising Contrast Agents for T1-weighted Magnetic Resonance Imaging, Adv. Funct. Mater., 2017, 27(27), 1606821 CrossRef.
- V. Patsula, M. Moskvin, S. Dutz and D. Horák, Journal of Physics and Chemistry of Solids Size-Dependent Magnetic Properties of Iron Oxide Nanoparticles, J. Phys. Chem. Solids, 2016, 88, 24–30 CrossRef CAS.
- S. Kumar, M. Kumar and A. Singh, Synthesis and Characterization of Iron Oxide Nanoparticles (Fe2O3, Fe3O4): A Brief Review, Contemp. Phys., 2021, 62(3), 144–164, DOI:10.1080/00107514.2022.2080910.
- A. S. Teja and P.-Y. Koh, Synthesis, Properties, and Applications of Magnetic Iron Oxide Nanoparticles, Prog. Cryst. Growth Charact. Mater., 2009, 55(1–2), 22–45 CrossRef CAS.
- E. A. Campos, D. V. B. S. Pinto, J. I. S. Oliveira, E. C. Mattos and R. C. L. Dutra, Synthesis, Characterization and Applications of Iron Oxide Nanoparticles-a Short Review, J. Aerosp. Technol. Manage., 2015, 7, 267–276 CrossRef CAS.
- A. Dey, M. Zubko, J. Kusz, V. R. Reddy, A. Banerjee and A. Bhattacharjee, Solventless Synthesis and Characterization of α-Fe, γ-Fe, Magnetite and Hematite Using Iron(III)Citrate, Solid State Sci., 2019, 95, 105932–105939, DOI:10.1016/j.solidstatesciences.2019.105932.
- M. Chakraborty, S. Kundu, B. Das and A. Bhattacharjee, Thermal Transformation of 1-(Ferrocenyl) Ethanol to Iron Oxide Nanoparticles Based on Reaction Atmosphere: Analysis of the Decomposition Reaction Using Non-Isothermal Thermogravimetry, J. Therm. Anal. Calorim., 2023, 1–17 Search PubMed.
- S. Kundu, T. Sarkar, G. Ghorai, M. Zubko, P. K. Sahoo, M. Weselski, V. R. Reddy and A. Bhattacharjee, Study on Co-Precursor Driven Solid State Thermal Conversion of Iron (III) Citrate to Iron Oxide Nanomaterials, Appl. Phys. A, 2023, 129(4), 264 CrossRef CAS.
- S. Kundu, M. Chakraborty and A. Bhattacharjee, Solid-State Reaction of Ferrocene Controlled by Co-Precursor and Reaction Atmosphere Leading to Hematite and Cohenite Nanomaterials: A Reaction Kinetic Study, J. Phys. Chem. C, 2023, 127(37), 18397–18408, DOI:10.1021/acs.jpcc.3c04772.
- M. E. Brown, Introduction to Thermal Analysis: Techniques and Applications, Springer Science & Business Media, 2001, vol 1 Search PubMed.
- M. Chakraborty, A. Dey and A. Bhattacharjee, Insights into the Thermal Decomposition of Organometallic Compound Ferrocene Carboxaldehyde as Precursor for Hematite Nanoparticles Synthesis, Z. fur Phys. Chem., 2022, 236(9), 1137–1161, DOI:10.1515/zpch-2021-3175.
- X. Zhang, Applications of Kinetic Methods in Thermal Analysis: A Review, Eng. Sci., 2020, 14(2), 1–13 Search PubMed.
- J. Farjas and P. Roura, Isoconversional Analysis of Solid State Transformations: A Critical Review. Part I. Single Step Transformations with Constant Activation Energy, J. Therm. Anal. Calorim., 2011, 105(3), 757–766, DOI:10.1007/s10973-011-1446-4.
- S. Vyazovkin, Modification of the Integral Isoconversional Method to Account for Variation in the Activation Energy, J. Comput. Chem., 2001, 22(2), 178–183, DOI:10.1002/1096-987X(20010130)22:2<178::AID-JCC5>3.0.CO;2-#.
- M. J. Starink, The Determination of Activation Energy from Linear Heating Rate Experiments: A Comparison of the Accuracy of Isoconversion Methods, Thermochim. Acta, 2003, 404(1–2), 163–176, DOI:10.1016/S0040-6031(03)00144-8.
- J. Farjas and P. Roura, Modification of the Kolmogorov–Johnson–Mehl–Avrami Rate Equation for Non-Isothermal Experiments and Its Analytical Solution, Acta Mater., 2006, 54(20), 5573–5579 CrossRef CAS.
- J. Zheng, J. Huang, L. Tao, Z. Li and Q. Wang, A Multifaceted Kinetic Model for the Thermal Decomposition of Calcium Carbonate, Crystals, 2020, 10(9), 849 CrossRef CAS.
- G. Tan, Q. Wang, H. Zheng, W. Zhao, S. Zhang and Z. Liu, Concept of Variable Activation Energy and Its Validity in Nonisothermal Kinetics, J. Phys. Chem. A, 2011, 115(22), 5517–5524 CrossRef CAS PubMed.
- B. Das, J. Kusz, V. R. Reddy, M. Zubko and A. Bhattacharjee, Solventless Synthesis, Morphology, Structure and Magnetic Properties of Iron Oxide Nanoparticles, Solid State Sci., 2017, 74, 62–69, DOI:10.1016/j.solidstatesciences.2017.10.010.
- W. Qian, Q. Chen, F. Cao and C. Chen, Synthesis and Characterization of Polyhedral Graphite Particles, Open Mater. Sci. J., 2008, 2(1), 19–22, DOI:10.2174/1874088x00802010019.
- A. Bhattacharjee, A. Rooj, M. Roy, J. Kusz and P. Gütlich, Solventless Synthesis of Hematite Nanoparticles Using Ferrocene, J. Mater. Sci., 2013, 48(7), 2961–2968, DOI:10.1007/s10853-012-7067-x.
- D. Amara, J. Grinblat and S. Margel, Solventless Thermal Decomposition of Ferrocene as a New Approach for One-Step
Synthesis of Magnetite Nanocubes and Nanospheres, J. Mater. Chem., 2012, 22(5), 2188–2195 RSC.
- A. Dey, M. Zubko, J. Kusz, V. Raghavendra Reddy, A. Banerjee and A. Bhattacharjee, Thermal Synthesis of Hematite Nanoparticles: Structural, Magnetic and Morphological Characterizations, Int. J. Nano Dimens., 2020, 11(2), 188–198 CAS.
- A. Bhattacharjee, A. Rooj, D. Roy and M. Roy, Thermal Decomposition Study of Ferrocene [(C5H5)2Fe], J. Exp. Phys., 2014, 2014, 1–8, DOI:10.1155/2014/513268.
- A. Rooj, M. Roy and A. Bhattacharjee, No Title, Int. J. Chem. Kinet., 2017, 49(5), 319–332, DOI:10.1002/kin.21077.
- B. Das and A. Bhattacharjee, Effect of Co-Precursor Maliec Anhydride on the Thermal Decomposition of Acetyl Ferrocene: A Reaction Kinetic Analysis, Curr. Phys. Chem., 2019, 9(1), 22–35 CrossRef CAS.
- S. Vyazovkin, A. K. Burnham, J. M. Criado, L. A. Pérez-Maqueda, C. Popescu and N. Sbirrazzuoli, ICTAC Kinetics Committee Recommendations for Performing Kinetic Computations on Thermal Analysis Data, Thermochim. Acta, 2011, 520(1–2), 1–19, DOI:10.1016/j.tca.2011.03.034.
- N. Koga, S. Vyazovkin, A. K. Burnham, L. Favergeon, N. V. Muravyev, L. A. Perez-Maqueda, C. Saggese and P. E. Sánchez-Jiménez, ICTAC Kinetics Committee Recommendations for Analysis of Thermal Decomposition Kinetics, Thermochim. Acta, 2022, 179384, DOI:10.1016/j.tca.2022.179384.
- T. A. Ozawa, New Method of Analyzing Thermogravimetric Data, Bull. Chem. Soc. Jpn., 1965, 38(11), 1881–1886, DOI:10.1246/bcsj.38.1881.
- J. H. Flynn and L. A. A. Q. Wall, Direct Method for the Determination of Activation Energy from Thermogravimetric Data, J. Polym. Sci., Part B: Polym. Lett., 1966, 4(5), 323–328 CrossRef CAS.
- H. E. Kissinger, Reaction Kinetics in Differential Thermal Analysis, Anal. Chem., 1957, 29(11), 1702–1706 CrossRef CAS.
- T. Akahira and T. Sunose, Method of Determining Activation Deterioration Constant of Electrical Insulating Materials, Res. Rep. Chiba. Inst. Technol., 1971, 16(1971), 22–31 Search PubMed.
- T. Wanjun and C. Donghua, An Integral Method to Determine Variation in Activation Energy with Extent of Conversion, Thermochim. Acta, 2005, 433(1–2), 72–76 CrossRef.
- H. Akbi, A. Mekki and S. Rafai, New Linear Integral Kinetic Parameters Assessment Method Based on an Accurate Approximate Formula of Temperature Integral, Int. J. Chem. Kinet., 2022, 54(1), 28–41 CrossRef CAS.
- L. Vlaev, N. Nedelchev, K. Gyurova and M. Zagorcheva, A Comparative Study of Non-Isothermal Kinetics of Decomposition of Calcium Oxalate Monohydrate, J. Anal. Appl. Pyrolysis, 2008, 81(2), 253–262, DOI:10.1016/j.jaap.2007.12.003.
- F. J. Gotor, J. M. Criado, J. Malek and N. Koga, Kinetic Analysis of Solid-State Reactions: The Universality of Master Plots for Analyzing Isothermal and Nonisothermal Experiments, J. Phys. Chem. A, 2000, 104(46), 10777–10782, DOI:10.1021/jp0022205.
- G. I. Senum and R. T. Yang, Rational Approximations of the Integral of the Arrhenius Function, J. Therm. Anal., 1977, 11(3), 445–447 CrossRef.
- B. Das and A. Bhattacharjee, Kinetic Analysis of Nonisothermal Decomposition of Acetyl Ferrocene, Int. J. Chem. Kinet., 2018, 50(4), 259–272 CrossRef CAS.
- Y. Yamamoto and N. Koga, Thermal Decomposition of Maya Blue: Extraction of Indigo Thermal Decomposition Steps from a Multistep Heterogeneous Reaction Using a Kinetic Deconvolution Analysis, Molecules, 2019, 24(13), 2515, DOI:10.3390/molecules24132515.
- A. Perejón, P. E. Sánchez-Jiménez, J. M. Criado and L. A. Pérez-Maqueda, Kinetic Analysis of Complex Solid-State Reactions. A New Deconvolution Procedure, J. Phys. Chem. B, 2011, 115(8), 1780–1791, DOI:10.1021/jp110895z.
- L. W. Finger and R. M. Hazen, Crystal Structure and Isothermal Compression of Fe2O3, Cr2O3, and V2O3 to 50 Kbars, J. Appl. Phys., 1980, 51(10), 5362–5367 CrossRef CAS.
- U. Holzwarth and N. Gibson, The Scherrer Equation versus the'Debye-Scherrer Equation, Nat. Nanotechnol., 2011, 6(9), 534 CrossRef CAS PubMed.
- S. Saadi, I. Ahmed, A. Khazaal, K. Jasim and K. Harbbi, Application the Halder – Wagner to Calculation Crystal Size and Micro Strain by X-Ray Diffraction Peaks Analysis, NeuroQuantology, 2022, 20, 199–204, DOI:10.14704/nq.2022.20.1.NQ22074.
- N. Mitrovic, M. Milosevic, G. Mladenovic and N. Investigations, Computational and Experimental Approaches in Materials Science and Engineering, ed. G. M. Nenad Mitrovic, Springer, 2019 Search PubMed.
- S. Vyazovkin and W. Linert, Kinetic Analysis of Reversible Thermal Decomposition of Solids, Int. J. Chem. Kinet., 1995, 27(1), 73–84, DOI:10.1002/kin.550270109.
- G. Wang, X. Min, N. Peng and Z. Wang, The Isothermal Kinetics of Zinc Ferrite Reduction with Carbon Monoxide, J. Therm. Anal. Calorim., 2021, 146(5), 2253–2260, DOI:10.1007/s10973-020-10542-z.
- V. Dhyani, J. Kumar and T. Bhaskar, Thermal Decomposition Kinetics of Sorghum Straw via Thermogravimetric Analysis, Bioresour. Technol., 2017, 245, 1122–1129, DOI:10.1016/j.biortech.2017.08.189.
- D. Wu, G. Ding, W. Chi and L. Jiang, Research on the Pyrolysis Kinetics of Resin Powder on Waste Printed Circuit Board with Different Particle Sizes at Different Heating Rates: Inspiration for the Pyrolysis Mechanism, J. Therm. Anal. Calorim., 2021, 147, 8047–8059, DOI:10.1007/s10973-021-11102-9.
- R. K. Mishra, A. Sahoo and K. Mohanty, Pyrolysis Kinetics and Synergistic Effect in Co-Pyrolysis of Samanea Saman Seeds and Polyethylene Terephthalate Using Thermogravimetric Analyser, Bioresour. Technol., 2019, 289, 121608 CrossRef CAS PubMed.
- A. Leonhardt, S. Hampel, C. Mueller, I. Moench, R. Koseva, M. Ritschel, D. Elefant, K. Biedermann and B. Buechner, Synthesis, Properties, and Applications of Ferromagnetic-Filled Carbon Nanotubes, Chem. Vap. Deposition, 2006, 12(6), 380–387 CrossRef CAS.
- G. Lapidus, D. Barton and P. E. Yankwich, Kinetics and Stoichiometry of the Gas-Phase Decomposition of Oxalic Acid, J. Phys. Chem., 1964, 68(7), 1863–1865 CrossRef CAS.
- S. Yamamoto and R. A. Back, The Gas-Phase Photochemistry of Oxalic Acid, J. Phys. Chem., 1985, 89(4), 622–625 CrossRef CAS.
- T. Kakumoto, K. Saito and A. Imamura, Unimolecular Decomposition of Oxalic Acid, J. Phys. Chem., 1987, 91(9), 2366–2371 CrossRef CAS.
- J. Higgins, X. Zhou, R. Liu and T. T.-S. Huang, Theoretical Study of Thermal Decomposition Mechanism of Oxalic Acid, J. Phys. Chem. A, 1997, 101(14), 2702–2708, DOI:10.1021/jp9638191.
- B. Das, J. Kusz, V. R. R. Reddy, M. Zubko and A. Bhattacharjee, Solventless Synthesis, Morphology, Structure and Magnetic Properties of Iron Oxide Nanoparticles, Solid State Sci., 2017, 74, 62–69, DOI:10.1016/j.solidstatesciences.2017.10.010.
- Q. Duan, H. Cao, X. Li and J. Sun, Effects of Urea on the Thermal Decomposition Behavior of Ammonium Nitrate: A Reliable Thermal Safety Performance Enhancer, Process Saf. Environ. Prot., 2023, 171, 482–492, DOI:10.1016/j.psep.2023.01.032.
|
This journal is © The Royal Society of Chemistry 2023 |