Next Article in Journal
Impact of Isomeric Dicarboxylate Ligands on the Formation of One-Dimensional Coordination Polymers and Metallocycles: A Novel cis→trans Isomerization
Previous Article in Journal
Influence of Calcium Binding on Conformations and Motions of Anionic Polyamino Acids. Effect of Side Chain Length
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interpretation and Physical Meaning of Kinetic Parameters Obtained from Isoconversional Kinetic Analysis of Polymers

by
Nicolas Sbirrazzuoli
Institut de Chimie de Nice, Université Côte d’Azur, CNRS, UMR 7272, 06108 Nice, France
Polymers 2020, 12(6), 1280; https://doi.org/10.3390/polym12061280
Submission received: 7 May 2020 / Revised: 31 May 2020 / Accepted: 1 June 2020 / Published: 3 June 2020
(This article belongs to the Section Biobased and Biodegradable Polymers)

Abstract

:
Several successful examples—where physically sounded kinetic information was obtained from thermoanalytical data in different application fields, such as polymerization of thermosetting resins, biobased polymers and nanocomposites, crystallization and glass transition of semi-crystalline polymers and their nanocomposites—are here presented and discussed. It is explained how the kinetic parameters obtained from advanced isoconversional methods can be interpreted in terms of reaction mechanisms or changes in the rate-limiting step of the overall process, in the case of complex chemical reactions or complex physical transitions, and how these parameters can be used to extract model-fitting parameters.

Graphical Abstract

1. Introduction

Obtaining kinetic parameters with a real physical meaning in the case of complex reactions is not straightforward. Isoconversional methods are amongst the more reliable kinetic methods for the treatment of thermoanalytical data [1,2]. These methods are very effective when the reaction mechanism is not known and/or in the case of complex reactions involving several parallel, consecutive steps, or a combination of these, as well as large variations in viscosity. The main advantages of isoconversional methods are that they afford an evaluation of the effective activation energy, Eα, without assuming any form of the reaction model, and can detect changes in the rate-limiting steps of the overall reaction rate as measured with thermoanalytical techniques. Changes in the reaction mechanism are obtained by analyzing the changes in the Eα variation, called the Eα-dependency. For this reason, isoconversional methods have largely been used for studying the polymerization kinetics of numerous thermoset resins, but have also been found to be very efficient in the elucidation of physical transition kinetics, such as the crystallization and glass transition of polymers, the gelation of biopolymers, and melting [1,3]. This article aims at presenting several successful examples where physically sounded kinetic information was obtained from thermoanalytical data in different application fields, and explaining how the kinetic parameters obtained can be interpreted in terms of reaction mechanisms or changes in the rate-limiting step of the overall process. These examples show that advanced isoconversional methods provide for the obtaining of kinetic parameters with a real physical meaning and, in some cases, can be used to extract model-fitting parameters. These methods were initially applied with success to petro-based polymers, and are now frequently used in this field, thus the present study aims at promoting their use in the field of biobased and biodegradable polymer chemistry.

2. Theoretical Part

2.1. Advanced Isoconversional Methods

The isoconversional principle states that the reaction rate at a constant extent of conversion is only a function of temperature. Thus, it is possible in this case to compute a value of Eα for each value of α without any assumption regarding the reaction model.
The general form of the basic rate equation (dα/dt) is usually written as [1]:
d α d t = k ( T ) f ( α ) = A exp ( E R T ) f ( α )
where k(T) is the rate constant, T is the temperature, f(α) is the differential form of the mathematical function that describe the reaction model that represents the reaction mechanism, E is the activation energy, and A is the pre-exponential factor. Linearization of Equation (1) leads to Equation (2) which, if computations are performed at a constant α value, is known as the Friedman’s method:
ln ( d α d t ) α , i = ln [ A α f ( α ) ] E α R T α , i
To use this Equation, one should determine the reaction rate (dα/dt)α,i and the temperature Tα,i corresponding to a given extent of conversion α for the i temperature programs used. For a given extent of conversion, the term [Aα f(α)] is constant, and the plot on the left hand side of Equation (2) versus the reciprocal of temperature gives a straight line, which allows us to obtain a value of E for each value of α, named Eα, without any assumption regarding the reaction model f(α). Note that using this method it is not possible to independently evaluate Aα and f(α) without introducing the reaction model in the form of the function f(α). The advantages of differential methods, such as Friedman’s method, are that they use no approximations and can be applied to any temperature program, i.e., isothermal, nonisothermal, or even more complex temperature programs. In such a way, the sample temperature, if available, can be used for the kinetic computations. The main drawback is that differential isoconversional methods can sometimes reveal numerical instabilities [4]. To overcome some drawbacks of differential and usual integral methods, advanced isoconversional methods have been developed [1,2,3]. One of these methods is the advanced non-linear isoconversional method (NLN), or Vyazovkin’s method, described by Equations (3) and (4) [5,6,7].
Φ ( E α ) = i = 1 n j i n J [ E α , T i ( t α ) ] J [ E α , T j ( t α ) ]
J [ E α , T ( t α ) ] = t α Δ α t α exp [ E α R T ( t ) ] d t
where Eα is called the effective activation energy. The Eα value is determined as the value that minimizes the function Φ(Eα). A value of Eα is computed for each value of α lying in between 0.02 and 0.98 with a step of 0.02. For each i-th temperature program, the time tα,i and temperature Tα,i related to selected values of α are determined by an accurate interpolation, using a Lagrangian algorithm. This nonlinear kinetic method (NLN) allows one to handle a set of n experiments carried out under different arbitrary temperature programs Ti(t), and uses a numerical integration performed using the trapezoidal rule of the integral with respect to time. The method developed by N. Sbirrazzuoli and implemented in his software can treat any kind of isothermal or non-isothermal data from differential thermal analyses (DTA), differential scanning calorimetry (DSC), calorimetry (Calvet-type and reaction calorimetry), fast scanning calorimetry (FSC), thermogravimetric analysis (TGA), dynamic mechanical analysis (DMA), temperature modulated DSC (TMDSC) rheometry, heat capacity variations, or other techniques [4,6,8,9,10].
Application of Equations (1)–(4) requires evaluation of the extent of conversion at time t, αt. This parameter can be extracted from data provided by different techniques by assuming that the measured property is proportional to the reaction progress. If one uses DSC or calorimetric data, the extent of conversion is computed according to Equation (5)
α t = t i t ( d Q / d t ) d t t i t f ( d Q / d t ) d t = Q t Q t o t
where ti represents the time at first integration bound of the DSC signal, and tf the time at last integration bound selected when the reaction is finished. dQ/dt is the heat flow measured by DSC at time t, Qtot is the total heat released or absorbed by the reaction and Qt is the current heat change.
For integral data, such as TGA, DMA, rheometry or heat capacity data, the extent of conversion is computed according to Equation (6)
α t = X t X t i X t f X t i
where X represents any physical property measured at time t.

2.2. Interpretation of the Eα-Dependency

It has taken time for model-free methods, mainly represented by isoconversional methods, to be accepted by the kinetic community. Indeed, the proposed approach constitutes a real paradigm-change. The classical kinetic approach starts by writing the supposed mechanism of the reaction and then the kinetic parameters are deduced, most often by fitting the data. In the case of complex reactions, an accurate fit of the data can be obtained with kinetic parameters without physical meaning. The model-free approach is completely different, since it involves analyzing the kinetic data obtained without formulating hypotheses on the reaction mechanisms, and trying to deduce information about a possible mechanism or a change in the rate-limiting step of the overall reaction rate that may occur during the reaction progress. This is accomplished by analyzing the so-called Eα-dependency. In the case of a single-step process, Eα does not change with α. In this case, evaluation of kinetic parameters is quite easy. However, most reactions, and especially those occurring with high viscosity variations, do not obey this simple kinetic scheme. Significant variations in Eα with α indicate a kinetically complex process (i.e., multi-steps) and Eα variations can be associated with a change in the reaction mechanism. Recently, Sbirrazzuoli has refined this approach by pointing out the interest in using variations of Aα, f(α) and k(T), evaluated in a model-free way, to identify rate-limiting steps [11].
The way to evaluate independently Aα and f(α) in a model-free way is explained in detail in [10,11]. The study highlighted correlations between the reaction rate, activation energy dependency, rate constants for the chemically controlled part of the reaction and the diffusion-controlled part, activation energy and pre-exponential factors of the individual steps, and changes in the rate-limiting steps. It was shown how some parameters computed using Friedman’s and NLN methods can help to identify change in the rate-limiting steps of the overall polymerization mechanism as measured by thermoanalytical techniques. The principle is based on the idea that the overall reaction rate (dα /dt)α, calculated for a given extent of conversion (α), can be put in the form of a product of two terms, [Aα f(α) and exp(−Eα/RT)]. If the rate decreases and one of the two terms increases while the other decreases, it means that this term has become dominant, which indicates a change in the limiting step of the reaction. In addition, calculating Aα, f(α) and kα(T) separately (but in a “model-free” way) using the method of compensation effect allows one to distinguish a change in the rate due to a change in the number of effective collisions or “entropic” effects (i.e., term Aα), or due to a change in the reaction mechanism (i.e., term f(α)) [10].

3. Discussion

3.1. Polymerization and Crosslinking

As for any chemical reaction, the overall reaction rate of polymerization and crosslinking reactions depends on several factors—such as reactivity of chemical species, activities of reactants, temperature, number and efficient collisions between molecules, and the energetic barrier to overcome for reactants to be transformed into products—but diffusion becomes the rate-limiting factor when molecular motion becomes hindered by the growing polymer chains. After the gel point, the mobility is hindered in the polymer network, and molecules can only participate in a reaction after a series of diffusion jumps. Crosslinking of the epoxy–amine system has demonstrated a decrease of Eα, from high values to very low values [12]. The initial high values correspond to the activation energy of the uncatalyzed reaction (70–90 kJ mol·1), while the energy of the catalyzed reaction corresponds to lower values (45–55 kJ mol·1) [12,13]. In these pioneering works, the two systems under study were respectively an epoxy novolac (Ciba-Geigy; araldite LY 5052) cured with isophorone diamine, and diglycidyl ether of bisphenol A (DGEBA) cured with 1,3-phenylenediamine. Although these systems were not biobased, it will be shown in the following that the methodology developed several years ago can be applied to the development of biobased polymers. In addition, as the initial methods have been refined since their first applications, more information on mechanisms can be obtained nowdays. The very low values reported at the end of the reaction (10–20 kJ mol·1) are related to the diffusion of small monomers. The overall rate can be parametrized in term of an autocatalytic process during the chemically controlled part of the reaction, followed by a diffusion model for the late stages of the reaction. This result has been confirmed by numerical simulations [11].

3.1.1. Sourour and Kamal Model

The initial stages of the reaction can be described by the autocatalytic model proposed by Sourour and Kamal, where the rate constant k(T) of Equation (1) is expressed as Equation (7) [14].
k C ( T , α ) = A 1 exp ( E 1 R T ) + A 2 exp ( E 2 R T ) α m
A1 being the pre-exponential factor for the non-catalyzed reaction and A2 the pre-exponential factor for the catalyzed reaction. E1 and E2 respectively represent the activation energy of the non-catalyzed reaction and the activation energy of the catalyzed reaction, and m and n are kinetic exponents.

3.1.2. Diffusion Model

The Sourour and Kamal model is not suitable to describe the end of the polymerization process when the reaction is diffusion-controlled, thus a rate constant that takes into account the diffusion processes has been proposed [12,15];
k D ( T , α ) = D 0 exp ( E D R T + K α )
where D0 represents the pre-exponential factor of the diffusion-controlled reaction, K is a constant accounting for the effect of the chemical reaction on the change in diffusivity, and ED represents the activation energy of the diffusion-controlled reaction.

3.1.3. Effective Activation Energy

Application of the isoconversional principle to the autocatalytic model of Sourour and Kamal leads to Equations (9)–(11) [11,12]:
E α = k 1 ( T ) E 1 + k 2 ( T ) E 2 α m k 1 ( T ) + k 2 ( T ) α m
E α = A 1 exp ( E 1 / R T ) E 1 + A 2 exp ( E 2 / R T ) E 2 α m A 1 exp ( E 1 / R T ) + A 2 exp ( E 2 / R T ) α m
E α = ( A 1 / A 2 ) exp ( E 1 / R T ) E 1 + exp ( E 2 / R T ) E 2 α m ( A 1 / A 2 ) exp ( E 1 / R T ) + exp ( E 2 / R T ) α m
These equations show that the effective activation energy depends on both temperature and the extent of conversion. The activation energies of the non-catalyzed (E1) and catalyzed (E2) reactions are obtained by non-linear fitting of Equation (11) to the effective activation energy dependency with temperature. The fit is realized in the interval where the reaction is chemically controlled (generally for α around 0.3–0.5), and leads to the activation energies E1 and E2 [11,12]. Note that at a low degree of conversion (α → 0), the effective activation energy Eα yields an estimate for the activation energy of the non-catalyzed reaction E1.
Application of the isoconversional principle to the diffusion model leads to Equations (12)–(14) [11,12]:
E α = k ( T ) E D + k D ( T , α ) E 2 k ( T ) + k D ( T , α )
E α = A 2 exp ( E 2 / R T ) E D + D 0 exp ( E D / R T + K α ) E 2 A 2 exp ( E 2 / R T ) + D 0 exp ( E D / R T + K α )
E α = ( A 2 / D 0 ) exp ( E 2 / R T ) E D + exp ( E D / R T + K α ) E 2 ( A 2 / D 0 ) exp ( E 2 / R T ) + exp ( E D / R T + K α )
A non-linear fitting of Equation (14) is realized in the interval where the reaction is diffusion-controlled, and leads to the activation energies E2 and ED.
Estimation of the activation energy of multi-step processes can be simplified by using a large excess of one reactant. As an example, a large excess of amine (1,3-phenylenediamine) was used to estimate the activation energy of primary amine addition (non-catalyzed) independently of that of the secondary amine addition to epoxy. This was possible because the reactivity of the former is much higher than that of the latter. Indeed, the use of a large excess of amine has led to a quasi-constant value of the activation energy, which attests to the control of the overall rate by a single-step process in this case [16]. If k1(T) ≫ k2(T) then EαE1 in Equation (9). The obtained value was assigned to the values of the activation energy of the primary amine addition (55 kJ mol·1). For the stoichiometric counterpart of this system, it was shown that the reaction shift from chemical to diffusion control for α > 0.30. This was first deduced from the sole Eα–dependence analysis, obtained from DSC measurements and confirmed later by TMDSC and rheological measurements [12,16]. Figure 1 illustrates the correlation between the decrease of Eα and of the complex heat capacity (Cp*) attributed to vitrification. Vitrification during isothermal crosslinking can be detected by rheometry or TMDSC measurements [16,17]. After this stage, the molecular mobility is considerably reduced, and the reaction is diffusion controlled.
In the nonisothermal mode, the Eα values obtained at the end of the reaction are very low (10–20 kJ mol·1), and do not match with the activation energy of any chemical reaction, such as primary and secondary amine addition or etherification (Figure 2). In fact, these values correspond to the diffusion of small molecules, such as the unreacted chain segments of the polymer chain [19]. At the gel point the growing individual macromolecules become connected into a single network so that the epoxy system loses the ability to flow, however, it is still capable of a long-range motion of polymer chain segments if the temperature is increased, as in nonisothermal conditions. This results, in this case, in the increasing Eα-dependencies at the end of the reaction, instead of the decreasing tendency observed for reactions controlled by motions of small molecules [20,21]. This is explained by the devitrification and/or reactivation of the chemical reaction due to an increase in molecular mobility. Devitrification can be detected by the increase of the complex heat capacity Cp* measured by TMDSC. An example is reported on Figure 2, where devitrification occurs for a temperature higher than 140 °C. It is worth noting that this temperature domain is in perfect agreement with the increase of Eα attributed to the reactivation of chemical reactions (Figure 2). The agreement between Eα-dependencies and Cp* variations presented in Figure 1 and Figure 2 is remarkable, and confirms that the Eα-dependence reflects a change in the overall polymerization mechanism.
Combining DSC and oscillatory rheometry, it was reported for this system that gelation occurs at α around 0.55 at a temperature of 124 °C, and vitrification at α around 0.90 at a temperature of 146 °C [16]. These values are in perfect agreement with the data of Figure 2. The temperature of 124 °C corresponds to the beginning of the Cp*’s decrease, while the temperature of 146 °C corresponds to the minimum of Eα (note that for nonisothermal conditions the heating rates used for the evaluation of Eα are higher than the rate used for the measurement of Cp*, which induces a slight shift of the Eα values to a higher temperature). For a temperature higher than 140–150 °C devitrification occurs, as attested by the Cp* increase, and chemical reactions are reactivated.
The control via diffusion of monomers at the end of the reaction is considerably reduced when nanoclays are added to the epoxy–amine system (DGEBA/1,3-phenylenediamine) [22]. The molecular mobility is higher at the end of the reaction for mixtures containing clay independent of the nature of the alkylammonium ions present in the galleries of the clay, as reflected by the lower decrease of Eα, characteristic of diffusion control. Using the method proposed by Sbirrazzuoli, based on Equations (11) and (14), activation energies consistent for the uncatalyzed and for the catalyzed reaction were obtained [13]. In addition, it was shown that addition of modified clay induces less ineffective shocks, and the reactivity is increased due to an enhancement of the efficiency of collisions [13]. For the chemically controlled part of the reaction, the high increase in the efficiency of collisions of the amine addition, initiated by a proton donor, results in an enhancement of the ring-opening. In the diffusion-controlled part of the reaction, the pre-exponential factor, computed according to the model-free method described, is also increased by addition of modified clay, and was explained by a high increase of the frequency of diffusion jumps (1000 times) [10,13]. A similar result was reported for Furfuryl alcohol (FA) polymerization in presence of organically modified montmorillonite [23].
Because the fit of the Eα-dependence alone may lead to several sets of correlated pre-exponential factors (A1, A2, D0), it is preferable to use first the Eα-dependence to get initial values and to use the ratios (A1/A2, A2/D0) as constraints [13]. Thus, Equations (11) and (14) were proposed by Sbirrazzuoli [13] to reduce the degree of freedom and avoid obtaining local minima when applying non-linear minimization. Only ratios are obtained. Then, the values A1, A2, D0 are obtained by fitting experimental rates and Equations (10) and (13) as additional constraints [24]. The use of the reaction rate at a constant value of the extent of conversion [dα/dt]α,i (see Equation (6) of reference [11]) facilitates the computations because a lower number of points is involved in the minimization [11]. These data ([dα/dt]α,i) are automatically computed when one use the Friedman method.
For the chemically controlled part of the reaction studied, an apparent contradiction was reported [13]. The reaction with nanoparticles started at a lower temperature, while the values of Eα were higher. This was explained by the higher values found for the pre-exponential factor when nanoparticles were introduced. Indeed, a higher value of A resulted in an increase of the reaction rate, which compensated for the increase of E (Equation (1)). This experimental observation was confirmed later using simulated data [11]. These simulations have confirmed the relevance of the use of isoconversional analysis to identify the rate-limiting steps in complex (multi-step) reactions. These examples illustrate the significance of computing the pre-exponential factor dependence (lnAα-dependence) in a model-free way using compensation effect, in addition to Eα. Then, once Eα and Aα dependencies have been determined, the effective rate constant (kef) and mathematical function related to the mechanism f(α) can be evaluated for each extent of conversion. These parameters are then used to identify changes in rate-limiting steps and to gain insights into the reaction mechanism.
Mechanistic insights were obtained for the polymerization of cyanate esters with different bridging fragments ((CH3)2C, S, and O) using isoconversional kinetic analysis. A higher reactivity of a sulfur-bridged monomer compared to its carbon-bridged analog was found. This was explained by the larger value of the pre-exponential factor for the metal-catalyzed step of the reaction, due to a higher polarizability that induces a stronger intermolecular interaction between molecules and promotes preferred orientation of the reacting monomers [24].
Granando et al. reported a quasi-constant value of Eα for the crosslinking of a Phenol–Formaldehyde (PF) resole system, which was attributed to a one-step mechanism of condensation of the methylol group [25]. They concluded that the diffusion-controlled regime is much less limiting than in epoxy systems, where the curing kinetics drastically slows down. This system was compared to the crosslinking of a phenol–terephthalaldehyde (PTPA) resole, which exhibits a more complex mechanism. The first plateau of Eα was assigned to the condensation of the secondary alcohol onto phenols, followed by a slight decrease associated with the shift in the polymerization mechanism, from the first reaction to the addition onto phenol of the second aldehyde function for the one-fold reacted TPA. The slight increase of Eα during the later stage of curing PF and PTFA was assigned to diffusional restrictions due to the forming of a three-dimensional crosslinked network.
The polymerization kinetics of Furfuryl alcohol (FA) into PolyFurfuryl alcohol (PFA) with organically modified Montmorillonite (Nanomer I30E) under acidic catalysis were studied by both DSC and rheometry [26]. DSC can detect residual heat after curing, but it is generally not easy to use these data to compute kinetic parameters, especially at high degrees of crosslinking. In this case, rheometry can be an efficient technique following the end of the reactions occurring in a highly viscous media. Polymerization of FA into PFA generates water, which makes the rheological measurements difficult. Thus, for low temperatures DSC is favored, and at high temperatures rheology is more suitable when the objective is to use these data to make complex computations which require high-quality data. The perfect continuity between the Eα values obtained from DSC and the rheometric data presented in Figure 3 is additional proof of the physical meaning of the Eα values obtained. Indeed, the Eα-dependence cannot be determined by DSC at a high temperature because no signal is detected at the end of the reaction when the material is highly crosslinked, however, the remaining crosslinking of the system under heating can be adequately determined by rheometric measurements and the Eα-dependence computed using this technique. In this case, the parameter X of Equation (6), used to follow the residual reaction occurring in the rubbery state, is the elastic modulus G’.
Polymerization of FA into PFA with maleic anhydride (MA) as acidic initiator was studied by DSC and DMA [27]. DSC was used to monitor the reaction in the liquid state, while DMA was performed on cured materials, via heating samples from 20 to 350 °C following the residual reactions occurring in the solid state. An increase of the elastic modulus was observed for all the systems between 170 and 290 °C, while no thermal event was recorded at a temperature higher than 220 °C on the first DSC heating curves. Thus, this increase was attributed to residual crosslinks occurring in the rubbery state after intensive post-curing. The effective activation energy dependency (Eα), as a function of temperature computed from DMA measurements, shows increasing values, which corresponds to a kinetic control by the diffusion of long polymeric chains (Figure 4). As the PFA material progressively cross-links in the solid state, the molecular mobility becomes more restricted, which leads to an increase in the activation energy. The values are lower in the presence of solvents, attesting that chains are less constraining in this case. More remarkably, there is a perfect continuity between the Eα-dependencies obtained for polymerization from the liquid state (DSC measurements) and those obtained for additional cross-links occurring in the rubbery state (DMA measurements), the final values of the reaction starting from the liquid state being in good agreement with the first values of the reaction starting from the rubbery state. The combination of DSC and DMA data allows one to study the polymerization kinetics and get mechanistic information over a wide temperature range, starting from the liquid state and ending in the solid state. This perfect continuity is an additional indication of the physical meaning of the kinetic parameters obtained (Figure 4).
An increase of the Eα-dependence attributed to kinetic control via the diffusion of long polymeric chains was also reported for the crosslinking of DGEBA with diamine, and for the polymerization of 2-hydroxyethyl methacrylate [20,21]. As already mentioned, the increasing dependence is expected when the process becomes determined by the diffusion of large molecules, such as polymer chains or their long segments. The difference between decreasing to very low values (5 kJ mol−1) obtained for the crosslinking of the FA/I30E system (Figure 3) and the increasing to high values (220 kJ mol−1) obtained for the FA/MA system (already cured) (Figure 4) also illustrate this statement [23,26,27]. In the case of the FA/I30E system, the elastic modulus reaches a plateau and the system is in the glassy state [23,26], while for the FA/MA system the elastic modulus starts to increase at high temperatures, when chemical reactions are re-activated. The Eα-dependence decreases when the rate becomes limited by the diffusion of small molecules or short segments of the polymer chain, while the Eα-dependence increase when the rate is limited by the mobility of long segments of the polymer chains [3]. Stochastically modulated temperature DSC has highlighted the vitrification and devitrification that may occur during crosslinking [16,28]. Vitrification induces a shift from chemical to diffusion control, which results in decreasing Eα-dependency. In nonisothermal conditions, the system may devitrify after vitrification, leading to the reactivation of residual chemical reactions, which leads to increasing Eα-dependency (Figure 2) [20].
Similar conclusions were drawn by analyzing the Eα-dependencies resulting from the polymerization of epoxidized linseed oil (ELO) with three different biobased dicarboxylic acids, i.e., sebacic acid, suberic acid and succinic acid, used as crosslinkers [29]. The viscosity measured by nonisothermal dynamic rheometry of ELO cured with suberic or sebacic acids show a plateau at the end of the polymerization, while this was not the case when using succinic acid. In this case, the viscosity curve is steadily increasing, and the Eα-dependency shows an increase at the end of the reaction. In contrast, when ELO is cured with suberic or sebacic acids, the Eα-dependencies decrease to low values, as a result of the shift to diffusion control of small unreacted molecules, and chemical reactions are not reactivated at high temperatures in this case. In this study, the determination of Eα, lnAα and f(α) was done without assumption in the reaction mechanism. Indeed, evaluation of the kinetic triplet in a fully model-free way may allow new insights into the reaction mechanism to be obtained. This was achieved by coupling kinetic analysis, FTIR and chemo-rheology investigations. The consistency of the conclusions drawn from these different techniques for analyzing the complex polymerization mechanism is remarkable. It was shown that the early stage of the reaction is controlled by an autocatalytic step, which results from an increase of the amount of hydroxyl groups in the reaction medium for the three systems under studies, i.e., ELO/Succinic, ELO/Suberic and ELO/Sebacic. The formation of a sufficient amount of β-hydroxy-esters plays the major role while at low viscosity. For ELO/Suberic and ELO/Sebacic, the transport of unreacted molecules and/or their concentration in the rubbery state determine and explain the deceleration of the reaction rate at the end of the reaction. With succinic acid, the reaction starts at a lower temperature and the reaction rate is higher for the early stages of the reaction, while the apparent activation energy is higher. Again, this seems to be contradictory. Computations have shown that this phenomenon is due to a higher frequency of diffusion jumps, expressed by higher values of the pre-exponential term. For ELO/Succinic acid system, the viscosity is higher at the end of the process, and only an increase of the temperature permits an increase of the diffusion. At the end of the process, the viscosity is high and the number and efficiency of collisions between molecules are weak, and only an increase of the temperature permits an increase of the diffusion. The pre-exponential factor is much lower for longer chains (i.e., sebacic and suberic), leading to a much more pronounced difference in the k(T) values between ELO/Succinic and the two other systems. Such a difference, between the efficiency of collision for long chains and short chains, is somehow compensated by the influence of the reaction mechanism f(α). FTIR spectrum of the ELO/Succinic system did not highlight remaining carboxylic acid after cross-linking. This agrees with the lower f(α) values for this system, which attests to a progressively lower concentration of carboxylic group at the end of the reaction. The high decrease of the term f(α) explains the lower values of the reaction rate for succinic acid. The viscosity increase after the vitrification plateau for ELO/Succinic suggests the continuation of additional cross-links. As the temperature is continuously increasing, some chemical reactions can occur. This phenomenon involves a cooperative motion of long polymer chains formed during the crosslinking, and the energy barrier is higher, as compared to the low values obtained for the diffusion of small unreacted molecules in the case of sebacic and suberic acids. In contrast, the viscosity reaches a plateau for the ELO/Suberic and ELO/Sebacic systems at the end of the reaction where remaining carboxylic acids are still present, as shown by FTIR measurements. Nevertheless, these species are unable to react due to the poor transport of these longer chains, in agreement with a diffusion control evidenced by the decrease of Eα and Aα to low values at the end of the reaction. This is in perfect agreement with the shift from chemical to diffusion control of small molecules for curing with suberic and sebacic acids, and the diffusion control of long chain segment motions for curing with succinic acid.
Humins are a biomass-derived material, a co-product of the acid-catalyzed conversion of cellulose and hemicellulose to platform chemicals. The Eα values computed from non-isothermal rheometric data on humins, and humins with acid initiator (p-Toluenesulfonic acid monohydrate; pTSA), show decreasing dependencies (Figure 5). The first steep decrease was attributed to an autocatalytic step. The first value of 110 kJ mol·1 represents the activation energy for the non-catalyzed crosslinking of crude humins. A second step is identified from α > 0.35, and at this stage the viscosity has strongly increased. The decrease to 25 kJ.mol·1 is in good agreement with a diffusion control at the end of the reaction. For α > 0.8, the high temperature reached might increase molecular mobility, promoting chemical reactions to be reactivated, and the rate is controlled by motions of longer chain segments, thus Eα increases [30]. For humins with the pTSA system, the initiation reactions occur at significantly lower temperatures and with lower value of Eα, which is in good agreement with the accelerating effect of pTSA.
In addition, the resulting Eα decrease is in good agreement with the data obtained from isothermal DSC and complex heat capacity (Cp*), obtained from quasi-isothermal Stochastically Temperature Modulated DSC (TOPEM©) [30].

3.2. Crystallization of Semi-Crystalline Polymers

At temperatures close to the melting temperature, the nucleation rate is low, and becomes the rate-limiting step for the crystallization rate of polymers from the melt. If the crystallization temperature is decreased, the crystal growth rate increases. This is typical of an anti-arrhenian behavior, which is explained by the control of the crystal growth rate by nucleation [3]. The opposite situation occurs when a sample is crystallized upon heating from the solid or glassy state. At a temperature above the glass transition temperature (Tg) and well below the melting temperature (Tm), the nucleation rate is high, and the diffusion is the rate-limiting step for the crystal growth rate. In this temperature domain, the crystallization rate obeys the classical arrhenius behavior. Application of isoconversional analysis to the non-isothermal crystallization of semicristalline polymers resulted in decreasing Eα values for crystallization upon heating from the solid state, and negative increasing Eα values for crystallization upon cooling from the melt [31,32,33]. It was shown that this result is in perfect agreement with the Hoffman–Lauritzen (HL) theory of secondary crystallization, which describes the dependence of the crystal growth rate on temperature. A consequence of this is that the effective activation energy Eα, computed with an isoconversional method, should be negative and increasing in the melt crystallization region, and positive and decreasing in the glass crystallization region. A value of Eα = 0 should be theoretically obtained at the temperature of maximum growth rate, where the crystal growth rate reaches the maximum and the effective activation energy Eα changes from a negative to a positive value.
Then, an entirely novel approach was proposed that allows for the evaluating of the parameters of the Hoffman–Lauritzen theory of polymer crystallization as derived via isoconversional analysis. This new method represents a unique way to extract important parameters regarding crystallization as well as additional information on the mechanisms, such as the anti-arrhenian dependence of the crystallization rate, with respect to the temperature and the mechanisms of chain folding. The dependence obtained gives real information regarding the changes in the mechanism during crystallization.
The Hoffman–Lauritzen theory describes the temperature dependence of the growth rate G measured microscopically, and is presented in Equation (15):
G = G 0 exp ( U * R ( T T ) ) exp ( K g T Δ T f )
where G0 is the pre-exponential factor, U* is the activation energy of the segmental jump (associated with diffusion process), T m 0 is the equilibrium melting temperature, ∆T = T m 0 T is the undercooling, f = 2T/( T m 0 + T) is the correction factor, and T is a hypothetical temperature where motion associated with viscous flow ceases, which is taken 30 K below the glass transition temperature Tg (i.e., T = Tg − 30 K).
Application of the isoconversional principle results in:
E α ( T ) = U * T 2 ( T T ) 2 + K g R ( T m 0 ) 2 T 2 T m 0 T ( T m 0 T ) 2 T
where U* corresponds to the activation energy of molecular diffusion at the interfacial boundary between melt and crystals, and Kg corresponds to the activation energy for the nucleation of a crystal with critical size. According to the Hoffman–Lauritzen theory, the crystallization rate reaches the maximum for a given temperature, denoted by Tmax. If the crystallization temperature Tc is in the region TmaxTm (Tm: melting temperature), the sample will follow an anti-Arrhenian behavior, which is characterized by the negative value of the temperature coefficient (represented by the apparent activation energy in our calculations), which increases when the temperature decreases (crystallization on cooling). In this region, the crystallization rate is controlled by nucleation, which is characterized by a negative temperature coefficient. If, on the other hand, the crystallization temperature Tc is in the region TgTmax, the sample will follow the classical Arrhenian behavior, which will result in positive values for the apparent activation energy as calculated with advanced isoconversional methods. In this temperature domain, the growth rate becomes controlled by diffusion and leads to a positive temperature coefficient, which decreases as the temperature increases because diffusion is facilitated (crystallization on heating). This new approach was applied with success to many different thermoplastic polymers [34,35,36,37,38,39,40,41,42,43,44].
As an additional proof of the physical meaning of this approach, in all these examples the overall shape of the Eα-dependence computed with various isoconversional methods is always the same. Changes in the slope of this dependence were used to identify a change in the crystallization regime, while deviation from the variations predicted by the Hoffman–Lauritzen theory have allowed us to identify more complex phenomena not taken into account in this theory, such as primary nucleation, chain folding, lamellar thickening, etc. [38]. In contrast, the Avrami model is limited to a single-step process, which is generally not the case for polymers, as evidenced by the strong Eα-dependencies frequently reported [31,32,33,35,36,37,38,39,40,41,42,43,44]. The physical meaning of the Avrami exponent remains questionable in the case of crystallization of polymers [38]. It was shown that application of the isoconversional principle to the Avrami equation leads to a constant value of Eα, in contradiction with experimental observation [45]. An abrupt increase in the effective activation energy was observed for poly(ethylene therephtalate) (PET) at α ~ 0.85 (around 202 °C), and was assigned to a change in the crystallization mechanism, in agreement with microscopic measurements [31,32]. This result shows that the negative values of the experimental activation energy, as well as its variation, appear to be physically sound.
The method was successfully applied to the study of the crystallization of a new and promising biobased polymer, the poly(ethylene 2,5-furandicarboxylate) (PEF). It was shown that crystallization upon heating can be described by the classical HL rate equation, while melt crystallization shows an unexpected rate depression compared to the HL predictions (Figure 6). For the melt crystallization, it was shown that nucleation governs the process, which is reflected by the negative increasing values of the apparent activation energy as computed by an advanced isoconversional method. In the cold crystallization region (crystallization from the glassy state), it was shown that diffusion governs the process and is the rate-limiting step for crystallization. Analysis of the Eα-dependence has shown the existence of a change in the crystallization mechanism from the liquid at 170 °C. This result was later confirmed by other authors, i.e., Stoclet et al. and Tsanaktsis et al., using isothermal DSC, small angle X-ray diffraction (SAXS) and wide-angle X-ray Scattering (WAXS) measurements, which found that a defective crystalline form was induced upon isothermal crystallization at Tc ≤ 170 °C, while a more ordered one was induced above this critical Tc [46,47]. An atypical depression of the crystallization rate for crystallization during cooling from the melt at T < ~170 °C has been demonstrated, and was attributed to a change in the crystallization regime. The unexpected depression of the melt crystallization rate compared to the Hoffman–Lauritzen predictions was explained by transition to a regime where sporadic secondary nucleation limits the lateral growth rate. The faster crystallization rate, found when PEF was crystallized from the glass, suggests the more pronounced role of the other mechanisms, such as primary nucleation, at lower temperatures. The Hoffman–Lauritzen parameters U * and Kg were used to calculate the temperature corresponding to the maximum crystal growth rate (167 °C), which was also confirmed experimentally a year later. These parameters give information on the phenomena of nucleation and diffusion. Thus, crystallization on cooling from the melt is controlled by nucleation, and crystallization on heating (from the solid) by diffusion. The PEF crystals formed during heating or cooling exhibit similar structures (WAXD), but the crystal growth dynamics differ.
The non-isothermal Polytetrafluoroethylene (PTFE) ultrafast crystallization under tremendously fast cooling rates was achieved using fast scanning calorimetry (FSC) and ultra-fast scanning calorimetry (UFSC). The dependence of the effective activation energy (Eα) on the relative extent of crystallization reveals that PTFE crystallization follows a transition from regime II to regime III at around 315–312 °C, which highlights a progressive transition from nucleation to a crystal growth-controlled process. The Eα-dependencies computed from DSC, FSC and UFSC, for a wide range of cooling rates that cover a wide range of temperatures, were in perfect continuity, showing negative increasing values that tend to zero at the temperature of maximum growth rate, in agreement with the HL theory [36]. For nanocomposites based on PTFE and nanosilica, it appears from the HL parameters computed using the Eα-dependencies that the nanosilica promotes the crystallization at slow cooling rates by inducing a nucleating effect, but appears to hinder the crystallization at fast cooling rates (i.e., for temperatures located further from the melting), and that the nanosilica promotes the formation of stable and more perfect crystals even when high cooling rates are employed [42].
The study of polydimethylsiloxane (PDMS) silica nanocomposites also confirmed the anti-arrhenius behavior in the melt crystallization region, and the arrhenius behavior in the glass crystallization region [37]. The study shows that glass and melt crystallizations share the same dynamics. Nucleation is promoted in the presence of silica and crystal growth is favored. Secondary crystallization is enhanced in the presence of silica nanoclusters, and higher crystal perfection is obtained. Similar conclusions were drawn for poly(butylene succinate) (PBS)/Graphene nanocomposites. The secondary crystallization upon heating is promoted with the insertion of graphene, and higher crystal perfection is obtained. A decrease of the nucleation barrier and an increase of the crystal growth rate are observed in the presence of graphene, as reflected by the values of the HL parameters U* and Kg evaluated from the Eα-dependencies. This resulted in an increase of the crystal growth rate, explained by the simultaneous enhancement of nucleation and diffusion contributions in the presence of graphene nanosheets [41].

3.3. Relaxation and Glass Transition Kinetics

The Eα-dependency is computed according to the following Equations (3) and (4). The extent of conversion during the glass transition α can be evaluated from DSC data with the normalized heat capacity, as follows:
C p N = ( C p C p g ) | T ( C p e C p g ) | T
where CP is the observed heat capacity, and Cpg and Cpe are respectively the glassy and equilibrium (liquid) heat capacity. Because the values of Cpg and Cpe are temperature-dependent, they must be extrapolated into the glass transition region. Then, an advanced isoconversional method is applied. Analysis of the Eα-dependence has shown a decrease in the activation energy of the relaxation, which reflects an increase in the molecular mobility when the system changes from the glassy state to the rubbery state upon heating [48]. The application of isoconversional analysis to the glass transition in various glassy systems has shown that the variation in the activation energy correlates with the dynamic fragility of the glass-forming liquids, being especially great for polymers [49]. This decrease is in agreement with the cooperative character of this transition, and has been found for several different polymers [37,48,49,50].
The variation in the effective activation energy (Eα) during the glass transition, for both poly(ethylene terephthalate) (PET) and PEF amorphous and semi-crystalline samples, as a function of temperature, exhibits a decrease in the effective activation energy for all samples (Figure 7). This decrease can be interpreted in terms of the cooperative motion of the chain segments involved in the glass transition. In the glassy state, only local motions of the chain segments may occur because of the very low molecular mobility. As the temperature rises, the molecular motion increases, and the translational motion of the segments and the whole chain becomes possible. This results in a decrease of the energetic constraints, which is reflected in a decrease of the effective activation energy. The effective activation energies (Eα) are lower for PET, compared to the values obtained for PEF. The values obtained for amorphous PET are higher and correspond to a lower temperature domain, as compared with the values of the semi-crystalline samples. In addition, these values are in perfect continuity with the values obtained for amorphous samples. The situation is very different for PEF. The values are very close for both amorphous and semi-crystalline PEF. These results confirm the hypothesis of a reduced coupling between amorphous and crystal phases for PEF, in agreement with results obtained from the evaluation of Tg for the Mobile Amorphous Fraction (MAF) of PET, and of PET by TMDSC and DMA. Thus, analysis of the Eα-dependencies has confirmed the conclusion that the relaxation of MAF is less influenced by the presence of crystals in PEF than in PET [50].

4. Conclusions

Advanced isoconversional methods can treat isothermal and nonisothermal data, as well as data provided from different techniques (i.e., DSC, rheometry, DMA), using the same equation. This is very attractive with regard to getting complementary information, especially in the case of complex reactions that are very dependent on the temperature domain investigated and that often need to be analyzed from the liquid state (DSC, calorimetry, FTIR, rheometry) to the solid state (TGA, DMA). This is very informative because some techniques are more sensitive at the beginning of the reaction, when the material is liquid, and less sensitive at the end when the material is highly viscous or solid, as is the case in polymerization. In this case, rheometry is a very complementary tool to DSC for extracting kinetic information on the viscous system, and DMA allows the obtaining of kinetic parameters in the solid state, when the material is cured but reactions are still in progress. In addition, FTIR and NMR data give additional information at the molecular level, used to confirm the hypotheses drawn regarding the proposed mechanism. The use of Fast Scanning Calorimetry (FSC) permits one to considerably enlarge the temperature domain for the computations. The examples presented in this study show that physically sounded kinetic information can be obtained from isoconversional kinetic analysis, that the parameters obtained can be used as starting values for nonlinear model-fitting methods, and that the exploitation of the complementarity of experimental technics and kinetic methods can be very valuable for the elucidation of complex reaction mechanisms.

Author Contributions

Conceptualization, N.S.; methodology, N.S.; software, N.S.; validation, N.S.; formal analysis, N.S.; investigation, N.S.; resources, N.S.; data curation, N.S.; writing—original draft preparation, N.S.; writing—review and editing, N.S.; visualization, N.S.; supervision, N.S.; project administration, N.S.; funding acquisition, N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received not external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vyazovkin, S.; Sbirrazzuoli, N. Isoconversional kinetic analysis of thermally stimulated processes in polymers. Macromol. Rapid Commun. 2006, 27, 1515–1532. [Google Scholar] [CrossRef]
  2. Vyazovkin, S.; Burnham, A.K.; Criado, J.M.; Pérez-Maqueda, L.A.; Popescu, C.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim. Acta 2011, 520, 1–19. [Google Scholar] [CrossRef]
  3. Vyazovkin, S. Isoconversional Kinetics of Thermally Stimulated Processes; Springer: Berlin, Germany, 2015. [Google Scholar]
  4. Sbirrazzuoli, N. Is the Friedman method applicable to transformations with temperature dependent reaction heat? Macromol. Chem. Phys. 2007, 208, 1592–1597. [Google Scholar] [CrossRef]
  5. Vyazovkin, S. Evaluation of activation energy of thermally stimulated solid-state reactions under arbitrary variation of temperature. J. Comput. Chem. 1997, 18, 393–402. [Google Scholar] [CrossRef]
  6. Sbirrazzuoli, N.; Vincent, L.; Vyazovkin, S. Comparison of several computational procedures for evaluating the kinetics of thermally stimulated condensed phase reactions. Chemometr. Intell. Lab. 2000, 54, 53–60. [Google Scholar] [CrossRef]
  7. Vyazovkin, S. Modification of the integral isoconversional method to account for variation in the activation energy. J. Comput. Chem. 2001, 22, 178–183. [Google Scholar] [CrossRef]
  8. Sbirrazzuoli, N.; Brunel, D.; Elegant, L. Different kinetic equations analysis. J. Therm. Anal. 1992, 38, 1509–1524. [Google Scholar] [CrossRef]
  9. Sbirrazzuoli, N.; Girault, Y.; Elegant, 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]
  10. Sbirrazzuoli, N. Determination of pre-exponential factors and of the mathematical functions f(α) or G(α) that describe the reaction mechanism in a model-free way. Thermochim. Acta 2013, 564, 59–69. [Google Scholar] [CrossRef]
  11. Sbirrazzuoli, N. Advanced Isoconversional kinetic analysis for the elucidation of complex reaction mechanisms: A new method for the identification of rate-limiting steps. Molecules 2019, 24, 1683–1699. [Google Scholar] [CrossRef] [Green Version]
  12. Vyazovkin, S.; Sbirrazzuoli, N. Mechanism and kinetics of epoxy-amine cure studied by differential scanning calorimetry. Macromolecules 1996, 29, 1867–1873. [Google Scholar] [CrossRef]
  13. Alzina, C.; Sbirrazzuoli, N.; Mija, A. Hybrid nanocomposites: Advanced nonlinear method for calculating key kinetic parameters of complex cure kinetics. J. Phys. Chem. B 2010, 114, 12480–12487. [Google Scholar] [CrossRef]
  14. Sourour, S.; Kamal, M.R. Differential scanning calorimetry of epoxy cure: Isothermal cure kinetics. Thermochim. Acta 1976, 14, 41–59. [Google Scholar] [CrossRef]
  15. Stolin, A.M.; Merzhanov, A.G.; Malkin, A. YA. Polym. Eng. Sci. 1979, 19, 1074–1080. [Google Scholar] [CrossRef]
  16. Sbirrazzuoli, N.; Vyazovkin, S.; Mititelu, A.; Sladic, C.; Vincent, L. A study of epoxy-amine cure kinetics by combining isoconversional analysis with temperature modulated DSC and dynamic rheometry. Macromol. Chem. Phys. 2003, 204, 1815–1821. [Google Scholar] [CrossRef]
  17. Van Assche, G.; Van Hemelrijck, A.; Rahier, H.; Van Mele, B. Modulated differential scanning calorimetry: Isothermal cure and vitrification of thermosetting systems. Thermochim. Acta 1995, 268, 121–142. [Google Scholar] [CrossRef]
  18. Sbirrazzuoli, N.; Vyazovkin, S. Learning about epoxy cure mechanisms from isoconversional analysis of DSC data. Thermochim. Acta 2002, 388, 289–298. [Google Scholar] [CrossRef]
  19. Glasstone, S.; Laidler, K.J.; Eyring, H. The Theory of Rate Processes; McGrow-Hill: New York, NY, USA; London, UK, 1941. [Google Scholar]
  20. Sbirrazzuoli, N.; Mititelu-Mija, A.; Vincent, L.; Alzina, C. Isoconversional kinetic analysis of stoichiometric and off-stoichiometric epoxy-amine cures. Thermochim. Acta 2006, 447, 167–177. [Google Scholar] [CrossRef]
  21. Achilias, D.S. Investigation of the radical polymerization kinetics using DSC and mechanistic or isoconversional methods. J. Therm. Anal. Calorim. 2014, 116, 1379–1386. [Google Scholar] [CrossRef]
  22. Alzina, C.; Mija, A.; Vincent, L.; Sbirrazzuoli, N. Effects of incorporation of organically modified montmorillonite on the reaction mechanism of epoxy/amine cure. J. Phys. Chem. B 2012, 116, 5786–5794. [Google Scholar] [CrossRef]
  23. Zavaglia, R.; Guigo, N.; Sbirrazzuoli, N.; Mija, A.; Vincent, L. Complex kinetic pathway of furfuryl alcohol polymerization catalyzed by green montmorillonite clays. J. Phys. Chem. B 2012, 116, 8259–8268. [Google Scholar] [CrossRef]
  24. Galukhin, A.; Liavitskaya, T.; Vyazovkin, S. Kinetic and mechanistic insights into thermally initiated polymerization of cyanate esters with different bridging groups. Macromol. Chem. Phys. 2019, 220, 1900141. [Google Scholar] [CrossRef]
  25. Granado, L.; Tavernier, R.; Foyer, G.; David, G.; Caillol, S. Comparative curing kinetics study of high char yield formaldehyde- and terephthalaldehyde-phenolic thermosets. Thermochim. Acta 2018, 667, 42–49. [Google Scholar] [CrossRef]
  26. Zavaglia, R. Nouveaux Matériaux Biosourcés. Étude et élaboration de nanocomposites à partir de PolyAlcool Furfurylique. Ph.D. Thesis, University of Nice Sophia Antipolis, Nice, France, 20 July 2011. [Google Scholar]
  27. Falco, G.; Guigo, N.; Vincent, L.; Sbirrazzuoli, N. FA polymerization disruption by protic polar solvent. Polymers 2018, 10, 529. [Google Scholar] [CrossRef] [Green Version]
  28. Van Assche, G.; Van Hemelrijck, A.; Rahier, H.; Van Mele, B. Modulated differential scanning calorimetry: Non-isothermal cure, vitrification, and devitrification of thermosetting systems. Thermochim. Acta 1996, 286, 209–224. [Google Scholar] [CrossRef]
  29. Menager, C.; Guigo, N.; Vincent, L.; Sbirrazzuoli, N. Polymerization kinetic pathways of epoxidized linseed oil with aliphatic biobased dicarboxylic acids. J. Polym. Sci. 2020, 58, 1–11. [Google Scholar]
  30. Sangregorio, A.; Guigo, N.; de Jong, E.; Sbirrazzuoli, N. Kinetics and chemorheological analysis of cross-linking reactions in humins. Polymers 2019, 11, 1804. [Google Scholar] [CrossRef] [Green Version]
  31. Vyazovkin, S.; Sbirrazzuoli, N. Isoconversional analysis of the nonisothermal crystallization of a polymer melt. Macromol. Rapid Commun. 2002, 23, 766–770. [Google Scholar] [CrossRef]
  32. Vyazovkin, S.; Sbirrazzuoli, N. Isoconversional analysis of calorimetric data on nonisothermal crystallization of a polymer melt. J. Phys. Chem. B 2003, 107, 882–888. [Google Scholar] [CrossRef]
  33. Vyazovkin, S.; Sbirrazzuoli, N. Estimating the activation energy for nonisothermal crystallization of polymer melts. J. Therm. Anal. Calorim. 2003, 72, 681–686. [Google Scholar] [CrossRef]
  34. Vyazovkin, S.; Sbirrazzuoli, N. Isoconversional approach to evaluating the Hoffman-Lauritzen parameters (U* and Kg) from the overall rates of nonisothermal melt crystallization. Macromol. Rapid Commun. 2004, 25, 733–738. [Google Scholar] [CrossRef]
  35. Vyazovkin, S.; Stone, J.; Sbirrazzuoli, N. Hoffman-Lauritzen parameters for nonisothermal crystallization of poly(ethylene terephthalate) and poly(ethylene oxide) melts. J. Therm. Anal. Calorim. 2005, 80, 177–180. [Google Scholar] [CrossRef]
  36. Bosq, N.; Guigo, N.; Zhuravlev, E.; Sbirrazzuoli, N. Non-isothermal crystallization of polytetrafluoroethylene in wide range of cooling rates. J. Phys. Chem. B 2013, 117, 3407–3415. [Google Scholar] [CrossRef] [PubMed]
  37. Bosq, N.; Guigo, N.; Persello, J.; Sbirrazzuoli, N. Melt and glass crystallization of PDMS and PDMS silica nanocomposites. Phys. Chem. Chem. Phys. 2014, 16, 7830–7840. [Google Scholar] [CrossRef] [PubMed]
  38. Codou, A.; Guigo, N.; van Berkel, J.; de Jong, E.; Sbirrazzuoli, N. Nonisothermal crystallization kinetics of biobased poly(ethylene 2,5-furandicarboxylate) synthesized via direct esterification process. Macromol. Chem. Phys. 2014, 215, 2065–2074. [Google Scholar] [CrossRef]
  39. Tsanaktsis, V.; Bikiaris, D.N.; Guigo, N.; Exarhopoulos, S.; Papageorgiou, D.G.; Sbirrazzuoli, N.; Papageorgiou, G.Z. Synthesis, properties and thermal behavior of poly(decylene 2,5-furanoate): A biobased polyester from 2,5-furan dicarboxylic acid. RSC Adv. 2015, 5, 74592–74604. [Google Scholar] [CrossRef]
  40. Papageorgiou, D.G.; Guigo, N.; Tsanaktsis, V.; Exarhopoulos, S.; Bikiaris, D.N.; Sbirrazzuoli, N.; Papageorgiou, G.Z. Fast crystallization and melting behavior of a long-spaced aliphatic furandicarboxylate biobased polyester, poly(dodecylene 2,5-furanoate). Ind. Eng. Chem. Res. 2016, 55, 5315–5326. [Google Scholar] [CrossRef]
  41. Bosq, N.; Guigo, N.; Aht-Ong, D.; Sbirrazzuoli, N. Crystallization of poly(butylene succinate) on rapid cooling and heating: Towards enhanced nucleation by graphene nanosheets. J. Phys. Chem. C 2017, 121, 11915–11925. [Google Scholar] [CrossRef]
  42. Bosq, N.; Guigo, N.; Sbirrazzuoli, N. Crystallization of polytetrafluoroethylene in a wide range of cooling rates: Nucleation and diffusion in the presence of nanosilica clusters. Molecules 2019, 24, 1797–1816. [Google Scholar] [CrossRef] [Green Version]
  43. Papageorgiou, G.Z.; Achilias, D.S.; Bikiaris, D.N. Crystallization kinetics of biodegradable poly(butylene succinate) under isothermal and non-isothermal conditions. Macromol. Chem. Phys. 2007, 208, 1250–1264. [Google Scholar] [CrossRef]
  44. Papageorgiou, G.Z.; Achilias, D.S.; Karayannidis, G.P. Estimation of thermal transitions in poly(ethylene naphthalate): Experiments and modeling using isoconversional methods. Polymer 2010, 51, 2565–2575. [Google Scholar] [CrossRef]
  45. Guigo, N.; van Berkel, J.; de Jong, E.; Sbirrazzuoli, N. Modelling the non-isothermal crystallization of polymers: Application to poly(ethylene 2,5-furandicarboxylate). Thermochim. Acta 2017, 650, 66–75. [Google Scholar] [CrossRef]
  46. Stoclet, G.; Gobius du Sart, G.; Yeniad, B.; Lefebvre, J.M. Isothermal crystallization and structural characterization of poly(ethylene-2,5-furanoate). Polymer 2015, 72, 165–176. [Google Scholar] [CrossRef]
  47. Tsanaktsis, V.; Papageorgiou, D.G.; Exarhopoulos, S.; Bikiaris, D.N.; Papageorgiou, G.Z. Crystallization and polymorphism of poly(ehylene furanoate). Cryst. Growth Des. 2015, 15, 5505–5512. [Google Scholar] [CrossRef]
  48. Vyazovkin, S.; Sbirrazzuoli, N.; Dranca, I. Variation of the effective activation energy throughout the glass transition. Macromol. Rapid Commun. 2004, 25, 1708–1713. [Google Scholar] [CrossRef]
  49. Vyazovkin, S.; Sbirrazzuoli, N.; Dranca, I. Variation in activation energy of the glass transition for polymers of different dynamic fragility. Macromol. Chem. Phys. 2006, 207, 1126–1130. [Google Scholar] [CrossRef]
  50. Codou, A.; Moncel, M.; van Berkel, J.G.; Guigo, N.; Sbirrazzuoli, N. Glass transition dynamics and cooperativity length of poly(ethylene 2,5-furandicarboxylate) compared to poly(ethylene terephthalate). Phys. Chem. Chem. Phys. 2016, 18, 16647–16658. [Google Scholar] [CrossRef]
Figure 1. Effective activation energy (Eα) (left) and complex heat capacity (Cp*) obtained from quasi-isothermal TMDSC measurements (right) vs extent conversion (α) during isothermal curing. Inset: Total TMDSC heat flow vs time. Adapted with permission from Copyright© 2002 Elsevier Ltd. All rights reserved [18].
Figure 1. Effective activation energy (Eα) (left) and complex heat capacity (Cp*) obtained from quasi-isothermal TMDSC measurements (right) vs extent conversion (α) during isothermal curing. Inset: Total TMDSC heat flow vs time. Adapted with permission from Copyright© 2002 Elsevier Ltd. All rights reserved [18].
Polymers 12 01280 g001
Figure 2. Effective activation energy (Eα) vs temperature during nonisothermal curing at 1, 2 and 4 °C min·1 (left) and complex heat capacity (Cp*) obtained from nonisothermal TMDSC measurement at 1 °C min·1 (right). The dot lines indicate the temperature at gelation and vitrification measured by oscillatory rheometry; the subsequent extent of conversions was obtained using DSC and rheometry. Lozanges: last integration bound between 149 and 205 °C; circles: last integration bound between 228 and 297 °C. Adapted with permission from Copyright© 2020, John Wiley and Sons and adapted with permission from Copyright© 2003 Elsevier Ltd. All rights reserved [16,20].
Figure 2. Effective activation energy (Eα) vs temperature during nonisothermal curing at 1, 2 and 4 °C min·1 (left) and complex heat capacity (Cp*) obtained from nonisothermal TMDSC measurement at 1 °C min·1 (right). The dot lines indicate the temperature at gelation and vitrification measured by oscillatory rheometry; the subsequent extent of conversions was obtained using DSC and rheometry. Lozanges: last integration bound between 149 and 205 °C; circles: last integration bound between 228 and 297 °C. Adapted with permission from Copyright© 2020, John Wiley and Sons and adapted with permission from Copyright© 2003 Elsevier Ltd. All rights reserved [16,20].
Polymers 12 01280 g002
Figure 3. Effective activation energy dependency (Eα) as a function of temperature for FA/I30E system, by DSC (triangle) and by rheometry (circle) [26].
Figure 3. Effective activation energy dependency (Eα) as a function of temperature for FA/I30E system, by DSC (triangle) and by rheometry (circle) [26].
Polymers 12 01280 g003
Figure 4. Effective activation energy dependency (Eα) as a function of temperature obtained from DSC (triangle) and DMA (circle) measurements for FA/MA system [27].
Figure 4. Effective activation energy dependency (Eα) as a function of temperature obtained from DSC (triangle) and DMA (circle) measurements for FA/MA system [27].
Polymers 12 01280 g004
Figure 5. Left: evolution of the effective activation energy (Eα) as function of the temperature for crude humins (black squares) and humins/pTSA (blue circles); and right: complex viscosity variation with temperature [30].
Figure 5. Left: evolution of the effective activation energy (Eα) as function of the temperature for crude humins (black squares) and humins/pTSA (blue circles); and right: complex viscosity variation with temperature [30].
Polymers 12 01280 g005
Figure 6. Eα-dependency vs (a) relative degree of crystallinity (α), (b) temperature. Dash line: fit of Equation (16) for glass and melt crystallization. Adapted with permission from Copyright© 2014, John Wiley and Sons [38].
Figure 6. Eα-dependency vs (a) relative degree of crystallinity (α), (b) temperature. Dash line: fit of Equation (16) for glass and melt crystallization. Adapted with permission from Copyright© 2014, John Wiley and Sons [38].
Polymers 12 01280 g006
Figure 7. Dependence of the effective activation energy Eα upon temperature, for amorphous PEF (solid triangle), semi-crystalline PEF (open triangle) (~31% crystallinity), amorphous PET (solid circle) and semi-crystalline PET (open circles) (~34% crystallinity). Adapted and reproduced from with permission from the Royal Society of Chemistry [50].
Figure 7. Dependence of the effective activation energy Eα upon temperature, for amorphous PEF (solid triangle), semi-crystalline PEF (open triangle) (~31% crystallinity), amorphous PET (solid circle) and semi-crystalline PET (open circles) (~34% crystallinity). Adapted and reproduced from with permission from the Royal Society of Chemistry [50].
Polymers 12 01280 g007

Share and Cite

MDPI and ACS Style

Sbirrazzuoli, N. Interpretation and Physical Meaning of Kinetic Parameters Obtained from Isoconversional Kinetic Analysis of Polymers. Polymers 2020, 12, 1280. https://doi.org/10.3390/polym12061280

AMA Style

Sbirrazzuoli N. Interpretation and Physical Meaning of Kinetic Parameters Obtained from Isoconversional Kinetic Analysis of Polymers. Polymers. 2020; 12(6):1280. https://doi.org/10.3390/polym12061280

Chicago/Turabian Style

Sbirrazzuoli, Nicolas. 2020. "Interpretation and Physical Meaning of Kinetic Parameters Obtained from Isoconversional Kinetic Analysis of Polymers" Polymers 12, no. 6: 1280. https://doi.org/10.3390/polym12061280

APA Style

Sbirrazzuoli, N. (2020). Interpretation and Physical Meaning of Kinetic Parameters Obtained from Isoconversional Kinetic Analysis of Polymers. Polymers, 12(6), 1280. https://doi.org/10.3390/polym12061280

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