Modelling nucleation from solution with the string method in the osmotic ensemble

Abstract

ABSTRACTDirect molecular simulation of nucleation from solution is a challenging task, requiring a combination of ‘rare events’ techniques and methods to control the chemical potential. Rare event methods usually keep the total number of molecules fixed, resulting in artificial free energy profiles due to the depletion of solute from the solution phase. In order to address this issue, we present a new approach that uses the string method in collective variables in the osmotic ensemble to obtain minimum free energy pathways for nucleation at constant supersaturation. Our method does not require using an explicit reservoir of solute molecules, or making additional assumptions about the activity coefficients in the solution. We apply the new method to the crystallisation of sulfamerazine from acetonitrile and methanol solutions, and compare the resulting potential of mean force profiles to those obtained using analytical corrections previously employed in the literature.