Related Experiment Video
Updated: Mar 28, 2026

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry
Published on: April 8, 2020
Water droplet excess free energy determined by cluster mitosis using guided molecular dynamics
Gabriel V Lau1, Patricia A Hunt2, Erich A Müller1
1Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom.
Abstract:
Atmospheric aerosols play a vital role in affecting climate by influencing the properties and lifetimes of clouds and precipitation. Understanding the underlying microscopic mechanisms involved in the nucleation of aerosol droplets from the vapour phase is therefore of great interest. One key thermodynamic quantity in nucleation is the excess free energy of cluster formation relative to that of the saturated vapour. In our current study, the excess free energy is extracted for clusters of pure water modelled with the TIP4P/2005 intermolecular potential using a method based on nonequilibrium molecular dynamics and the Jarzynski relation. The change in free energy associated with the "mitosis" or division of a cluster of N water molecules into two N/2 sub-clusters is evaluated. This methodology is an extension of the disassembly procedure used recently to calculate the excess free energy of argon clusters [H. Y. Tang and I. J. Ford, Phys. Rev. E 91, 023308 (2015)]. Our findings are compared to the corresponding excess free energies obtained from classical nucleation theory (CNT) as well as internally consistent classical theory (ICCT). The values of the excess free energy that we obtain with the mitosis method are consistent with CNT for large cluster sizes but for the smallest clusters, the results tend towards ICCT; for intermediate sized clusters, we obtain values between the ICCT and CNT predictions. Furthermore, the curvature-dependent surface tension which can be obtained by regarding the clusters as spherical droplets of bulk density is found to be a monotonically increasing function of cluster size for the studied range. The data are compared to other values reported in the literature, agreeing qualitatively with some but disagreeing with the values determined by Joswiak et al. [J. Phys. Chem. Lett. 4, 4267 (2013)] using a biased mitosis approach; an assessment of the differences is the main motivation for our current study.
More Related Videos
Related Concept Videos
Excess Pressure Inside a Drop and a Bubble
Phase Transitions: Vaporization and Condensation
Intermolecular Forces
Intermolecular Forces
States of Water
Water freezes when the intermolecular forces are greater than the kinetic energy. Unlike most other substances, water is less dense in its solid state than in its liquid state. This is because each water molecule can form...
Surface Tension of Fluid
Surface tension varies...

