Identifying Slow Molecular Motions in Complex Chemical Reactions

We have studied the cyclization reaction of deprotonated 4-chloro-1-butanethiol to tetrahydrothiophene by means of well-tempered metadynamics. To properly select the collective variables, we used the recently proposed variational approach to conformational dynamics within the framework of metadyanmics. This allowed us to select the appropriate linear combinations from a set of collective variables representing the slow degrees of freedom that best describe the slow modes of the reaction. We performed our calculations at three different temperatures, namely, 300, 350, and 400 K. We show that the choice of such collective variables allows one to easily interpret the complex free-energy surface of such a reaction by univocal identification of the conformers belonging to reactants and product states playing a fundamental role in the reaction mechanism.