Modeling the crystallization kinetics of polymers displaying high levels of secondary crystallization

The isothermal crystallization of poly(3-hydroxybutyrate-co-3-hydroxyvalerate) was evaluated using a range of models, namely, Avrami, simplified Hillier, Tobin, Malkin, Urbanovici–Segal, Velisaris–Seferis, and Hay. Two methods of model evaluation were used: determination of the parameters through traditional double log plots and curve fitting via nonlinear, multivariable regression. Visual inspection of the cumulative crystallization curves, calculation of the R2 value and standard error of the regression, and evaluation of the returned parameters were used to assess which model best describes the experimental data. The Hay model was found to generate the best fit, closely followed by the Velisaris–Seferis parallel model, suggesting that primary and secondary crystallization occur concurrently. The Avrami, Malkin, and Tobin models were found to perform well when the data is restricted to the region where primary crystallization dominates; however, they could not be used to successfully model the entire crystallization process. This work highlights the importance of selecting the most appropriate model for analyzing kinetics, especially when high levels of lamellar thickening and infilling occur during crystallization. The isothermal crystallization of poly(3-hydroxybutyrate-co-3-hydroxyvalerate) was evaluated using a range of models namely, Avrami, simplified Hillier, Tobin, Malkin, Urbanovici-Segal, Velisaris-Seferis and Hay. The Hay model was found to generate the best fit, closely followed by the Velisaris-Seferis parallel model, suggesting that primary and secondary crystallization occur concurrently. This work highlights the importance of selecting the most appropriate model for analyzing isothermal crystallization kinetics, especially when high levels of secondary crystallization occur.


Introduction
The kinetics of crystallization are an important consideration in polymer processing because they significantly affect the final thermal and mechanical properties of a material. If primary crystallization occurs slowly, when cooling from the melt phase, only a few nuclei are generated that grow into large spherulites, rendering the material brittle. As the rate of crystallization is dependent on the degree of undercooling from the melt, knowledge of the kinetics allows precise control of material morphology during processing.
Polymer crystallization consists of two processes: primary and secondary. Primary crystallization is initiated in the amorphous polymer through the formation of nuclei from small regions of chain entanglements. Growth of crystalline lamella proceeds from these nuclei to form spherulites that expand radially until they begin to impinge onto one another, marking the end of the primary process. The secondary crystallization process relates to thickening of the lamella due to reorganization of the crystal structure and the formation of small crystalline structures within the amorphous regions, known as infilling. This secondary process has been reported to follow a logarithmic time scale [1]; more recently, a square root time dependency has been found, consistent with a diffusion-controlled process [2].
The Avrami model is commonly used to determine the crystallization kinetics of polymers, where the Avrami exponent (n a ) provides information for the nucleation and growth mechanisms [3][4][5]. Limitations have been found to this model, including only being applicable to primary crystallization and meaningless fractional n a values [6][7][8]. Alternatives, namely, the simplified Hillier [9,10], Tobin [11][12][13], Malkin [14], and Urbanovici-Segal [15,16] models, many of which are extensions or adaptations of the standard Avrami analysis, have been proposed with the aim of accounting for the presence of secondary crystallization. Supaphol has previously analyzed the effectiveness of these models for the isothermal crystallization of syndiotactic polypropylenes [8,17] with others comparing combinations of these models for pure polymers [10,14,16,[18][19][20], blends [21][22][23] and nanocomposites [24][25][26]. Following an in-depth review of these papers, no singular model has been found to be consistently the best over a wide range of polymers. In addition, none of these analyses have been performed for polymers with substantial secondary crystallization, and where such crystallization is present, the data is limited to considering the primary crystallization region in isolation [25]. Kelly et al. have previously shown the Avrami model to be unsuitable for poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHB-co-HV) as a result of the high degree of secondary crystallization that follows the primary mechanism [27]. The researchers proposed an alternative method, the Hay model, giving meaningful integer n a and R 2 values in excess of 0.995 for PHB-co-HV. In addition, Velisaris and Seferis developed a model specifically to account for the primary and secondary crystallization processes occurring in either parallel or series [28]. This theory has been shown to be more effective than both the Avrami and Malkin models in the modeling of the isothermal crystallization of polypropylene and polypropylene composites containing talc [29].
Poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBco-HV) is a bacterially synthesized polymer that is formed via an energy storage mechanism when specific bacteria are fed a sugar source in the absence of key nutrients [30,31]. This sustainable production method, coupled with nontoxic, biocompatible, and biodegradable properties, has led to possible applications in the packaging industry [32]. However, a relatively low glass transition temperature of approximately 4°C results in secondary crystallization at room temperature, resulting in progressive embrittlement and hindering widespread use [33,34]. A thorough understanding of all aspects of the crystallization process for PHB-co-HV is paramount to be able to create suitable packaging materials; therefore, it is important to ensure that the most suitable models are utilized.
In this paper, the isothermal crystallization of PHB-co-HV is analyzed using the Avrami, simplified Hillier, Tobin, Malkin, Urbanovici-Segal, Velisaris-Seferis, and Hay models. Nonlinear regression is utilized to determine the fit of the models, and the significance of the parameters obtained and the suitability of each model are discussed. Finally, the most suitable model for the analysis of PHB-co-HV is proposed.

Theoretical background
Differential scanning calorimetry (DSC) is often used to study isothermal crystallization. As the evolution of crystallinity is directly proportional to the heat released during the process, the relative crystallinity as a function of time, θ(t) can be obtained through integration of the crystallization exotherms [17]. It is important, especially in cases of prolonged isothermal crystallization experiments with heat-labile polymers, that the temperatures selected do not lead to degradation of the material during the analysis. Many models have been proposed that can be used to analyze the development of crystallinity during primary and secondary processes.
The Avrami equation (Eq. (1)) is the most widely used equation to analyze the crystallization kinetics of polymers; however, it is only suitable for the primary mechanism, resulting in noninteger values being regularly determined for the Avrami exponent [8,14]. The parameters of the Avrami equation can be determined from the gradient (n a ) and antilogarithm of the y-intercept (k a ) of a double log plot of log(−ln(1 − X t )) versus log time.
where X t is the fractional crystallinity at time t; k a is the Avrami rate constant; and n a , which should be an integer between one and four, is the Avrami exponent representing the mechanism of nucleation and geometry of growth. The simplified Hillier model (Eq. (2)) aims to limit the data to the primary process alone [9,10]. This is achieved by incorporating a value, X p,∞ for the degree of crystallinity at the conclusion of the primary process. The value of X p,∞ is taken from the intersection of tangents drawn to the primary and secondary regions of fractional crystallinity versus time plots [35]. The parameters of the Hillier equation can be determined from the gradient (n a ) and antilogarithm of the y-intercept (k a ) of a double log plot of log (−ln(1 − X t / X p,∞ )) versus log time.
where X t is the fractional crystallinity at time t; X p,∞ is the fractional crystallinity at the end of the primary process; k a is the Avrami rate constant; and n a is the Avrami exponent. Tobin modified the Avrami equation in an attempt to improve the modeling for the latter stages of crystallization by considering impingement of the spherulites. The resulting equation (Eq. (3)) also contains both a rate constant (k t ) and exponent (n t ), which have similar physical meanings to the Avrami model with the exception that n t does not need to be an integer [11][12][13]. The parameters of the Tobin equation can be determined from the gradient (n t ) and antilogarithm of the y-intercept (k t ) in a double log plot of log(X t /(1 − X t )) versus log time.
where X t is the fractional crystallinity at time t; k t is the Tobin rate constant; and n t is the Tobin exponent.
Malkin derived his model (Eq. (4)) based on the dependence of the overall rate of crystallization on the sum of the rates of primary nuclei generation and crystal growth [14]: where C 0 is the Malkin exponent and is proportional to the ratio of the primary crystal growth rate to the primary nucleation rate and C 1 is the Malkin crystallization rate constant, and is related to the addition of various fractions of these two parameters [20]. Unlike previous models, these parameters can only be determined through curve fitting procedures; however, they can be estimated from the Avrami parameters, as illustrated in Eqs. (5) and (6) [8]: Urbanovici and Segal [15] derived a model (Eq. (7)) to describe the crystallization process, which simplifies to the general Avrami equation when r tends to one. The parameters of this model can only be determined via curve fitting.
where k us and n us are the Urbanovici-Segal rate constant and the Urbanovici-Segal exponent, respectively, and r is a parameter that satisfies the condition r > 0. In addition to these equations, researchers have also devised models to specifically account for the secondary crystallization process. The Velisaris-Seferis model is based on two distinct regions observed in the Avrami double log plot, which are attributed to primary and secondary crystallization [28]. It describes how these processes can occur in either parallel or series, giving rise to two separate equations (Equations (8) and (9), respectively) based on the standard Avrami model and apportioned to the fraction of each process that has occurred at time t. The parameters of these equations can only be determined via curve fitting.
Parallel : Series : where X t and X ∞ are the fractional crystallinities at time t and infinite time, respectively; w p and w s are the weight factors for the primary and secondary processes, respectively, where the sum of these two factors equals 1; k p and k s are the rate constants of the primary and secondary processes, respectively; and n vs and m vs are the exponents for both processes.
Recently, a new model, the Hay model, has been proposed that assumes that both processes occur concurrently and that the total crystallinity at time t is the sum of the two individual contributions. The primary crystallization mechanism follows the standard Avrami equation, whereas the secondary crystallization is based on a square root dependence of fractional crystallinity, as reported previously [2]. This has led to the following Hay model (Eq. (10)) [36,37] and, more recently, a new method to enable this theory to be applied to isothermal studies performed by DSC [27]. The parameters in this equation can be determined through a series of steps involving separation of the primary and secondary crystallization processes and evaluation of the primary data by the Avrami model in isolation [27].
where X t and X p,∞ are the fractional crystallinities at time t and at the end of the primary process, respectively; k p and k s are the Avrami rate constants of the primary and secondary processes, respectively; and n a is the Avrami exponent.
In this paper, each of the models discussed is applied to PHB-co-HV, which is known to display significant secondary crystallization. The possible parameters are determined through traditional double log plots or calculated from the Avrami parameters, and curve fitting procedures are applied. A range of methods are applied to compare the models and determine which are best suited to polymers displaying significant degrees of secondary crystallization.

Experimental procedure
This paper applies the data obtained by Kelly et al. [27] to a range of crystallization kinetic models. For convenience, the experimental procedure used to obtain these data is reported briefly here.
Plaques of PHB-co-HV were generated via compression molding using a Moore E1127 hydraulic hot press (George E. Moore & Sons Ltd, Birmingham). PHB-co-HV (8 g) was placed into a mold (152 × 158 × 0.266 mm) and inserted into the press at 190°C. The polymer was left to melt for 5 min before a load of 10 tons was applied for an additional 3 min. The plaques were allowed to cool in the press for 20 min.
A Mettler Toledo DSC1 differential scanning calorimeter (Mettler Toledo, Schwerzenbach, Switzerland) was used to determine the crystallization kinetics of PHB-co-HV across a range of temperatures (138-146°C) under nitrogen. This temperature range was limited by the rapid crystallization kinetics at lower temperatures and the propensity for degradation during prolonged periods at higher temperatures [34]. Discs (Ø 5 mm) were punched out from the plaques, and two were analyzed for each test to give a mass of approximately 10 mg. The discs were placed into 40 µL aluminum pans and capped with aluminum lids (Mettler Toledo) before being sealed with a press. Samples were held at 210°C for 2 min under nitrogen to melt the polymer before being cooled to the crystallization temperature at a rate of 30°Cmin −1 . The samples were held at this temperature for approximately six half-lives, as estimated from the minima in the crystallization exotherm. This timescale was chosen to obtain enough data from the secondary crystallization-only region while minimizing the degradation of the polymer due to prolonged heating at elevated temperatures. Initial studies indicated the presence of degradation when PHB-co-HV was heated at 146°C for in excess of 500 min (data not shown). Three samples were analyzed for each temperature (labeled A, B, and C). A heat of fusion of 140 Jg −1 was used to calculate the fractional crystallinity, as in the absence of a published value for PHB-co-HV, this value has been adopted by many others when analyzing copolymers containing low hydroxyvalerate contents [39][40][41].

Evaluation of crystallization kinetic models
The data, following baseline and induction time corrections, obtained by Kelly et al. [27] was used to analyze each of the crystallization kinetic models. For the Avrami, simplified Hillier, Tobin, and Hay models, the parameters were determined via the traditional method of creating double log plots, as discussed in more detail above. In addition, the Makin parameters were calculated from the Avrami results (Eqs. (5) and (6)). These parameters were then applied to the respective models to generate fraction crystallinity versus time plots, and the results were compared to the experimental data using SPSS software (IBM, Portsmouth, UK). The visual fit for the curves was evaluated along with the coefficient of determination (R 2 ) values and standard error of regression.
SPSS software was also used to perform direct curve fittings for each of the models using nonlinear regression methods and to calculate the respective kinetic parameters. These parameters were again used to generate cumulative crystallization curves, and the results were evaluated in the same way.
To enable visual comparisons between the models, all cumulative data were converted to the actual fractional crystallinity (Equation (11)) instead of the traditional method of normalizing data to the final crystallinity.
Results and discussion The Avrami, simplified Hillier, Tobin, and Hay parameters were determined through the traditional method of generating double log plots, and SPSS software was used to determine the fit of these parameters to the experimental data. In addition, the Malkin parameters were calculated directly from the Avrami results (Eqs. (5) and (6)). Figure 1 illustrates the fit of these models to the crystallization kinetics for PHB-co-HV at 138°C. Similar results were obtained per repeat measurement (A, B, and C) and temperature (data not shown). The results highlight that although the models are able to sufficiently describe the primary crystallization, there is an inability to accurately model the region dominated by the secondary crystallization process (X t > 0.6, as determined by the Hay model). The exception to this is the Hay model, which is designed to include contributions from both primary and secondary crystallization throughout the entire curve.
Additionally, SPSS software was used to directly curve fit the data to the respective models using nonlinear, multivariable regression analysis (Fig. 2). The results show that the Hay and Velisaris-Seferis models are most effective in describing the crystallization kinetics of PHB-co-HV. This is not surprising, as these theories have been developed to consider both primary and secondary crystallization, occurring concurrently or in series, and, therefore, more accurately represent the experimental data. Again, similar observations were found for each temperature analyzed (data not shown). It is also important to note that when curve fitting is used, the accuracy in the modeling of the primary crystallization region is reduced for those models that do not consider the secondary process, most notably the Avrami, simplified Hillier, and Malkin models. This is because the curve fitting process aims to minimize the level of error throughout the entire curve.
A more sophisticated approach to compare these models is to consider the coefficient of determination (R 2 ) values (Supplementary information S1-8). Each of the models generated R 2 values in excess of 0.9, showing a good fit to the experimental data. In addition, the results obtained from direct curve fitting gave higher R 2 values due to the whole curve being considered; however, this also leads to a deterioration in the quality of the fit for the primary crystallization region (Figs. 1 and 2). This is especially prevalent for the Avrami and Malkin models, which are based on primary crystallization alone. The Urbanovici-Segal, Velisaris-Seferis, and Hay models all generated high R 2 values in excess of 0.997, with the Hay model proving the most effective returning R 2 values of 1 for every sample when applying the curve fitting procedure.
As shown in Fig. 2, a high R 2 value does not necessarily result in a more accurate representation of the experimental data, especially for the case of models based on primary crystallization alone. It is, therefore, important to also consider individual regions of the traces. Four areas of interest from the cumulative curve were analyzed using the parameters obtained from the double log plot and curve fitting procedure, and the standard error of the regression, s, was calculated (Eq. (12)). These regions correspond to the initial stages of crystallization, the primary dominant region, the transition from primary to secondary crystallization, and the secondary crystallization region.
where s is the standard error of the regression, X t is the experimental fractional crystallinity, X t ′ is the modeled fractional crystallinity and n is the number of data points. For the purpose of comparing these models, the region corresponding to the initial stages of crystallinity was defined to be below 10% crystallinity (0 < X t < 0.1). The parameters derived from the double log plots generated a   (Table 1, Supplementary information S9). When applying the curve fitting procedure, SPSS software takes into consideration the entire plot, and is, therefore, unable to effectively process those models that focus mainly on primary crystallization. This is evident when comparing the difference in the errors between the curves calculated from the double log parameters and those from the curve fitting procedure for the Avrami and Hay models. The Avrami, simplified Hillier and Hay models, using parameters calculated from the double log plots, were found to provide the closest fit to the experimental data.
The region where primary crystallization dominates (0.25 < X t < 0.45) is taken around the inflection point of the crystallization curve, as this is where the rate of primary crystallization is greatest. In this region, the curves derived from the double log plots again proved to be the most effective (Table 1, Supplementary information S10). The simplified Hillier, using the calculated parameters, and both versions of the Hay method were found to be most suitable for modeling the primary crystallization with the standard error for the regression generally calculated to be less than 0.005. Both of these methods dictate the completion of primary crystallization and analyze this portion independently, which is perhaps why these models are the most successful.
The majority of the models were unable to successfully describe the transition from primary to secondary crystallization, with greater standard errors observed compared to the other regions (Table 1, Supplementary information S11). In general, results obtained using parameters calculated from the double log plots gave greater standard errors of regression than curve fitting due to the transitional region being rarely considered within models. As the Tobin, Urbanovici-Segal, Velisaris-Seferis, and Hay models account for the secondary crystallization process, they were more successful than those solely based on primary crystallization, with the SPSS curve fitting of the Hay model proving to be significantly superior.
Only the parallel form of the Velisaris-Seferis and the Hay equations could be used to successfully model the secondary crystallization portion of the curves (Table 1, Supplementary information S12) with standard errors of regression of 0.005 or less obtained. This is to be expected, as these are the only models that account for the presence of both primary and secondary crystallization.
Interestingly, the parallel form of the Velisaris-Seferis model outperformed the series method over the majority of  DLP curve produced from parameters calculated through double log plots, CF Curve produced via curve fitting by SPSS Fig. 3 Primary and secondary crystallization plots extracted from the cumulative curve. "Total" denotes the experimental data, and "Primary" and "Secondary" denote the calculated curves the selected areas. This confirms the Hay theory that both primary and secondary crystallization occur concurrently [27] and provides an explanation as to why the Hay model significantly improves on the other models when studying across the entire crystallization curve. Manipulation of the Hay equation can also enable the primary and secondary aspects of the cumulative crystallization curve to be plotted in isolation, highlighting the parallel nature of the two processes (Fig. 3).

Assessment of the parameters
One assumption in the Avrami equation is that the n value should be an integer that directly relates to the nucleation and growth mechanisms occurring during primary crystallization. It has been well documented that integer values are rarely observed [8,14], and this has been found in this work with both the Avrami and simplified Hillier approaches (Supplementary information S1 and S2). When curve fitting was used to calculate the Avrami exponent, the values deviated even further from an integer number as a result of the inability to consider the secondary crystallization region. The Hay model generated an integer value of 3.0 for all of the samples analyzed, but again, the n value was reduced by up to 0.2 when direct curve fitting was used (Supplementary information S8). Unlike the Avrami exponent, the Tobin (n t ) (Supplementary information S3) and Urbanovici-Segal (n us ) (Supplementary Information S5) exponents do not need to be integer values. The Malkin parameter C 0 can be calculated from the Avrami exponent, and, therefore, should also be an integer; in this case, 60. Although the results are more precise, an integer value was not obtained (Supplementary information S4). In addition, the curve fitting procedure for the Malkin model gives considerably lower C 0 values, proving that this method is unsuccessful (Supplementary information S4). As the Velisaris-Seferis method is based on the Avrami equations, n vs should also be an integer, which is again not observed from the curve fitting procedure (Supplementary information S6 and S7). In addition, a wide variation in the exponential values is found between the repeat measurements and different temperatures, highlighting that high R 2 values in the curve fitting procedure do not guarantee the calculation of the correct values for the parameters.

Analysis of the primary region alone
Evaluation of the models thus far has focused on a polymer that shows significant secondary crystallization. As most of the models assessed have been developed for primary crystallization, it is not surprising that they do not adequately describe the entire crystallization process effectively. The same crystallization-time curves were all reassessed using those models that emphasize primary crystallization (Avrami, simplified Hillier, Tobin, Malkin, and Urbanovici-Segal), with the data limited to this region alone. Modeling of the Avrami equation for a nominal range of time points shows that for an 'n' value of three, the rate of primary crystallization minimizes at a time equivalent to double the half-life of crystallization [27]. If the 'n' value is two, this increases to triple the half-life. Therefore, if the data are restricted to a time period of triple the half-life (as estimated from the minima of the exotherm), one can be confident that it mainly contains primary crystallization. The majority of the results from this analysis are presented in the supplementary information (S13-20). As shown by the cumulative curves created by the models, this restriction of the data improves the effectiveness of the models (Fig. 4). This is reflected in an increase in the R 2 values for each of the models (Supplementary information S13-16), with the exception of the Urbanovici-Segal model (Supplementary information S17). In addition, the curve fitting procedure produced R 2 values ranging between 0.999 and 1 for the simplified Hillier and Tobin models across the majority of samples analyzed. In most cases, the Avrami exponent was also found to be approximately equal to an integer value of three. Despite the improvement in the models when only considering the primary portion of the curve, the Hay model still gave the highest R 2 values using both calculated parameters and curve fitting. This again highlights the concurrent nature of primary and secondary crystallization and the requirement to analyze both processes in isolation.
Comparison of the standard error of regression for defined areas of the crystallization curve also pointed to an improvement when only the primary region was considered (Supplementary information S18-20). A marginal improvement in the region corresponding to the initial stages of crystallization was found for those curves produced from the parameters calculated by the double log plot; however, a large reduction in the standard error of the regression was found from the curve fitting procedure. This is because the fitting procedure is no longer distorted by a large degree of secondary crystallization that is not accounted for by the models. The greatest improvement was found where the rate of primary crystallization is at a maximum. With the exception of the Urbanovici-Segal model, large reductions to the standard error of regression were observed from both the calculated and curve fitting procedures compared to the full crystallization curve. The improvement in the transition region is inconsistent across each of these models, as this process is not considered by those models focusing mainly on primary crystallization.
Restricting the data to consider only the region where the primary mechanism dominates improves the ability of Avrami, simplified Hillier, Tobin, and Malkin models to model the crystallization kinetics of PHB-co-HV with high R 2 values, and, where applicable, n a values close to three are obtained (Supplementary information S13-16). However, as the above results show, the models are still unable to account for the entire curve; in particular, the transition from primary to secondary crystallization. Comparison across all the studies performed shows that the Hay model is the most successful at describing the crystallization kinetics of PHB-co-HV and most likely any polymer containing a high degree of secondary crystallization.

Conclusions
This work highlights the importance of selecting the optimum model when analyzing the crystallization kinetics of polymers, especially those displaying high levels of secondary crystallization.
The isothermal crystallization of PHB-co-HV was evaluated using a range of models, namely, Avrami, simplified Hillier, Tobin, Malkin, Urbanovici-Segal, Velisaris-Seferis, and Hay. Two methods of model evaluation were used: determination of the parameters through traditional double log plots and nonlinear, multivariable regression to perform curve fittings. The results show that although the curve fitting method produces higher R 2 values, modeling of the primary region is less accurate than when the parameters are calculated from double log plots. The exception to this is the Hay model, which is designed to account for both primary and secondary crystallization occurring simultaneously.
The Hay and Velisaris-Seferis models were found to be the most successful at modeling the entire crystallization curve as a result of incorporating the secondary crystallization process. In addition, the parallel form of the Velisaris-Seferis model outperforms the series model, confirming the Hay theory that primary and secondary crystallization occur concurrently. When assessing the parameters returned by the models, the Avrami exponent was found to have a noninteger value for the Velisaris-Seferis equations and varied vastly between the samples. The Avrami exponent was calculated via a double log plot of the Hay model to be 3. This work shows that the use of the Hay model, and in particular calculation of model parameters through the double log plot, to be the optimum method for analyzing the crystallization kinetics of polymers displaying significant secondary crystallization.

Compliance with ethical standards
Conflict of interest The authors declare no competing interests.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/. Fig. 4 Comparison of the Avrami, Simplified Hillier, Tobin, Malkin, and Urbanovici-Segal models with experimental data in the primary crystallization region: a Curves generated by imputing the crystallization kinetic parameters calculated from double log plots (b) Nonlinear multivariable curve fitting of the kinetic models. Solid black line corresponds to the experimental data and the model curves are as follows: Avrami -dark blue dots; Simplified Hillier -red short dashes; Tobin -green medium dashes; Malkin -pink short dash and dot; Urbanovici-Segal -light green long dashes