Critical radius explained

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]

\Delta G _{T} = \Delta G _{V} + \Delta G _{S}

It has two components, the volume energy

\Delta G _{V} and the surface energy \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

\Delta G _{V} , considering spherical particles, is given by:\Delta G _{V} = (4)/(3) \pi r ^{3} \Delta g _{v}

where

\Delta g _{v} is the Gibbs free energy per volume and obeys -inf ty < \Delta g _{v} < inf ty . 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: \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

\Delta G _{S} and considering spherical particles, its mathematical expression is given by:\Delta G _{S} = 4 \pi r ^{2} \gamma > 0

where

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

The total Gibbs free energy is therefore:

\Delta G _{T =- (4 \pi)/(3)} r ^{3} \Delta g _{v} + 4 \pi r ^{2} \gamma

The critical radius

r _{c} is found by optimization, setting the derivative of \Delta G _{T} equal to zero.(d\Delta G _{T})/(dr) =-4 \pi r ^{2} _{c} \Delta g _{v+} 8 \pi r _{c} \gamma = 0

yielding

r _{c} = (2 \gamma)/(|\Delta g _{v|}) ,

where

\gamma is the surface tension and |\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.

\Delta G _{c} = (16 \pi\gamma ^{3})/(3 (\Delta g ^{2} _{v)})

Interpretation

When the Gibbs free energy change is positive, the nucleation process will not be prosperous. The nanoparticle radius is small, the surface term prevails the volume term

\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 volume term overcomes the superficial term \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.

Methods for reducing the critical radius

Supercooling

In order to decrease the value of the critical radius

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

\Delta T = T _{f} - T be the temperature difference, where T _{f} is the phase transition temperature. Let \Delta g _{v} = \Delta h _{v} - T \Delta s _{v} be the volume Gibbs free energy, enthalpy and entropy respectively.

When

T = T _{f} , the system has null Gibbs free energy, so:\Delta g _{f,v} = 0 \Leftrightarrow \Delta h _{f,v} = T _{f} \Delta s _{f,v}

In general, the following approximations can be done:

\Delta h _{v} → \Delta h _{f,v} and \Delta s _{v} → \Delta s _{f,v}

Consequently:

\Delta g _{v} \simeq \Delta h _{f,v} - T\Delta s _{f,v} = \Delta h _{f,v} - (T\Delta h _{f,v})/(T _{f}) = \Delta h _{f,v} (T _{f} - T)/(T _{f})

So:

\Delta g _{v} = \Delta h _{f,v} (\Delta T)/(T _{f})

Substituting this result on the expressions for

r _{c} and \Delta G _{c} , the following equations are obtained:r _{c} = (2 \gamma T _{f})/(\Delta h _{f,v}) (1)/(\Delta T) \Delta G _{c} = (16 \pi \gamma ^{3} T ^{2} _{f})/(3 (\Delta h _{f,v} ) ^{2}) (1)/((\Delta T) ^{2})

Notice that

r _{c} and \Delta G _{c} diminish with an increasing supercooling. Analogously, a mathematical derivation for the superheating can be done.

Supersaturation

Supersaturation is a phenomenon where the concentration of a solute exceeds the value of the equilibrium concentration.

\Delta \mu = - k _{B} T ln \left( (c _{0})/(c _{eq}) \right) , where k _{B} is the Boltzmann constant, c _{0} is the solute concentration and c _{eq} is the equilibrium concentration. For a stoichiometric compound and considering \mu = (\partial G)/(\partial N) and N = (V)/(v _{a}) , where v _{a} is the atomic volume:\Delta g _{v} = (\Delta \mu)/(v _{a}) = - (k _{B} T)/(v _{a}) ln \left( (c _{0})/(c _{eq}) \right).

Defining the supersaturation as

S = (c _{0-c} _{eq})/(c _{eq}) , this can be rewritten as\Delta g _{v} = - (k _{B} T)/(v _{a}) ln \left( 1+ S\right).

Finally, the critical radius

r _{c} and the Gibbs free energy of nuclear formation \Delta G _{c} can be obtained asr _{c} = (2 \gamma v _{a})/(k _{B} T ln \left( 1+ S\right)) ,\Delta G _{c} = (16 \pi \gamma ^{3} V ^{2} _{M})/(3 (RT ln \left( 1+ S\right)) ^{2}) ,

where

V _{M} is the molar volume and R is the molar gas constant.

See also

• Nucleation
• Ostwald ripening
• Supercooling
• Superheating

References

• 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.

Notes and References

1. Web site: Crystallization Kinetics. 16 August 2018. 13 August 2018. https://web.archive.org/web/20180813184136/http://community.dur.ac.uk/sharon.cooper/lectures/cryskinetics/handoutsalla.html. dead.