Next Article in Journal
Privacy-Constrained Biometric System for Non-Cooperative Users
Next Article in Special Issue
Estimating the Nucleation Ability of Various Surfaces Towards Isotactic Polypropylene via Light Intensity Induction Time Measurements
Previous Article in Journal
Multivariate Shortfall and Divergence Risk Statistics
Previous Article in Special Issue
The Criterion for the Crystallization Ability Assessment as Applied to Borate Glass Powders and Monoliths
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Critical Size of Secondary Nuclei Determined via Nucleation Theorem Reveals Selective Nucleation in Three-Component Co-Crystals

Department of Chemical Engineering, Advanced Materials Laboratory of Ministry of Education (MOE), Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(11), 1032; https://doi.org/10.3390/e21111032
Submission received: 7 September 2019 / Revised: 20 October 2019 / Accepted: 21 October 2019 / Published: 24 October 2019
(This article belongs to the Special Issue Crystallization Thermodynamics)

Abstract

:
The critical size of the secondary nuclei plays an important role in determining the crystal growth rate. In the past, the Nucleation Theorem has been applied to determine the number of molecules in the critical nuclei of a single-component crystal via variation of the crystal growth rate with dilution by the non-crystallizable component. In this work, we extend the method to the three-component co-crystal poly (ethylene oxide)/urea/thiourea inclusion compound. The theoretical crystal growth kinetics were deduced and the dependence of the radial growth rate of the inclusion compound spherulites on the mass fraction of urea in urea/thiourea was measured. The results reveal that the secondary nuclei of the poly (ethylene oxide)/urea/thiourea inclusion compound consist mainly of ethylene oxide repeating units and urea molecules. We propose that only urea molecules and ethylene oxide repeating units are selected to form the secondary nuclei while co-crystallization of the three components happens at the lateral spreading stage. As a result, the composition of the critical secondary nuclei is different from that of the bulk inclusion compound crystals. The work is expected to deepen our understanding of the nucleation of multi-component co-crystals.

1. Introduction

Nucleation is a central problem in the research field of crystallization [1,2,3]. According to whether the newly formed crystal nucleus is from melt or on the pre-existing surface of the same substance, nucleation can be classified as primary nucleation or secondary nucleation [4]. The latter usually occurs during the crystal growth process and determines the growth rate. If the crystal growth process adopts the surface nucleation mechanism, formation of the critical secondary nuclei with the highest Gibbs free energy acts as the rate-determining step during the crystal growth process [5,6,7,8,9,10,11].
The Nucleation Theorem based on the classical nucleation theory (CNT), proposed by Kashchiev in 1982, has been applied to determine the critical size of nuclei. Nucleation Theorem states that the derivative of nucleation work, with respect to the difference of the bulk Gibbs free energy of a crystallizable unit, equals the critical nucleus size with a negligible bias for the excess Gibbs free energy [12,13,14]. Its statistical thermodynamic foundation was proposed by Ford later [15,16]. The difference of the bulk Gibbs free energy of a crystallizable unit during the crystallization can be changed by different ratios of dilution with non-crystallizable molecules, which leads to the change of the entropy barrier of nucleation.
Following the above idea, we have adopted the Nucleation Theorem to measure the critical secondary nucleus size in a two-component co-crystal, poly (ethylene oxide) (PEO)/urea inclusion compound (IC) during its crystal growth process [17]. The PEO/urea ICs consist of guest PEO chains residing in the cavities formed by the host urea molecules [18,19,20,21,22,23,24,25,26]. By replacing a small amount of urea with N,N′-dimethyl urea with the PEO concentration unchanged, we could change the rate of secondary nucleation and determine the critical size of the secondary nuclei of the PEO/urea IC β crystal. Here, dimethylurea acted as a dilute agent, which was completely precluded from the crystal lattice. In the PEO/urea IC β crystal, the fixed stoichiometric ratio of 3/2 exists between ethylene oxide (EO) repeating units and urea molecules [22,23,24,27]; the crystallization unit is an entity containing 1 urea and 1.5 EO repeating units. Taking an entity as a unit of crystallization, the inclusion compound can be simplified as one-component crystal. A power-law relationship between the radial growth rate G and the mass fraction x of urea in the total mass of urea and the diluent was observed, and the exponent gave the critical size of the secondary nucleus. Our results revealed that a critical secondary nucleus of the PEO/urea IC β crystal consisted of 4 to 9 entities in the studied temperature range [17].
In this work, we attempt to extend the Nucleation Theorem to determine the critical size of the secondary nuclei in a three-component co-crystal system poly (ethylene oxide)/urea/thiourea inclusion compound. Thiourea is an analogue molecule of urea, and as a result, the three components co-crystallize. PEO/urea/thiourea IC prepared by electrospinning was reported by Ye et al. [28], but the crystal growth kinetics were not examined. Here we find that the spherulites of PEO/urea/thiourea IC can form by isothermal crystallization from the homogeneous PEO/urea/thiourea mixed melt, which makes the investigation of crystal growth kinetics possible.
This article is organized as follows. First, based on the Nucleation Theorem, we deduce the theoretical variation of the crystal growth kinetics of two-component co-crystal with the changed ratio of the two components. PEO/urea/thiourea IC with varied ratios of urea to thiourea molecules and fixed concentration of PEO can be considered as two-component co-crystals, so the PEO/urea/thiourea IC system can be treated with the deductions. Second, we report the crystal structure and the radial growth rates (G) of PEO/urea/thiourea IC spherulites with different content of urea to check the theory, and we obtain the critical nuclei size. To our surprise, the critical size of the secondary nucleus during crystal growth of the PEO/urea/thiourea IC is very similar to that of the PEO/urea IC. Finally, we will discuss the results of the crystal structure and the critical size of secondary nuclei. We propose that the critical secondary nuclei of the three-component crystal system consist of only two components, ethylene oxide segment and urea molecules, while all the three components are incorporated into the crystal lattices during the following crystal growth. Namely, the segregation of small molecules occurs only during the nucleation stage.

2. Theory

Considering a two-component co-crystal system consisting of crystallization units A and B, we try to deduce the crystal growth rate G varying with the composition. When the crystallization temperature is above the glass transition temperature and high enough, the secondary nucleation is the rate-determining step during the crystal growth. From the classical nucleation theory [5,6,7,29,30,31,32], the crystal growth rate G is determined by the following equation,
G = G 0 exp ( W * / k B T c )
In Equation (1), G0 is the pre-exponential factor, kB is the Boltzmann constant, Tc is the isothermal crystallization temperature, and W* is the formation work of the critical secondary nuclei.
We consider the critical secondary nuclei to be composed of components A and B. Δg1 and Δg2 are defined as the absolute value of the bulk Gibbs free energy difference caused by the crystallization of one A and B crystallization unit, respectively. The number of A and B crystallization unit in the critical secondary nucleus is termed as n1 and n2. The Gibbs free energy difference W(n1, Δg1, n2, Δg2) during the formation of one cluster containing n1 A units and n2 B units is expressed as follows,
W ( n 1 , Δ g 1 , n 2 , Δ g 2 ) = n 1 Δ g 1 n 2 Δ g 2 + F s ( n 1 , Δ g 1 , n 2 , Δ g 2 )
where Fs represents the excess Gibbs free energy during the cluster formation. If the value of Δg1 or Δg2 is invariable with cluster size, W will be a binary function of n1 and n2, and the work to form a critical secondary nucleus is equal to the maximum point W(n1*, Δg1, n2*, Δg2) of this binary function (we rewrite it as W*(Δg1, Δg2)), where n1* and n2* represents the number of the crystallization units A and B in the critical secondary nucleus. So, the following equation establishes,
W ( n 1 , Δ g 1 , n 2 , Δ g 2 ) n i n 1 = n 1 * ( Δ g 1 , Δ g 2 ) , n 2 = n 2 * ( Δ g 1 , Δ g 2 ) = 0 , i = 1 , 2 .
Therefore, combining the Equations (2) and (3), we have,
W * Δ g i = W Δ g i n 1 = n 1 * , n 2 = n 2 * = n i * + F s Δ g i n 1 = n 1 * , n 2 = n 2 * n i * , i = 1 , 2 .
The capillarity approximation is adopted, namely, the interfacial free energy (per area) does not vary with the change of Δg and the cluster size. This assumption in the classical nucleation theories has been generally accepted. Thus, the independence of Fs on Δgi leads to the first Nucleation Theorem described in Equation (4). The details of the deduction can be found in Kashchiev’s papers [12,13,14] and in the comments by Schmelzer [2]. Schmelzer has deduced more detailed deductions of the Nucleation Theorem and it generalized formulation, and in addition, the valid range and limits of the Nucleation Theorem are given in the comments [33]. The extension of the Nucleation Theorem to the multi-component system has also been derived by Schmelzer. In fact, Equation (4) can be derived from the Equation (27) in Schmelzer’s paper [33], assuming that the densities of the two components are the same.
The next step is to find the quantitative relation between Δgi and the composition in the system. We define x’ as the volume fraction of component A in the critical secondary nucleus consisting of n1* A molecules and n2* B molecules, and x’ is approximately equal to the molar fraction n1*/(n1* + n2*). Considering that the mixing of A and B in the melt is homogeneous at the molecular level, and A molecules and B molecules randomly distribute in the critical secondary nucleus, which brings the entropy of mixing to the critical secondary nucleus, and assuming that the enthalpy of mixing A and B is zero in both the melt and co-crystal nucleus (due to the similarity between the urea and thiourea molecules when they form co-crystals in this work), the equivalent thermodynamic process is shown in Scheme 1 to calculate the absolute value of Δg1, the Gibbs free energy change of one A molecule (or entity). Equation (5) is deduced from the equation ΔG1 + ΔG2 = ΔG3 + ΔG4 in Scheme 1, and a similar analysis can be used to deduce Equation (6).
Δ g 1 = Δ g 10 + k B T c ln ( x / x )
Δ g 2 = Δ g 20 + k B T c ln [ ( 1 x ) / ( 1 x ) ] .
In the above equations, Δg10 and Δg20 are the absolute values of the bulk Gibbs free energy difference caused by the crystallization of one A and B unit from the pure A and B melts, respectively. x represents the volume fraction of the major component A in the mixed melt of A and B. Equations (5) and (6) can be applied when x’ is neither equal to 0 nor equal to 1; otherwise, the equations are simplified to those for the one-component crystal system.
Combining Equations (1) and (4)–(6), and n1*/x’ = n2*/(1 − x’), we can calculate the fitting slope (termed as K) of lnG versus lnx as follows,
K = ln G ln x = ln G 0 ln x 1 k B T c W * ln x 1 k B T c W * ln x = 1 k B T c ( W * Δ g 1 Δ g 1 ln x + W * Δ g 2 Δ g 2 ln x ) = n 1 * ( ln x ln x ln x ln x ) + n 2 * ( ln ( 1 x ) ln x ln ( 1 x ) ln x ) = n 1 * n 1 * x x ln x x 1 x n 2 * + n 2 * 1 x x ln x = n 1 * x 1 x n 2 *
Equation (7) gives the relationship between lnG versus lnx expressed by the parameters of the critical secondary nucleus. The prefactor G0 may vary with the composition x. For the two-component system, we can first consider the following two extreme cases: if only one component (with the volume concentration noted as x) is involved in the secondary nucleation, G0 is proportional to x. If the two components are simultaneously involved in nucleation and the nucleation rates of the two components are the same, G0 does not change with x. The other cases lie between the above two extreme cases, so the partial derivative of lnG0 with respect to lnx ranges from 0 to 1. Thus, considering the possible variation of G0 with x, the difference between the slope K and the number of molecules (or entities) in the critical nuclei is less than 1, and we can omit this term. When n2* equals zero, Equation (7) is simplified to the case of the nucleation of a single-component crystal.
To apply Equation (7) to the three-component PEO/urea/thiourea IC system, we regard an entity consisting of 1 EO with 1.5 urea and that of 1 EO with 1.5 thiourea as the two types of crystallization units A and B, respectively. The stoichiometric ratio of EO units to small molecules in the PEO/urea inclusion compound is 3/2, and each unit cell averagely contains 1.5 EO and 1 urea molecule [22,23,24,27]. We use 1 small molecule and 1.5 ethylene oxide as an entity of crystal nucleation. Consequently, the three-component co-crystal PEO/urea/thiourea IC can be treated as a two-component A/B co-crystal. For simplicity, we approximate the volume ratio by mass ratio since urea and thiourea have similar densities. So, in the later section, the mass ratio of the PEO/small molecules is fixed to 1/0.91, which is the same as the EO/urea mass ratio in PEO/urea IC β crystal, and x represent the mass fraction of urea in the total mass of urea/thiourea.

3. Experimental Section

3.1. Sample Preparation

Poly (ethylene oxide) (PEO) with a viscosity-average molecular weight (Mv) of 100,000 g/moL and urea were both purchased from the Sigma-Aldrich Company. Thiourea was purchased from the Aladdin Company. All the materials were used without further purification.
The mass ratio of EO to the sum of urea and thiourea was fixed at 1/0.91, which is equal to the mass ratio of EO to urea in the PEO/urea IC β phase. The mass ratio of urea to thiourea was set from 10/0 to 0/10. Complete dissolution of PEO (0.2 g), urea, and thiourea (0.182 g of the latter two in total) in the deionized water (3 mL) was accelerated by heating and stirring. Then the water solution of PEO/urea/thiourea was rapidly cooled down and crystallized in liquid nitrogen to ensure the homogeneous composition in the obtained ice block. After being freeze-dried for 48 h to completely remove the water solvent, a white spongy solid block of homogenous PEO/urea/thiourea mixture was obtained.
The crystals of the PEO/urea/thiourea inclusion compounds with different compositions were prepared by isothermal crystallization at 70 °C, after being melted at 155 °C for 15 s, for the relevant characterizations.

3.2. Characterization

Each film of PEO/urea/thiourea IC crystal with a thickness of about 200 μm was prepared by sandwiching milligrams of the as-prepared bulk sample between a cover glass and a piece of polyimide film. The isothermally crystallized samples were obtained by quenching to 70 °C, and were held at the same temperature for enough time (usually less than 30 min) after melting at 155 °C for 15 s on a hot-stage. After cooling down to the room temperature, the cover polyimide film was carefully removed to obtain a bare, flat surface of the PEO/urea/thiourea IC crystals for one-dimensional X-ray diffraction (XRD) measurement, which was carried out on a diffractometer (Riguka SmartLab, Rigaku, Tokyo, Japan) with a scanning rate of 10°/min with 2θ step of 0.01° from 5 to 50° under the Ni-filtered Cu Kα radiation (40 kV, 30 mA).
The thermal behaviors of the PEO/urea/thiourea IC crystals were studied by differential scanning calorimeter (Shimadzu DSC-60, Shimadzu, Kyoto, Japan) under nitrogen atmosphere. An indium standard was used for calibration before the measurement. All the samples were melted at 150 °C for 1 min, and then quenched to 70 °C and isothermally crystallized completely for 30 min, followed by heating to 160 °C with a rate of 10 °C/min to measure the melting points and enthalpies of fusion.
To observe the spherulite morphology and to measure the radial growth rates, a polarized optical microscope (POM) (Olympus BX41P, Olympus, Tokyo, Japan) equipped with a CCD digital camera (Motic Moticam Pro 282A, Motic, Xiamen, China) was used. A sample sandwiched between two glass slides of 10 mm × 10 mm was first melted at 155 °C for 15 s to eliminate the thermal history, and then pressed to a film with the thickness around 10 μm. Each growth rate of spherulites was measured for three different samples and the averaged value of growth rate and standard deviation were calculated.

4. Results and Discussion

4.1. Crystalline Structure of the Peo/Urea/Thiourea Ternary Inclusion Compounds

The wide-angle X-ray diffractograms of a series of PEO/urea/thiourea IC samples are presented in Figure 1a. The mass ratio of PEO and small molecules (urea plus thiourea) is fixed to 1/0.91, with the mass ratio of urea to thiourea varying from 10/0 to 1/9. All the diffractograms appear almost the same, except that the positions of the main diffraction peaks shift to the low-angle direction steadily with the increasing content of thiourea molecules (Figure 1b). The results reveal that the crystal structure is similar in all the samples, namely, the crystal structure of all the samples are the same as the orthorhombic system of the PEO/urea IC β phase. The main diffraction peaks are indexed in Figure 1.
From the position shift of the diffraction peaks, we can estimate the expansion rate of the spacing of the diffraction planes in PEO/urea/thiourea IC crystals. The interplanar spacings of the different (hkl) planes, dhkl(x), with x indicating the mass fraction of urea molecules in urea/thiourea, are calculated based on the Bragg’s law in Equation (8). The percentage of expansion αhkl is defined by dividing the slope of the linear fitting of the expansion distance dhkl(x) versus x by the interplanar spacing of PEO/urea IC crystal without thiourea molecules dhkl(x = 1) (Figure 1c and Equation (9)). The expansion percentages of the different crystal faces range from 3% to 7% (Figure 1d). We conclude that the continuous change of the interplanar spacings is due to the gradual addition of the thiourea molecule, which suggests that the thiourea molecules replace the urea molecules stochastically and uniformly. In other words, PEO/urea/thiourea IC is a ternary isomorphic system.
d h k l ( x ) = λ 2 sin θ h k l ( x )
α h k l = 1 d h k l ( 1 ) d d x d h k l ( x ) .

4.2. Thermodynamic Stability of the PEO/Urea/Thiourea Ternary Inclusion Compounds

The DSC heating curves of the ternary isomorphic PEO/urea/thiourea IC completely crystallized at 70 °C are shown in Figure 2a. With the addition of the thiourea component, the melting point Tm of the PEO/urea/thiourea IC decreased first, and then increased (Figure 2b), which reflects the same tendency of the thermodynamic stability. The enthalpy of fusion of the PEO/urea/thiourea ternary isomorphic IC decreased first, and then it leveled off to about 100 J/g (Figure 2b). The obvious decrease of the enthalpy of fusion in the sample with the urea mass fraction in urea/thiourea of 0.1 is mainly caused by the separate crystallization of some thiourea molecules.
We can conclude that the TU molecules participate in the crystallization of IC. If the TU molecules did not participate in the crystallization of IC, the theoretical enthalpy of fusion of the three-component samples should be x times of that of PEO/U IC (141.4 J/g), which is much lower than the values in Figure 2b obtained from the DSC thermograms. These support that the PEO, urea, and thiourea have formed the ternary isomorphic inclusion compound system.
In a strict isomorphic system, the melting point and enthalpy of fusion should change linearly with crystal composition [34,35,36], but our PEO/urea/thiourea IC system does not follow this law. There are two possible reasons: (1) The two types of host small molecules do not distribute homogeneously at the molecular level in the IC crystals, but exhibit the homogeneous mixing for the XRD characterization; (2) when urea and thiourea molecules crystallize together, some defects in the IC crystals might form. The density of defects would be the largest at the mediate ratio of urea/thiourea, so the melting point of IC shows a V-shape tendency with the content of urea in urea/thiourea and the enthalpy of fusion decreases with the increasing of thiourea content.

4.3. Spherulite Morphologies and Radial Growth Rate of PEO/Urea/Thiourea IC

Here we focus on the crystal growth kinetics of PEO/urea/thiourea IC so as to explore the effect of the urea/thiourea composition on the radial growth rate of PEO/urea/thiourea IC spherulites.
Figure 3 presents the POM images of the spherulites of PEO/urea/thiourea IC with different compositions isothermally crystallized at 70 °C. No phase separation in the mixed melt suggests that PEO chains, urea, and thiourea molecules form a mixed melt at the molecular level, which has also been proved by Ye et al. via infrared spectroscopy [28]. The hydrogen-bond network between the three components prevents the occurrence of phase separation in the mixed melt. Therefore, we believe that the concentration of urea molecules at the crystal growth front of spherulite is the same as the initial global concentration.
Figure 4 plots the measured radial growth rate of PEO/urea/thiourea IC spherulites varying with the content of urea at different isothermal crystallization temperatures. The urea mass fractions in the urea/thiourea host small molecules are from 100% to 60%. The results of samples with the urea mass fractions in urea/thiourea lower than 60% are not shown here because of the large uncertainties of the data. The radial growth rate decreases with the increasing temperature in each sample, and its special tendency supports the surface nucleation mechanism for spherulite growth [4,17].
Figure 5 illustrates the relationship between the spherulite radial growth rates G of IC and the mass fractions x of urea in urea/thiourea. Figure 6 shows the plot of lnG and lnx. Deduced from the theory based on the Nucleation Theorem, lnG is proportional to lnx of the crystallizable component in the binary blend system if the other component does not participate in the crystallization, as confirmed by our previous work on the PEO/urea/dimethylurea blend system [17], and Equation (7) dominates the crystallization kinetics of the double-component crystallization system. However, to our surprise, the preferable linear dependence of lnG on lnx still occurs in the PEO/urea/thiourea ternary isomorphic IC system. The linear fitting of lnG versus lnx gives the slopes ranging from 5 to 12 at isothermal crystallization temperatures from 60 °C to 85 °C with the square of linear correlation coefficients all above 0.97, which inspires us to consider the reason behind the crystal growth kinetics and the detailed secondary nucleation mechanism.

4.4. Critical Size of the Secondary Nucleus of PEO/Urea/Thiourea IC

According to Equation (7), if the number ratio of urea/thiourea in the critical secondary nucleus is equal to the volume ratio (approximately equal to the mass ratio) of urea/thiourea in the three-component melt before crystallization, i.e., n1*/n2* ≈ x/(1 − x), the value of K should be zero. However, our experimental results demonstrate that the fitting slope K is a non-zero value, which does not change with x when x ranges from 60% to 100%. If n2* adopts a non-zero value, K will vary with x obviously according to Equation (7). Consequently, according to our experimental results that lnG shows linear dependence on lnx, i.e., K is almost constant for all x at each crystallization temperature, the simplest explanation is that n2* is zero (namely, x’ = 1) and the slope K equals n1* for all x. Specifically, the critical secondary nucleus does not contain thiourea molecules. We propose that during the crystal growth process of PEO/urea/thiourea IC, the critical secondary nuclei consist mainly of urea molecules and EO units. Namely, the nucleation event selects urea and EO units from the three-component PEO/urea/thiourea melt, implying a stronger nucleating ability of urea with EO units compared to thiourea molecules. The reasonability of our proposition will be discussed in the following text.
According to our previous research results, the critical secondary nucleus of PEO/urea IC consists of 3 to 9 urea molecules when the crystallization temperature ranges from 65 to 85 °C [17]. The number of urea molecules in the critical secondary nucleus of the three-component PEO/urea/thiourea IC is similar to that of the two-component IC of PEO/urea, confirming that only the urea molecules are selected to form the critical secondary nuclei with EO repeating units. In other words, the nucleation ability of urea molecules to form an inclusion compound with ethylene oxide units is stronger than thiourea molecules, so the selection of urea molecules occurs at the secondary nucleation stage. At the later lateral spreading growth stage used to complete the growth of each layer on the growth front, EO, urea, and thiourea can co-crystallize together. The diagrammatic sketch of the above microscopic process is drawn in Figure 7. Because the secondary nuclei only take a tiny percentage in the samples, X-ray diffraction could not detect the existence of the segregated PEO/urea IC nuclei. A larger size in the critical secondary nucleus compared with the PEO/urea system is probably caused by the more difficult diffusion of urea molecules in the melt in the presence of thiourea.
Another piece of evidence supporting the stronger nucleating ability of urea with EO units comes from the fact that introducing a small amount of urea molecules (e.g., 10%) can make the crystal structure of PEO/urea/thiourea IC adopt the orthorhombic crystal system of the PEO/urea IC β phase, as shown in the above Figure 1a, rather than the monoclinic crystal system of the PEO/thiourea IC α phase [28,37,38].

5. Conclusions

Via theoretical deductions from the Nucleation Theorem based on the entropy effect due to dilution, we establish a method to determine the critical size of secondary nuclei in an isomorphic co-crystal system (the composition of each component in the crystal is the same as that in the melt) from the variation of the radial growth rates of spherulites in the composition. The method is further applied to the three-component system of inclusion compounds formed between small molecules and polymer chains. In the work presented here, we adopt a method to determine the critical size of the secondary nuclei of three-component co-crystal PEO/urea/thiourea inclusion compounds. After the gradual replacement of urea by thiourea molecules with fixed PEO concentration, the crystal structure of PEO/urea/thiourea IC maintains the orthorhombic system of PEO/urea IC β crystals, though the cell lattice parameters expand a bit, with less than 3% on average. This result indicates that the co-crystal of the three components adopts an isomorphic crystal structure. The thermodynamic stability of PEO/urea/thiourea IC and the radial growth rate of spherulites first decrease and then increase with the increasing content of thiourea molecules. Furthermore, the addition of a small amount of urea (such as 10%) can change the crystal structure of the three-component IC from α form of PEO/thiourea IC to the crystal structure close to that of PEO/urea β IC.
A good power-law relationship is observed between the radial growth rate G of the PEO/urea/thiourea IC spherulites and the mass fraction x of urea in urea/thiourea when the isothermal crystallization temperature Tc ranges from 60 to 85 °C. The slope of the linear fitting curve increases from 5 to 12 with increasing temperature, but does not vary with the urea mass fraction. The results reveal that during the crystal growth process of PEO/urea/thiourea IC, the critical secondary nuclei consist mainly of urea molecules and EO units, but in the lateral spreading to complete the growth of each layer, the three components can co-crystallize together to form the three-component inclusion compound. Our finding that the selection of a component occurs at the secondary nucleation stage rather than the following lateral spreading growth of the co-crystal system deserves further study.

Author Contributions

Y.G., B.G., J.X. conceived and designed the experiments; Y.G. performed the experiments; Y.G., B.G., J.X. analyzed the data; Y.G. and J.X. wrote the paper, all the authors discussed the results.

Funding

This research was funded by National Natural Science Foundation of China (NSFC, Grant No. 21873054) and NSFC-DFG joint project (Grant No. 21861132018).

Acknowledgments

The authors are deeply grateful to the anonymous reviewers and the academic editor, Juern W.P. Schmelzer for the valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schick, C.; Androsch, R.; Schmelzer, J.W.P. Homogeneous crystal nucleation in polymers. J. Phys. Condens. Matter 2017, 29, 453002. [Google Scholar] [CrossRef]
  2. Schmelzer, J.; Röpke, G.; Priezzhev, V.B. Nucleation Theory and Applications; Wiley Online Library: Hoboken, NJ, USA, 2005; Volume 76. [Google Scholar]
  3. Xu, J.; Zhang, S.; Guo, B. Insights from polymer crystallization: Chirality, recognition and competition. Chin. Chem. Lett. 2017, 28, 2092–2098. [Google Scholar] [CrossRef]
  4. Mutaftschiev, B. The Atomistic Nature of Crystal Growth; Springer Science & Business Media: Berlin, Germany, 2001; Volume 43. [Google Scholar]
  5. Volmer, M.; Weber, A. Keimbildung in übersättigten Gebilden. Z. Phys. Chem. 1926, 119, 277–301. [Google Scholar] [CrossRef]
  6. Becker, R.; Döring, W. Kinetische Behandlung der Keimbildung in übersättigten Dämpfen. Ann. Phys. 1935, 416, 719–752. [Google Scholar] [CrossRef]
  7. Turnbull, D.; Fisher, J.C. Rate of Nucleation in Condensed Systems. J. Chem. Phys. 1949, 17, 71–73. [Google Scholar] [CrossRef]
  8. Armitstead, K.; Goldbeck-Wood, G.; Keller, A. Polymer crystallization theories. In Macromolecules: Synthesis, Order and Advanced Properties; Springer: Berlin/Heidelberg, Germany, 1992; pp. 219–312. [Google Scholar]
  9. Zhang, M.C.; Guo, B.; Xu, J. A Review on Polymer Crystallization Theories. Crystals 2017, 7, 4. [Google Scholar] [CrossRef]
  10. Xu, K.-L.; Guo, B.-H.; Reiter, R.; Reiter, G.; Xu, J. Simulation of secondary nucleation of polymer crystallization via a model of microscopic kinetics. Chin. Chem. Lett. 2015, 26, 1105–1108. [Google Scholar] [CrossRef]
  11. Šarić, A.; Michaels, T.C.T.; Zaccone, A.; Knowles, T.P.J.; Frenkel, D. Kinetics of spontaneous filament nucleation via oligomers: Insights from theory and simulation. J. Chem. Phys. 2016, 145, 211926. [Google Scholar] [CrossRef]
  12. Kashchiev, D. On the relation between nucleation work, nucleus size, and nucleation rate. J. Chem. Phys. 1982, 76, 5098–5102. [Google Scholar] [CrossRef]
  13. Oxtoby, D.W.; Kashchiev, D. A general relation between the nucleation work and the size of the nucleus in multicomponent nucleation. J. Chem. Phys. 1994, 100, 7665–7671. [Google Scholar] [CrossRef]
  14. Kashchiev, D. Forms and applications of the nucleation theorem. J. Chem. Phys. 2006, 125, 014502. [Google Scholar] [CrossRef] [PubMed]
  15. Ford, I.J. Nucleation theorems, the statistical mechanics of molecular clusters, and a revision of classical nucleation theory. Phys. Rev. E 1997, 56, 5615–5629. [Google Scholar] [CrossRef] [Green Version]
  16. Ford, I.J. Statistical mechanics of nucleation: A review. Proc. Inst. Mech. Eng. Part C: J. Mech. Eng. Sci. 2004, 218, 883–899. [Google Scholar] [CrossRef]
  17. Gao, Y.; Zhang, S.; Shi, J.; Guo, B.; Xu, J. Study of the Crystal Growth Mechanism and Critical Secondary Nucleus Size of Poly (Ethylene oxide)/Urea Inclusion Compound. Cryst. Growth Des. 2019, 19, 3834–3842. [Google Scholar] [CrossRef]
  18. Tadokoro, H.; Yoshihara, T.; Chatani, Y.; Murahashi, S. A preliminary report of structural studies on polyethylene oxide-urea complex. J. Polym. Sci. Part B Polym. Lett. 1964, 2, 363–368. [Google Scholar] [CrossRef]
  19. Suehiro, K.; Nagano, Y. Structural studies on molecular complexes of polyethers, 1. Urea-ethylene oxide oligomer complexes. Makromol. Chem. 1983, 184, 669–674. [Google Scholar] [CrossRef]
  20. Chenite, A.; Brisse, F. Structure and conformation of poly (ethylene oxide), PEO, in the trigonal form of the PEO-urea complex at 173 K. Macromolecules 1991, 24, 2221–2225. [Google Scholar] [CrossRef]
  21. Ye, H.-M.; Peng, M.; Xu, J.; Guo, B.-H.; Chen, Q.; Yun, T.-L.; Ma, H. Conformation transition and molecular mobility of isolated poly (Ethylene oxide) chains confined in urea nanochannels. Polymer 2007, 48, 7364–7373. [Google Scholar] [CrossRef]
  22. Liu, Y.; Pellerin, C. Highly oriented electrospun fibers of self-assembled inclusion complexes of poly (Ethylene oxide) and urea. Macromolecules 2006, 39, 8886–8888. [Google Scholar] [CrossRef]
  23. Liu, Y.; Antaya, H.; Pellerin, C. Characterization of the stable and metastable poly (Ethylene oxide)-urea complexes in electrospun fibers. J. Polym. Sci. Part B Polym. Phys. 2008, 46, 1903–1913. [Google Scholar] [CrossRef]
  24. Gao, Y.; Yao, S.; Ye, H.; Guo, B.; Xu, J. Orientation of polymer chains in spherulites of poly (Ethylene oxide)-urea inclusion compounds. Polymer 2017, 130, 209–217. [Google Scholar] [CrossRef]
  25. Lin, D.-W.; Zhong, Z.; Tang, Y.-R.; Gao, Y.; He, Y.-N.; Guo, B.-H.; Xu, J. Morphology and crystalline structure of inclusion compounds formed between poly (Ethylene glycol) and urea. Chin. J. Polym. Sci. 2014, 32, 1234–1242. [Google Scholar] [CrossRef]
  26. Zhong, Z.; Yang, X.; Wang, B.-H.; Yao, Y.-F.; Guo, B.; Yu, L.; Huang, Y.; Xu, J. Solvent-polymer guest exchange in a carbamazepine inclusion complex: Structure, kinetics and implication for guest selection. CrystEngComm 2019, 21, 2164–2173. [Google Scholar] [CrossRef]
  27. Liu, Y.; Pellerin, C. Stability and phase behavior of the poly (Ethylene oxide)–urea complexes prepared by electrospinning. Polymer 2009, 50, 2601–2607. [Google Scholar] [CrossRef]
  28. Ye, H.-M.; Hong, L.-T.; Gao, Y.; Xu, J. Isomorphism in ternary complex: Poly (Ethylene oxide), urea and thiourea. Chin. Chem. Lett. 2017, 28, 888–892. [Google Scholar] [CrossRef]
  29. Oxtoby, D.W. Nucleation of First-Order Phase Transitions. Acc. Chem. Res. 1998, 31, 91–97. [Google Scholar] [CrossRef]
  30. Sear, R.P. Nucleation: Theory and applications to protein solutions and colloidal suspensions. J. Phys. Condens. Matter 2007, 19, 033101. [Google Scholar] [CrossRef]
  31. Erdemir, D.; Lee, A.Y.; Myerson, A.S. Nucleation of Crystals from Solution: Classical and Two-Step Models. Acc. Chem. Res. 2009, 42, 621–629. [Google Scholar] [CrossRef]
  32. Vekilov, P.G. Nucleation. Cryst. Growth Des. 2010, 10, 5007–5019. [Google Scholar] [CrossRef]
  33. Schmelzer, J.W.P. Comments on the Nucleation Theorem. J. Colloid Interface Sci. 2001, 242, 354–372. [Google Scholar] [CrossRef]
  34. Ye, H.-M.; Wang, R.-D.; Liu, J.; Xu, J.; Guo, B.-H. Isomorphism in Poly (butylene succinate-co-butylene fumarate) and Its Application as Polymeric Nucleating Agent for Poly (Butylene succinate). Macromolecules 2012, 45, 5667–5675. [Google Scholar] [CrossRef]
  35. Pérez-Camargo, R.A.; Arandia, I.; Safari, M.; Cavallo, D.; Lotti, N.; Soccio, M.; Müller, A.J. Crystallization of isodimorphic aliphatic random copolyesters: Pseudo-eutectic behavior and double-crystalline materials. Eur. Polym. J. 2018, 101, 233–247. [Google Scholar] [CrossRef]
  36. Ye, H.-M.; Wang, J.; Wang, C.-S.; Li, H.-F. Unique Isodimorphism of Poly (decamethylene succinate-ran-decamethylene fumarate): Large Pseudoeutectic Region and Fantastic Crystallization/Melting Behavior. Macromolecules 2019, 52, 1447–1457. [Google Scholar] [CrossRef]
  37. Liu, Y.; Antaya, H.; Pellerin, C. Structure and Phase Behavior of the Poly (Ethylene oxide)—Thiourea Complex Prepared by Electrospinning. J. Phys. Chem. B 2010, 114, 2373–2378. [Google Scholar] [CrossRef]
  38. Ye, H.-M.; Liu, K.-S.; Zhou, Q. Inclusion complex formed between poly (Ethylene oxide) and thiourea. Colloids Surf. Physicochem. Eng. Asp. 2015, 467, 251–258. [Google Scholar] [CrossRef]
Scheme 1. The equivalent thermodynamic process used to calculate the absolute value Δg1 of the Gibbs free energy difference caused by the crystallization of one A molecule. 1 moL A molecules are assigned as the research system, and we have the equation ΔG1 + ΔG2 = ΔG3 + ΔG4.
Scheme 1. The equivalent thermodynamic process used to calculate the absolute value Δg1 of the Gibbs free energy difference caused by the crystallization of one A molecule. 1 moL A molecules are assigned as the research system, and we have the equation ΔG1 + ΔG2 = ΔG3 + ΔG4.
Entropy 21 01032 sch001
Figure 1. (a) X-ray diffractograms of the poly (ethylene oxide)/urea/thiourea inclusion compound samples with various urea/thiourea mass ratios, the variation of (b) the main diffraction peak positions and (c) the corresponding interplanar spacings belonging to different crystal faces with the urea mass fraction in urea/thiourea, and (d) the expansion percentages of the interplanar spacings belonging to crystal faces (400), (410), (220), (420)/(321), (610) and (032) in the PEO/urea/thiourea IC crystals.
Figure 1. (a) X-ray diffractograms of the poly (ethylene oxide)/urea/thiourea inclusion compound samples with various urea/thiourea mass ratios, the variation of (b) the main diffraction peak positions and (c) the corresponding interplanar spacings belonging to different crystal faces with the urea mass fraction in urea/thiourea, and (d) the expansion percentages of the interplanar spacings belonging to crystal faces (400), (410), (220), (420)/(321), (610) and (032) in the PEO/urea/thiourea IC crystals.
Entropy 21 01032 g001
Figure 2. (a) Differential scanning calometry heating curves of the PEO/urea/thiourea IC samples with various urea/thiourea mass ratios prepared by the isothermal crystallization at 70 °C from the PEO/urea/thiourea mixed melt. (b) The melting points and the enthalpies of fusion of the same samples obtained from the DSC curves.
Figure 2. (a) Differential scanning calometry heating curves of the PEO/urea/thiourea IC samples with various urea/thiourea mass ratios prepared by the isothermal crystallization at 70 °C from the PEO/urea/thiourea mixed melt. (b) The melting points and the enthalpies of fusion of the same samples obtained from the DSC curves.
Entropy 21 01032 g002
Figure 3. Polarized optical microscopy images of PEO/urea/thiourea IC spherulites crystallized from the melt at the isothermal crystallization temperature of 70 °C with gradual mass fraction of urea in urea/thiourea from 100% to 10%. The scale bar is 500 μm in (ae), and 200 μm in (fj).
Figure 3. Polarized optical microscopy images of PEO/urea/thiourea IC spherulites crystallized from the melt at the isothermal crystallization temperature of 70 °C with gradual mass fraction of urea in urea/thiourea from 100% to 10%. The scale bar is 500 μm in (ae), and 200 μm in (fj).
Entropy 21 01032 g003
Figure 4. Radial growth rates of PEO/urea/thiourea IC spherulite with different urea/thiourea mass fraction (from 10/0 to 6/4) varying with the isothermal crystallization temperatures from 60 °C to 85 °C.
Figure 4. Radial growth rates of PEO/urea/thiourea IC spherulite with different urea/thiourea mass fraction (from 10/0 to 6/4) varying with the isothermal crystallization temperatures from 60 °C to 85 °C.
Entropy 21 01032 g004
Figure 5. Radial growth rates of PEO/urea/thiourea IC spherulites measured by POM varying with the urea mass fraction in urea/thiourea at the isothermal crystallization temperature ranging from 60 °C to 85 °C.
Figure 5. Radial growth rates of PEO/urea/thiourea IC spherulites measured by POM varying with the urea mass fraction in urea/thiourea at the isothermal crystallization temperature ranging from 60 °C to 85 °C.
Entropy 21 01032 g005
Figure 6. Plot of the double logarithmic coordinates of the above Figure 5. The perfect linear dependence of lnG on lnx appears in the plot.
Figure 6. Plot of the double logarithmic coordinates of the above Figure 5. The perfect linear dependence of lnG on lnx appears in the plot.
Entropy 21 01032 g006
Figure 7. The diagrammatic sketch of the secondary nucleation and the following lateral spreading growth process. Only the urea molecules and PEO chains can participate in the critical secondary nucleus.
Figure 7. The diagrammatic sketch of the secondary nucleation and the following lateral spreading growth process. Only the urea molecules and PEO chains can participate in the critical secondary nucleus.
Entropy 21 01032 g007

Share and Cite

MDPI and ACS Style

Gao, Y.; Guo, B.; Xu, J. Critical Size of Secondary Nuclei Determined via Nucleation Theorem Reveals Selective Nucleation in Three-Component Co-Crystals. Entropy 2019, 21, 1032. https://doi.org/10.3390/e21111032

AMA Style

Gao Y, Guo B, Xu J. Critical Size of Secondary Nuclei Determined via Nucleation Theorem Reveals Selective Nucleation in Three-Component Co-Crystals. Entropy. 2019; 21(11):1032. https://doi.org/10.3390/e21111032

Chicago/Turabian Style

Gao, Yang, Baohua Guo, and Jun Xu. 2019. "Critical Size of Secondary Nuclei Determined via Nucleation Theorem Reveals Selective Nucleation in Three-Component Co-Crystals" Entropy 21, no. 11: 1032. https://doi.org/10.3390/e21111032

APA Style

Gao, Y., Guo, B., & Xu, J. (2019). Critical Size of Secondary Nuclei Determined via Nucleation Theorem Reveals Selective Nucleation in Three-Component Co-Crystals. Entropy, 21(11), 1032. https://doi.org/10.3390/e21111032

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