Critical radius

Critical radius is the minimum particle size from which an aggregate is thermodynamically stable. In other words, it is the lowest radius formed by atoms or molecules clustering together (in a gas, liquid or solid matrix) before a new phase inclusion (a bubble, a droplet or a solid particle) is viable and begins to grow. Formation of such stable nuclei is called nucleation.

At the beginning of the nucleation process, the system finds itself in an initial phase. Afterwards, the formation of aggregates or clusters from the new phase occurs gradually and randomly at the nanoscale. Subsequently, if the process is feasible, the nucleus is formed. Notice that the formation of aggregates is conceivable under specific conditions. When these conditions are not satisfied, a rapid creation-annihilation of aggregates takes place and the nucleation and posterior crystal growth process does not happen.

In precipitation models, nucleation is generally a prelude to models of the crystal growth process. Sometimes precipitation is rate-limited by the nucleation process. An example would be when someone takes a cup of superheated water from a microwave and, when jiggling it with a spoon or against the wall of the cup, heterogeneous nucleation occurs and most of water particles convert into steam.

If the change in phase forms a crystalline solid in a liquid matrix, the atoms might then form a dendrite. The crystal growth continues in three dimensions, the atoms attaching themselves in certain preferred directions, usually along the axes of a crystal, forming a characteristic tree-like structure of a dendrite.

Mathematical derivation

The critical radius of a system can be determined from its Gibbs free energy.[1]

${\displaystyle \Delta G_{T}=\Delta G_{V}+\Delta G_{S}}$

It has two components, the volume energy ${\displaystyle \Delta G_{V}}$ and the surface energy ${\displaystyle \Delta G_{S}}$. The first one describes how probable it is to have a phase change and the second one is the amount of energy needed to create an interface.

The mathematical expression of ${\displaystyle \Delta G_{V}}$ , considering spherical particles, is given by:

${\displaystyle \Delta G_{V}={\frac {4}{3}}\pi r^{3}\Delta g_{v}}$

where ${\displaystyle \Delta g_{v}}$ is the Gibbs free energy per volume and obeys ${\displaystyle -\infty <\Delta g_{v}<\infty }$. It is defined as the energy difference between one system at a certain temperature and the same system at the fusion temperature and it depends on pressure, the number of particles and temperature: ${\displaystyle \Delta g_{v}(T,p,N)}$. For a low temperature, far from the fusion point, this energy is big (it is more difficult to change the phase) and for a temperature close to the fusion point it is small (the system will tend to change its phase).

Regarding ${\displaystyle \Delta G_{S}}$ and considering spherical particles, its mathematical expression is given by:

Free energy change versus the nanoparticle radius. Below the critical radius, the clusters are not big enough to start the nucleation process. The Gibbs free energy change is positive and the process is not prosperous. This critical radius corresponds to the minimum size at which a particle can survive in solution without being redissolved. Above the critical radius, the particles will form and grow as it is thermodinamically favourable.

${\displaystyle \Delta G_{S}=4\pi r^{2}\gamma >0}$

where ${\displaystyle \gamma }$ is the surface tension we need to break to create a nucleus. The value of the ${\displaystyle \Delta G_{S}}$ is never negative as it always takes energy to create an interface.

The total Gibbs free energy is therefore:

${\displaystyle \Delta G_{T}=-{\frac {4\pi }{3}}r^{3}\Delta g_{v}+4\pi r^{2}\gamma }$

The critical radius ${\displaystyle r_{c}}$ is found by optimization, setting the derivative of ${\displaystyle \Delta G_{T}}$ equal to zero.

${\displaystyle {\frac {d\Delta G_{T}}{dr}}=-4\pi r_{c}^{2}\Delta g_{v}+8\pi r_{c}\gamma =0}$

yielding

${\displaystyle r_{c}={\frac {2\gamma }{|\Delta g_{v}|}}}$,

where ${\displaystyle \gamma }$ is the surface tension and ${\displaystyle |\Delta g_{v}|}$ is the absolute value of the Gibbs free energy per volume.

The Gibbs free energy of nuclear formation is found replacing the critical radius expression in the general formula.

${\displaystyle \Delta G_{c}={\frac {16\pi \gamma ^{3}}{3(\Delta g_{v})^{2}}}}$

Interpretation

When the Gibbs free energy change is positive, the nucleation process will not be prosperous. The nanoparticle radius is small, the superficial term prevails the volum term ${\displaystyle \Delta G_{S}>\Delta G_{V}}$. Contrary, if the variation rate is negative, it will be thermodynamically stable. The size of the cluster surpasses the critical radius. In this occasion, the volum term overcomes the superficial term ${\displaystyle \Delta G_{S}<\Delta G_{V}}$.

From the expression of the critical radius, as the Gibbs volume energy increases, the critical radius will decrease and hence, it will be easier achieving the formation of nuclei and begin the crystallization process.

Example of reducing the critical radius

In order to decrease the value of the critical radius ${\displaystyle r_{c}}$ and promote nucleation, a supercooling or superheating process may be used.

Supercooling is a phenomenon in which the system's temperature is lowered under the phase transition temperature without the creation of the new phase. Let ${\displaystyle \Delta T=T_{f}-T}$ be the temperature difference, where ${\displaystyle T_{f}}$ is the phase transition temperature. Let ${\displaystyle \Delta g_{v}=\Delta h_{v}-T\Delta s_{v}}$ be the volume Gibbs free energy, enthalpy and entropy respectively.

When ${\displaystyle T=T_{f}}$, the system has null Gibbs free energy, so:

${\displaystyle \Delta g_{f,v}=0\Leftrightarrow \Delta h_{f,v}=T_{f}\Delta s_{f,v}}$

In general, the following approximations can be done:

${\displaystyle \Delta h_{v}\rightarrow \Delta h_{f,v}}$ and ${\displaystyle \Delta s_{v}\rightarrow \Delta s_{f,v}}$

Consequently:

${\displaystyle \Delta g_{v}\simeq \Delta h_{f,v}-T\Delta s_{f,v}=\Delta h_{f,v}-{\frac {T\Delta h_{f,v}}{T_{f}}}=\Delta h_{f,v}{\frac {T_{f}-T}{T_{f}}}}$

So:

${\displaystyle \Delta g_{v}=\Delta h_{f,v}{\frac {\Delta T}{T_{f}}}}$

Subtituting this result on the expressions for ${\displaystyle r_{c}}$ and ${\displaystyle \Delta G_{c}}$, the following equations are obtained:

${\displaystyle r_{c}={\frac {2\gamma T_{f}}{\Delta h_{f,v}}}{\frac {1}{\Delta T}}}$

${\displaystyle \Delta G_{c}={\frac {16\pi \gamma ^{3}T_{f}^{2}}{3(\Delta h_{f,v})^{2}}}{\frac {1}{(\Delta T)^{2}}}}$

Notice that ${\displaystyle r_{c}}$ and ${\displaystyle \Delta G_{c}}$ diminish with an increasing supercooling. Analogously, a mathematical derivation for the superheating can be done.

See also

• Nucleation
• Homogeneous nucleation
• Heterogeneous nucleation
• Ostwald ripening
• Supercooling
• Superheating

References

1. "Crystallization Kinetics". Retrieved 16 August 2018.

• N.H.Fletcher, Size Effect in Heterogeneous Nucleation, J.Chem.Phys.29, 1958, 572.
• Nguyen T. K. Thanh,* N. Maclean, and S. Mahiddine, Mechanisms of Nucleation and Growth of Nanoparticles in Solution, Chem. Rev. 2014, 114, 15, 7610-7630.