Search for: All records

Abstract Dynamic motion often controls selectivity in reactions featuring two consecutive potential‐energy transition states. Here we report density functional theory (DFT)‐based direct dynamics trajectories and machine learning classification analysis for cyclopentadienone dimerization and a N 2 extrusion reaction leading to semibullvalene. These reactions have consecutive transition states, and there is dynamic selectivity that determines which of two possible C‐C bonds is formed after the first transition state. For cyclopentadienone dimerization with a bispericyclic first transition state, machine learning analysis using transition‐state based features provided >90% trajectory classification accuracy, but only using AdaBoost and random forest algorithms. Many other relatively sophisticated machine learning algorithms showed poor accuracy despite the obvious motion responsible for selectivity. Feature importance analysis confirmed that the sigmatropic rearrangement vibrational motion in the bispericyclic transition state provides prediction of which of the second C‐C bonds is dynamically formed. For the reaction leading to semibullvalene, machine learning analysis provides solid accuracy for classifying trajectories and predicting which C‐C bond is formed and which C‐C bond is broken immediately after N 2 ejection. Like the cyclopentadienone dimerization reaction, machine learning feature importance analysis showed that the sigmatropic rearrangement vibrational motion in the N 2 extrusion transition state determines which C‐C bond is formed and which is broken. Surprisingly, machine learning struggles to predict which trajectories undergo a subsequent [3,3] sigmatropic rearrangement process, which isomerizes equivalent forms of semibullvalene. 

Quasiclassical trajectory analysis is now a standard tool to analyze non-minimum energy pathway motion of organic reactions. However, due to the large amount of information associated with trajectories, quantitative analysis of the dynamic origin of reaction selectivity is complex. For the electrocyclic ring opening of cyclopropyl radical, more than 4000 trajectories were run showing that allyl radicals are formed through a mixture of disrotatory intrinsic reaction coordinate (IRC) motion as well as conrotatory non-IRC motion. Geometric, vibrational mode, and atomic velocity transition-state features from these trajectories were used for supervised machine learning analysis with classification algorithms. Accuracy >80% with a random forest model enabled quantitative and qualitative assessment of transition-state trajectory features controlling disrotatory IRC versus conrotatory non-IRC motion. This analysis revealed that there are two key vibrational modes where their directional combination provides prediction of IRC versus non-IRC motion. 

The transfer of a β-hydrogen from a metal-alkyl group to ethylene is a fundamental organometallic transformation. Previously proposed mechanisms for this transformation involve either a two-step β-hydrogen elimination and migratory insertion sequence with a metal hydride intermediate or a one-step concerted pathway. Here, we report density functional theory (DFT) quasiclassical direct dynamics trajectories that reveal new dynamical mechanisms for the β-hydrogen transfer of [Cp*Rh III (Et)(ethylene)] + . Despite the DFT energy landscape showing a two-step mechanism with a Rh–H intermediate, quasiclassical trajectories commencing from the β-hydrogen elimination transition state revealed complete dynamical skipping of this intermediate. The skipping occurred either extremely fast (typically <100 femtoseconds (fs)) through a dynamically ballistic mechanism or slower through a dynamically unrelaxed mechanism. Consistent with trajectories begun at the transition state, all trajectories initiated at the Rh–H intermediate show continuation along the reaction coordinate. All of these trajectory outcomes are consistent with the Rh–H intermediate <1 kcal mol −1 stabilized relative to the β-hydrogen elimination and migratory insertion transition states. For Co, which on the energy landscape is a one-step concerted mechanism, trajectories showed extremely fast traversing of the transition-state zone (<50 fs), and this concerted mechanism is dynamically different than the Rh ballistic mechanism. In contrast to Rh, for Ir, in addition to dynamically ballistic and unrelaxed mechanisms, trajectories also stopped at the Ir–H intermediate. This is consistent with an Ir–H intermediate that is stabilized by ∼3 kcal mol −1 relative to the β-hydrogen elimination and migratory insertion transition states. Overall, comparison of Rh to Co and Ir provides understanding of the relationship between the energy surface shape and resulting dynamical mechanisms of an organometallic transformation.