471740 Polymorphic Transformations at Ambient Conditions Elucidated with Molecular Dynamics
In this work, we present the prediction of finite temperature solid-solid transformations using temperature and pressure dependent free energies derived from a novel molecular dynamics approach. Molecular dynamics simulations implicitly include the thermal motions present at finite temperatures and can therefore predict the emergence of a more stable crystal structure at high temperatures and pressures. We examine the five organic molecules benzene, formamide, glycine, acetic acid, and imidazole which are known experimentally to have different globally stable structures at low and high temperatures. The relative free energies for these solids in the point-charge OPLS potential correctly identify the high temperature / high pressure polymorphs as being more stable than the low temperature form.
We also compute the temperature and pressure dependent free energies in the more expensive polarizable AMOEBA potential to determine the added benefit of using a more accurate energy function. We use a novel Boltzmann reweighting approach to circumvent the additional expense of the polarizable potential which has kept the more accurate model from being used in previous MD studies. The polarizable AMOEBA potential is shown to have a higher accuracy in predicting the low temperature crystal structure in these systems. However, both OPLS and AMOEBA have equal fidelity at predicting the high temperature and high pressure structure. These preliminary results suggest that high level energy functions are more accurate at predicting low temperature structures, but have a less significant advantage at high temperatures relative to cheaper levels of theory. Finally, we show that these short ambient temperature MD simulations can allow metastable lattice minima to restructure into more stable minima.
See more of this Group/Topical: Computational Molecular Science and Engineering Forum