Crystallization --> Nucleation

Formulas for numbers of molecules on the surface of a large nucleus.

by Leonid Sakharov

Numerical simulation of nucleation phenomena during crystallization demands a definition several parameters of the clusters of crystalline molecules in initial amorphous phase. Among of such key parameters are a numbers of molecules on the surface of cluster these can be transformed from one phase to other changing number of molecules in the cluster.

Let’s cluster has i molecules with volume vm each.

Effective size of d is defined by equation vm = d3

The volume of the cluster is:

V[i] = (4π/3) * r3 = i * vm

where r is effective radius of the cluster.

The volumes of surface areas where molecules can move from crystal into feeding environment Vs- or into cluster Vs-  from surrounding the nucleus phase are defined by formulas:

Vs- = (4π/3) * (r3 – (r – g*d)3)


Vs+ = (4π/3) * ((r + g*d)3 - r3)

where g is geometrical coefficient that has meaning of effective thickness of surface layer measured in size of molecule. Let also introduce a constant p = 4π/3 = 4.188790205 that will help to simplify an equations.

Number of molecules on surface of the cluster can be calculated by expressions:

Ns- = Vs-/vm

Ns+ = Vs+/vm

Following are series of mathematical transformations to present Ns- and Ns+  in most convenient for calculation form:


Ns- = (p/vm) * ( r3 – r3 + 3 r2*gd – 3 r*g2d2 +g3d3)

Ns+ = (p/vm) * (r3 + 3 r2*gd + 3 r*g2d2 +g3d3  - r3)

Defining radius of nucleus from number of molecules:

r = p-1/3 * i1/3*vm1/3

and substituting it in equation for number molecules on surface of nucleus the following formula appears:

Ns- = (p/vm) * (3 p-2/3 * i2/3*vm2/3*g*vm1/3 – 3 p-1/3 * i1/3*vm2/3*g2 vm2/3  +g3 vm)

Ns+ = (p/vm) * (3 p-2/3 * i2/3*vm2/3*g*vm1/3 + 3 p-1/3 * i1/3*vm2/3*g2 vm2/3  +g3 vm)

And finally getting equations:

Ns- = 3 p1/3 * i2/3*g – 3 p2/3 * i1/3*g2  +g3

Ns+ = 3 p1/3 * i2/3*g + 3 p2/3 * i1/3*g2  +g3

Taking that the value of geometrical coefficient could be safely estimated having value 0.85 formulas above can be presented as


Ns- = 4.1*i2/3 – 5.6*i1/3  + 0.6

Ns+ = 4.1*i2/3 + 5.6*i1/3  + 0.6

It is important to point out on the remarkable fact that formulas above do not include a volume of molecule.

Next figure shows numbers of molecules on surface of cluster calculated by formulas above.


There is significant difference between number molecules in the side of cluster and outside it these belong to amorphous phase. This difference is specially is pronounced for small clusters creating pure geometrical advantage for agglomeration of molecules.

Submitted at Aug. 18, 2010; 15:16


Post your comment:

 Name:           Terms and policy.

Home   Our software   Software archive   Expert database   Forum   Photo album   Partners   About  
Copyright © by LeoKrut