Crystallization of Supercooled Liquids: Self-Consistency Correction of the Steady-State Nucleation Rate

Crystal nucleation can be described by a set of kinetic equations that appropriately account for both the thermodynamic and kinetic factors governing this process. The mathematical analysis of this set of equations allows one to formulate analytical expressions for the basic characteristics of nucleation, i.e., the steady-state nucleation rate and the steady-state cluster-size distribution. These two quantities depend on the work of formation, [Formula: see text] , of crystal clusters of size n and, in particular, on the work of critical cluster formation, [Formula: see text]. The first term in the expression for [Formula: see text] describes changes in the bulk contributions (expressed by the chemical potential difference, [Formula: see text]) to the Gibbs free energy caused by cluster formation, whereas the second one reflects surface contributions (expressed by the surface tension, [Formula: see text]: [Formula: see text] , [Formula: see text] , where [Formula: see text] is a parameter describing the size of the particles in the liquid undergoing crystallization), n is the number of particles (atoms or molecules) in a crystallite, and [Formula: see text] defines the size of the critical crystallite, corresponding to the maximum (in general, a saddle point) of the Gibbs free energy, G. The work of cluster formation is commonly identified with the difference between the Gibbs free energy of a system containing a cluster with n particles and the homogeneous initial state. For the formation of a “cluster” of size [Formula: see text] , no work is required. However, the commonly used relation for [Formula: see text] given above leads to a finite value for [Formula: see text]. By this reason, for a correct determination of the work of cluster formation, a self-consistency correction should be introduced employing instead of [Formula: see text] an expression of the form [Formula: see text]. Such self-consistency correction is usually omitted assuming that the inequality [Formula: see text] holds. In the present paper, we show that: (i) This inequality is frequently not fulfilled in crystal nucleation processes. (ii) The form and the results of the numerical solution of the set of kinetic equations are not affected by self-consistency corrections. However, (iii) the predictions of the analytical relations for the steady-state nucleation rate and the steady-state cluster-size distribution differ considerably in dependence of whether such correction is introduced or not. In particular, neglecting the self-consistency correction overestimates the work of critical cluster formation and leads, consequently, to far too low theoretical values for the steady-state nucleation rates. For the system studied here as a typical example (lithium disilicate, [Formula: see text]), the resulting deviations from the correct values may reach 20 orders of magnitude. Consequently, neglecting self-consistency corrections may result in severe errors in the interpretation of experimental data if, as it is usually done, the analytical relations for the steady-state nucleation rate or the steady-state cluster-size distribution are employed for their determination.