Free energies of molecular clusters determined by guided mechanical disassembly

Calculating free energies of organic molecules on insulating substrates

The challenges and limitations in calculating free energies and entropies of adsorption and interaction of organic molecules on an insulating substrate are discussed. The adhesion of 1,3,5-tri(4'-cyano-[1,1'-biphenyl]-4-yl)benzene (TCB) and 1,4-bis(4-cyanophenyl)-2,5-bis(decyloxy)benzene (CDB) molecules to step edges on the KCl(001) surface and the formation of molecular dimers were studied using classical molecular dynamics. Both molecules contain the same anchoring groups and benzene ring structures, yet differ in their flexibility. Therefore, the entropic contributions to their free energy differ, which affects surface processes. Using potential of mean force and thermodynamic integration techniques, free energy profiles and entropy changes were calculated for step adhesion and dimer formation of these molecules. However, converging these calculations is nontrivial and comes at large computational cost. We illustrate the difficulties as well as the possibilities of applying these methods towards understanding dynamic processes of organic molecules on insulating substrates.

Water droplet excess free energy determined by cluster mitosis using guided molecular dynamics

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.