Predicting transition state geometries is one of the most challenging tasks in computational chemistry, which often requires expert-based knowledge and permanent human intervention. This short communication reports technical details and preliminary results of a python-based tool (AMADAR) designed to generate any Diels-Alder (DA) transition state geometry (TS) and analyze determined IRC paths in a (quasi-)automated fashion, given the product SMILES. Two modules of the package are devoted to performing, from IRC paths, reaction force analyses (RFA) and atomic (fragment) decompositions of the reaction force F and reaction force constant [Formula: see text]. The performance of the protocol has been assessed using a dataset of 2000 DA cycloadducts retrieved from the ZINC database. The sequential location of the corresponding TSs was achieved with a success rate of 95%. RFA plots confirmed the reaction force constant [Formula: see text] to be a good indicator of the (non)synchronicity of the associated DA reactions. Moreover, the atomic decomposition of [Formula: see text] allows for the rationalization of the (a)synchronicity of each DA reaction in terms of contributions stemming from pairs of interacting atoms. The source code of the AMADAR tool is available on GitHub [ CMCDD/AMADAR(github.com) ] and can be used directly with minor customizations, mostly regarding the local working environment of the user.
Keywords: Diels–Alder reaction; Finite difference approach; Hellman–Feynman force; RDKit; Reaction force analysis.
© 2022. The Author(s).