Next Article in Journal
Stochastic Schrödinger Equations and Conditional States: A General Non-Markovian Quantum Electron Transport Simulator for THz Electronics
Next Article in Special Issue
Ice-Crystal Nucleation in Water: Thermodynamic Driving Force and Surface Tension. Part I: Theoretical Foundation
Previous Article in Journal
Robust Frame Synchronization Scheme for Continuous-Variable Quantum Key Distribution with Simple Process
Previous Article in Special Issue
Estimating the Nucleation Ability of Various Surfaces Towards Isotactic Polypropylene via Light Intensity Induction Time Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of the Nucleation Theorem to Crystallization of Liquids: Some General Theoretical Results

by
Jürn W. P. Schmelzer
Institute of Physics, University of Rostock, Albert-Einstein-Strasse 23-25, 18059 Rostock, Germany
Entropy 2019, 21(12), 1147; https://doi.org/10.3390/e21121147
Submission received: 4 November 2019 / Revised: 20 November 2019 / Accepted: 22 November 2019 / Published: 25 November 2019
(This article belongs to the Special Issue Crystallization Thermodynamics)

Abstract

:
Different aspects in applying the nucleation theorem to the description of crystallization of liquids are analyzed. It is shown that, by employing the classical Gibbs’ approach in the thermodynamic description of heterogeneous systems, a general form of the nucleation theorem can be formulated that is valid not only for one-component but generally for multi-component systems. In this analysis, one basic assumption of classical nucleation theory is utilized. In addition, commonly employed in application to crystallization, it is supposed that the bulk properties of the critical clusters are widely identical to the properties of the newly evolving crystal phase. It is shown that the formulation of the nucleation theorem as proposed by Kashchiev [J. Chem. Phys. 76, 5098-5102 (1982)], also relying widely on the standard classical approach in the description of crystal nucleation, holds for multi-component systems as well. The general form of the nucleation theorem derived by us is taken then as the starting point for the derivation of particular forms of this theorem for the cases that the deviation from equilibrium is caused by variations of either composition of the liquid phase, temperature, or pressure. In this procedure, expressions recently developed by us for the curvature dependence of the surface tension, respectively, its dependence on pressure and/or temperature are employed. The basic assumption of classical nucleation theory mentioned above is, however, in general, not true. The bulk and surface properties of the critical crystal clusters may differ considerably from the properties of the evolving macroscopic phases. Such effects can be incorporated into the theoretical description by the application of the generalized Gibbs approach for the specification of the dependence of the properties of critical crystal clusters on the degree of metastability of the liquid phase. Applying this method, it is demonstrated that a similar formulation of the nucleation theorem, as derived based on classical nucleation theory, holds true also in cases when a dependence of the state parameters of the critical clusters on the degree of deviation from equilibrium is appropriately accounted for.
PACS:
64.60.Bd General theory of phase transitions; 64.60.Q Nucleation; 81.10.Aj Theory and models of crystal growth; 82.60.Nh Thermodynamics of nucleation in physical chemistry and chemical physics

1. Introduction

Crystallization processes play a decisive role in a broad spectrum of phenomena occurring in nature and technology. They are governed by nucleation and growth of crystallites evolving in the ambient liquid phase. Due to its importance, intensive research is devoted to the analysis of crystallization (see, e.g., [1,2,3,4]) attempting to exhibit the major factors either catalyzing or inhibiting crystal nucleation and growth.
The interpretation of experimental data on crystal nucleation and growth processes is performed up to now predominantly via the classical nucleation theory (CNT) [5,6,7,8,9]. It has been developed based on the work of Farkas [10], Stranski and Kaischew [11], Becker and Döring [12], Volmer [13], Frenkel [14], Fisher and Turnbull [15], and others. Thermodynamic aspects of CNT are described by the fundamental approach to the description of heterogeneous systems developed by J. W. Gibbs [16,17] originally in application to condensation and boiling. Clusters of the newly evolving phase are assumed to be of spherical shape and surface contributions to the thermodynamic quantities are described by the surface tension. The theoretical foundation of the possibility to treat critical crystal nucleation (critical crystals may have, in general, quite different shapes and the relations describing the surface contributions to the thermodynamic potentials are, in general, much more complex) in terms of such simplified model is described in detail in [18,19]. This model we will use in the present study.
Utilizing CNT, the expression for the steady-state nucleation rate, J ( T , p ) , is commonly written in the form [4]
J ( T , p ) = J 0 exp W c k B T , J 0 = c σ k B T D d 0 c σ k B T η d 0 2 ,
W c = 16 π 3 σ 3 ( Δ g ) 2 , R c = 2 σ Δ g .
Here, W c is the work of critical cluster formation, k B is the Boltzmann constant, and T is the absolute temperature. Assuming a spherical shape (and choosing, as it is commonly done in applications to nucleation, the surface of tension as the dividing surface [16,17]), the surface area of the critical cluster is given by A c = 4 π R c 2 , where R c is the critical cluster radius and σ is the surface tension both taken for this particular dividing surface.
The basic concepts of Gibbs’ thermodynamic theory of heterogeneous systems, used in derivation of Equation (2), are also employed in the generalization of Gibbs’ classical treatment performed by us [20,21,22,23]. Being aware of possible different opinions concerning the correct form of Equation (2) (denoted commonly as “Renninger–Wilemski problem” [24]), we would like to underline that Equation (2) hold—as far as the surface of tension is chosen as the dividing surface as we do—both in cases when the surface tension is taken as a constant involving the capillarity approximation and similarly also when a curvature dependence of the surface tension is accounted for [16,17,25,26,27,28,29]. The critical cluster size is determined in such approach by the ratio of surface tension, σ , and the thermodynamic driving force of crystallization, Δ g , the change of the bulk contributions to the Gibbs’ free energy per unit volume of the crystal phase both for one-component and multi-component systems [4,30].
In applications of CNT to crystal nucleation, the thermodynamic driving force of crystallization is specified in accordance with Gibbs’ theory of surface phenomena. It has the same value for both nucleation and growth processes, i.e., the possibility that critical clusters have properties different from the properties of the evolving crystalline phase is, as a rule, not accounted for. Corrections in the theoretical description of crystal nucleation rates are introduced in CNT into the description exclusively via curvature corrections to the surface tension [18,19]. In contrast, employing the generalized Gibbs approach [4,21,22,23], the bulk properties of the critical clusters are also shown to deviate, as a rule, significantly from the properties of the newly evolving macroscopic phases and can be appropriately accounted for.
While the exponential term in Equation (1) is determined by thermodynamic aspects of nucleation, the pre-exponential factor, J 0 , in this relation is governed by the kinetic mechanisms underlying the aggregation processes. The expressions given in Equation (1) have been advanced originally for the case of homogeneous crystallization of one-component liquids. These liquids are characterized by a particle density, c = 1 / d 0 3 , where d 0 is a measure of the size of the ambient phase particles, and D is their diffusion coefficient in the liquid. Since the relevant for nucleation and growth diffusion coefficients are frequently not known, they are replaced in applications of CNT to crystallization by the Newtonian viscosity, η . This replacement is performed in terms of the Stokes–Einstein–Eyring relation, D k B T / η d 0 . Although formulated originally in application to one-component systems, Equation (1) can be employed also in more general and more realistic cases, as discussed in detail in [4,31,32]. In particular, for multi-component systems nucleation may be limited not kinetically but by diffusion processes of the different components or by some intermediate form of aggregation kinetics (see also [33]). Consequently, in such cases the parameter d 0 may become dependent on the size of the critical clusters.
For one-component systems we can perform the substitution Δ g = Δ μ / v α , where Δ μ is the difference between the chemical potential per particle in the liquid and the crystal phases and v α d 0 3 is the volume per particle in the crystal phase. Utilizing these notations, we obtain approximately
d W c d ( Δ μ ) = n c .
Here, n c is the number of particles in the critical cluster. The assumptions employed in the derivation of Equation (3) will be specified below. One of them consists in the application of the capillarity approximation, i.e., the surface tension is assumed to be equal to the respective value for an equilibrium coexistence of liquid and crystal at planar interfaces.
Relations of the form of Equation (3) are widely discussed in the literature [3,4,16,21,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51], starting immediately with the work by Gibbs. They are denoted, following Kashchiev [36], as the first nucleation theorem. The interest in such relations is caused by the fact that they allow one to derive conclusions on the properties of critical clusters in nucleation based on measurements of the nucleation rate. Indeed, in line with Equation (3), variations of the nucleation rate are intrinsically connected with changes of the work of critical cluster formation. Consequently, from measurements of the nucleation rate one can make conclusions on the change of the work of critical cluster formation caused by variations of the external control parameters and, utilizing relations of the form of Equation (3), on the properties of critical clusters. Of course, in such procedures the variations of the kinetic coefficients like diffusion coefficient or viscosity have to also be appropriately accounted for (see, e.g., [3,46]). Most studies of applications of the nucleation theorem are devoted to condensation. Here we concentrate our attention on crystal nucleation utilizing appropriate approximations valid for this particular case of phase formation (for more general approaches and the discussion of different forms of the nucleation theorem, see, e.g., [40,42,43,44]).
It was stressed by Kashchiev and Oxtoby [36,40] that relations of the form of Equation (3) may hold independently on the method of how surface correction terms are introduced into the expression for the change of the Gibbs free energy. However, once such surface corrections terms may be of significance, then their proper specification is of huge importance for a correct application of the nucleation theorem to the determination of the properties of critical clusters (see, e.g., [42,43]). As noted in another context by Albert Einstein in conversation with Werner Heisenberg: “Only the theory decides what one can observe” [52]. In our case, this statement refers to the question how one can determine the properties of critical clusters from the derivatives of the work of critical cluster formation by the thermodynamic driving force or other appropriate control parameters specifying theoretically correctly both the bulk and, in particular, the surface contributions to critical cluster formation.
In line with the above given considerations, in the present paper, we derive several versions of Equation (3) utilizing different approaches in the specification of both the bulk and surface contributions to the work of critical cluster formation. We start the analysis first taking as granted the assumptions commonly employed in CNT (Section 2.1) assuming, in particular, that the expressions for the thermodynamic driving force employed are widely correct. We perform here the analysis immediately for the general case of crystal nucleation in multi-component systems going beyond the case of one-component systems as studied in [36] and extending the results obtained by us earlier in [42]. Once the expressions for the thermodynamic driving force are considered to be widely correct, the only method which can be used in CNT in order to reconcile experimental data with theoretical predictions consists in the introduction of a size dependence of the surface tension [16,21,53,54,55,56] of the critical clusters (or, equivalently, a dependence of the surface tension on pressure and temperature determining the degree of metastability of the melt, respectively, the thermodynamic driving force). Based on the Stefan–Skapski–Turnbull rule, in recent papers [18,19,30,57,58] expressions for the dependence of the surface tension on both external control parameters have been derived. These relations are utilized in the present analysis to arrive at expressions for the derivatives of the work of critical cluster formation with respect to the thermodynamic driving force. In experiments, the thermodynamic driving force can be varied by changes of composition, pressure, and/or temperature. By this reason, we will formulate dependencies of the form of Equation (3) also for these particular cases (Section 2.2). Our approach is compared with the method employed by Kashchiev [36] in order to arrive at expressions of the form of Equation (3). It is shown that, provided the underlying assumptions hold, the approach developed by Kashchiev is valid also for multi-component systems (Section 2.3).
In a second step of the present analysis, an extension of Equation (3) is derived for the case that the work of critical cluster formation is described in terms of the generalized Gibbs approach not utilizing a variety of assumptions inherent in CNT and, in particular, accounting appropriately for the dependence of the bulk properties of critical crystallites on the degree of deviation from thermodynamic equilibrium. As will be shown, also for such a general situation, a formulation of the nucleation theorem in application to crystal nucleation can be given (Section 3) which is widely similar to the result obtained by us in the present paper in terms of CNT. A summary of the results and their brief discussion completes the paper (Section 4).

2. Nucleation Theorem: Analysis in Terms of Classical Nucleation Theory

2.1. Basic Equations and Results

In CNT, the thermodynamic driving force of crystallization, Δ g , is expressed commonly in a first approximation via the difference of the bulk contributions to the Gibbs free energy between liquid and crystal per unit volume of the crystal phase [3,4,19,30] both taken at the pressure, p, and temperature, T, of the liquid phase, i.e., Δ g Δ g C N T ( T , p ) with
Δ g C N T ( T , p ) = i = 1 k ρ i α μ i β T , p , { x i β } μ i α T , p , { x i α } .
Here, μ i α and μ i β are the chemical potentials per particle of the different components in the crystal (specified by the subscript α ) and the liquid (specified by the subscript β ) phases, x i α / β and ρ i α / β are the molar fractions and particle densities of the different components in the respective phases. In the derivation of this relation it is assumed that the composition of the critical crystallites is widely identical to the composition of the newly evolving macroscopic phases [4,30]. Consequently, the thermodynamic driving force is essentially a function of pressure and temperature.
Instead of Δ g , frequently the chemical potential difference, Δ μ , is also utilized for the description of the thermodynamic driving force being for multi-component systems of the form [4,30]
Δ μ ( T , p ) = i = 1 k x i α μ i β T , p , { x i β } μ i α T , p , { x i α } ,
Δ g C N T ( T , p ) = 1 v α Δ μ ( T , p ) , x i α = n i α n , v α = V α n .
Here, n is the total number of particles and n i α are the numbers of particles of the i-th component, x i α are its molar fractions and v α the (average) volume per particle all of them referring to the crystal phase, { x i β } is the set of independent molar fractions of the different components in the liquid phase. Employing these notations, Equation (2) can be reformulated as
W c = 16 π 3 σ 3 v α 2 ( Δ μ ) 2 , R c = 2 σ v α Δ μ ,
and supplemented by a relation for the number of particles, n c , in the critical cluster
n c = 32 π 3 σ Δ μ 3 v α 2 .
By taking the derivative of W c with respect to Δ μ , we obtain then
d W c d ( Δ μ ) = n c 1 3 2 d ln σ d ln ( Δ μ ) .
Utilizing the capillarity approximation (assuming the surface tension to be equal to the respective value for a planar equilibrium coexistence of liquid and crystal phases at a planar interface, i.e., that it does not change in response on variations of the degree of deviation from equilibrium, Δ μ ), Equation (9) is reduced to the particular form given by Equation (3) but here referring to nucleation in multi-component systems.

2.2. Some Alternative Forms of the Nucleation Theorem

Equation (9) connects variations of the work of critical cluster formation in dependence on the thermodynamic driving force of the phase transformation, expressed here via Δ μ ( T , p ) , with changes of the surface tension caused by such variations of the thermodynamic driving force. In order to derive particular expressions for the application of this relation to experiment, both the expressions for the thermodynamic driving force and the surface tension have to be specified expressing them via those control parameters, which can be varied in experiment. Such specification will be performed here again for the general case of crystal nucleation in multi-component systems.
Equation (4) is the most general expression for the description of the thermodynamic driving force of crystallization in CNT. However, detailed knowledge of the chemical potentials of the different components in both phases is commonly not available. By this reasoning, this general relation is transformed into expressions allowing one to compute Δ g ( T , p ) directly by knowing the parameters ( T m , p m ) of some particular macroscopic equilibrium state liquid-crystal and the respective deviations of temperature, ( T T m ) , and pressure, ( p p m ) , from it. Following such ideas, the dependence of thermodynamic driving force of crystallization on temperature and pressure can be described in a good approximation in the form (for details see [4,19,30,58])
Δ g C N T ( T , p ) Δ h m 1 T T m 1 Δ c p ( T m , p m ) 2 Δ s m 1 T T m
+ p m Δ v m p p m 1 1 p m Δ κ T ( T m , p m ) 2 Δ v m p p m 1 .
Here, Δ h m and Δ s m are the melting enthalpy and melting entropy per unit volume of the crystal at the temperature, T m , and the pressure, p m , Δ v m = ( V l i q u i d V c r y s t a l ) / V c r y s t a l , Δ c p and Δ κ T are the differences of specific volume, specific heat, and isothermal compressibility also taken at ( T m , p m ) . These relations can be employed for both one- and multi-component (considered here) systems and both stoichiometric and non-stoichiometric crystallization. As a rule, the inequalities Δ h m > 0 and Δ v m > 0 hold. However, for some substances like water, crystallization may be accompanied by an increase of the volume resulting in Δ v m < 0 . In such cases, a decrease of pressure is required in order to initiate crystal nucleation. These phenomena are discussed in detail in an accompanying paper [59] (see also [60]).
Similarly to the thermodynamic driving force, the dependence of the surface tension on temperature and pressure has to be known for the determination of the work of critical cluster formation, Equation (2). This dependence can be written as [19,57,58]
σ ( T , p ) σ ( T m , p m ) T T m 1 Δ c p Δ s m 1 T T m p m Δ α T Δ s m p p m 1 .
Here, Δ α T is the difference in the thermal expansion coefficients between liquid and crystal. A combination of Equations (9)–(11) allows us to determine the change of the work of critical cluster formation in dependence on either pressure or temperature. The respective computations are straightforward and by this reason not given here in the general form. Anyway, to have some impression below two special cases are considered employing truncated expressions for the dependencies of the thermodynamic driving force and the surface tension.
In case the chemical potential difference and surface tension is caused by temperature variations described in such first-order approximation as
Δ μ ( T , p m ) v α Δ h m 1 T T m , σ ( T , p m ) σ ( T m , p m ) T T m ,
Equation (9) yields
d W c d T n c 3 v α Δ h m 2 T 1 T 3 T m .
For pressure-induced nucleation, using the approximations,
Δ μ ( T m , p ) v α p m Δ v m p p m 1 ,
σ ( T m , p ) σ ( T m , p m ) 1 p m Δ α T Δ s m p p m 1 ,
we obtain instead
d W c d p n c v α Δ v m 2 2 + p m Δ α T Δ s m p p m 1 1 p m Δ α T Δ s m p p m 1
both as zeroth-order estimates. As evident, in both cases the partial derivatives of the work of critical cluster formation taken with respect to the external control parameters are proportional to the number of particles in the critical clusters. However, in these cases the coefficient of proportionality is generally quite different from one.
As noted above, in CNT the bulk state parameters of the critical crystallites are assumed to be widely identical to the properties of the newly evolving macroscopic phases. Consequently, in line with Equation (5), the thermodynamic driving force can be considered as a function of the set of variations of the chemical potential differences, { Δ μ i } , of the different components. Provided it would be possible to vary only one of the chemical potential differences, say Δ μ j , then (see, e.g., also [50,51])
W c ( Δ μ j ) = d W c d ( Δ μ ) Δ μ ( Δ μ j ) = n j c 1 3 2 d ln σ d ln ( Δ μ ) , n j c = n c x j α .
Here, n j c is the number of particles of the j-th component in the critical cluster. However, in any attempts to realize such variations in experiment, one has to take care about possible limitations of the applicability of these equations caused by necessity to fulfil the Gibbs–Duhem relations [42,43,61].
In contrast to Δ μ j , the independent molar fractions, x j β , j = 1 , 2 , , ( k 1 ) , of the liquid phase can be varied separately. Assuming that such process of variation of only one independent molar fraction is performed at constant pressure and temperature, we obtain
W c x j β = d W c d ( Δ μ ) Δ μ x j β
or
W c x j β = n c 1 3 2 d ln σ d ln ( Δ μ ) x j β i = 1 k x i α μ i β T , p , { x i β } μ i α T , p , { x i α } .
A change of the composition of the liquid may lead to variations of the composition of the crystal phase coexisting in equilibrium with the liquid. If such variations of the composition of the crystalline phase can be neglected, we obtain as a special case
W c x j β = n c 1 3 2 d ln σ d ln ( Δ μ ) i = 1 k x i α μ i β T , p , { x i β } x j β .
Again, the partial derivative ( W c / x j β ) is proportional to the critical cluster size, n c , but the evaluation of the coefficient of proportionality requires a detailed knowledge both on the bulk properties of the coexisting phases and the change of the surface tension with respect to variations of the thermodynamic driving force for this particular method of variation of the state of the liquid.

2.3. Comparison with the Approach Employed by Kashchiev

In his first paper ([36], see also [3,40]), Kashchiev considered exclusively one-component systems. He formulated an equation for the work of critical cluster formation of the form of Equation (7) for this particular case. He assumed then validity of the capillarity approximation and arrived at Equation (3). As shown above, this method is similarly applicable to the description of crystal nucleation for the general case of multi-component systems and accounting for a dependence of the surface tension of critical crystal nuclei immediately yields the generalization given by Equation (9).
The account of surface free energy contributions was performed by Kashchiev in a different way as compared to the standard one used here. In line with Gibbs’ thermodynamic approach and the approximations usually employed in the description of crystal nucleation, the bulk contribution to the change of the Gibbs’ free energy was written as n Δ μ supplemented but surface energy terms given by F ( n , Δ μ ) . Accounting for both bulk and surface contributions, the work of formation of a cluster containing n particles is expressed then as [36]
W ( n , Δ μ ) = n Δ μ + F ( n , Δ μ ) .
It is stated that the derivative of the work of critical cluster formation with respect to Δ μ reads
d W c d ( Δ μ ) = n c + F c ( Δ μ ) ,
denoting this expression as nucleation theorem. In the second term on the right hand side of Equation (21), the partial derivative with respect to Δ μ implies that the derivative has to be taken at constant values of the cluster size.
Utilizing Gibbs’ method and certain approximations, one obtains as a special case F ( n , Δ μ ) σ A then. Here, A = 4 π R 2 is the surface area of the crystallite. A substitution of this relation into Equation (21) leads immediately to Equations (7) and (9), again, when Δ μ is expressed via Equation (5). Consequently, in contrast to the respective statements in [36,40], Equations (20) and (21) hold not only for one- but also for multi-component systems as far as the assumptions commonly utilized in CNT are fulfilled. Such extension of the approach formulated by Kashchiev [36] to multi-component systems was also performed already by Yang Gao et al. in [50,51].
Although both methods lead for this particular case to the same results, we consider our approach to be preferable. It involves exclusively the analysis of parameters of critical clusters and does not require the specification of the bulk and surface contributions to cluster formation for arbitrary cluster sizes as it has to be done in the formulation of Equation (20). If this procedure is realized in a general form, very detailed considerations have to be performed, as discussed in detail in [20,21,22,62,63]. Moreover, it is supposed in the formulation of Equation (20) that the surface contributions are a function of the number of particles in the cluster, n, and of the driving force, Δ μ . In general, this is as a rule not the case as evident from Gibbs’ theory of heterogeneous systems [16,17]. According to Gibbs’ adsorption equation, the surface tension of critical clusters is a function of either the set of state parameters of the newly evolving or of the ambient phase. For example, as evident from Equation (11), the surface tension of critical clusters can be expressed as a function of pressure and temperature. In line with Equation (10), one and the same values of Δ μ may yield, in principle, different values of the surface tension. Moreover, utilizing exclusively Gibbs’ theory (referring to thermodynamic equilibrium states) nothing can be said concerning the value of the surface tension for clusters not being in equilibrium with the ambient liquid.
In our approach, it is merely required that infinitesimal variations of the state parameters of the ambient phase like ( T , p , { x i β } ) lead to infinitesimal variations of W c , Δ μ , and σ (referred also exclusively to the critical cluster) and that the differential quotient, ( d σ / d ( Δ μ ) ) , has well-defined values. This difference in the interpretation is even more important if the expressions for the thermodynamic driving force and the surface tension are more complex as compared with the relations formulated based on the assumptions commonly employed in CNT and underlying also Equation (20). The respective generalizations and their consequences will be analyzed in the next section.

3. Thermodynamics of Cluster Formation: Beyond the Classical Gibbs’ Approach

3.1. One Main Deficiency of Classical Nucleation Theory

As noted, for example, by John von Neumann [64], a theory has to fulfil two demands, it has “correctly to describe phenomena from a reasonably wide area. Furthermore, it must satisfy certain esthetic criteria—that is, in relation to how much it describes, it must be rather simple”. In CNT, the possibility to describe in a relatively simple way experimental data on crystal nucleation is related to the assumption that the bulk properties of the critical clusters are widely identical to the properties of the newly evolving macroscopic phases. Once the critical clusters have such properties, also sub- and super-critical clusters can be considered to have the same properties.
However, already in the course of developing the basic concepts of CNT it has been understood clearly that the bulk properties of the critical (and, consequently, of sub- and super-critical) crystallites may also be quite different as compared to the properties of the newly evolving macroscopic phases (see, e.g., [4,7,65,66,67,68]). As a consequence, alternative continuums approaches to the description of heterogeneous systems have been developed and applied to the description of phase formation processes [69,70,71], repeating and advancing the work performed much earlier by J. D. van der Waals [72,73]. These theoretical methods and their generalization, as well as computer simulations, confirm in line with an increasing number of experimental studies the point of view that the parameters of the critical cluster may differ significantly from the properties of the evolving macroscopic phases (e.g., [74,75,76,77,78,79,80,81,82]). Consequently, the computation of the thermodynamic driving force performed under the assumption of nearly constant macroscopic values of the bulk state parameters of the crystal clusters can be only an approximation. In some cases, this assumption may appropriately reflect the kinetics of crystallization for some selected systems. As a rule, however, it may lead to severe inconsistencies in the description of the experimental data.
The mentioned basic assumption of CNT is supported by Gibbs’ classical method of description of heterogeneous systems. Consequently, the question arises whether one has to abandon this approach at all. However, as shown by us these problems can be resolved by developing a generalization of the classical Gibbs approach [4,21,22,23,83,84,85,86]. As a first step in such a generalization, we developed a theory of heterogeneous systems for non-equilibrium states. This is a straightforward generalization of Gibbs’ classical approach, who restricted his considerations exclusively to equilibrium states. It accounts, in particular, appropriately for the fact that in non-equilibrium states the surface tension has to depend, in general, on the state parameters of both the ambient and newly evolving phases.
In order to determine the thermodynamic potentials for non-equilibrium states, one of the problems one has to solve consists in the determination of the values of the chemical potentials, μ j σ , and the temperature, T σ , referring to the surface phase. In accordance with an intensive analysis of such problems by Defay et al. [87], Prigogine and Bellemans [88] (“a surface phase has no real autonomy, in general”), and Rowlinson and Widom [89] (“we cannot measure or define unambiguously the thermodynamic properties of the surface phase”) we postulated [20,21,22] that the superficial surface parameters, T σ and μ j σ , have to be set equal to the respective quantities of the ambient phase
T σ = T β , μ j σ = μ j β .
As one consequence, in above equations, the superficial particle numbers, n i σ , and the superficial entropy, S σ , which enter the description in terms of Gibbs theory as correction terms describing effects of the interface, are included into the number of particles of the newly evolving phase, n i α = ρ i α V α , and the volume density of the entropy, s α , of this phase, correspondingly.
Based on such generalization, more complex equations for the description of the properties of the critical clusters governing nucleation are obtained as compared to Gibbs’ classical treatment. Accounting for such effects and the described method of their description, we will now advance a reformulation of the nucleation theorem.

3.2. Basic Equations and Results

Both in the classical and generalized Gibbs’ approaches [16,20,21,22,23,23], the work of critical cluster formation, W c , and the critical cluster radius, R c (referred to the surface of tension), can be expressed as
W c = 16 π 3 σ 3 ( p α , T α , { x i α } ; p β , T β , { x i β } ) Δ g ( p α , T α , { x i α } ; p β , T β , { x i β } ) 2 ,
R c = 2 σ ( p α , T α , { x i α } ; p β , T β , { x i β } ) Δ g ( p α , T α , { x i α } ; p β , T β , { x i β } ) .
The thermodynamic driving force of cluster formation is given generally by [29,30]
Δ g ( p α , T α , { x i α } ; p β , T β , { x i β } ) = s α ( T β T α ) + ( p α p β )
+ i = 1 k ρ i α μ i β p β , T β , { x i β } μ i α p α , T α , { x i α } .
In the derivation of these relations, the fundamental equation of Gibbs’ approach was generalized to allow the surface tension to be a function of the state parameters of both ambient and newly evolving phases.
Equation (4), employed commonly in CNT, is obtained as a limiting case from Equation (25) provided that temperature and pressure of the critical cluster and the ambient liquid coincide. The first of these assumptions ( T α T β ) is in agreement with the equilibrium conditions
μ i β p β , T β , { x j β } μ i α p α , T α , { x j α } = 0 for i = 1 , 2 , , k ,
T β T α = 0 ,
p β p α + 2 σ R c = 0 ,
as obtained in the classical Gibbs approach [16,17], while the second ( p α p β ), in general, is not. Approximately, Equation (4) can be obtained from Equation (25) by a truncated Taylor expansion [4,30]. In this procedure it is assumed that the bulk state parameters of the critical clusters depend neither on pressure nor on temperature. The conditions T α T β and p α p β are also employed in the derivation of Equations (8) and (11) for the dependence of the specific surface energy on temperature and pressure.
In the framework of the generalized Gibbs approach, the critical cluster size is determined by Equation (24) in combination with Equation (25), again. Moreover, instead of Equations (26) and (27) we get the additional equilibrium conditions now in the form [4,21,22,23,29],
μ i β μ i α = 3 R c σ ρ i α { ρ i β } , s β for i = 1 , 2 , , k ,
T β T α = 3 R c σ s α { ρ i β } , s β .
Both the equilibrium conditions for pressure and also for thermal and diffusion equilibrium are modified. In particular, we conclude that the first of the assumptions employed in CNT ( T α = T β ) for the determination of the thermodynamic driving force also also lacks any theoretical foundation (see also [90,91]). Note further that in this more general situation (described by Equations (29) and (30)), Equation (24) is not reduced to Equation (28) as this is the case in Gibbs’ classical approach.
It follows as an additional consequence that the bulk parameters of the critical clusters may vary considerably in dependence on the degree of metastability. This effect is determined by the derivatives of the surface tension with respect to the state parameters of the cluster phase. In the generalized Gibbs’ approach, the surface tension is a function of the state parameters of both ambient and newly evolving phases. By this reason, here derivatives of the surface tension occur in the equilibrium conditions. Indeed, utilizing Equations (29) and (30), the correct expression, Equation (25), for the thermodynamic driving force of nucleation can be expressed as
Δ g = ( p α p β ) + 3 R c s α σ s α { ρ j β } , s β + i = 1 k ρ i α σ ρ i α { ρ j β } , s β .
Variations of the bulk state parameters of the cluster phase in response to variations of the degree of metastability may be, consequently, of considerable effect concerning the value of the thermodynamic driving force.
In order to employ above derived equations, we have to know the surface tension (and here in an even more complex form as compared to the classical approach) as a function of the state parameters of both coexisting phases. An expression for the surface tension can be obtained, again, via the Stefan–Skapski–Turnbull relation. We may write [57,58]
σ σ ( T m , p m ) = 1 T m Δ s m T β s β T α s α + i = 1 k ρ i α μ i β μ i α .
This expression allows us to determine the dependence of the surface tension on pressure and temperature, again. However, in contrast to Equation (11), here the dependence of the bulk state parameters of the critical clusters on pressure and temperature has to also be accounted for. These state parameters have to be determined by the equilibrium conditions Equations (29) and (30), again. Substituting the equilibrium conditions into Equation (32), we get instead of Equation (8)
σ σ ( T m , p m ) = 1 T m Δ s m T β ( s β s α ) + 3 R c s α σ s α { ρ j β } , s β + i = 1 k ρ i α σ ρ i α { ρ j β } , s β
(for details see [57]).
The specification of the dependence of the thermodynamic driving force and the surface tension, accounting for possible variations of both bulk and surface state parameters of the critical clusters, is, consequently, a highly complex task (see also [23]). For any set of state parameters of the ambient liquid phase ( p β , T β , { x i β } ) , first the state parameters of the critical clusters have to be determined. This procedure can be performed via the equilibrium conditions, Equations (24), (29) and (30), supplying us with ( k + 2 ) relations for the determination of the ( k + 1 ) bulk state parameters and the size of the critical cluster. Having at one’s disposal these parameters, both the surface tension and the thermodynamic driving force can be computed. Infinitesimal variations of the state parameters of the ambient liquid phase will result, consequently, in similar variations of both the thermodynamic driving force, d Δ g , and the surface tension, d σ . For the formulation of the general form of the nucleation theorem employing the generalized Gibbs’ approach only latter mentioned conclusion is of importance, i.e., that the ratio ( d σ / d ( Δ g ) ) is well-defined.
Indeed, as noted, the work of critical cluster formation is given also in the generalized Gibbs’ approach by the simple relation, Equation (23). Taking the derivative with respect to Δ g , we obtain similarly to Equations (8) and (9)
V c = 32 π 3 σ Δ g 3 ,
d W c d ( Δ g ) = V c 1 3 2 d ln σ d ln ( Δ g ) .
Here, V c is the volume of the critical cluster. The general form of the nucleation theorem remains the same as derived in employing basic assumptions of CNT. Of course, introducing the assumptions of CNT at this stage, Equation (9) is reproduced as a special limiting case, again.
Note, however, that also another expression derived in terms of Gibbs classical approach can be obtained based on Equation (35) (but not via Equation (9) involving additional approximations sketched here earlier) as a special case. Utilizing Gibbs’ equilibrium conditions, Equations (26) and (27), Δ g can be replaced by ( p α p β ) . In such a limiting case, we obtain
d W c d ( p α p β ) = V c 1 3 R c d σ d ( p α p β )
derived by us earlier in [42] utilizing Gibbs’ classical approach [16].
Although the nucleation theorem is of very similar form to the one derived in terms of CNT, the specification of the conclusions of what its consequences are when, e.g., either composition, pressure, or temperature are varied becomes a much more complex task. Consequently, the predictions concerning the parameters of the critical clusters as obtained via the generalized Gibbs approach may be quite different from the conclusions obtained in terms of CNT (see also [42]).

4. Results and Discussion

In the first part of the present analysis of the nucleation theorem, the main ideas of CNT are utilized. The results can be summarized as follows: (i) By employing the classical Gibbs’ approach in the thermodynamic description of heterogeneous systems and the basic assumptions of classical nucleation theory, a general form of the nucleation theorem, Equation (9), can be formulated valid not only for one-component but generally for multi-component systems. (ii) This result is taken then as the starting point for the derivation of particular forms, Equations (13), (15), (16) and (18) of this theorem for the cases that the degree of deviation from equilibrium is caused by variations of either composition of the liquid phase, temperature, and/or pressure. In this procedure, the expressions mentioned above for the curvature dependence of the surface tension, respectively, on pressure and/or temperature are utilized. (iii) It is shown that the nucleation theorem in the form of Equation (21) as expressed by Kashchiev [36] can be extended in application to crystallization to multi-component systems, provided the basic assumptions commonly employed in classical nucleation theory are fulfilled.
In the second part of the analysis, effects caused by changes of the properties of critical clusters are incorporated via the generalized Gibbs’ approach. It is demonstrated that (iv) a similar formulation of the nucleation theorem, Equation (35), as obtained via the assumptions employed in CNT, Equation (9), also holds true in such cases. However, (v) specific forms of the nucleation theorem, describing the change of the work of critical cluster formation with respect to variations of temperature, pressure, and/or composition, will have a considerably more complex shape. Similarly to the specification of the parameters of critical clusters in crystal nucleation, they have to be appropriately accounted for in order to arrive at a correct method of determination of critical cluster properties via the different forms of the nucleation theorem.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Christian, J.W. The Theory of Transformations in Metals and Alloys; Pergamon Press: Oxford, UK, 1981. [Google Scholar]
  2. Skripov, V.P.; Koverda, V.P. Spontaneous Crystallization of Undercooled Liquids; Nauka: Moscow, Russia, 1984. (In Russian) [Google Scholar]
  3. Kelton, K.F.; Greer, A.L. Nucleation in Condensed Matter: Applications in Materials and Biology; Pergamon: Amsterdam, The Netherlands, 2010. [Google Scholar]
  4. Gutzow, I.S.; Schmelzer, J.W.P. The Vitreous State: Thermodynamics, Structure, Rheology, and Crystallization, 2nd ed.; (First Edition, 1995); Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  5. Herlach, D.; Galenko, P.; Holland-Moritz, D. Metastable Solids from Undercooled Melts; Cahn, R.W., Ed.; Pergamon Materials Series; Elsevier: Amsterdam, The Netherlands, 2007; Volume 10. [Google Scholar]
  6. Debenedetti, P.G. Metastable liquids: Concepts and Principles; Princeton University Press: Princeton, NJ, USA, 1996. [Google Scholar]
  7. Fokin, V.M.; Yuritsyn, N.S.; Zanotto, E.D.; Schmelzer, J.W.P. Homogeneous crystal nucleation in silicate glasses: A 40 years perspective. J. Non-Cryst. Solids 2006, 352, 2681–2714. [Google Scholar] [CrossRef]
  8. Schmelzer, J.W.P. (Ed.) Glass: Selected Properties and Crystallization; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2014. [Google Scholar]
  9. Komatsu, T. Design and control of crystallization in oxide glasses. J. Non-Cryst. Solids 2015, 428, 156–175. [Google Scholar] [CrossRef]
  10. Farkas, L. Keimbildungsgeschwindigkeit in übersättigten Dämpfen (English Translation: Nucleation rate in supersaturated vapors). Z. Phys. Chem. 1927, 125, 236–242. [Google Scholar]
  11. Kaischew, R.; Stranski, I.N. Zur kinetischen Ableitung der Keimbildungsgeschwindigkeit (English Translation: On the Kinetic Derivation of the Nucleation Rate). Z. Phys. Chem. 1934, B 26, 317–326. [Google Scholar]
  12. Becker, R.; Döring, W. Kinetische Behandlung der Keimbildung in übersättigten Dämpfen (English Translation: Kinetic Treatment of Nucleation in Supersaturated Vapors). Ann. Phys. 1935, 24, 719–752. [Google Scholar] [CrossRef]
  13. Volmer, M. Kinetik der Phasenbildung (English Translation: Kinetics of Phase Formation); Th. Steinkopff: Dresden, Germany, 1939. [Google Scholar]
  14. Frenkel, Y.I. The Kinetic Theory of Liquids; Oxford University Press: Oxford, UK, 1946. [Google Scholar]
  15. Turnbull, D.; Fisher, J.C. Rate of Nucleation in Condensed Systems. J. Chem. Phys. 1949, 17, 71–73. [Google Scholar] [CrossRef]
  16. Gibbs, J.W. On the Equilibrium of Heterogeneous Substances. Trans. Connect. Acad. Sci. 1879, 3, 108, 343. [Google Scholar] [CrossRef]
  17. Gibbs, J.W. The Collected Works. In Thermodynamics; Longmans & Green: New York, NY, USA; London, UK; Toronto, ON, Canada, 1928; Volume 1. [Google Scholar]
  18. Schmelzer, J.W.P.; Abyzov, A.S.; Ferreira, E.B.; Fokin, V.M. Curvature dependence of the surface tension and crystal nucleation in liquids. Int. J. Appl. Glass Sci. 2019, 10, 57–60. [Google Scholar] [CrossRef]
  19. Schmelzer, J.W.P.; Abyzov, A.S.; Baidakov, V.G. Entropy and the Tolman Parameter in Nucleation Theory. Entropy 2019, 21, 670. [Google Scholar] [CrossRef]
  20. Ulbricht, H.; Schmelzer, J.W.P.; Mahnke, R.; Schweitzer, F. Thermodynamics of Finite Systems and the Kinetics of First-Order Phase Transitions. In Teubner-Texte zur Physik; Teubner-Verlag: Leipzig, Germany, 1988; Volume 17. [Google Scholar]
  21. Schmelzer, J.W.P.; Boltachev, G.S.; Baidakov, V.G. Is Gibbs’ thermodynamic theory of heterogeneous systems really perfect? In Nucleation Theory and Applications; Schmelzer, J.W.P., Ed.; WILEY-VCH: Berlin, Germany; Weinheim, Germany, 2005; pp. 418–446. [Google Scholar]
  22. Schmelzer, J.W.P.; Boltachev, G.S.; Baidakov, V.G. Classical and generalized Gibbs approaches and the work of critical cluster formation in nucleation theory. J. Chem. Phys. 2006, 124, 194503. [Google Scholar] [CrossRef]
  23. Schmelzer, J.W.P.; Abyzov, A.S. Crystallization of glass-forming melts: New answers to old questions. J. Non-Cryst. Solids 2018, 501, 11–20. [Google Scholar] [CrossRef]
  24. Reguera, D.; Reiss, H. Nucleation in confined ideal binary mixtures: The Renninger–Wilemski problem revisited. J. Chem. Phys. 2003, 119, 1533–1546. [Google Scholar] [CrossRef]
  25. Ono, S.; Kondo, S. Molecular theory of surface tension in liquids. In Handbuch der Physik (Encyclopedia of Physics); Flügge, S., Ed.; Springer: Berlin/Heidelberg, Germany; Göttingen, Germany, 1960; Volume 10, pp. 134–280. [Google Scholar]
  26. Renninger, R.G.; Hiller, F.C.; Bone, R.C. Comment on “Self-nucleation in the sulfuric acid-water system”. J. Chem. Phys. 1981, 75, 1584–1585. [Google Scholar] [CrossRef]
  27. Wilemski, G. Composition of the critical nucleus in multi-component vapor nucleation. J. Chem. Phys. 1984, 80, 1370–1372. [Google Scholar] [CrossRef]
  28. Schmelzer, J.W.P. Curvature Dependent Surface Tension and Nucleation Theory. In Nucleation Theory and Applications; Schmelzer, J.W.P., Röpke, G., Priezzhev, V.B., Eds.; Joint Institute for Nuclear Research Publishing House: Dubna, Russia, 1999; Chapter 9; pp. 268–289. ISBN 5-851165-539-9. [Google Scholar]
  29. Schmelzer, J.W.P.; Baidakov, V.G.; Boltachev, G.S. Kinetics of boiling in binary liquid-gas solutions: Comparison of different approaches. J. Chem. Phys. 2003, 119, 6166–6183. [Google Scholar] [CrossRef]
  30. Schmelzer, J.W.P.; Abyzov, A.S. Crystallization of glass-forming liquids: Thermodynamic driving force. J. Non-Cryst. Solids 2016, 449, 41–49. [Google Scholar] [CrossRef]
  31. Schmelzer, J.W.P. On the determination of the kinetic pre-factor in classical nucleation theory. J. Non-Cryst. Solids 2010, 356, 2901–2907. [Google Scholar] [CrossRef]
  32. Fokin, V.M.; Abyzov, A.S.; Zanotto, E.D.; Cassar, D.R.; Rodrigues, A.M.; Schmelzer, J.W.P. Crystal nucleation in glass-forming liquids: Variation of the size of the structural units with temperature. J. Non-Cryst. Solids 2016, 447, 35–44. [Google Scholar] [CrossRef]
  33. Slezov, V.V. Kinetics of First-Order Phase Transitions; Wiley-VCH: Berlin, Germany; Weinheim, Germany, 2009. [Google Scholar]
  34. Hill, T.L. Thermodynamics of Small Systems. J. Chem. Phys. 1962, 36, 3182–3197. [Google Scholar] [CrossRef]
  35. Nielsen, E. Kinetics of Precipitation; Pergamon: Oxford, UK, 1964. [Google Scholar]
  36. Kashchiev, D. On the relation between nucleation work, nucleus size, and nucleation rate. J. Chem. Phys. 1982, 76, 5098–5102. [Google Scholar] [CrossRef]
  37. Anisimov, M.P.; Vershinin, S.N. Spontaneous nucleation rates, sizes and compositions of criticval nuclei in many-component supersaturated vapor. In Atmospheric Aerosols and Nucleation; Wagner, P.E., Vali, G., Eds.; Springer: Berlin, Germany, 1988; pp. 393–396. [Google Scholar]
  38. Baidakov, V.G.; Kaverin, A.M.; Sulla, I.I. Accessible superheating of liquid ethane. Teplofizika Vysokikh Temperatur 1989, 27, 410–412. (In Russian) [Google Scholar]
  39. Viisanen, Y.; Strey, R.; Reiss, H. Homogeneous nucleation rates for water. J. Chem. Phys. 1993, 99, 4680–4692. [Google Scholar] [CrossRef]
  40. Oxtoby, D.; 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]
  41. Bowles, R.K.; McGraw, R.; Schaaf, P.; Senger, B.; Voegel, J.-C.; Reiss, H. A molecular based derivation of the nucleation theorem. J. Chem. Phys. 2000, 113, 4524–4532. [Google Scholar] [CrossRef]
  42. Schmelzer, J.W.P. Comments on the nucleation theorem. J. Colloid Interface Sci. 2001, 242, 354–372. [Google Scholar] [CrossRef]
  43. Schmelzer, J.W.P. Some additional comments on the nucleation theorem. Russ. J. Phys. Chem. 2003, 77 (Suppl. 1), S143–S145. [Google Scholar]
  44. Kashchiev, D. Forms and applications of the nucleation theorem. J. Chem. Phys. 2006, 125, 014502. [Google Scholar] [CrossRef]
  45. Wilemski, G. Volumes of critical bubbles from the nucleation theorem. J. Chem. Phys. 2006, 125, 114507. [Google Scholar] [CrossRef]
  46. Baidakov, V.G.; Tipeev, A.O.; Bobrov, K.S.; Ionov, G.V. Crystal nucleation rate isotherms in Lennard-Jones liquids. J. Chem. Phys. 2010, 132, 234505. [Google Scholar] [CrossRef]
  47. Vehkamäki, H.; McGrath, M.J.; Kurten, T.; Julin, J.; Lehtinen, K.E.J.; Kulmala, M. Rethinking the application of the first nucleation theorem to particle formation. J. Chem. Phys. 2012, 136, 094107. [Google Scholar] [CrossRef]
  48. Malila, J.; McGraw, R.; Laaksonen, A.; Lehtinen, K.E.J. Communication: Kinetics of scavenging of small, nucleating clusters: First nucleation theorem and sum rules. J. Chem. Phys. 2015, 142, 011102. [Google Scholar] [CrossRef]
  49. Perepezko, J.H.; Wilde, G. Melt undercooling and nucleation kinetics. Curr. Opin. Solid State Mater. Sci. 2016, 20, 3–12. [Google Scholar] [CrossRef]
  50. Gao, Y.; Guo, B.; Xu, J. Critical Size of Secondary Nuclei Determined via the Nucleation Theorem Reveals Selective Nucleation in Three-Component Co-Crystals. Entropy 2019, 21, 1032. [Google Scholar] [CrossRef]
  51. 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]
  52. Holton, G. Werner Heisenberg and Albert Einstein. Phys. Today 2000, 53, 38–42. [Google Scholar] [CrossRef] [Green Version]
  53. Tolman, R. Consideration of the Gibbs Theory of Surface Tension. J. Chem. Phys. 1948, 16, 758–774. [Google Scholar] [CrossRef] [Green Version]
  54. Tolman, R. The Effect of Droplet Size on Surface Tension. J. Chem. Phys. 1949, 17, 333–337. [Google Scholar] [CrossRef] [Green Version]
  55. Fokin, V.M.; Zanotto, E.D. Crystal nucleation in silicate glasses: The temperature and size dependence of crystal-liquid surface energy. J. Non-Cryst. Solids 2000, 265, 105–112. [Google Scholar] [CrossRef] [Green Version]
  56. Abyzov, A.S.; Fokin, V.M.; Zanotto, D.; Rodrigues, A.M.; Schmelzer, J.W.P. The effect of elastic stresses on the thermodynamic barrier for crystal nucleation. J. Non-Cryst. Solids 2016, 432, 325–333. [Google Scholar] [CrossRef]
  57. Schmelzer, J.W.P.; Abyzov, A.S. Crystallization of glass-forming liquids: Specific surface energy. J. Chem. Phys. 2016, 145, 064512. [Google Scholar] [CrossRef]
  58. Schmelzer, J.W.P.; Abyzov, A.S.; Fokin, V.M. Thermodynamic Aspects of Pressure-Induced Crystallization: Kauzmann Pressure. Int. J. Appl. Glass Sci. 2016, 7, 474–485. [Google Scholar] [CrossRef]
  59. Hellmuth, O.; Schmelzer, J.W.P.; Feistel, R. Ice-crystal nucleation in water: Thermodynamic driving force and surface tension. Entropy 2019. submitted. [Google Scholar]
  60. Marcolli, C. Technical note: Fundamental aspects of ice nucleation via pore condensation and freezing including Laplace pressure and growth into macroscopic ice. Atmos. Chem. Phys. 2019. [Google Scholar] [CrossRef]
  61. Kubo, R. Thermodynamics; North-Holland Publishing Company: Amsterdam, The Netherlands, 1968. [Google Scholar]
  62. Nishioka, K.; Kusaka, I. Thermodynamic formulas of liquid phase nucleation from vapor in multicomponent systems. J. Chem. Phys. 1992, 96, 5370–5376. [Google Scholar] [CrossRef]
  63. Debenedetti, P.G.; Reiss, H. Reversible work of formation of an embryo of a new phase within a uniform macroscopic mother phase. J. Chem. Phys. 1998, 108, 5498–5505. [Google Scholar] [CrossRef] [Green Version]
  64. John von Neumann, Quotes. Available online: https://www.goodreads.com/author/quotes/205578.John_von_Neuman (accessed on 15 November 2019).
  65. Cahn, J.W. Reflections on Diffusive Interfaces and Spinodal Decomposition. In The Selected Works of J. W. Cahn; Carter, W.C., Johnson, W.C., Eds.; The Minerals, Metals, and Materials Society: Warrendale, PA, USA, 1998; pp. 1–8. [Google Scholar]
  66. Schmelzer, J.W.P.; Schmelzer, J., Jr.; Gutzow, I.S. Reconciling Gibbs and van der Waals: A new approach to nucleation theory. J. Chem. Phys. 2000, 112, 3820–3831. [Google Scholar] [CrossRef]
  67. Schmelzer, J.W.P.; Abyzov, A.S. How Do Crystals Nucleate and Grow: Ostwalds Rule of Stages and Beyond. In Thermal Physics and Thermal Analysis, Hot Topics in Thermal Analysis and Calorimetry; Sestak, J., Hubik, P., Mares, J.J., Eds.; Springer International Publishing: Cham, Switzerland, 2017; Volume 11, pp. 195–211. [Google Scholar]
  68. Abyzov, A.S.; Davydov, L.N.; Schmelzer, J.W.P. Heterogeneous Nucleation in Solutions on Rough Solid Surfaces: Generalized Gibbs Approach. Entropy 2019, 21, 782. [Google Scholar] [CrossRef] [Green Version]
  69. Hillert, M. A Theory of Nucleation of Solid Metallic Solutions. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 1956. [Google Scholar]
  70. Cahn, J.W.; Hilliard, J.E. Free energy of a non-uniform system: I. Interfacial free energy. J. Chem. Phys. 1958, 28, 258–267. [Google Scholar] [CrossRef]
  71. Cahn, J.W.; Hilliard, J.E. Nucleation in a Two-Component Incompressible Fluid. J. Chem. Phys. 1959, 31, 688–699. [Google Scholar] [CrossRef]
  72. Van der Waals, J.D.; Kohnstamm, P. Lehrbuch der Thermodynamik (English Translation: Textbook of Thermodynamics); Johann Ambrosius Barth: Leipzig, Germany; Amsterdam, The Netherlands, 1908. [Google Scholar]
  73. Rowlinson, J.S.; van der Waals, J.D. The thermodynamic theory of capillarity under the hypothesis of a continuous variation of density. J. Stat. Phys. 1979, 20, 197–200. [Google Scholar] [CrossRef]
  74. Zanotto, E.D.; James, P.F. Experimental tests of the classical nucleation theory for glasses. J. Non-Cryst. Solids 1985, 74, 373–394. [Google Scholar] [CrossRef]
  75. Granasy, L.; James, P.F. Nucleation and growth in cluster dynamics: A quantitative test of the classical kinetic approach. J. Chem. Phys. 2000, 113, 9810–9821. [Google Scholar] [CrossRef]
  76. Fokin, V.M.; Zanotto, E.D.; Schmelzer, J.W.P. Method to estimate crystal/liquid surface energy by dissolution of subcritical nuclei. J. Non-Cryst. Solids 2000, 278, 24–34. [Google Scholar] [CrossRef]
  77. Fokin, V.M.; Potapov, O.V.; Zanotto, E.D.; Spiandorello, F.M.; Ugolkov, V.L.; Pevzner, B.Z. Mutant crystals in Na2O·2CaO·3SiO2-glasses. J. Non-Cryst. Solids 2003, 331, 240–253. [Google Scholar] [CrossRef]
  78. Schmelzer, J.W.P.; Gokhman, A.R.; Fokin, V.M. On the Dynamics of First-Order Phase Transitions in Multi-Component Systems. J. Colloid Interface Sci. 2004, 272, 109–133. [Google Scholar] [CrossRef]
  79. Tatchev, D.; Hoell, A.; Kranold, R.; Armyanov, S. Size distribution and composition of magnetic precipitates in amorphous Ni-P Alloys. Physica 2005, B 369, 8–19. [Google Scholar] [CrossRef]
  80. Tatchev, D.; Goerigk, G.; Valova, E.; Dille, J.; Kranold, R.; Armyanov, S.; Delplancke, J.-L. Investigation of the primary crystallization of Ni-17 at.% P alloy by ASAXS. J. Appl. Crystallogr. 2005, 38, 787–794. [Google Scholar] [CrossRef] [Green Version]
  81. Roskosz, M.; Toplis, M.J.; Richet, P. Kinetic vs. Thermodynamic Control of Crystal Nucleation and Growth in Molten Silicates. J. Non-Cryst. Solids 2005, 352, 180–184. [Google Scholar] [CrossRef]
  82. Fokin, V.M.; Zanotto, E.D. Continuous compositional changes of crystal and liquid during crystallization of a sodium calcium silicate glass. J. Non-Cryst. Solids 2007, 353, 2459–2468. [Google Scholar] [CrossRef]
  83. Schmelzer, J.W.P.; Abyzov, A.S. Generalized Gibbs’ Approach to the Thermodynamics of Heterogeneous Systems and the Kinetics of First-Order Phase Transitions. J. Eng. Thermophys. 2007, 16, 119–129. [Google Scholar] [CrossRef]
  84. Schmelzer, J.W.P.; Fokin, V.M.; Abyzov, A.S.; Zanotto, E.D.; Gutzow, I. How Do Crystals Form and Grow in Glass-Forming Liquids: Ostwalds Rule of Stages and Beyond. Int. J. Appl. Glass Sci. 2010, 1, 16–26. [Google Scholar] [CrossRef]
  85. Schmelzer, J.W.P.; Abyzov, A.S.; Fokin, V.M. Crystallization of glass: What we know, what we need to know. Int. J. Appl. Glass Sci. 2016, 7, 253–261. [Google Scholar] [CrossRef]
  86. Johari, G.P.; Schmelzer, J.W.P. Crystal Nucleation and Growth in Glass-forming Systems: Some New Results and Open Problems. In Glass: Selected Properties and Crystallization; Schmelzer, J.W.P., Ed.; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2014; pp. 521–585. [Google Scholar]
  87. Defay, R.; Prigogine, I.; Bellemans, R.; Everett, D.H. Surface Tension and Adsorption; Longmans, Green & Co.: London, UK, 1966. [Google Scholar]
  88. Prigogine, I.; Bellemans, A. Statistical Mechanics of Surface Tension and Adsorption. In Adhesion and Adsorption of Polymers Part A; Lee, L.H., Ed.; Plenum Press: New York, NY, USA, 1980. [Google Scholar]
  89. Rowlinson, J.S.; Widom, B. Molecular Theory of Capillarity; Clarendon Press: Oxford, UK, 1982. [Google Scholar]
  90. Boltachev, G.S.; Schmelzer, J.W.P. On the definition of temperature and its fluctuations in small systems. J. Chem. Phys. 2010, 133, 134509. [Google Scholar] [CrossRef]
  91. Schmelzer, J.W.P.; Boltachev, G.S.; Abyzov, A.S. On the temperature of the critical clusters in nucleation theory. J. Chem. Phys. 2013, 139, 034702. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Schmelzer, J.W.P. Application of the Nucleation Theorem to Crystallization of Liquids: Some General Theoretical Results. Entropy 2019, 21, 1147. https://doi.org/10.3390/e21121147

AMA Style

Schmelzer JWP. Application of the Nucleation Theorem to Crystallization of Liquids: Some General Theoretical Results. Entropy. 2019; 21(12):1147. https://doi.org/10.3390/e21121147

Chicago/Turabian Style

Schmelzer, Jürn W. P. 2019. "Application of the Nucleation Theorem to Crystallization of Liquids: Some General Theoretical Results" Entropy 21, no. 12: 1147. https://doi.org/10.3390/e21121147

APA Style

Schmelzer, J. W. P. (2019). Application of the Nucleation Theorem to Crystallization of Liquids: Some General Theoretical Results. Entropy, 21(12), 1147. https://doi.org/10.3390/e21121147

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