#### Abstract

The method to get the glass media with necessary optical characteristics is proposed. This method is based on inserting a necessary number of heterogeneous centers of given activity into the system. The theoretical description of the nucleation process in such situation is given and the resulting formulas allow choosing the characteristics of heterogeneous centers necessary to get the required optical characteristics of the media.

#### 1. Introduction

The theory of a vitreous state inevitably includes the nucleation phenomena [1]. Every vistreous medium contains the embryos of a crystal phase which are growing extremely slowly. So one can say that they are practically stable objects. Such embryos are the centers of scattering and they affect the optical properties of glass. The start of systematic scattering theory was given by Lord Rayleigh [2]; the scattering on fluctuations was considered by Smoluchowsky and Einstein [3, 4]. Essential contribution in the consideration of the scattering with account of magnetic properties was made by Cabannes [5, 6]. A bright example of the theory of scattering on particles of essential sizes is given by the famous Mie scattering theory [7]. This approach was essentially generalized by Jobst [8] and Debye [9].

In every scattering theory the distribution of the embryos sizes is reflected in the optical coefficients of material. Ordinary it is necessary to know only several first momentums of the distribution of the embryos over sizes. Sometimes it is sufficient to know only the mean size of the embryos.

Impurities are inevitably presented in the media and they are like some heterogeneous centers. The nucleation theory on heterogeneous centers was investigated already by Frenkel in [10] and this paper completes the classical theory of nucleation created by Becker and Doering [11], Volmer [12], Kramers [13], and Zeldovich [14]. The full analysis of the stationary nucleation rate in the case of heterogeneous condensation was given by Kuni [15]. Even without the solution of kinetic equation one can evidently see that the main factor in the rate of nucleation is the Boltzmann distribution , where is the free energy of the critical embryo counted in the thermal units. This value is the central object of the theory and the main efforts were paid to find a correct expression for . Unfortunately it is rather difficult to do and even in the simplest case of the charged nuclei it is necessary to build a rather complex theory with decompositions [16–18]. These decompositions were constructed in frames of the general Gibbs dividing surfaces theory [19]. In the formalism of these decompositions there appear several coefficients with very specific physical meaning.

During the process of the glass preparation it is necessary to have the glass with the necessary optical properties. So it is necessary to have impurities in the necessary quantity and of the given size.

Around the embryos there appear specific profiles of density of the substance. The role of the density profiles in the light propagation has been studied in [20]. The analogous property in frames of the embryos in nucleation theory was investigated in [21].

The classical theory of nucleation [22, 23] is like an elegant necklace which has to be around the free energy of critical embryo. The problem here is not only in the necessity to know what the surface tension has to be put in these constructions [24] and to describe the surface layer [25], but also to know the influence of the heterogeneous center on these objects. Even in the case without the heterogeneous center the question of the free energy profile is one of the most actual questions [26]. Considering the possibility of the direct measurement of the free energy of embryos one has to state that except special situations (see [27]) there is no clear way to do it. The direct application of the statistical mechanical approach [28] also causes certain difficulties. The same can be said about the comparison with the 2D Ising model as it is proposed in [29]. The density functional approach [30–33] also can not give the absolutely precise result. But the accuracy of the last approach is rather satisfactory. One has to mention also the possibility of using the approximate scaling formulas [34–36] based on some characteristic basic points for decompositions as parameters. One can also insert here the generalizations associated with the dependence of the surface tension on the size [37, 38], the effects in constrained systems [39], and the deformation of heterogeneous centers [40]. Here the possibility of fluctuations [41] and the Tolman-like corrections [42] is not considered. The refinements of [43, 44] can be introduced here directly.

The most evident way to get the necessary properties is to govern the law of creation of the ideal metastability , that is, the imaginary metastability which would be in the system when the process of the new phase formation is forbidden. The value of the ideal metastability is fully governed externally and it seems that this function can be regulated to have the necessary properties. Unfortunately, the preparation of the glass requires the absolutely prescribed conditions to have the glass with optimal properties and it is impossible to change the conditions of the glass preparation. So as a function of time is supposed to be fixed.

But one can propose another rather simple and effective method. The main idea is to inject a given number of heterogeneous centers of the given activity to change the effective conditions of the nucleation of the main quantity of impurities which occurs pseudo-homogeneously (i.e., without exhaustion of the potential sites). The theoretical description of this method will be the matter of the current paper. This paper is organized as follows: (i)At first the* properties of the nucleation rate* will be recalled. The exponential approximation for the nucleation rate as a function of a supersaturation will be introduced. (ii)To construct the description of the multistage process which is the nucleation with the presence of the active heterogeneous impurities at first the* nucleation without heterogeneous impurities* will be studied. This process is more simple but the extraction of the characteristic features will help in construction of the more complex process. (iii)Having started the description of the nucleation process with the presence of the active heterogeneous impurities one can consider the separate process of exhaustion of these impurities or the* formation of the embryos on the active heterogeneous impurities*. This part of consideration will give the parameters of the influence of the embryos on heterogeneous centers on the process of formation of the pseudo-homogeneous embryos of a crystal phase. (iv)Under the influence of heterogeneously formed embryos the* formation of the pseudo-homogeneous embryos* will be studied. Namely, this part of consideration will give the necessary optical characteristics of the media. (v)To get the necessary optical characteristics it is necessary to formulate the* conditions on parameters of the active heterogeneous centers*. This is the matter of the next section. Here the ways to solve the system of algebraic equations of parameters are also discussed. (vi)The final section of the paper is Conclusions where the results are summarized and the limitations and restrictions of the theory are described.

#### 2. Properties of the Nucleation Rate

The rate of nucleation is proportional to exponent of the free energy of the critical embryo. This factor is the main factor and that is why the free energy of the critical embryo is the main object of interest in the investigations devoted to the nucleation.

The general formula for the rate of nucleation is the following:In (1) is the number of sites (or the molecules) which can be the starting point for the nucleation (crystallization) formation of the embryo and the factor is the Zeldovich factor. The value of is given by the numerous considerations but the approach appearing in the classical theory of nucleation is not yet radically reconsidered. Although there are some attempts to refine the Zeldovich factor these refinements do not lead to the change of this value in the order of the magnitude. Some interesting ideas for appear in the multidimensional consideration [45, 46] but this is not the subject of the current paper.

It is reasonable to extract in the kinetic factor (the number of adsorbed molecules in the time unit) and to writespeaking about the truncated Zeldovich factor because contrary to the liquid-vapor nucleation or the solid-liquid nucleation in the situation considered here (i.e., in the process of vitrification) the kinetic factor has specific behavior: at some moment of time it becomes so small that even the supercritical embryos practically stop growing. Otherwise it is impossible to explain why the crystal embryos do not gradually transform the whole volume of the glass into the set of crystals.

In the traditional approach mentioned above it is supposed that the intensity of absorption of new molecules by the embryo iswhere is the density of the molecules in the vicinity of the embryo, is the mean thermal velocity of the molecule, is the condensation coefficient, and is the surface area of the embryo.

The last formula (3) is supposed to be valid both for the critical and for the supercritical embryos. The evolution of undercritical embryos is not too interesting because in the region of undercritical embryos there is the quasi-equilibrium state.

Really, formula (3) has a rather general range of application. Then it follows that the peculiarities of are caused by the peculiarities of . This explains the factorization (2). The condensation coefficient contains the main dependence specific for the vitrification process. Ordinary it is presented asimplying that in the process of adsorption there is a necessity to overcome some barrier. The height of barrier is denoted by . This barrier can have a very specific origin including the barriers for displacement of the molecules in the environment of the embryo in order to give enough space for the appropriate installation of the new molecule in the embryo (this seems to be the main reason for annulation of ). So only some general remarks can be made concerning (4).

Instead of the previous formula (4) it is better to writeimplying that even when the thermodynamic conditions correspond to the absence of imaginary barrier there remain other factors (like thermal relaxation, noncorrespondence of the position of the molecule and the profile of the adsorbing region of the embryo, etc.) which affects the probability of adsorption. Note that barrier in (5) is an abstract effective barrier.

In the liquid-vapor nucleation it is reasonable to suppose that the explicit dependence of on the time is rather smooth. In the vitrification process has to go to zero when grows. It is rather easy to explain why it is so. It is supposed that ordinary there are no thermal effects of nucleation and the creation of is attained by the cooling of the system. Then the vapor consumption leads to the decrease of the number of the molecules but the temperature is not affected by embryos. The decrease of ordinary leads to increase of the activation barriers. If it is supposed that , where is the activation barrier of accommodation, then we need the increase of with decrease of temperature. Really, it is quite natural because the activation barrier is measured in thermal units. Then it can be seen that goes to zero.

If is linearized as a function of then it is possible to come towith the linearization coefficient . The linearization (6) invokes a question of the dependence of on time. Here it is also supposed that the linearization can be made and one can come towith parameters and . Parameter can be put to zero by the appropriate shift of the time scale.

Certainly, it is necessary to have some concrete dependence in order to present some concrete formulas. Analogous theory can be constructed for some other concrete dependence of on .

Certainly can not be greater than . The exponential approximation does not correspond to this property. A more simple and effective approximation is the following: for and for . Namely, this approximation will be used. Fortunately it is possible to determine , at least to estimate in order to ensure the correct final size of embryos.

One has to stress that appears both in expression for and in expression for the rate of growth. So after there will be no appearance of the supercritical droplets nor growth of the already existing droplets. The process is therefore stopped.

Since it is required that the system is metastable it means that it is supposed that without active heterogeneous centers there will be some embryos of the new phase. This means thatwhere is the moment of intensive formation of the pseudo-homogeneous embryos. Certainly the moment in (8) of characteristic time of appearance of embryos on the active heterogeneous centers has to be smaller than . Then

One can prove that the moment of the pseudo-homogeneous formation in the presence of active heterogeneous centers satisfies inequalityThenInequalities (8), (9), (10), and (11) presented above are rather important because they allow stating that the process of nucleation can be considered at the practically constant value of .

A special question which has to be analyzed here is whether the rate of nucleation is really the stationary one. Justifications of the stationarity are quite similar to the liquid-vapor case.

The mechanism of the embryos growth can be chosen as the free molecular one which corresponds to the case of crystallization. This leads to , where is the mean concentration of the molecules in the noncrystalline phase.

The rate of nucleation is the function of the power of metastability which is called the supersaturation. Ordinary is presented aswhere is the concentration of molecules in the noncrystalline phase at the state of the phase coexistence. Certainly thenwhere is the number of the molecules in the crystalline phase calculated in units of . Because of the growth of the crystals the value grows also. At the value stops to grow but the process of the new droplets formation stops also.

One can see that according to (12), (13) the supersaturation is not governed purely externally but depends on the process of nucleation and on the process of the substance consumption and the media heating by the regular growth of supercritical embryos. So there appears the self-consistent problem which is the subject of investigation in the theories of the global nucleation.

To see how depends on it is reasonable to notice that the main dependence occurs through the factor . The free energy is the smooth regular function of thermodynamic parameters and of and it is reasonable to linearize it on the deviation in . It leads towhere is some base for decomposition and is a parameter of linearization. For in (14) one can getwhere is the chemical potential and is the number of the molecules in the critical embryo. The function is rather smooth and takes moderate values. For of the mother phase system like the ideal gas (i.e., for every mother phase where the mean field approximation works) one can easily come toSo according to (15), (16) one can estimate as and it is possible to see that gets the big values ( to ensure the thermodynamic description of the critical embryo). It means that the dependence of on is rather sharp.

Later the mentioned estimate leads to the estimateduring the period of intensive nucleation. Here is the value of at some characteristic moment which belongs to the nucleation period.

Also one can get the estimate analogous to (17) for and on the base of (18) one can see that alsoThe last estimate shows that the nucleation period is relatively short in time. Namely, estimate (19) allows speaking above about , , and as some moments of time for nucleation.

To see that the exponential approximation really works one can present the characteristic behavior of the Zeldovich factor for the pseudo-homogeneous case. It is shown in Figure 1. One can see here two lines: the practically constant line which is the reduced Zeldovich factor and the exponent which is the main Boltzmann factor. So one can see that the preexponential factor (Zeldovich factor) does not really change. Here the worst value for the parameter is and the worst value for the renormalized surface tension is (the biggest possible value in the normal conditions). Namely, these values are taken for this example. Certainly, one can not include here the explicit time dependence of the kinetic coefficient of adsorption.

Really, the relative change of the preexponential factor here has the order . It can be hardly seen in Figure 1 because the variation is very small.

#### 3. Nucleation without Heterogeneous Impurities

To calculate one has to know not only the rate of nucleation but also the rate of the embryos growth. Fortunately the free molecular regime leads to the following rather simple rate of growth of the number of the molecules inside the supercritical embryo:with characteristic parameter . Certainly,Namely, (21) shows that the time of freezing stops the growth.

One can prove that in calculation of only the input from the supercritical embryos has to be taken into account.

Fortunately according to (20) is rather unsensitive to the supersaturation and there is no need to observe the influence of the nucleation on the rate of the embryos growth. Then for the embryo appeared at the number of molecules will be approximatelyand with account of (17) one can transform (22) into

On the base of (23) the expression for will be

Then balance equation (13) on the base of (24) can be written as

Having inserted the mentioned linearizations and exponential approximation (14), (15) for it is possible to come in (25) toExpression (26) has enough high accuracy.

One can prove thatEstimate (27) allows taking instead of and writing (26) asHaving shifted the time to have the zero point at one can come instead of (28) to

The last equation (29) is typical for nucleation and is formulated, for example, in [47] as one of the limit cases. It can be solved by iterationsThe lower index in brackets indicates the number of iteration. The initial approximation in the iteration procedure (30) is

The first approximation calculated on the base of (31) isSinceit is easy to come in (32) on the base of (33) toEquation (34) has enough high accuracy.

One can prove that already the first approximation is rather good and the size spectrumis rather accurate.

The rate of nucleation in the first (see (35)) and the second approximations of the iteration procedure is drawn in Figure 2. One can see here two lines: the lower line is the first iteration and the upper line is the second iteration for the nucleation rate as a function of the special shifted and scaled time.

One can see that the first iteration is rather close to the second one. One can see the similarity of the forms of the first and the second iteration. If it is necessary to increase the accuracy one can make a simple shift.

One can choose as the moment of the maximal supersaturation. Thenis the algebraic equation on parameters (actually on ) and the good approximation to solve (36) isApproximation (37) is rather accurate.

All what has been done is necessary to show that at , that is, , can be estimated asAlthough (38) is rather simple it will be absolutely sufficient for further constructions.

Now the further evolution will be considered. One can see there two periods. The first period is the rapid consumption of the surplus substanceIt is finished at some moment of time when the whole (or practically all) surplus substance is in the embryos. Now one has to estimate . It is possible to see on the base of (39) thatSo according to (40) this period is rather rapid. It allows supposing that the condition also ensuresThe probability that (41) is not observed is lower than .

Later the supersaturation is fallen and there will be an asymptotic period where . Here the mean effective linear coordinate of the spectrum grows very slow like . At it stops to grow. It is reasonable to putand (42) gives the natural condition for the mean linear size of the embryos.

Here one has to make one important notation. In the liquid-vapor transition the volume per molecule in a vapor phase is a thousand times greater than the volume per particle in a liquid phase. This ensures that the relative volume of the new phase in vapor-liquid transition is small.

In the opposite liquid-vapor transition it seems that the situation will be opposite and the whole volume begins to be the vapor phase. But here the main effect of the bubble formation is to take away the surplus stretching. It is necessary to keep in mind.

In the situation of the liquid-crystal transition the volume per the molecule in a crystal phase is comparable with and it is possible to see the crystallization of the whole volume (ordinary this does not take place because of the heat release effects). But in the process of vitrification the mother phase is so viscous that the main effect of the appearance of the crystal embryo is to take away the stretching and the tension in the mother media. Then the relative volume of the new crystal phase will be small at .

Later without the crystal phase would cover practically the whole volume. But the cut-off of the nucleation and condensation process at will stop the consumption of the volume by the crystal phase. Unfortunately one can not radically change and . That is why the method of heterogeneous centers insertion is suggested.

#### 4. Formation of the Embryos on the Active Heterogeneous Impurities

Now the process of the heterogeneous condensation will be studied. It is supposed that the centers are relatively active. This means that all of them are exhausted in the process of nucleation. This property simplifies the theoretical description of kinetics of the phase formation.

One can note that the attribution of the centers (sites) to the group of the active ones depends on the intensity of the change of the external conditions. So this property is not absolute.

It is supposed also that the embryos on the heterogeneous centers can be described thermodynamically. It means that the number of molecules in the critical embryo is big enough. As for the number of the molecules inside the equilibrium center one can avoid the restriction necessary for the thermodynamic description of the equilibrium embryo. It can be done since only the behavior of the height of the activation barrieris interesting and it determines the results of nucleation. Here the variation of because of uncertainty of deviation of the embryos description from the thermodynamic one will be small in comparison with .

The behavior of the heterogeneous embryo free energy together with the behavior of the homogeneous free energy is drawn in Figure 3. The solid line is the homogeneous free energy and the dashed line is the heterogeneous free energy. Here the capillary modelfor is used for the homogeneous case. As in Figure 1 the values and are used. In the heterogeneous caseHere is used as a value comparable with . Certainly, this example (namely, (44), and (45)) is no more than a model.

Since one can prove thatwhere is a chemical potential, one can see that according to (43), (46) only when is the value essential in . So when is small one can not describe the equilibrium embryo adequately but the dependence of on (which means the dependence on ) is not important.

For heterogeneous nucleation the nucleation rate is given bywhere is the number of the free heterogeneous centers (unoccupied by the supercritical embryos) and comes from the equilibrium distributionSo the structure of the rate of nucleation (47), (48) remains the same; it is the equilibrium distribution formally prolonged until the critical embryo and multiplied by the kinetic factor and by the truncated Zeldovich factor .

Instead of the general theory (see [47]) one can simply note that all centers will be exhausted when the ideal number of the formed droplets which is calculated as where is the total number of heterogeneous centers, is equal to . The requirement of essential exhaustion of heterogeneous centers together with (49) gives an equationon .

To solve (50) more effectively it is reasonable to introduce an approximationwhereThe accuracy of (51), (52) is rather high.

One can easily estimate as

Then (50) can be rewritten with the help of (53) asand givesorEquation (56) does not contain . Certainly (54) and (54) do not contain it also. So the power of exhaustion and the level of exhaustion are purely determined by the nature of centers and do not depend on their quantity. This allows speaking about the active centers. It is possible to see that having changed the number of centers one can change the numbers of the embryos which begin to consume metastability rather actively after .

The detailed justification of the presented method is rather long but it exists at the level of precise analytical derivation.

#### 5. Formation of the Pseudo-Homogeneous Embryos

Now the formation of the pseudo-homogeneous embryos will be studied. The difference from the consideration of the homogeneous nucleation presented above will be the change of conditions for the new process because there are the heterogeneous consumers of metastability.

Certainly, the new forming embryos consume metastability and this leads to the perturbation of the rate of growth of the supercritical embryos formed on heterogeneous centers. But one can prove that this influence is very small and can be neglected. So one can see two separate problems here. The first problem is to determine the supersaturation which appears after the heterogeneous embryos have been formed. The second problem is to determine the characteristics of the embryos size spectrum when the homogeneous process takes place.

The balance equation with account of only heterogeneous centers is written aswhere is the coordinate of the linear size of embryos and it can be found asSo (57), (58) lead to the first-order differential equation which can be solved. One can reduce (57) with the help of (58) towith parameters and and get the solution of (59):with the arbitrary constant which has to be chosen to satisfy the initial condition

So then (61) leads to the concrete value of the constant in (60) which will be the following:Expression (62) has to be inserted into (60).

Unfortunately it is impossible to inverse the dependence given by (60) and get the explicit solution. The dependence for , is drawn as an example in Figure 4.

One has to note that here it is impossible to use an approximationanalogous to (7) but it is necessary instead of (63) to integrate explicitly.

Now the behavior of will be analyzed qualitatively. The difference at first increases in time very rapidly (faster than ) and this leads to the maximum of even without homogeneous appearance of embryos. After attaining the maximum of the supersaturation it decreases and the conditions for the intensive homogeneous formation of new embryos disappear. So there exist two possibilities for the further homogeneous nucleation.

The first possibility is to see the nucleation in conditions of the growing which isHere in (64) the homogeneously formed embryos take part in the metastability consumption rather actively.

Alternative possibility is the following: the supersaturation fully coincides with and homogeneously formed embryos do not play any essential role in evolution. Here the number of embryos isand the embryos linear size spectrum iswherefor andfor . The accuracy of expressions (65), (66), (67), and (68) is sufficient.

One can use the exponential approximation based on to calculate the number of droplets and the size spectrum. Also one can use the power decomposition of at the point of maximum.

The second possibility is not interesting also from the practical point of view because it requires essential number of heterogeneous centers (sites).

The only necessary fact which remains to be clarified is the existence of the time of the end of the rapid growth of the embryos sizes. Really, the rapid increase of can be stopped only when comes to zero. But it has to be stopped because can not be greater than . So exists.

How can one calculate and the linear size at the end of intensive growth (here these values are marked as and )? There is no sharp sensitivity of these values to the supersaturation and one can get these values by the simple iteration procedureExpressions (69) present the first step of this procedure.

The time is determined as to satisfy the following equation:The value isExpressions (70) and (71) have the necessary accuracy.

It is also necessary to determine the mean size for the homogeneously formed embryos. It iswhere is the point of the maximum of . The last integral can be easily taken analytically since the integrand is the power function. Then (72) becomes a simple algebraic formula.

#### 6. Conditions on Parameters of the Active Heterogeneous Centers

Now the first possibility will be considered. Here consideration is similar to the case of the pure homogeneous nucleation but it is necessary to consider instead of .

The function is determined and it is possible to approximately linearize it at the point and to getwhereFormula (73) (with account of (74)) is quite analogous to (7).

The balance equation isand is counted from .

The last equation (75) can be solved by iterationsThe lower index in brackets which appeared in (76) indicates the number of iteration. The initial approximation will be

The first approximation on the base of (77) isOne can prove that already the first approximation (78) is rather good and gives the size spectrumThe form of the size spectrum given by (79) is rather accurate.

One can choose as the moment of the maximal supersaturation. Thenis the algebraic equation on . The solution of (80) is approximatelyFormula (81) solves the problem to find parameters of the spectrum.

Now the further evolution will be considered. Since already ensured , here one can also see which will be marked as . One can show thatif it is supposed that is equal to . So here alsoand one can see on the base of (82), (83) that practically ensures .

The final values will beEquations (84) demonstrate the necessary accuracy.

The value is calculated aswhereApproximation of (85), (86) has enough high accuracy.

The value is calculated aswhere is the total number of embryos formed homogeneously in units of . The value appearing in (87) is already determined. If one takes an approximationone can get the following expression for :Expressions (88), (89) complete the recipe of calculations.

#### 7. Conclusions

The known value of allows getting the total number of the homogeneously formed embryos as The integration can be performed analytically which brings (90) toSimple formula (91) allows calculating .

Changing one can choose . Changing one can choose the force of the action of heterogeneous centers and then choose . Then one can get necessary and .

The second possibility when the supersaturation is determined by the heterogeneously formed embryos is not too attractive because here the number of heterogeneous centers is greater than in the first variant and can be even greater than the number of homogeneously formed embryos. The only positive property is that in the majority of situationsand one can hope thatSo the process of nucleation is frozen only after the end of essential stages of nucleation. The probability of violation of (92), (93) is practically negligible.

In the first variant the number is many times smaller than :Then one can simply say that the optical properties are determined according to (94) by and . Now the method of how to regulate these values is presented.

The most serious limitation of the theory is the free molecular regime of growth. It leads to the collective consumption of metastability. But when the quantity of the heterogeneous centers is small the region of influence of every heterogeneous center is large enough and it is impossible to keep the free molecular regime here. Fortunately there is a certain analogy between the situation of the collective consumption and the situation with the zones of depletion. This analogy was established in [48]. The mentioned fact was seen for the period of nucleation. Here the situation is more complex because there is no appearance of new embryos after . So here there is a nucleation in the system with some very curved boundaries (the boundaries of exhausted regions) and these boundaries move rather slowly. But although the situation is another it is possible to see that the maximum of can be taken in the first iteration, that is, on the base of the metastability consumption which occurs by the ideal law, that is, in approximation . This approximation is good for all with , that is, for convex functions. But for there appears another stabilizing factor: the rate of the supersaturation growth rapidly increases and the consumption of metastability increases. It means that the boundaries begin to move rather rapidly and the error in the position of these boundaries is not too important. So the form of the free volume can be smoothed and it is possible to come to the absence of the effects of the noncollective consumption. In any case this question forms a matter of a separate publication.

#### List of Variables

Surface area of the embryo | |

: | Parameter in the equation for the evolution after the heterogeneous nucleation stops |

: | Surplus chemical potential |

: | Parameter in the equation for the evolution after the heterogeneous nucleation stops |

: | Parameter in the model for the free energy of the embryo on heterogeneous center |

: | Index to mark the critical embryo |

: | Index to mark to equilibrium embryo |

: | Free energy of the embryo measured in the thermal units ( is the Boltzmann constant and is the absolute temperature) |

: | Size spectrum of the embryos |

: | Number of the molecules in the crystalline phase calculated in units of |

: | Rate of nucleation |

: | Condensation coefficient |

: | Basis for the approximation for |

: | Coefficient in linearization of |

: | Number of molecules in a new phase at the end of intensive growth |

: | Mean size for the homogeneously formed embryos |

: | Number of heterogeneously formed droplets calculated in approximation of the total number of the free heterogeneous centers |

: | Mean concentration of the molecules in the noncrystalline phase |

: | Number of sites (or the molecules) which can be the starting point for the nucleation (crystallization) formation of the embryo |

: | Equilibrium distribution |

: | Number of the free heterogeneous centers (unoccupied by the supercritical embryos) |

: | Density of the molecules in the vicinity of the embryo |

: | Total number of heterogeneous centers |

: | Concentration of molecules in the noncrystalline phase at the state of the phase coexistence |

: | Linear size of the embryo |

: | Absolute temperature |

: | Time |

: | Time of the “freezing” of the system, which is a parameter in approximation for kinetic coefficient |

: | Moment of intensive formation of the pseudo-homogeneous embryos |

: | Time of appearance of embryos on the active heterogeneous centers |

: | Moment of the pseudo-homogeneous formation in the presence of active heterogeneous centers |

: | Time of the end of the rapid growth of the embryos sizes |

: | Mean thermal velocity of the molecule |

: | Kinetic factor (intensity of adsorption of the molecules by the embryo) |

: | Zeldovich factor |

: | Truncated Zeldovich factor |

: | Coordinate of the spectrum of the heterogeneously formed embryos |

: | Height of the activation barrier of nucleation |

: | Width of the equilibrium zone |

: | Parameter of decompositions for the free energy, which ordinary equals the total number of the molecules inside the homogeneous critical embryo or the difference between the number of molecules in the critical embryo and the number of molecules in the equilibrium embryo |

: | Parameter for heterogeneous nucleation |

: | Chemical potential |

: | Number of molecules in the embryo |

: | Number of molecules in the critical embryo |

: | Number of molecules in the equilibrium embryo |

: | Activation barrier of accommodation |

: | Imaginary metastability which would be in the system when the process of the new phase formation is forbidden |

: | Basis for linearization of |

: | Coefficient in linearization of |

: | Supersaturation that would be in a system when the heterogeneous formation is allowed and the homogeneous mechanism is forbidden |

: | Parameter of linearization of |

: | A cubic root of |

: | Renormalized surface tension |

: | Characteristic time, parameter in the law of growth, which ordinary equals the mean time between collisions of some given molecule of the saturated disordered phase with other molecules in this phase |

: | Power of metastability or the supersaturation |

: | Basis for decompositions for |

: | Index denoting the values at the period of nucleation. |

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.