Next Article in Journal
Penetration Routes of Oxygen and Moisture into the Insulation of FR-EPDM Cables for Nuclear Power Plants
Next Article in Special Issue
Microfluidic-Assisted Formulation of ε-Polycaprolactone Nanoparticles and Evaluation of Their Properties and In Vitro Cell Uptake
Previous Article in Journal
Advanced Optical Wavefront Technologies to Improve Patient Quality of Vision and Meet Clinical Requests
Previous Article in Special Issue
Impact of Crosslinking on the Characteristics of Pectin Monolith Cryogels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Curing Kinetics of Bioderived Furan-Based Epoxy Resins: Study on the Effect of the Epoxy Monomer/Hardener Ratio

1
Department of Chemical, Materials and Production Engineering, University of Naples Federico II, Piazzale Vincenzo Tecchio, 80, 80125 Naples, Italy
2
Institute for Polymers, Composites and Biomaterials (IPCB)—CNR, Via Campi Flegrei 34, 80078 Pozzuoli, Italy
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(23), 5322; https://doi.org/10.3390/polym14235322
Submission received: 1 November 2022 / Revised: 27 November 2022 / Accepted: 27 November 2022 / Published: 5 December 2022

Abstract

:
The potential of furan-based epoxy thermosets as a greener alternative to diglycidyl ether of Bisphenol A (DGEBA)-based resins has been demonstrated in recent literature. Therefore, a deep investigation of the curing behaviour of these systems may allow their use for industrial applications. In this work, the curing mechanism of 2,5-bis[(oxiran-2-ylmethoxy)methyl]furan (BOMF) with methyl nadic anhydride (MNA) in the presence of 2-methylimidazole as a catalyst is analyzed. In particular, three systems characterized by different epoxy/anhydride molar ratios are investigated. The curing kinetics are studied through differential scanning calorimetry, both in isothermal and non-isothermal modes. The total heat of reaction of the epoxy resin as well as its activation energy are estimated by the non-isothermal measurements, while the fitting of isothermal data with Kamal’s autocatalytic model provides the kinetic parameters. The results are discussed as a function of the resin composition. The global activation energy for the curing process of BOMF/MNA resins is in the range 72–79 kJ/mol, depending on both the model used and the sample composition; higher values are experienced by the system with balanced stoichiometry. By the fitting of the isothermal analysis, it emerged that the order of reaction is not only dependent on the temperature, but also on the composition, even though the values range between 0.31 and 1.24.

Graphical Abstract

1. Introduction

Epoxy thermoset resins are polymeric materials widely used in the industrial field thanks to their excellent mechanical and thermal properties, as well as their chemical resistance [1]. Specific applications require remarkable structural characteristics of epoxy resins, which in turn are dictated by the nature of the epoxy monomer and hardener, as well as the network feature developed upon curing [2]. The curing reaction proceeds through different mechanisms depending on the nature of the crosslinking agent and the catalyst, as well as on the curing conditions [2].
In particular, the latter have a significant impact on the network structure, and as a result, the deeper the understanding of the curing mechanism and kinetics of an epoxy resin, the greater the control on the final properties of the material [3]. Several experimental techniques and methods are used to assess chemical reaction kinetics, such as Fourier transform infrared spectroscopy (FTIR) [4,5], nuclear magnetic resonance (NMR) [6,7], and dynamic mechanical analysis [8,9]. In this sense, one of the most widely used techniques to identify the kinetic parameters, such as the rate constant, activation energy, and reaction order, associated with a curing reaction is differential scanning calorimetry (DSC), used [10] under both isothermal and dynamic conditions [11]. The widespread use of the DSC technique is due to its higher feasibility compared to FTIR and NMR. In fact, NMR analysis requires highly expensive and difficult-to-use apparatus. FTIR instead, despite the easiness of use, is a less accurate technique due to the possible overlapping of the characteristic peaks associated with the reacting groups [12]. In fact, Shnawa [13] investigated the structure of a tannin-based epoxy resin (TER) by means of FTIR. Although this equipment provided some insights about the chemical architecture of this compound, it was possible to investigate the formation of the crosslinked network by DSC only. In particular, it was observed that the TER was successfully cured with an amine-based curing agent. Moreover, DSC analysis confirmed that TER blended to commercial epoxy resin could act as a curing accelerator by monitoring the exothermicity of the curing peaks and the completion of the curing temperatures. In another work, Yang and co-authors [14] demonstrated that the addition of a phosphorous/imidazole-containing compound to epoxy resins significantly decreases the curing time.
In recent studies carried out by the authors, the curing behavior of 2,5-bis[(oxiran-2-ylmethoxy)methyl]furan (BOMF) with methyl nadic anhydride (MNA) has been already reported [12]. This resin has been applied as a tinplate coating, showing excellent adhesion and chemical resistance, thus being a good substitute for DGEBA-based products in the field of metal packaging [15]. An interesting aspect of such systems relies on the possibility of tuning their properties by changing the epoxy/anhydride ratio, the nature of the catalyst, or by including an opportunely selected inorganic filler.
Considering the already well-established use of carbohydrate derivatives concerning the category of thermosetting epoxy resins, furanic compounds attracted an increased interest as a potential alternative in the replacement of oil-based ones. Under this perspective, furan-based epoxy resins have been studied in the last years but have not been fully exploited yet. To assess the potential use of these compounds, a study on the reaction kinetics by FTIR [12] and the coating properties [15] was already performed. To our knowledge, a deeper study on the reaction mechanisms and the kinetic parameters of BOMF/MNA has not been performed yet.
Taking all the above considerations into account, the goal of this work is not only to provide insight into the curing kinetics of the biobased BOMF/MNA resin, but also to evaluate the effects of different epoxy/anhydride combinations in the reacting mixture. To this aim, DSC has been selected as a tool to carry out these studies. In particular, non-isothermal and isothermal DSC analyses were conducted, and different thermokinetic methods were used to calculate the characteristic kinetic parameters of the BOMF/MNA curing process as function of the resin composition.

2. Materials and Methods

2.1. Materials

In this work 2,5-bis[(oxiran-2-ylmethoxy)methyl]furan (BOMF) was used as a biobased epoxy monomer and methyl-5-norbornene-2,3-dicarboxylic anhydride (methyl nadic anhydride, MNA, 90%), purchased from Sigma-Aldrich (Darmstadt, Germany), was used as a curing agent. To ensure the obtainment of a crosslinked material, 2-methylimidazole (2-MI, 99%), purchased from Acros Organics (Geel, Belgium), was used as an initiator. MNA and 2-MI were used as received without further purification, while BOMF was synthesized following the procedure present in the literature [12].

2.2. Sample Preparation

Epoxy/anhydride mixtures were prepared by mixing at room temperature the epoxy monomer and the curing agent in the proper amounts and the initiator to the extent of 0.5 wt.% of the total amount of the epoxy monomer and anhydride. The weights of BOMF (the epoxy monomer) and MNA (the anhydride curing agent) were calculated to obtain different molar epoxy-to-anhydride ratios, namely 1/1 (stoichiometric ratio), 1/0.8 (epoxy excess), and 0.8/1 (anhydride ratio).

2.3. Differential Scanning Calorimetry

Differential scanning calorimetric analyses (DSC) were performed by means of a DSC Q2000 (TA Instrument, New Castle, DE, USA) equipped with a refrigerator cooling system (RCS). The analyses were performed on samples of about 7–8 mg sealed in aluminum pans immediately after the preparation of the mixtures. The samples were stored at −30 °C for a maximum of 2 days. Measurements under inert nitrogen flux (50 mL/min) in both non-isothermal conditions, heating from 25 to 200 °C at heating rates of 1, 1.5, 3, 5, and 10 °C/min, and isothermal conditions, at temperatures of 90, 100, 110, 120, and 130 °C, were performed.

3. Modeling

The degree of cure (α) corresponds to the integrated area of a specific DSC peak and can range from 0 to 1. At the current time (t), α can be measured as reported [16,17,18]:
α ( t ) = 1 Δ H t o t 0 t d H = Δ H ( t ) Δ H t o t
where ∆H(t) = 0 t d H is the reaction enthalpy at the time t and ∆Htot is the total reaction enthalpy [19,20]. As commonly performed in kinetic analyses conducted by DSC [16], it is assumed that there is a direct proportionality among the heat flow (dQ/dt) detected during the analysis and the conversion rate (/dt) [21] according to the following equation:
d Q d t = Q t o t d α d t
In the case of the isothermal mode, the functional dependencies of the conversion rate can be regarded as the product of two basic functions, respectively, related to the dependency on the temperature and the degree of cure:
d α d t = k ( T )   f ( α )
where k(T) is the rate constant and f(α) is the reaction model. According to the transition-state theory [22], the activation energy represents the difference in the energy amount between the reactant molecule at its initial configuration and the corresponding transition-state molecule in its activated configuration. In the context of a chemical reaction, this quantity can be formally regarded as the minimum amount of energy required to activate the reactant molecules so that they can undergo chemical transformations [23]. Under this perspective, the activation energy can be enclosed in a rate constant that is dependent on the temperature, not the concentration of the reactants, and can be formally expressed by the Arrhenius equation:
k ( T ) = A   e ( E a R T )
where A is a preexponential factor, R is the gas constant, T is the temperature, and Ea is the activation energy. The reaction model f(α) reflects the reaction mechanism and, in the specific case of an epoxy curing reaction, it generally takes the form of a nth-order kinetic model (Equation (5)) or an autocatalytic process (Kamal’s model, Equation (6)) [23,24,25]:
f ( α ) = ( 1 α ) n
f ( α ) = α m ( 1 α ) n
In general, the nth-order and the autocatalytic model provide a suitable description of the behavior of the thermosetting materials [26,27]. In particular, the nth-order kinetics report a reaction rate proportional to the n exponent of the reactant concentration. In this case, the model assumes that only one reaction occurs during the curing process and, therefore, can lead to some limitations since more simultaneous reactions can take place during the whole process. With reference to the autocatalytic model, reaction products also participate in the reaction, leading to a reaction rate increase during the initial stage, reaching a maximum, and then decreasing [16,28].
The abovementioned models imply that a system achieves a full extent of cure, that is, α = 1. As a complete cure is not always observed in epoxy curing systems, Equations (5) and (6) can be manipulated so that the kinetic models are also applicable to incomplete cure, and Equations (7) and (8) are, respectively, obtained:
f ( α ) = ( α m a x α ) n
f ( α ) = α m ( α m a x α ) n
where αmax is the maximum degree of conversion (experimentally measured) and must be implemented to cope with the incomplete curing reaction. As these concepts are generally valid in the simple case of one-step chemical reactions, the aforementioned considerations and the formalism introduced for the activation energy can be extended to complex multistep reactions, taking into account the rate-controlling step. In our case, for multistep reaction mechanisms, the reaction rate can be expressed as a function of more than one rate constant. Indeed, Kamal’s model can be rearranged and expressed as reported in Equation (9) [29]:
d α d t = ( k 1 + k 2 α m ) ( α m a x α ) n
where k1 is the rate constant of a catalyzed n-order reaction and k2 is the rate constant of an autocatalytic m-order reaction.
When dealing with the non-isothermal mode, the formalism implied for the expression of the conversion rate in the case of the isothermal analysis (Equation (3)) can be adopted in the same way to express the functional dependencies. As heating experiments are performed at constant heating rates β = dT/dt > 0, the conversion rate can be expressed as the following:
d α d t = 1 β k ( T )   f ( α )
Typically, the rate of reaction at a constant heating rate as function of temperature exhibits a relative maximum within the temperature interval under study. Moreover, the rate constant can be expressed by the Arrhenius equation also in the case of the non-isothermal analysis. In our case, the Kissinger and the Flynn–Wall–Ozawa models have been applied for the determination of the energy of activation of the overall reaction.
The Kissinger model [30,31] retrieved the Ea at several peak temperatures (Tp) due to the different applied heating rates β, assuming that the reaction rate in a constant conversion is solely dependent on the temperature. In this case, the following linear relationship is reported:
ln ( β T p 2 ) = E a R T p + l n ( A R E a )
Under these conditions, Kissinger [30,31,32] reported that, for a series of non-isothermal tests, the following dependency of the heating rate on the inverse of the peak temperature holds true:
d   ln ( β / T p 2 ) d ( 1 / T p ) = E a R
In the case of the Flynn–Wall–Ozawa model, the Arrhenius rate law was integrated and then the Doyle approximation [33,34] was applied to obtain Equation (13). According to ASTM E698, [35] also in this case, the analysis was performed on several heating rates, with Tp being independent of the heating rate:
ln ( β ) = 1.052   E a R T p c o n s t . + l n ( A   E a R   f ( α ) )
where const. stands for a constant term depending on the integration method applied to the Arrhenius law, while f(α) is the reaction model. As well in this case, for a series of non-isothermal tests, the following dependency of the heating rate on the inverse of the peak temperature holds true: [36,37]
d   ln ( β ) d ( 1 / T p ) = 1.052     E a R
Both Equations (12) and (14) suggest two alternative ways to estimate the activation energy. Through the linear fitting of data derived by dynamic DSC measurements, from a plot of ln(β/ T p 2 ) [or ln(β)] versus 1/Tp, it is possible to determine the apparent activation energy value with further insight on the preexponential factor A in the case of the Kissinger method.

4. Results

Differential scanning calorimetry was used to investigate the curing behavior of epoxy/anhydride mixtures with different molar ratios in both dynamic and isothermal conditions. DSC thermograms collected at different heating rates (from 1 to 10 °C/min) for the three formulations (stoichiometric, epoxy excess, and anhydride excess) are respectively reported in Figure 1a,c,e with their respective peak temperature (Tp) values listed in Table 1. The absence of any residual post-cure reactions was certified by the absence of a further exotherm peak during the second heating ramp performed on all the studied samples. As no further peak was detected, samples were assumed as fully cured during the dynamic run. Therefore, the total heat of reaction ΔHtot was evaluated using the DSC analysis conducted at 10 °C/min for each sample. The calculated values of ΔHtot are 353.7, 301.1, and 377.2 J/g for the samples at a stoichiometric epoxy/anhydride ratio, epoxy excess, and anhydride excess, respectively. The thermal stability of the same three epoxy/anhydride systems at different molar ratios was investigated by means of thermogravimetric analysis (TGA) and reported in the supplementary information of a previous work [12]. The systems did not experience any thermal degradation in the range of temperatures utilized in the DSC analyses.
All the systems exhibit a single exothermic peak with a Tp, which shifts to higher temperatures by increasing the heating rate regardless of the epoxy/anhydride ratio. Furthermore, when an epoxy excess is used, at each heating rate, Tp is about 5 °C lower than that of other samples. It has been reported that a lower Tp appears to be an indicator of higher curing reactivity [38]; however, this occurrence will be confuted by considerations regarding the activation energy of the curing reaction, reported further in this section. This trend can be, instead, explained by the increased mobility of the growing macromolecular units associated with the decrease in the mixture viscosity when the epoxy monomer is added in stoichiometric excess [12].
As both the Kissinger and Ozawa approaches are based on fitting procedures applied at various heating rates [30,31,32,35,36,37] and rely on the assumption that the extent of reaction α at the peak temperature (Tp) is constant and independent of such heating rates, these methods were applied to calculate the activation energy as solely a function of the epoxy/anhydride ratio and provide a better elucidation on the stoichiometry effects. Figure 2a reports the values of ln(β/ T p 2 ) as function of 1/Tp; these data were linearly fitted, and the obtained parameters were used to evaluate the activation energy as well as the preexponential factor A for the overall reaction, according to Equation (11). Similarly, Figure 2b reports the linear plot of ln(β) versus 1/Tp according to the Ozawa model (Equation (13)), and its linear fitting can be observed. The calculated values of the kinetic parameters are summarized in Table 2.
Both the Kissinger and Ozawa models return similar values of activation energy for systems with a stoichiometric ratio and an excess of epoxide, indicating the same global reactivity of the two systems. In particular, it can be observed that the lower peak temperatures when epoxides are added in larger quantities, compared to those of mixtures at a stoichiometric ratio, may suggest an enhanced reactivity of the samples at epoxy excess. However, we remark that this observation may lead to an erroneous evaluation of the tendency of reaction of the system. On the other hand, for the system with an excess of anhydride, the Ea is lower, and this occurrence indicates an inhibited reactivity of the system despite the fact that Tp is comparable with the stoichiometric system. The Ea values evaluated with the Kissinger method are systematically lower than those evaluated with the Ozawa method, as expected by the comparison with similar data reported in the literature [30,31,32,35,36,37,39,40,41].
In the case of the isothermal DSC tests conducted on the BOMF/MNA resins at different epoxy/anhydride ratios, to prevent vitrification as much as possible, five temperature values lower than the onset temperature value of the cure peak detected in the dynamic DSC analyses were chosen. In Figure 1b,d,f, the thermograms of the isothermal analyses are reported in which the heat flows as a function of time (directly related to the rate of conversion according to Equation (2)) are depicted. As can be seen, when decreasing the curing temperature, a longer reaction time is needed for the system to complete the cure and to reach the maximum conversion rate. This behavior is characteristic of autocatalytic reaction kinetics [18,42]. As already observed, reactions characterized by n-order kinetics show a peak in isothermal scans, eventually approaching zero for longer times regardless of the temperature. In our case, all systems cured at 90 °C did not fully react during the analysis time (two hours). It was proved by performing a heating scan on the samples that underwent the isothermal analysis in which a residual cure was detected. In addition, lower peak temperatures and higher values of the conversion rate are associated with an increase in the curing temperature. For all the samples cured above 90 °C, no post-cure reactions were detected by dynamic scans performed on the samples that underwent isothermal tests.
The kinetic parameters of the curing reaction were calculated using the analyses of thermograms obtained from the isothermal DSC experiments. To this regard, an accurate baseline construction was required to obtain consistent results relevant to the kinetic analysis. In this work, for all the samples, a baseline was constructed, according to the procedure reported by Barton et al. [10], by performing an isothermal scan on the post-cured samples. The degree of cure as a function of time for each temperature (Figure 3a,c,e) was determined relating the reaction enthalpy obtained by isothermal analysis to the ΔHtot value derived from the dynamic measurements previously reported, following Equation (1).
The sigmoidal shape of the α–time curves, especially for low temperatures, is characteristic of autocatalytic processes for which the reaction rate is relatively low at the beginning of the process, and the conversion degree only slowly increase due to the initial reaction steps. As new functional groups, which can catalyze the curing process, are created by the first reaction steps, the cure degree rapidly increases until a high degree of cure is achieved. At this time, the reaction rate drastically decreases, and α reaches a plateau. Nevertheless, no sample reaches the maximum degree of cure (Table 3). In particular, a linear increase in αmax with temperature can be noticed for the system with an equimolar epoxy/anhydride ratio, while systems with unbalanced stoichiometry exhibited a non-linear behavior.
Reaction rate data were fitted to the Kamal model (Equation (9)) by minimizing the least square error function to evaluate the k1, k2, m, and n kinetic parameters with no prior assumption for the total reaction order m + n. The best-fitting values obtained for the kinetic parameters at various curing temperatures are listed in Table 4. In Figure 3b,d,f the experimental data (symbols) are compared to the fitted data (black solid lines). Overall, the curing degree rate is adequately described by autocatalytic behavior.
The total reaction order m + n as well as the rate constants show an increasing trend with the curing temperature at each epoxy/anhydride ratio.
Further analysis can be carried out on the obtained values for the rate constants k1 and k2, which express, respectively, the influence of the initiation step on the overall reaction and the dependence of the reaction rate on the species formed during the reaction steps (that also act as catalysts in autocatalytic processes) [18,42]. As the rate constants comply with the Arrhenius law (Equation (4)), by plotting their logarithm against the inverse of the cure temperatures (Figure 4), it is possible to gather values of the preexponential factor A and the activation energy Ea for the single-reaction processes (Table 5).

5. Discussion

One of the key factors influencing the morphological changes in the physical structure during thermosetting processes is the curing kinetics. In particular, the epoxy/anhydride crosslinking mechanism is extremely complex, as it involves the simultaneous occurrence of competitive reactions, such as etherification and esterification reactions, which globally lead to the network formation [39,43,44].
When anhydrides are used as curing agents in presence of imidazole as an initiator, the first reaction step is the initiation, a nucleophilic attack of the imidazole nitrogen to the anhydride group (Scheme 1a) that generates a zwitterion intermediate containing a quaternary nitrogen cation and an active oxyanion. The latter reacts with the anhydride or the epoxy group, giving rise to a carboxylate anion (Scheme 1b), which initiates chain-wise polymerization [43,44,45]. The formed anion can further react both with an epoxy group (polyetherification, Scheme 1c) or with another anhydride (polyesterification, Scheme 1d), propagating the chain growth through an alternating copolymerization anionic mechanism.
Due to the complexity of the mechanism underlying the generation of the thermosetting network of epoxy/anhydride resins, a single-step reaction mechanism is unlikely to correctly fit the experimental data. Moreover, the continuous generation of reactive groups throughout the propagation steps suggests the autocatalytic behavior of the epoxy/anhydride reaction. In fact, the autocatalytic nature of the BOMF/MNA curing reaction clearly is suggested by the shape of the thermograms derived from isothermal DSC analysis: the presence of a maximum at different times, depending on the cure temperature, is particular to autocatalytic processes [46,47,48].
Non-isothermal DSC analyses were elaborated by means of two different methods to derive information about the activation energy of the overall curing process. As expected, the Kissinger and Ozawa methods did not return the same value of Ea, [40,41] but the latter are reasonable values for an epoxy/anhydride cure, in trend with the reported literature [39,49]. In the specific case of BOMF/MNA, if an excess of anhydride is used, even though the maximum reaction rate is reached at a temperature comparable with the systems with a stoichiometric epoxy/anhydride molar ratio, the system is more reactive (lower activation energy, Table 2). An accelerating effect is indicated by a decrease in the activation energy since less energy from the reacting components is required to complete the reaction [40]. This finding could be related to the higher number of initiation sites formed in the presence of higher amounts of anhydride [12]. Consequently, the number of propagating chains is increased, and so is the number of reactive sites, with an autocatalytic effect responsible for a globally faster reaction (reflected also by the lower value of Ak, Table 2). We remark that the activation energy only provides data on the reaction rate, not about the reaction mechanism, and that both the Kissinger and Ozawa methods provide a single value of the activation energy for the whole process despite the fact that, in complex systems such as epoxy/anhydride resins, this parameter fluctuates with the curing period.
The isothermal analysis of the curing process can, instead, give a major degree of detail about the reaction mechanisms. Isothermal DSC data were fitted to an autocatalytic model developed by Kamal, which well characterizes systems wherein the reaction rate increases due to the catalytic impact of the groups formed during the reaction itself, to define the relationship among conversion, time, and temperature. First of all, the conversion α as a function of time was evaluated at different temperatures for the three studied systems (Figure 3a,c,e). The lack of a linear dependance of α with temperature for the systems with epoxy excess and anhydride excess, together with the partial extent of cure experienced by all the samples, can be explained by the occurrence, to different extents, of vitrification. The reactivity of the molecules dominates the rate of reaction at the beginning of a curing reaction. As the reaction proceeds, the degree of cure increases and the system first undergoes gelation (i.e., the development of molecular branching), which is related to an increase in viscosity and a subsequent worsening of the molecular mobility. When diffusion-controlled processes overstep the mobility of the reacting groups, the reaction becomes diffusion-controlled and vitrification occurs [50]. By correlating this phenomenon to the conversion degree, we can observe that at low conversion values, the network is far from being formed and the reacting species, still present in high concentrations, can easily diffuse, and the reaction proceeds fast. Once high conversions are reached, long chains are formed, resulting in a viscosity increase which limits the chain mobility; at this point, the reaction becomes diffusion-controlled, and a plateau value for the conversion is reached. When the system vitrifies, increasing the temperature above the glass transition temperature (Tg), the viscosity decreases, and subsequently, the chain mobility is somehow improved, enabling the reaction of the remaining unreacted functional groups [51]. For the systems studied in this work, all the temperatures selected for the isothermal analysis are above Tg [12]. However, due to the narrowness of the cure peak (Figure 1a,c,e), even though the temperature of isothermal analysis is higher than Tg, this latter is still too close to the peak temperature, so vitrification is unavoidable. It is noteworthy to highlight the behavior of the system with an epoxy excess, for which strong vitrification is verified at temperatures of 120 and 130 °C. Using the non-isothermal DSC analysis, it emerged that when BOMF is present in stoichiometric excess, the curing reaction begins at a lower temperature (Tonset) compared to the other systems. In particular, when this system was heated at 10 °C/min to the temperature of 120 °C, the degree of conversion reached was higher than the α experienced by the other samples at the same temperature. This behavior was also reflected in the dependence of conversion on temperature resulting from the isothermal analysis.
The fitting of experimental data with the mechanistic Kamal method was also carried out. As can be seen from Figure 3b,d,f, the model fits well the experimental data, especially for the system with epoxy excess, regardless of the conversion degree. For both mixtures at a stoichiometric ratio and anhydride excess, instead, the fit is in good agreement with the model, particularly at low α values, and a lack of accuracy at higher conversions is noted, when the reaction is controlled by diffusion phenomena rather than by kinetic parameters [52].
The parameters obtained by the fitting of the experimental data, reported in Table 4, reflect, as far as k1 and m are concerned, the contribution of non-autocatalytic phenomena (in this specific case of study, the initiation step). On the other hand, k2 and n account for the contribution of autocatalytic phenomena (i.e., propagation in epoxy/anhydride cure). Regardless of the resin composition, all four parameters generally increase with temperature, even though with a different trend, and k2 and m are systematically higher than k1 and n. Specifically, when the stoichiometry is unbalanced, the order of reaction increases faster with temperature. For the stoichiometric system, instead, m increases much more than n with temperature, indicating a stronger temperature dependance of propagation compared to initiation. Noteworthy are the values of m and n for the system with an excess of anhydride: they are almost the same, with an exception made at 90 °C. This occurrence indicates that initiation and propagation proceed simultaneously. Generally speaking, the overall reaction order, m + n, is in the range of 0.35–1.23 regardless the composition.
Further analysis was conducted on the rate constants obtained by the fitting procedure, extrapolating the values of the activation energy (Ea1 and Ea2) and the preexponential factor (A1 and A2) for k1 and k2, respectively. In particular, while the trend of k2 with temperature is well fitted by a linear trend, the behavior of k1 is less prominent of a linear behavior, although it is still well described by a linear fit. This is due to the fact that k1 is only calculated using two or three of the first experimental datapoints, which might result in inaccurate estimated values [52]. As a general consideration, it emerges that the propagation step is much more kinetically disfavored.
All the considerations expressed here can be summarized into practical advice for the application of BOMF/MNA resins. As the curing degree of samples with a stoichiometric ratio among the reactants has an increasing trend with the curing temperature, this composition should be used for applications in which significant differences in curing temperatures might be experienced by different samples, or also at different points of the same sample (i.e., thick samples). If lower curing temperatures must be used, instead, a composition with an epoxy excess should be chosen, which ensures a higher αmax at 90 °C. However, for the cure of this sample, it is mandatory to control the curing temperature, as an uncontrolled increase in the temperature will have a detrimental effect on the degree of cure. On the other hand, when energy consumption is a limiting factor, the use of anhydride excess may be preferred due to the lower activation energy overall required for the curing process.

6. Conclusions

The cure kinetics of bioderived furan-based epoxy resins reacted with anhydride in the presence of 2-methylimidazole as an initiator were investigated by both non-isothermal and isothermal differential scanning calorimetry (DSC). The cure kinetics of BOMF/MNA resulted in being dependent on the epoxy/anhydride molar ratio. The main findings of this study can be summed up as follows:
When anhydride is in a molar excess, the promotion of the initiation step leads to an overall improved reactivity. In fact, it was observed that the activation energy of the curing process decreased from 76.2 to 72.4 kJ/mol (Kissinger method) and from 78.9 kJ/mol to 75.3 kJ/mol (Ozawa method) for stoichiometric and anhydride excess, respectively;
The high reactivity of the system, independent of the resin composition, leads to incomplete cure due to vitrification. This occurrence was observed even at the higher curing temperature (130 °C) from the maximum degree of conversion (0.79, 0.57, and 0.65 for stoichiometric, epoxy excess, and anhydride excess, respectively);
Both propagation and initiation are favored by increasing the temperature, but with a different dependence, as a function of the composition.
The autocatalytic Kamal model fits well the experimental data but did not succeed in fully describing the complex epoxy anhydride curing process. The global order of reaction (m + n) was in the range 0.35–1.24, 0.31–1.11, and 0.33–1.23 for stoichiometric, epoxy excess, and anhydride excess, respectively.
Overall, the observed results enlighten a promising use of these compounds, and a deeper understanding of the reaction kinetics and mechanisms of these systems could help in promoting the application of these biobased materials at an industrial level, addressing the global dependence on oil-derived polymers.

Author Contributions

Conceptualization, A.M.; methodology, A.M. and C.B.; investigation, A.M., N.F. and C.B.; data curation, A.M., N.F. and C.B.; writing—original draft preparation, A.M. and N.F.; writing—review and editing, A.M. and C.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, S.; Chen, Z.; Tusiime, R.; Cheng, C.; Sun, Z.; Xu, L.; Liu, Y.; Jiang, M.; Zhou, J.; Zhang, H.; et al. Highly improving the mechanical and thermal properties of epoxy resin via blending with polyetherketone cardo. Compos. Commun. 2019, 13, 80–84. [Google Scholar] [CrossRef]
  2. Vidil, T.; Tournilhac, F.; Musso, S.; Robisson, A.; Leibler, L. Control of reactions and network structures of epoxy thermosets. Prog. Polym. Sci. 2016, 62, 126–179. [Google Scholar] [CrossRef] [Green Version]
  3. Guo, Q. Use of Thermosets in the Building and Construction Industry. In Thermosets Structure, Properties, and Applications, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 279–300. [Google Scholar]
  4. Musto, P.; Martuscelli, E.; Ragosta, G.; Russo, P.; Villano, P. Tetrafunctional epoxy resins: Modeling the curing kinetics based on FTIR spectroscopy data. J. Appl. Polym. Sci. 1999, 74, 532–540. [Google Scholar] [CrossRef]
  5. Fernández-Francos, X.; Rybak, A.; Sekula, R.; Ramis, X.; Serra, A. Modification of epoxy–anhydride thermosets using a hyperbranched poly (ester-amide): I. Kinetic study. Polym. Int. 2012, 61, 1710–1725. [Google Scholar] [CrossRef]
  6. Martin-Gallego, M.; González-Jiménez, A.; Verdejo, R.; Lopez-Manchado, M.A.; Valentin, J.L. Epoxy resin curing reaction studied by proton multiple-quantum NMR. J. Polym. Sci. Part B Polym. Phys. 2015, 53, 1324–1332. [Google Scholar] [CrossRef]
  7. Buist, G.J.; Hamerton, I.; Howlin, B.J.; Jones, J.R.; Liu, S.; Barton, J.M. Comparative kinetic analyses for epoxy resins cured with imidazole–metal complexes. J. Mater. Chem. 1994, 4, 1793–1797. [Google Scholar] [CrossRef]
  8. Nuñez, L.; Fraga, F.; Castro, A.; Fraga, L. Elastic Moduli and Activation Energies for an Epoxy/m-XDA System by DMA and DSC. J. Therm. Anal. 1998, 52, 1013–1022. [Google Scholar] [CrossRef]
  9. Stark, W. Investigation of the curing behaviour of carbon fibre epoxy prepreg by Dynamic Mechanical Analysis DMA. Polym. Test. 2013, 32, 231–239. [Google Scholar] [CrossRef]
  10. Barton, J.M. The Application of Differential Scanning Calorimetry (DSC) to the Study of Epoxy Resin Curing Reactions; Springer: Berlin, Germany, 1985; pp. 111–154. [Google Scholar] [CrossRef]
  11. Vyazovkin, S.; Sbirrazzuoli, N. Kinetic methods to study isothermal and nonisothermal epoxy-anhydride cure. Macromol. Chem. Phys. 1999, 200, 2294–2303. [Google Scholar] [CrossRef]
  12. Marotta, A.; Faggio, N.; Ambrogi, V.; Cerruti, P.; Gentile, G.; Mija, A. Curing Behavior and Properties of Sustainable Furan-Based Epoxy/Anhydride Resins. Biomacromolecules 2019, 20, 3831–3841. [Google Scholar] [CrossRef]
  13. Shnawa, H.A. Curing and thermal properties of tannin-based epoxy and its blends with commercial epoxy resin. Polym. Bull. 2020, 78, 1925–1940. [Google Scholar] [CrossRef]
  14. Yang, S.; Huo, S.; Wang, J.; Zhang, B.; Wang, J.; Ran, S.; Fang, Z.; Song, P.; Wang, H. A highly fire-safe and smoke-suppressive single-component epoxy resin with switchable curing temperature and rapid curing rate. Compos. Part Eng. 2020, 207, 108601. [Google Scholar] [CrossRef]
  15. Marotta, A.; Faggio, N.; Ambrogi, V.; Mija, A.; Gentile, G.; Cerruti, P. Biobased furan-based epoxy/TiO2 nanocomposites for the preparation of coatings with improved chemical resistance. Chem. Eng. J. 2020, 406, 127107. [Google Scholar] [CrossRef]
  16. Yousefi, A.; LaFleur, P.G.; Gauvin, R. Kinetic studies of thermoset cure reactions: A review. Polym. Compos. 1997, 18, 157–168. [Google Scholar] [CrossRef]
  17. Hardis, R.; Jessop, J.L.; Peters, F.E.; Kessler, M.R. Cure kinetics characterization and monitoring of an epoxy resin using DSC, Raman spectroscopy, and DEA. Compos. Part A Appl. Sci. Manuf. 2013, 49, 100–108. [Google Scholar] [CrossRef] [Green Version]
  18. Kenny, J.M. Determination of autocatalytic kinetic model parameters describing thermoset cure. J. Appl. Polym. Sci. 1994, 51, 761–764. [Google Scholar] [CrossRef]
  19. Kessler, M.R.; White, S.R. Cure kinetics of the ring-opening metathesis polymerization of dicyclopentadiene. J. Polym. Sci. Part A Polym. Chem. 2002, 40, 2373–2383. [Google Scholar] [CrossRef]
  20. Roudsari, G.M.; Mohanty, A.K.; Misra, M. Study of the Curing Kinetics of Epoxy Resins with Biobased Hardener and Epoxidized Soybean Oil. ACS Sustain. Chem. Eng. 2014, 2, 2111–2116. [Google Scholar] [CrossRef]
  21. Vyazovkin, S.; Chrissafis, K.; Di Lorenzo, M.L.; Koga, N.; Pijolat, M.; Roduit, B.; Sbirrazzuoli, N.; Suñol, J.J. ICTAC Kinetics Committee recommendations for collecting experimental thermal analysis data for kinetic computations. Thermochim. Acta 2014, 590, 1–23. [Google Scholar] [CrossRef]
  22. Laidler, K.J.; King, M.C. Development of transition-state theory. J. Phys. Chem. 1983, 87, 2657–2664. [Google Scholar] [CrossRef]
  23. Anslyn, E.V.; Dougherty, D.A. Transition State Theory and Related Topics. In Modern Physical Organic Chemistry; University Science Books: Melville, NY, USA, 2006; pp. 365–373. [Google Scholar]
  24. Ghaemy, M.; Barghamadi, M.; Behmadi, H. Cure kinetics of epoxy resin and aromatic diamines. J. Appl. Polym. Sci. 2004, 94, 1049–1056. [Google Scholar] [CrossRef]
  25. Li, C.; Fan, H.; Hu, J.; Li, B. Novel silicone aliphatic amine curing agent for epoxy resin: 1,3-Bis(2-aminoethylaminomethyl) tetramethyldisiloxane. 2. Isothermal cure, and dynamic mechanical property. Thermochim. Acta 2012, 549, 132–139. [Google Scholar] [CrossRef]
  26. Ferdosian, F.; Zhang, Y.; Yuan, Z.; Anderson, M.; Xu, C. Curing kinetics and mechanical properties of bio-based epoxy composites comprising lignin-based epoxy resins. Eur. Polym. J. 2016, 82, 153–165. [Google Scholar] [CrossRef]
  27. Wan, J.; Gan, B.; Li, C.; Molina-Aldareguia, J.; Kalali, E.N.; Wang, X.; Wang, D.-Y. A sustainable, eugenol-derived epoxy resin with high biobased content, modulus, hardness and low flammability: Synthesis, curing kinetics and structure–property relationship. Chem. Eng. J. 2016, 284, 1080–1093. [Google Scholar] [CrossRef] [Green Version]
  28. Hu, J.; Shan, J.; Wen, D.; Liu, X.; Zhao, J.; Tong, Z. Flame retardant, mechanical properties and curing kinetics of DOPO-based epoxy resins. Polym. Degrad. Stab. 2014, 109, 218–225. [Google Scholar] [CrossRef]
  29. Kamal, M.R. Thermoset characterization for moldability analysis. Polym. Eng. Sci. 1974, 14, 231–239. [Google Scholar] [CrossRef]
  30. Kissinger, H.E. Reaction Kinetics in Differential Thermal Analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  31. Boey, F.; Qiang, W. Experimental modeling of the cure kinetics of an epoxy-hexaanhydro-4-methylphthalicanhydride (MHHPA) system. Polymer 2000, 41, 2081–2094. [Google Scholar] [CrossRef]
  32. Barros, J.J.P.; Silva, I.D.D.S.; Jaques, N.G.; Wellen, R.M.R. Approaches on the non-isothermal curing kinetics of epoxy/PCL blends. J. Mater. Res. Technol. 2020, 9, 13539–13554. [Google Scholar] [CrossRef]
  33. Doyle, C.D. Kinetic analysis of thermogravimetric data. J. Appl. Polym. Sci. 1961, 5, 285–292. [Google Scholar] [CrossRef]
  34. Doyle, C.D. Estimating isothermal life from thermogravimetric data. J. Appl. Polym. Sci. 1962, 6, 639–642. [Google Scholar] [CrossRef]
  35. ASTM E698; Standard Test Method for Kinetic Parameters for Thermally Unstable Materials Using Differential Scanning Calorimetry and the Flynn/Wall/Ozawa Method. ASTM International: West Conshohocken, PA, USA, 2016.
  36. Ozawa, T. A New Method of Analyzing Thermogravimetric Data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef] [Green Version]
  37. Flynn, J.H. A General Differential Technique for the Determination of Parameters for d (α)/dt= f (α) A exp (− E/RT) Energy of Activation, Preexponential Factor and Order of Reaction (when Applicable). J. Therm. Anal. Calorim. 1991, 37, 293–305. [Google Scholar] [CrossRef]
  38. Deng, J.; Liu, X.; Li, C.; Jiang, Y.; Zhu, J. Synthesis and properties of a bio-based epoxy resin from 2,5-furandicarboxylic acid (FDCA). RSC Adv. 2015, 5, 15930–15939. [Google Scholar] [CrossRef]
  39. Kumar, S.; Samal, S.K.; Mohanty, S.; Nayak, S.K. Study of curing kinetics of anhydride cured petroleum-based (DGEBA) epoxy resin and renewable resource based epoxidized soybean oil (ESO) systems catalyzed by 2-methylimidazole. Thermochim. Acta 2017, 654, 112–120. [Google Scholar] [CrossRef]
  40. Harsch, M.; Karger-Kocsis, J.; Holst, M. Influence of fillers and additives on the cure kinetics of an epoxy/anhydride resin. Eur. Polym. J. 2007, 43, 1168–1178. [Google Scholar] [CrossRef]
  41. Sbirrazzuoli, N.; Girault, Y.; Elégant, L. Simulations for evaluation of kinetic methods in differential scanning calorimetry. Part 3—Peak maximum evolution methods and isoconversional methods. Thermochim. Acta. 1997, 293, 25–37. [Google Scholar] [CrossRef]
  42. Barghamadi, M. Kinetics and Thermodynamics of Isothermal Curing Reaction of Epoxy-4, 40-Diaminoazobenzene Reinforced with Nanosilica and Nanoclay Particles. Polym. Compos. 2010, 31, 1442–1448. [Google Scholar] [CrossRef]
  43. Pin, J.-M.; Sbirrazzuoli, N.; Mija, A. From Epoxidized Linseed Oil to Bioresin: An Overall Approach of Epoxy/Anhydride Cross-Linking. ChemSusChem 2015, 8, 1232–1243. [Google Scholar] [CrossRef]
  44. Fernàndez-Francos, X.; Ramis, X.; Serra, À. From curing kinetics to network structure: A novel approach to the modeling of the network buildup of epoxy-anhydride thermosets. J. Polym. Sci. Part A Polym. Chem. 2014, 52, 61–75. [Google Scholar] [CrossRef]
  45. Kumar, S.; Samal, S.K.; Mohanty, S.; Nayak, S.K. Epoxidized Soybean Oil-Based Epoxy Blend Cured with Anhydride-Based Cross-Linker: Thermal and Mechanical Characterization. Ind. Eng. Chem. Res. 2017, 56, 687–698. [Google Scholar] [CrossRef]
  46. Paramarta, A.; Webster, D.C. Bio-based high performance epoxy-anhydride thermosets for structural composites: The effect of composition variables. React. Funct. Polym. 2016, 105, 140–149. [Google Scholar] [CrossRef] [Green Version]
  47. Pan, X.; Sengupta, P.; Webster, D.C. High Biobased Content Epoxy–Anhydride Thermosets from Epoxidized Sucrose Esters of Fatty Acids. Biomacromolecules 2011, 12, 2416–2428. [Google Scholar] [CrossRef] [PubMed]
  48. Mahendran, A.R.; Wuzella, G.; Kandelbauer, A.; Aust, N. Thermal cure kinetics of epoxidized linseed oil with anhydride hardener. J. Therm. Anal. 2011, 107, 989–998. [Google Scholar] [CrossRef]
  49. Sun, G.; Sun, H.; Liu, Y.; Zhao, B.; Zhu, N.; Hu, K. Comparative study on the curing kinetics and mechanism of a lignin-based-epoxy/anhydride resin system. Polymer 2006, 48, 330–337. [Google Scholar] [CrossRef]
  50. Atarsia, A.; Boukhili, R. Relationship between isothermal and dynamic cure of thermosets via the isoconversion representation. Polym. Eng. Sci. 2000, 40, 607–620. [Google Scholar] [CrossRef]
  51. Zhao, L.; Hu, X. Autocatalytic curing kinetics of thermosetting polymers: A new model based on temperature dependent reaction orders. Polymer 2010, 51, 3814–3820. [Google Scholar] [CrossRef]
  52. Thomas, R.; Durix, S.; Sinturel, C.; Omonov, T.; Goossens, S.; Groeninckx, G.; Moldenaers, P.; Thomas, S. Cure kinetics, morphology and miscibility of modified DGEBA-based epoxy resin—Effects of a liquid rubber inclusion. Polymers 2007, 48, 1695–1710. [Google Scholar] [CrossRef]
Figure 1. Non-isothermal DSC thermograms of BOMF/MNA resins with (a) a stoichiometric epoxy/anhydride ratio, (c) epoxy excess, and (e) anhydride excess. Dynamic DSC analyses were conducted at five heating rates (i.e., 1 °C/min, 1.5 °C/min, 3 °C/min, 5 °C/min, and 10 °C/min) for each sample. Isothermal DSC thermograms of BOMF/MNA resins with (b) a stoichiometric epoxy/anhydride ratio, (d) epoxy excess, and (f) anhydride excess at five different temperatures (i.e., 90, 100, 110, 120, and 130 °C) for each sample.
Figure 1. Non-isothermal DSC thermograms of BOMF/MNA resins with (a) a stoichiometric epoxy/anhydride ratio, (c) epoxy excess, and (e) anhydride excess. Dynamic DSC analyses were conducted at five heating rates (i.e., 1 °C/min, 1.5 °C/min, 3 °C/min, 5 °C/min, and 10 °C/min) for each sample. Isothermal DSC thermograms of BOMF/MNA resins with (b) a stoichiometric epoxy/anhydride ratio, (d) epoxy excess, and (f) anhydride excess at five different temperatures (i.e., 90, 100, 110, 120, and 130 °C) for each sample.
Polymers 14 05322 g001aPolymers 14 05322 g001b
Figure 2. Plots of the heating rate versus 1/Tp. Dash lines represent the linear fit according to (a) Kissinger (Equation (11)) and (b) Ozawa (Equation (13)).
Figure 2. Plots of the heating rate versus 1/Tp. Dash lines represent the linear fit according to (a) Kissinger (Equation (11)) and (b) Ozawa (Equation (13)).
Polymers 14 05322 g002
Figure 3. Plots of the degree of conversion (α) as a function of time and conversion rate (/dt) as function of conversion for epoxy/anhydride mixtures with (a,b) a stoichiometric ratio, (c,d) epoxy excess, and (e,f) anhydride excess, respectively. Solid lines in (b,d,f) represent the fitting curves of the experimental data by the Kamal model (Equation (9)).
Figure 3. Plots of the degree of conversion (α) as a function of time and conversion rate (/dt) as function of conversion for epoxy/anhydride mixtures with (a,b) a stoichiometric ratio, (c,d) epoxy excess, and (e,f) anhydride excess, respectively. Solid lines in (b,d,f) represent the fitting curves of the experimental data by the Kamal model (Equation (9)).
Polymers 14 05322 g003aPolymers 14 05322 g003b
Figure 4. Linear plots of the logarithm of the rate constants, (a) ln(k1) and (b) ln(k2)), obtained by fitting the isothermal data with the Kamal model (Equation (9)) as function of the inverse of the cure temperature (1/Tc) at different epoxy/anhydride ratios.
Figure 4. Linear plots of the logarithm of the rate constants, (a) ln(k1) and (b) ln(k2)), obtained by fitting the isothermal data with the Kamal model (Equation (9)) as function of the inverse of the cure temperature (1/Tc) at different epoxy/anhydride ratios.
Polymers 14 05322 g004
Scheme 1. Reaction steps of epoxy/anhydride cure: (a) initiation, (b) propagation, (c) polyetherification and (d) polyesterification.
Scheme 1. Reaction steps of epoxy/anhydride cure: (a) initiation, (b) propagation, (c) polyetherification and (d) polyesterification.
Polymers 14 05322 sch001
Table 1. Peak temperature (Tp) values for the non-isothermal DSC analyses at different heating rates for mixtures with balanced and unbalanced epoxy/anhydride ratios.
Table 1. Peak temperature (Tp) values for the non-isothermal DSC analyses at different heating rates for mixtures with balanced and unbalanced epoxy/anhydride ratios.
10 °C/min5 °C/min3 °C/min1.5 °C/min1 °C/min
Stoichiometric158.4 °C145.6 °C136.7 °C126.1 °C119.0 °C
Epoxy Excess153.1 °C140.7 °C130.8 °C120.9 °C114.5 °C
Anhydride Excess158.3 °C144.9 °C137.3 °C125.9 °C116.7 °C
Table 2. Values of the global activation energy (EaK) and the preexponential factor (AK) determined by the Kissinger method and the global activation energy (EaO) determined by the Ozawa method.
Table 2. Values of the global activation energy (EaK) and the preexponential factor (AK) determined by the Kissinger method and the global activation energy (EaO) determined by the Ozawa method.
KissingerOzawa
EaK [kJ/mol]AK [s−1]EaO [kJ/mol]
Stoichiometric76.21.37 × 10778.9
Epoxy Excess75.21.40 × 10777.9
Anhydride Excess72.44.51 × 10675.3
Table 3. Values of the maximum degree of cure αmax reached after 2 h of isothermal DSC analysis for samples with different epoxy/anhydride ratios.
Table 3. Values of the maximum degree of cure αmax reached after 2 h of isothermal DSC analysis for samples with different epoxy/anhydride ratios.
αmax
90 °C100 °C110 °C120 °C130 °C
Stoichiometric0.480.740.770.780.79
Epoxy Excess0.700.760.700.430.57
Anhydride Excess0.510.750.670.840.65
Table 4. Kinetic parameters obtained by least square fitting of the curing degree rate dα/dt versus α by the Kamal method for systems with different epoxy/anhydride ratios as function of curing temperature (Tc).
Table 4. Kinetic parameters obtained by least square fitting of the curing degree rate dα/dt versus α by the Kamal method for systems with different epoxy/anhydride ratios as function of curing temperature (Tc).
Tc [°C]k1 [s−1]k2 [s−1]mnm + n
Stoichiometric901.31 × 10−42.27 × 10−40.290.060.35
1001.40 × 10−45.87 × 10−40.400.280.68
1102.18 × 10−41.41 × 10−30.580.370.95
1204.29 × 10−43.49 × 10−30.860.501.36
1305.09 × 10−45.92 × 10−30.750.491.24
Epoxy Excess901.21 × 10−42.86 × 10−40.180.130.31
1001.60 × 10−47.72 × 10−40.510.360.87
1101.89 × 10−41.67 × 10−30.570.400.97
1202.21 × 10−43.23 × 10−30.560.471.03
1302.82 × 10−46.18 × 10−30.650.461.11
Anhydride Excess901.19 × 10−42.39 × 10−40.280.050.33
1001.45 × 10−43.76 × 10−40.270.250.52
1101.77 × 10−47.49 × 10−40.270.310.58
1203.02 × 10−43.35 × 10−40.790.781.57
1303.28 × 10−45.03 × 10−40.590.641.23
Table 5. Values of the preexponential factor (A1) and activation energy (Ea1) from the linear fitting of the rate constant k1 and values of the preexponential factor (A2) and activation energy (Ea2) from the linear fitting of the rate constant k2 for systems at different epoxy/anhydride ratios.
Table 5. Values of the preexponential factor (A1) and activation energy (Ea1) from the linear fitting of the rate constant k1 and values of the preexponential factor (A2) and activation energy (Ea2) from the linear fitting of the rate constant k2 for systems at different epoxy/anhydride ratios.
A1 [s−1]Ea1 [kJ/mol]A2 [s−1]Ea2 [kJ/mol]
Stoichiometric5.49 × 10246.58.83 × 1010101.3
Epoxy excess4.13 × 10−124.56.22 × 10992.4
Anhydride excess7.6433.62.85 × 101098.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marotta, A.; Faggio, N.; Brondi, C. Curing Kinetics of Bioderived Furan-Based Epoxy Resins: Study on the Effect of the Epoxy Monomer/Hardener Ratio. Polymers 2022, 14, 5322. https://doi.org/10.3390/polym14235322

AMA Style

Marotta A, Faggio N, Brondi C. Curing Kinetics of Bioderived Furan-Based Epoxy Resins: Study on the Effect of the Epoxy Monomer/Hardener Ratio. Polymers. 2022; 14(23):5322. https://doi.org/10.3390/polym14235322

Chicago/Turabian Style

Marotta, Angela, Noemi Faggio, and Cosimo Brondi. 2022. "Curing Kinetics of Bioderived Furan-Based Epoxy Resins: Study on the Effect of the Epoxy Monomer/Hardener Ratio" Polymers 14, no. 23: 5322. https://doi.org/10.3390/polym14235322

APA Style

Marotta, A., Faggio, N., & Brondi, C. (2022). Curing Kinetics of Bioderived Furan-Based Epoxy Resins: Study on the Effect of the Epoxy Monomer/Hardener Ratio. Polymers, 14(23), 5322. https://doi.org/10.3390/polym14235322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop