Evaluating the predictive character of the method of constrained geometries simulate external force with density functional theory

Mechanochemistry is a fast-developing field of interdisciplinary research with a growing number of applications. Therefore, many theoretical methods have been developed to quickly predict the outcome of mechanically induced reactions. Constrained geometries simulate External Force (CoGEF) is one of the earlier methods in this field. It is easily implemented and can be conducted with most DFT codes. However, recently, we observed totally different predictions for model systems of epoxy resins in different conformations and with different density functionals. To better understand the conformational and functional dependence in typical CoGEF calculations we present a systematic evaluation of the CoGEF method for different model systems covering homolytic and heterolytic bond cleavage reactions, electrocyclic ring opening reactions and scission of noncovalent interactions in hydrogen-bond complexes. From our calculations we observe that many mechanochemical descriptors strongly depend on the functional used, however, a systematic trend exists for the relative maximum Force. In general, we observe that the CoGEF procedure is forcing the system to high energetic regions on the molecular potential energy profiles, which can lead to unexpected and uncorrelated predictions of mechanochemical reactions. This is questioning the true predictive character of the method.

Automated summary

Mechanochemistry is a rapidly growing interdisciplinary research field that involves predicting the outcome of mechanically induced reactions. In this study, the CoGEF method was systematically evaluated for different model systems, revealing strong dependence on conformation and density functional. The findings question the predictive nature of the method.