Nucleation and superstabilization in small systems

Oiling-Out: Insights from Gibbsian Surface Thermodynamics by Stability Analysis

The crystallization process is a significant stage in the pharmaceutical industry. During the process of crystallization with cooling, it is possible for a secondary liquid phase to appear before the formation of crystals. This phenomenon is called "oiling out" or liquid-liquid phase separation (LLPS). In this article, we explore the oiling-out phenomenon in a binary system of water and vanillin using stability analysis based on Gibbsian surface thermodynamics. To obtain the full picture of oiling out, we investigated three cases: droplet-solute-lean liquid equilibrium (DLE), crystal-solute-rich liquid equilibrium (CL'E), and crystal-solute-lean liquid equilibrium (CLE). The phase diagram of the system is plotted using the NRTL model for activity coefficients, along with considering the effect of the interfacial curvature on the phase diagram. From the phase boundaries and free-energy diagram of each case, we showed that the occurrence of the oiling-out phenomenon is justified based on the lower energy barrier of the droplet formation compared to that of the crystal formation. However, the energy level of a stable crystal is significantly lower and hence more stable than that of a stable droplet. Finally, we have determined different regions for droplet and crystal formation in the metastable phase diagram based on their supersaturation and provide insight for the oiling-out phenomenon.

On the thermodynamics of curved interfaces and the nucleation of hard spheres in a finite system

We determine, for hard spheres, the Helmholtz free energy of a liquid that contains a solid cluster as a function of the size of the solid cluster by means of the formalism of the thermodynamics of curved interfaces. This is done at the constant total number of particles, volume, and temperature. We show that under certain conditions, one may have several local minima in the free energy profile, one for the homogeneous liquid and others for the spherical, cylindrical, and planar solid clusters surrounded by liquid. The variation of the interfacial free energy with the radius of the solid cluster and the distance between equimolar and tension surfaces are inputs from simulation results of nucleation studies. This is possible because stable solid clusters in the canonical ensemble become critical in the isothermal-isobaric ensemble. At each local minimum, we find no difference in chemical potential between the phases. At local maxima, we also find equal chemical potential, albeit in this case the nucleus is unstable. Moreover, the theory allows us to describe the stable solid clusters found in simulations. Therefore, we can use it for any combination of the total number of particles, volume, and global density as long as a minimum in the Helmholtz free energy occurs. We also study under which conditions the absolute minimum in the free energy corresponds to a homogeneous liquid or to a heterogeneous system having either spherical, cylindrical, or planar geometry. This work shows that the thermodynamics of curved interfaces at equilibrium can be used to describe nucleation.

Nucleation in a liquid droplet

Nucleation of multi-microclusters in a droplet.