1. Introduction
The nucleation behavior in a supersaturated solution is closely related to the induction time or metastable zone width (MSZW) measurements [
1,
2,
3]. As opposed to a lag time during the temperature-decreasing process for the prepared supersaturated solution at a higher temperature lowered to the desired constant temperature in the induction time measurements, there is no lag time for the supersaturated solution cooled at a constant cooling rate in the MSZW measurements. Thus, the MSZW data should be more reliable than the induction time data in determination of the nucleation rate for a system.
In classical nucleation theory (CNT) [
1,
2,
3], the nucleation rate of a crystallization system depends on both the pre-exponential factor and interfacial energy, which are usually determined using induction time data by assuming
, where
is the nucleation rate and
is the induction time [
4,
5,
6,
7,
8,
9,
10]. Due to the complicated data interpretation method, determination of the pre-exponential factor and interfacial energy using MSZW data has long been a challenging task.
Based on the Nyvlt’s approach [
11,
12], Sangwal [
13,
14,
15,
16] proposed a self-consistent Nyvlt-like equation using the power-law nucleation rate, leading to a linear relationship between
and
, where
is the initial saturated temperature,
is the MSZW, and
is the cooling rate. Similarly, using the nucleation rate based on CNT, Sangwal [
13,
14,
15,
16] developed a linear relationship of
versus
. Kashchiev et al. [
17] presented a general expression for the total volume and number of crystallites as functions of the cooling rate for progressive nucleation based on CNT. They verified that the linear dependence of
on
is the linear approximation to their presented equation. For instantaneous nucleation, a more simplified model is derived by Kashchiev et al. [
18] which also yields a linear relationship between
and
. Kubota [
19] proposed a model based on progressive nucleation to account for the MSZW limit corresponding to a point at which the number density of accumulated crystals has reached a fixed (but unknown) value during the cooling process. The simple power-law form of nucleation rate is adopted in the Kubota’s model, leading to a linear relationship of
versus
for the MSZW data.
The above-mentioned models have been widely applied in the literature to correlate MSZW with cooling rate in various crystallization systems [
20,
21,
22,
23,
24,
25,
26,
27]. Although the pre-exponential factor and interfacial energy can provide important information in understanding the nucleation behavior, the pre-exponential factor and interfacial energy of the crystallized substance are usually not determined in these studies. Recently, some researchers [
28,
29,
30] investigated the relationship between MSZW and process parameters in terms of the pre-exponential factor and interfacial energy using classical 3D nucleation theory approach based on the Sangwal’s theory [
13,
14,
15,
16]. Shiau and Lu [
31,
32,
33] developed an integral model to determine the pre-exponential factor and interfacial energy based on CNT using the MSZW data. As the nonlinear regression along with numerical integration involved in the integral model is complicated, the objective of this work is to present a linear regression model to determine the pre-exponential factor and interfacial energy of the crystallized substance based on CNT using the MSZW data.
2. Theory
2.1. Integral Method
The nucleation rate based on CNT is expressed as [
1,
2,
3]
where
is the nucleation pre-exponential factor,
is the interfacial energy,
is the Boltzmann constant, and
is the molecular volume.
To interpret the onset of nucleation based on the progressive nucleation observed in the induction time and MSZW measurements, the nucleation event is assumed corresponding to a point at which the total number density of the nuclei has reached a fixed (but unknown) value,
[
19,
30,
31,
32]. Thus, one can derive
at the induction time, which is consistent with
reported in the literature [
1]. As nuclei are continuously born during the cooling process, Shiau and Lu [
31,
32,
33] derived at the MSZW limit:
where
depends on the measurement device and on the substance.
If the solution is cooled at a constant rate defined as:
one obtains
during the cooling process. Substituting Equations (1) and (3) into Equation (2) yields [
31,
32,
33]:
where, as shown in
Figure 1,
is the initial saturated temperature at
,
is the maximum supercooling temperature at
,
is the MSZW,
is the saturated concentration at
,
is the temperature-dependent supersaturation during the cooling process, and
is the temperature-dependent solubility. Note that
remains unchanged before the onset of nucleation. As
usually decreases with decreasing temperature,
increases with decreasing temperature; and subsequently
defined in Equation (1) increases during the cooling process.
Shiau and Lu [
31,
32,
33] proposed the following nonlinear regression procedure to determine two parameters,
and
, in Equation (4) from the experimental data of
versus
: (a) Guess
; (b) determine
for each pair of
versus
data by numerical integration; (c) calculate
; (d) calculate the coefficient of variation
among all
. For the guessed
,
is defined as:
where
is the number of data points
. For the guessed
,
among all
is defined as:
The same procedure from (a) to (d) is repeated by guessing various values of
until the optimal
with the minimum
is found. Then, the corresponding
calculated from the optimal
is taken as the optimal
. If
is known,
can be determined.
2.2. Linearized Integral Method
As the nonlinear regression along with numerical integration involved in Equation (4) is complicated, a simplified linear regression model is proposed in this study to extend the applicability of integral model as follows. Based on the trapezoidal rule for the numerical integration, Equation (2) is approximated as:
Here, represents the total time required during the cooling process from to at a constant cooling rate . and represent the nucleation rate at and , respectively. Note that due to .
As shown in
Figure 1,
increases gradually from 1 at
during the cooling process. As defined in Equation (1),
starts from
at
and increases gradually as temperature decreases from
to
. The nucleation rate
at
can be expressed as:
with
where
is the saturated concentration at
and
is the supersaturation at
.
Substituting Equation (8) into Equation (7) yields:
Taking logarithm on both sides of Equation (10) gives:
By rearranging Equation (11), the linearized integral method I is expressed as:
A plot of versus at a given should give a straight line, the slope and intercept of which permit determination of and , respectively, without the knowledge of .
If the temperature-dependent solubility is described in terms of the van’t Hoff equation [
1], one obtains:
where
is the heat of dissolution and
is the gas constant. By substituting Equation (13) into Equation (12), linearized integral method II is expressed as:
Similarly, a plot of versus at a given should give a straight line, the slope and intercept of which permit determination of and , respectively, without the knowledge of .
It should be noted in the application of Equations (12) and (14) that
is not influenced by the chosen value of
although
needs to be determined based on
. Based on the study of 28 inorganic systems, Mersmann and Bartosch [
34] concluded that the minimum detectable volume fraction of nuclei in solution corresponds to
with the minimum detectable size of
. The intermediate value,
, was adopted at the detection of the nucleation point for the Lasentec focus beam reflectance measurements reported by Lindenberg and Mazzotti [
35] and for the turbidity measurements reported by Shiau et al. [
31,
32,
33,
36]. If the uniform-sized spherical nuclei of
with
are assumed for simplicity,
corresponds to
[
32].
2.3. Classical 3D Nucleation Theory Approach
Sangwal [
13,
14,
15,
16] related
with the rate of change of solution supersaturation as:
where
,
, and
is a constant defined as the number of clusters per unit volume. The value of
is governed by the aggregation and diffusion processes in the solution.
If
is slightly greater than 1, one obtains:
Combining Equations (13) and (16) yields:
Substituting Equation (9) into Equation (17) gives:
By combining Equations (8), (15) and (18), Sangwal [
13,
14,
15,
16] derived the classical 3D nucleation theory approach as:
where:
A plot of versus at a given should give a straight line, the slope and intercept of which permit determination of and , respectively. Subsequently, the values of and can be calculated without the knowledge of .
It should be noted in the application of Equation (19) that
is not influenced by the chosen value of
although
needs to be determined based on
. Sangwal [
14] proposed that the upper limit of
may be estimated from solute concentration in the saturated solution. For example,
, the saturated concentration for aqueous glycine solutions at 308.15 K (
), corresponds to
based on the number of solute molecules per unit solution volume, which is close to
proposed by Sangwal [
14].
3. Results and Discussion
The experimental MSZW data for three crystallization systems, including glutamic acid, glycine and, 3-nito-1,2,4-triazol-5-one (NTO), reported in the literature are analyzed as follows.
Figure 2a shows that the MSZW data of aqueous glutamic acid solutions fitted to linearized integral method I Equation (12) at various
, where the original experimental MSZW data listed in
Table 1 are taken from Shiau and Lu [
32] in a 200 mL vessel. The solubility of glutamic acid in water is taken as
;
[
37]. Note that
for glutamic acid.
Figure 2b shows the same MSZW data of aqueous glutamic acid solutions fitted to linearized integral method II Equation (14) at various
, where
is used for the van’t Hoff solubility equation. For comparison,
Figure 2c,d show the same MSZW data of aqueous glutamic acid solutions fitted to the classical 3D approach Equation (19) and integral method Equation (4), respectively, at various
.
Figure 3a shows that the MSZW data of aqueous glycine solutions fitted to linearized integral method I Equation (12) at various
, where the original experimental MSZW data listed in
Table 2 are taken from Shiau [
38] in a 200 mL vessel. The solubility of glycine in water is taken as
;
[
39]. Note that
for glycine.
Figure 3b shows the same MSZW data of aqueous glycine solutions fitted to linearized integral method II Equation (14) at various
, where
is used for the van’t Hoff solubility equation. For comparison,
Figure 3c,d show the same MSZW data of aqueous glycine solutions fitted to classical 3D approach Equation (19) and integral method Equation (4), respectively, at various
.
Figure 4a shows that the MSZW data of aqueous NTO solutions fitted to linearized integral method I Equation (12) at various
, where the original experimental MSZW data listed in
Table 3 are taken from Kim et al. [
40] in a 300 mL vessel. The solubility of NTO in water is taken as
;
[
40]. Note that
for NTO.
Figure 4b shows the same MSZW data of aqueous NTO solutions fitted to linearized integral method II Equation (14) at various
, where
is used for the van’t Hoff solubility equation. For comparison,
Figure 4c,d show the same MSZW data of aqueous NTO solutions fitted to the classical 3D approach Equation (19) and integral method Equation (4), respectively, at various
.
The fitted results for glutamic acid, glycine, and NTO are listed in
Table 4,
Table 5 and
Table 6. As integral method Equation (4) is numerically integrated without any approximations in this study,
and
obtained from integral method Equation (4) represent the exact solution to Equation (2) based on the nucleation event assumed corresponding to a point at which the total number density of the nuclei has reached
. For all three studied systems, as opposed to
obtained from classical 3D approach Equation (19) compared with that obtained from integral method Equation (4),
obtained from linearized integral methods Equations (12) and (14) is closer to that obtained from integral method Equation (4) at each condition. Furthermore,
obtained from linearized integral methods Equations (12) and (14) is also consistent with that obtained from integral method Equation (4) at each condition. As
and
are two different parameters,
determined from Equations (4), (12) and (14) based on
is not strictly comparable with that determined from Equation (19) based on
. It should be noted in
Table 4,
Table 5 and
Table 6 that
and
(or
) are determined first without the knowledge of
(or
). Consequently,
is not influenced by the chosen value of
(or
) although
needs to be determined based on
(or
). For example, if the chosen value of
(or
) is increased by ten times,
is also increased by ten times at each condition while
remains unchanged.
As compared in
Table 4,
Table 5 and
Table 6,
and
obtained from linearized integral method I Equation (12) are consistent with those obtained from linearized integral method II Equation (14) at each condition for all three studied systems. Thus, both equations can be applied to determine
and
of the crystallized substance using the MSZW data. As opposed to the temperature-dependent solubility required for linearized integral method I Equation (12), only the value of
is required in application of linearized integral method II Equation (14).
4. Conclusions
A linear regression method is proposed in this work to determine the pre-exponential factor and interfacial energy based on CNT using the MSZW data. Linearized integral method I Equation (12) is derived for any temperature-dependent functional form of the solubility while linearized integral method II Equation (14) is derived for the van’t Hoff temperature-dependent solubility. Only the value of is required in application of linearized integral method II Equation (14), as opposed to the temperature-dependent solubility required for linearized integral method I Equation (12). The experimental MSZW data for all three studied systems, including glutamic acid, glycine, and NTO, are fitted well to these two equations. The pre-exponential factor and interfacial energy obtained from linearized integral method I Equation (12) are consistent with those obtained from linearized integral method II Equation (14) for these systems.
As the integral method is numerically integrated without any approximations, the pre-exponential factor and interfacial energy obtained from the integral method represent the exact values based on the nucleation event assumed corresponding to a point at which the total number density of the nuclei has reached a fixed value. As opposed to the interfacial energy obtained from classical 3D nucleation theory approach compared with that from the integral method, the interfacial energy obtained from linearized integral methods Equations (12) and (14) is closer to that from the integral method at each condition. Furthermore, the pre-exponential factor obtained from linearized integral method Equations (12) and (14) is also consistent with that from the integral method at each condition.