1. Introduction
According to classical nucleation theory (CNT), only nuclei greater than a critical nucleus size are thermodynamically stable and can continue to grow to a detectable size [
1,
2,
3]. The formation of critical nuclei is closely related to the interfacial energy of the crystallized substance, which is usually calculated from the induction time data in the literature [
4,
5,
6,
7,
8,
9,
10].
The induction time is defined as the elapsed time between the creation of the supersaturation and the appearance of detectable nuclei at a constant temperature. Although the induction time can be detected by visual observation of the crystal’s appearance [
7,
11], turbidity measurements have been commonly adopted in recent years to determine the induction time by detecting the change in the intensity of transmitted light in solution at the onset of nucleation [
12,
13,
14,
15,
16,
17]. Traditionally, determination of the interfacial energy from the induction time data is often simplified by assuming
[
1,
4,
5,
6,
7,
8,
9,
10]. Thus, it is implicitly assumed that at the detection of the nucleation point, only the number of the nuclei is accounted for regarding the change in the intensity of transmitted light in solution.
The detection of nucleation point based on turbidity measurements should be influenced by both the number and the size of the nuclei [
18] as the change in the intensity of transmitted light in solution is proportional to the size distribution of the nuclei instead of the number of the nuclei. To incorporate the effect of the nuclei size distribution on the detection of nucleation, Shiau and coworkers [
18,
19] have developed a model to examine the turbidity induction time data of aqueous L-glutamic acid solutions using the L-glutamic acid growth kinetics reported by Scholl et al. [
20]. It is found that the obtained interfacial energy and growth activation energy of L-glutamic acid [
19] are consistent with the literature data. L-glycine is the simplest amino acid and is often used as a model compound in the study of solution nucleation [
21,
22,
23,
24,
25]. The objective of this work is to develop a model to study the nucleation and growth of aqueous L-glycine solutions based on the turbidity induction time data. To validate the obtained L-glycine growth kinetics from this model, the photomicroscopic growth experiments of aqueous L-glycine solutions are also conducted to determine the growth kinetics of nuclei at the same conditions for comparison.
2. Theory
The nucleation rate based on CNT [
1,
2,
3] is expressed as:
where
and
. For simplicity, the nucleation event is assumed to correspond to a point at which the total number density of the nuclei has reached a fixed (but unknown) value,
fN [
26,
27]. One obtains at the induction time
ti using:
Thus, it is implicitly assumed that the detection of the nucleation point is related to the number of nuclei. Substituting Equation (1) into Equation (2) yields:
This is consistent with the common method adopted in the literature to calculate
γ from induction time data [
1,
4,
5,
6,
7,
8,
9,
10].
The turbidity induction time measurements are based on the change in intensity of transmitted or scattered light along the detector direction, which should be related to the size distribution of the nuclei instead of the number of the nuclei at the detection point [
18,
19]. As nuclei are progressively generated during the induction time period (
t = 0–
ti), the nuclei born in the earlier stage will grow to a greater size than those born in the later stage at
. To incorporate the effect of crystal growth at the nucleation point, Shiau and Lu [
18] proposed a model to correlate the nucleation and growth with the turbidity induction time data using the predetermined growth kinetics. However, as it is often difficult to experimentally measure the growth kinetics of small nuclei, the application of this model is restricted.
As the turbidity induction time measurements are influenced by the size distribution of the nuclei, these data should provide important information on both nucleation and growth. A model is developed in the following to investigate the nucleation and growth of nuclei based on the turbidity induction time data without the predetermined growth kinetics. In the derivation, a simple empirical power-law growth rate is proposed as:
where the value of
g mostly falls between 1 and 2. Based on Burton–Cabrera–Frank (BCF) growth theory [
28,
29,
30], the value of
g is found close to 2 for low supersaturations [
31]. Mohan and Myerson [
32] indicated for aqueous L-glycine solutions at 293.15 K that Equation (4) with
g = 2 is consistent with the BCF growth kinetics reported by Li and Rodriguez-Hornedo [
33].
In the induction time study, nuclei born at any time
t (0 <
t <
ti) can grow from
t to
ti and their size at time
ti is:
The corresponding volume of nuclei with size
L at time
ti is
As nuclei are progressively generated from
t = 0 to
ti, the total volume of all the nuclei per unit solution volume at time
ti is given by:
Substituting Equation (6) into Equation (7) yields:
Note that J and G remain unchanged as S is kept at a particular supersaturation during each induction time experiment.
Substituting Equations (1) and (4) into Equation (8) with
g = 2 leads to:
Rearranging Equation (9) yields:
A plot of versus at a given temperature should give a straight line, the slope and intercept of which permit determination of γ and , respectively.
The temperature dependence of
kG can be expressed in terms of the Arrhenius equation as:
Once
is determined at different temperatures, substitution of Equation (11) yields:
Thus, a plot of versus should give a straight line, the slope and intercept of which permit determination of EG and , respectively. It should be noted that AJAG3 can be determined if fV is known.
4. Results and Discussion
The induction time data for aqueous L-glycine solutions were measured for various supersaturations from 293.15 K to 313.15 K. Each run was carried out at least three times to determine the average induction time under each condition. The average induction times and the corresponding standard deviations (SD) are listed in
Table 1. The equilibrium concentration of the α-form L-glycine in water was given by
Ceq(
T) = 5.4397 × 10
−3T2 − 3.2022 × 10
−1T − 188.2 (
Ceq in kg/m
3, and
T in K) [
37]. Note that
MW = 0.075 kg/mol,
ρC = 1607 kg/m
3, and
ν = 7.757 × 10
−29 m
3 for L-glycine.
Although L-glycine can be crystallized in different polymorphs, the α-form is usually achieved from nucleation of pure aqueous solutions [
21,
22,
23,
24,
25]. To validate the polymorphm of the L-glycine crystals, the final dried crystals at the end of the induction time experiments were analyzed using both optical microscopy and Raman spectroscopy (P/N LSI-DP2-785 Dimension-P2 System, 785 nm, manufactured by Lambda Solutions, INC., Seattle, WA, USA). The results all indicated that the needle-like α-form crystals were obtained from aqueous L-glycine solutions for various supersaturations from 293.15 K to 313.15 K.
Figure 4 shows the Raman spectra of pure α-form crystals and product crystals obtained at various supersaturations. As compared with the Raman spectra of pure α-form crystals reported by Murli et al. [
38], it was confirmed that α-form crystals were produced for various supersaturations at 303.15 K. The section of the Raman spectra of α-, β-, and γ-glycine used for characterization are also depicted by Bouchard et al. [
39].
Figure 5 shows the measured induction time data from 293.15 K to 313.15 K fitted to Equation (3) based on
fN. The fitted results are listed in
Table 2, where
γ was in the range 1.93–2.37 mJ/m
2 and
was in the range 4.35 × 10
−3–6.42 × 10
−3 s
−1. Although the exact value of
AJ could only be determined with a known value of
fN,
γ was not influenced by the chosen value of
fN.
As the induction time data are measured by the intensity change of the transmitted light,
Figure 6 shows the measured induction time data from 293.15 K to 313.15 K fitted to Equation (10) based on
fV. Note that
g was assumed to be 2 due to low supersaturations (
S = 1.4–2.4) in the induction time experiments. The fitted results are listed in
Table 3, where
γ was in the range 2.49–2.93 mJ/m
2 and
was in the range 2.78 × 10
−5–2.58 × 10
−3 s
−4. It should be noted that
γ was not influenced by the chosen value of
fV.
The turbidity induction time measurements in the current experiments were based on the intensity change of the transmitted light, which is related to
fV. Thus, compared to
γ based on
fN,
γ based on
fV should more accurately represent the actual interfacial energy of L-glycine. Shiau [
40] reported
γ = 1.35–2.02 mJ/m
2 for aqueous α-form L-glycine solutions using the turbidity metastable zone width measurements at the saturation temperature between 308.15 K and 328.15 K. Using the visual observation of the induction time data, Devi and Srinivasan [
25] reported
γ = 5 mJ/m
2 for aqueous α-form L-glycine solutions at 303.15 K.
Figure 7 shows the plot of
versus
fitted to Equation (12). The fitted results are listed in
Table 4, which indicates
EG = 58 kJ/mol and
= 2.30 × 10
26 s
−4. Although the exact value of
AJAG3 could only be determined with a known value of
fV,
EG was not influenced by the chosen value of
fV. Because activation energy is usually 10–20 kJ/mol for diffusion and 40–60 kJ/mol for surface integration [
1],
EG = 58 kJ/mol obtained for the growth of L-glycine in the induction time experiments should be integration controlled.
Based on the study of 28 inorganic systems, Mersmann and Bartosch [
41] estimated
fV = 10
−4–10
−3 with a detectable size of 10 μm. In the calculations here, the intermediate value,
fV = 4 × 10
−4, for spherical nuclei with
kV =
was assumed, leading to
fN = 7.64 × 10
11 m
−3. Consequently, as indicated in
Table 2,
AJ was in the range 3.32 × 10
9–4.91 × 10
9 m
−3·s
−1 based on
fN = 7.64 × 10
11 m
−3.
Table 4 indicates
AJAG3 = 9.18 × 10
22 s
−4 based on
fV = 4 × 10
−4. For simplicity, if
AJ obtained based on
fN was adopted to find
AG, one obtains
AG = 2.66 × 10
4–3.02 × 10
4 m/s.
The average growth rates and the corresponding standard deviations (SD) obtained from the photomicroscopic growth experiments are listed in
Table 5.
Figure 8 displays the growth rate data for various supersaturation from 303.15 K to 323.15 K. The growth rate obtained here is consistent with that reported by Han et al. [
23] for aqueous α-form L-glycine solutions at
S = 1.35 and
T = 303.15 K. Substituting Equation (11) into Equation (4) for
g = 2 yields:
Rearranging Equation (13) leads to:
As shown in
Figure 9, a plot of
versus
should give a straight line, leading to
EG = 57 kJ/mol and
AG = 6.05 × 10
4 m/s.
Figure 8 shows the growth rate data fitted well to Equation (14) using
g = 2 and the fitted values of
EG and
AG. Thus, it was reasonable to adopt the power-law growth rate of Equation (4) with
g = 2 in derivation of Equation (10). It should be noted that the turbidity induction time data were measured for nuclei of near-zero size (<10 μm, as assumed here for
fV = 4 × 10
−4 and
fN = 7.64 × 10
11 m
−3) in the 200 mL stirred solution while the photomicroscopic growth data were measured for nuclei of size
L = 50–100 μm in the 20 mL stagnant solution. Nevertheless,
EG obtained from the growth rate data was close to that obtained from the induction time data while
AG obtained from the growth rate data was still quite consistent with that obtained from the induction time data.