We introduce a method to obtain one-dimensional collective variables for studying rarely occurring transitions between two metastable states separated by a high free energy barrier. No previous information, not even approximated, on the path followed during the transition is needed. The only requirement is to know the fluctuations of the system while in the two metastable states. With this information in hand, we build the collective variable using a modified version of Fisher's linear discriminant analysis. The usefulness of this approach is tested on the metadynamics simulation of two representative systems. The first is the freezing of silver iodide into the superionic α-phase, and the second is the study of a classical Diels-Alder reaction. The collective variable works very well in these two diverse cases.

We introduce an extension of a recently published method [ Mendels et al. J. Phys. Chem. Lett. 2018 , 9 , 2776 - 2781 ] to obtain low-dimensional collective variables for studying multiple states free energy processes in chemical reactions. The only information needed is a collection of simple statistics of the equilibrium properties of the reactants and product states. No information on the reaction mechanism has to be given. The method allows studying a large variety of chemical reactivity problems including multiple reaction pathways, isomerization, and stereo- and regiospecificity. We applied the method to two fundamental organic chemical reactions. First we study the S2 nucleophilic substitution reaction of a Cl in CHCl leading to an understanding of the kinetic origin of the chirality inversion in such processes. Subsequently, we tackle the problem of regioselectivity in the hydrobromination of propene revealing that the nature of empirical observations such as Markovnikov's rule lies in the chemical kinetics rather than the thermodynamic stability of the products.

Several enhanced sampling methods such as umbrella sampling or metadynamics rely on the identification of an appropriate set of collective variables. Recently two methods have been proposed to alleviate the task of determining efficient collective variables. One is based on linear discriminant analysis, the other on a variational approach to conformational dynamics, and uses time-lagged independent component analysis. In this paper, we compare the performance of these two approaches in the study of the homogeneous crystallization of two simple metals. We focus on Na and Al and search for the most efficient collective variables that can be expressed as a linear combination of X-ray diffraction peak intensities. We find that the performances of the two methods are very similar. However, the method based on linear discriminant analysis, in its harmonic version, is to be preferred because it is simpler and much less computationally demanding.