Automated discovery of chemically reasonable elementary reaction steps

Computational framework for discovery of degradation mechanisms of organic flow battery electrolytes

The stability of organic redox-active molecules is a key challenge for the long-term viability of organic redox flow batteries (ORFBs). Electrolyte degradation leads to capacity fade, reducing the efficiency and lifespan of ORFBs. To systematically investigate degradation mechanisms, we present a computational framework that automates the exploration of degradation pathways. The approach integrates local reactivity descriptors to generate reactive complexes and employs a single-ended process search to discover elementary reaction steps, including transition states and intermediates. The resulting reaction network is iteratively refined with heuristics and human-guided validation. The framework is applied to study the degradation mechanisms of quinone- and quinoxaline-based electrolytes under acidic and basic aqueous conditions. The predicted reaction pathways and degradation products align with experimental observations, highlighting key degradation modes such as Michael addition, disproportionation, dimerization, and electrochemical transformation. The framework provides a valuable tool for in silico screening of stable electrolyte candidates and guiding the molecular design of next-generation ORFBs.

Facilitating Transition State Search with Minimal Conformational Sampling Using Reaction Graph

Elucidating transition states (TSs) is crucial for understanding chemical reactions. The reliability of traditional TS search approaches depends on input conformations that require significant effort to prepare. Previous automated methods for generating input reaction conformations typically involve extensive exploration of a large conformational space. Such exhaustive search can be complicated by the rapid growth of the conformational space, especially for reactions involving many rotatable bonds, multiple reacting molecules, and numerous bond formations and dissociations. To address this problem, we propose a new approach that generates reaction conformations for TS searches with minimal reliance on sampling. This method constructs a pseudo-TS structure based on a reaction graph containing bond formation and dissociation information and modifies it to produce reactant and product conformations. Tested on three different benchmarks, our method consistently generated suitable conformations without necessitating extensive sampling, demonstrating its potential to significantly improve the applicability of automated TS searches. This approach offers a valuable tool for a broad range of applications such as reaction mechanism analysis and network exploration.

Chemical Reaction Networks from Scratch with Reaction Prediction and Kinetics-Guided Exploration

Algorithmic reaction explorations based on transition state searches can now routinely predict relatively short reaction sequences involving small molecules. However, applying these algorithms to deeper chemical reaction network (CRN) exploration still requires the development of more efficient and accurate exploration policies. Here, an exploration algorithm, which we name yet another kinetic strategy (YAKS), is demonstrated that uses microkinetic simulations of the nascent network to achieve cost-effective, deep network exploration. Key features of the algorithm are the automatic incorporation of bimolecular reactions between network intermediates, compatibility with short-lived but kinetically important species, and incorporation of rate uncertainty into the exploration policy. In validation case studies of glucose pyrolysis, the algorithm rediscovers reaction pathways previously discovered by heuristic exploration policies and elucidates new reaction pathways for experimentally obtained products. The resulting CRN is the first to connect all major experimental pyrolysis products to glucose. Additional case studies are presented that investigate the role of reaction rules, rate uncertainty, and bimolecular reactions. These case studies show that naïve exponential growth estimates can vastly overestimate the actual number of kinetically relevant pathways in the physical reaction networks. In light of this, further improvements in exploration policies and reaction prediction algorithms make it feasible that CRNs might soon be routinely predictable in some contexts.

LASP to the Future of Atomic Simulation: Intelligence and Automation

Atomic simulations aim to understand and predict complex physical phenomena, the success of which relies largely on the accuracy of the potential energy surface description and the efficiency to capture important rare events. LASP software (large-scale atomic simulation with a Neural Network Potential), released in 2018, incorporates the key ingredients to fulfill the ultimate goal of atomic simulations by combining advanced neural network potentials with efficient global optimization methods. This review introduces the recent development of the software along two main streams, namely, higher intelligence and more automation, to solve complex material and reaction problems. The latest version of LASP (LASP 3.7) features the global many-body function corrected neural network (G-MBNN) to improve the PES accuracy with low cost, which achieves a linear scaling efficiency for large-scale atomic simulations. The key functionalities of LASP are updated to incorporate (i) the ASOP and ML-interface methods for finding complex surface and interface structures under grand canonic conditions; (ii) the ML-TS and MMLPS methods to identify the lowest energy reaction pathway. With these powerful functionalities, LASP now serves as an intelligent data generator to create computational databases for end users. We exemplify the recent LASP database construction in zeolite and the metal-ligand properties for a new catalyst design.

Discovery of New Synthetic Routes of Amino Acids in Prebiotic Chemistry

The origin of life on Earth remains one of the most perplexing challenges in biochemistry. While numerous bottom-up experiments under prebiotic conditions have provided valuable insights into the spontaneous chemical genesis of life, there remains a significant gap in the theoretical understanding of the complex reaction processes involved. In this study, we propose a novel approach using a roto-translationally invariant potential (RTIP) formulated with pristine Cartesian coordinates to facilitate the simulation of chemical reactions. By employing RTIP pathway sampling to explore the reactivity of primitive molecules, we identified several low-energy reaction mechanisms, such as two-hydrogen-transfer hydrogenation and HCOOH-catalyzed hydration and amination. This led to the construction of a comprehensive reaction network, illustrating the synthesis pathways for glycine, serine, and alanine. Further thermodynamic analysis highlights the pivotal role of formaldimine as a key precursor in amino acid synthesis, owing to its more favorable reactivity in coupling reactions compared to the traditionally recognized hydrogen cyanide. Our study demonstrates that the RTIP methodology, coupled with a divide-and-conquer strategy, provides new insights into the simulation of complex reaction processes, offering promising applications for advancing organic design and synthesis.

Reproducing Reaction Mechanisms with Machine-Learning Models Trained on a Large-Scale Mechanistic Dataset

Mechanistic understanding of organic reactions can facilitate reaction development, impurity prediction, and in principle, reaction discovery. While several machine learning models have sought to address the task of predicting reaction products, their extension to predicting reaction mechanisms has been impeded by the lack of a corresponding mechanistic dataset. In this study, we construct such a dataset by imputing intermediates between experimentally reported reactants and products using expert reaction templates and train several machine learning models on the resulting dataset of 5,184,184 elementary steps. We explore the performance and capabilities of these models, focusing on their ability to predict reaction pathways and recapitulate the roles of catalysts and reagents. Additionally, we demonstrate the potential of mechanistic models in predicting impurities, often overlooked by conventional models. We conclude by evaluating the generalizability of mechanistic models to new reaction types, revealing challenges related to dataset diversity, consecutive predictions, and violations of atom conservation.

TS-tools: Rapid and automated localization of transition states based on a textual reaction SMILES input

Here, TS-tools is presented, a Python package facilitating the automated localization of transition states (TS) based on a textual reaction SMILES input. TS searches can either be performed at xTB or DFT level of theory, with the former yielding guesses at marginal computational cost, and the latter directly yielding accurate structures at greater expense. On a benchmarking dataset of mono- and bimolecular reactions, TS-tools reaches an excellent success rate of 95% already at xTB level of theory. For tri- and multimolecular reaction pathways - which are typically not benchmarked when developing new automated TS search approaches, yet are relevant for various types of reactivity, cf. solvent- and autocatalysis and enzymatic reactivity - TS-tools retains its ability to identify TS geometries, though a DFT treatment becomes essential in many cases. Throughout the presented applications, a particular emphasis is placed on solvation-induced mechanistic changes, another issue that received limited attention in the automated TS search literature so far.

POMSimulator: An open-source tool for predicting the aqueous speciation and self-assembly mechanisms of polyoxometalates

Elucidating the speciation (in terms of concentration versus pH) and understanding the formation mechanisms of polyoxometalates remains a significant challenge, both in experimental and computational domains. POMSimulator is a new methodology that tackles this problem from a purely computational perspective. The methodology uses results from quantum mechanics based methods to automatically set up the chemical reaction network, and to build speciation models. As a result, it becomes possible to predict speciation and phase diagrams, as well as to derive new insights into the formation mechanisms of large molecular clusters. In this work we present the main features of the first open-source version of the software. Since the first report [Chem. Sci. 2020, 11, 8448-8456], POMSimulator has undergone several improvements to keep up with the growing challenges that were tackled. After four years of research, we recognize that the source code is sufficiently stable to share a polished and user-friendly version. The Python code, manual, examples, and install instructions can be found at https://github.com/petrusen/pomsimulator.

Synthesis of Macrolactone Core of ent-Formosalide A via Regioselective Ether Cyclization

Formosalide A is a cytotoxic macrolide isolated from the dinoflagellate Prorocentrum sp, whose structure is characterized by functionalized 5- and 6-membered ether rings embedded in the macrolactone and an all cis-tetraene side chain. Here, we report the synthesis of the macrolactone core of ent-formosalide A. Our approach is highlighted by the Au-mediated 6-exo-dig cyclization for the synthesis of the 6-membered ether ring, which proceeded in a highly regioselective manner. Control experiments demonstrated that the acyclic protecting group of the C9,C10-diol was crucial for the desired 6-exo-dig cyclization. Theoretical studies were performed focusing on structural component analysis, which suggested that the C8-C9-C10-C11 dihedral angle induced by the protecting group controlled the regioselectivity. An additional 6 steps including Shiina macrolactone formation from the 6-membered ether ring completed the synthesis of the macrolactone core of ent-formosalide A.

Graph to Activation Energy Models Easily Reach Irreducible Errors but Show Limited Transferability

Activation energy characterization of competing reactions is a costly but crucial step for understanding the kinetic relevance of distinct reaction pathways, product yields, and myriad other properties of reacting systems. The standard methodology for activation energy characterization has historically been a transition state search using the highest level of theory that can be afforded. However, recently, several groups have popularized the idea of predicting activation energies directly based on nothing more than the reactant and product graphs, a sufficiently complex neural network, and a broad enough data set. Here, we have revisited this task using the recently developed Reaction Graph Depth 1 (RGD1) transition state data set and several newly developed graph attention architectures. All of these new architectures achieve similar state-of-the-art results of ∼4 kcal/mol mean absolute error on withheld testing sets of reactions but poor performance on external testing sets composed of reactions with differing mechanisms, reaction molecularity, or reactant size distribution. Limited transferability is also shown to be shared by other contemporary graph to activation energy architectures through a series of case studies. We conclude that an array of standard graph architectures can already achieve results comparable to the irreducible error of available reaction data sets but that out-of-distribution performance remains poor.

Green Light Promoted Iridium(III)/Copper(I)-Catalyzed Addition of Alkynes to Aziridinoquinoxalines Through the Intermediacy of Azomethine Ylides

This manuscript describes the development of alkyne addition to the aziridine moiety of aziridinoquinoxalines using dual Ir(III)/Cu(I) catalytic system under green light-emitting diode (LED) photolysis (λmax =525 nm). This mild method features high levels of chemo- and regioselectivity and was used to generate 30 highly functionalized substituted dihydroquinoxalines in 36-98 % yield. This transformation was also carried asymmetrically using phthalazinamine-based chiral ligand to provide 9 chiral addition products in 96 : 4 to 86 : 14 e.r. The experimental and quantum chemical explorations of this reaction suggest a mechanism that involves Ir(III)-catalyzed triplet energy transfer followed by a ring-opening reaction ultimately leading to the formation of azomethine ylide intermediates. These azomethine intermediates undergo sequential protonation/copper(I) acetylide addition to provide the products.

Exploration of biochemical reactivity with a QM/MM growing string method

In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.

Toward De Novo Catalyst Discovery: Fast Identification of New Catalyst Candidates for Alcohol-Mediated Morita-Baylis-Hillman Reactions

Recently we have demonstrated how a genetic algorithm (GA) starting from random tertiary amines can be used to discover a new and efficient catalyst for the alcohol-mediated Morita-Baylis-Hillman (MBH) reaction. In particular, the discovered catalyst was shown experimentally to be eight times more active than DABCO, commonly used to catalyze the MBH reaction. This represents a breakthrough in using generative models for catalyst optimization. However, the GA procedure, and hence discovery, relied on two important pieces of information; 1) the knowledge that tertiary amines catalyze the reaction and 2) the mechanism and reaction profile for the catalyzed reaction, in particular the transition state structure of the rate-determining step. Thus, truly de novo catalyst discovery must include these steps. Here we present such a method for discovering catalyst candidates for a specific reaction while simultaneously proposing a mechanism for the catalyzed reaction. We show that tertiary amines and phosphines are potential catalysts for the MBH reaction by screening 11 molecular templates representing common functional groups. The method relies on an automated reaction discovery workflow using meta-dynamics calculations. Combining this method for catalyst candidate discovery with our GA-based catalyst optimization method results in an algorithm for truly de novo catalyst discovery.

Automated Exploration of Reaction Networks and Mechanisms Based on Metadynamics Nanoreactor Simulations

We developed an automated approach to construct a complex reaction network and explore the reaction mechanisms for numerous reactant molecules by integrating several theoretical approaches. Nanoreactor-type molecular dynamics was used to generate possible chemical reactions, in which the metadynamics was used to overcome the reaction barriers, and the semiempirical GFN2-xTB method was used to reduce the computational cost. Reaction events were identified from trajectories using the hidden Markov model based on the evolution of the molecular connectivity. This provided the starting points for further transition-state searches at the electronic structure levels of density functional theory to obtain the reaction mechanism. Finally, the entire reaction network containing multiple pathways was built. The feasibility and efficiency of the automated construction of the reaction network were investigated using the HCHO and NH3 biomolecular reaction and the reaction network for a multispecies system comprising dozens of HCN and H2O molecules. The results indicated that the proposed approach provides a valuable and effective tool for the automated exploration of the reaction networks.

Efficient Acceleration of Reaction Discovery in the Ab Initio Nanoreactor: Phenyl Radical Oxidation Chemistry

Over the years, many computational strategies have been employed to elucidate reaction networks. One of these methods is accelerated molecular dynamics, which can circumvent the expense required in dynamics to find all reactants and products (local minima) and transition states (first-order saddle points) on a potential energy surface (PES) by using fictitious forces that promote reaction events. The ab initio nanoreactor uses these accelerating forces to study large chemical reaction networks from first-principles quantum mechanics. In the initial nanoreactor studies, this acceleration was done through a piston periodic compression potential, which pushes molecules together to induce entropically unfavorable bimolecular reactions. However, the piston is not effective for discovering intramolecular and dissociative reactions, such as those integral to the decomposition channels of phenyl radical oxidation. In fact, the choice of accelerating forces dictates not only the rate of reaction discovery but also the types of reactions discovered; thus, it is critical to understand the biases and efficacies of these forces. In this study, we examine forces using metadynamics, attractive potentials, and local thermostats for accelerating reaction discovery. For each force, we construct a separate phenyl radical combustion reaction network using solely that force in discovery trajectories. We elucidate the enthalpic and entropic trends of each accelerating force and highlight their efficiency in reaction discovery. Comparing the nanoreactor-constructed reaction networks with literature renditions of the phenyl radical combustion PES shows that a combination of accelerating forces is best suited for reaction discovery.

Combined molecular dynamics and coordinate driving method for automatically searching complicated reaction pathways

The combined molecular dynamics and coordinate driving (MD/CD) method is updated and generalized in this work to broaden its applications in automatically searching reaction pathways for complicated reactions. In this updated version, MD simulations are performed with the GFN's family of methods to systematically sample conformers for almost any systems with atomic numbers Z ≤ 86. The improved CD procedure is greatly accelerated by applying a pre-screening stage at the semiempirical GFN2-xTB level. An automatic module based on the Marcus theory and its improved version (the Wolynes theory) is designed to include single electron transfer (SET) processes into reaction pathways. The capabilities of this method are demonstrated by exploring the most possible reaction pathways of three experimentally reported reactions: the organophosphine-catalyzed trans phosphinoboration, the Fe(II) complex-mediated C(sp2)-H borylation reaction, and the SET-triggered deaminative radical cross-coupling reaction. Comprehensive reaction networks are obtained for all three reactions with reasonable computational costs. Detailed mechanisms for these reactions can account for the reported experimental facts.

Multi-Time-Scale Simulation of Complex Reactive Mixtures: How Do Polyoxometalates Form?

Understanding the dynamics of reactive mixtures still challenges both experiments and theory. A relevant example can be found in the chemistry of molecular metal-oxide nanoclusters, also known as polyoxometalates. The high number of species potentially involved, the interconnectivity of the reaction network, and the precise control of the pH and concentrations needed in the synthesis of such species make the theoretical/computational treatment of such processes cumbersome. This work addresses this issue relying on a unique combination of recently developed computational methods that tackle the construction, kinetic simulation, and analysis of complex chemical reaction networks. By using the Bell-Evans-Polanyi approximation for estimating activation energies, and an accurate and robust linear scaling for correcting the computed pKa values, we report herein multi-time-scale kinetic simulations for the self-assembly processes of polyoxotungstates that comprise 22 orders of magnitude, from tens of femtoseconds to months of reaction time. This very large time span was required to reproduce very fast processes such as the acid/base equilibria (at 10-12 s), relatively slow reactions such as the formation of key clusters such as the metatungstate (at 103 s), and the very slow assembly of the decatungstate (at 106 s). Analysis of the kinetic data and of the reaction network topology shed light onto the details of the main reaction mechanisms, which explains the origin of kinetic and thermodynamic control followed by the reaction. Simulations at alkaline pH fully reproduce experimental evidence since clusters do not form under those conditions.

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

A Structure-Based Gaussian Expansion for Quantum Reaction Dynamics in Molecules: Application to Hydrogen Tunneling in Malonaldehyde

We developed an approximate method for quantum reaction dynamics simulations, namely, a structure-based Gaussian (SBG) expansion approach, where SBG bases for the expansion of the wave function Ψ, expressed by a product of single-atom Cartesian Gaussians centered at the positions of respective nuclei, are mainly placed around critical structures on reaction pathways such as on the intrinsic reaction coordinate (IRC) through a transition state. In the present approach, the "pseudo-lattice points" at which SBGs are deployed are selected in a perturbative manner so as to make moderate the expansion length. We first applied the SBG idea to a two-dimensional quadruple-well model and obtained accurate tunneling splitting values between the lowest four states. We then applied it to hydrogen tunneling in malonaldehyde and achieved a tunneling splitting of 27.1 cm^-1 with only 875 SBGs at the MP2/6-31G(d,p) level of theory, in good agreement with 25 cm^-1 by the more elaborate multiconfiguration time-dependent Hartree method. Reasonable results were also obtained for singly and doubly deuterated malonaldehyde. We analyzed the tunneling states by utilizing expansion coefficients of individual SBGs and found that 40-45% of the SBGs in Ψ are nonplanar structures and SBGs away from the IRC contribute a little to hydrogen transfer.

Thermally Accessible Prebiotic Pathways for Forming Ribonucleic Acid and Protein Precursors from Aqueous Hydrogen Cyanide

The search for prebiotic chemical pathways to biologically relevant molecules is a long-standing puzzle that has generated a menagerie of competing hypotheses with limited experimental prospects for falsification. However, the advent of computational network exploration methodologies has created the opportunity to compare the kinetic plausibility of various channels and even propose new pathways. Here, the space of organic molecules that can be formed within four polar or pericyclic reactions from water and hydrogen cyanide (HCN), two established prebiotic candidates for generating biological precursors, was comprehensively explored with a state-of-the-art exploration algorithm. A surprisingly diverse reactivity landscape was revealed within just a few steps of these simple molecules. Reaction pathways to several biologically relevant molecules were discovered involving lower activation energies and fewer reaction steps compared with recently proposed alternatives. Accounting for water-catalyzed reactions qualitatively affects the interpretation of the network kinetics. The case-study also highlights omissions of simpler and lower barrier reaction pathways to certain products by other algorithms that qualitatively affect the interpretation of HCN reactivity.

Exploring the Quantum Chemical Energy Landscape with GNN-Guided Artificial Force

Artificial force has been proven useful to get over energy barriers and quickly search a large portion of the energy landscape. This work proposes a method based on graph neural networks to optimize the choice of transformation patterns to examine and accelerate energy landscape exploration. In open search from glutathione, the search efficiency was largely improved in comparison to random selection. We also applied transfer learning from glutathione to tuftsin, resulting in further efficiency gains.

Predictive chemistry: machine learning for reaction deployment, reaction development, and reaction discovery

The field of predictive chemistry relates to the development of models able to describe how molecules interact and react. It encompasses the long-standing task of computer-aided retrosynthesis, but is far more reaching and ambitious in its goals. In this review, we summarize several areas where predictive chemistry models hold the potential to accelerate the deployment, development, and discovery of organic reactions and advance synthetic chemistry.

Prediction of High-Yielding Single-Step or Cascade Pericyclic Reactions for the Synthesis of Complex Synthetic Targets

Pericyclic reactions, which involve cyclic concerted transition states without ionic or radical intermediates, have been extensively studied since their definition in the 1960s, and the famous Woodward-Hoffmann rules predict their stereoselectivity and chemoselectivity. Here, we describe the application of a fully automated reaction-path search method, that is, the artificial force induced reaction (AFIR), to trace an input compound back to reasonable starting materials through thermally allowed pericyclic reactions via product-based quantum-chemistry-aided retrosynthetic analysis (QCaRA) without using any a priori experimental knowledge. All categories of pericyclic reactions, including cycloadditions, ene reactions, group-transfer, cheletropic, electrocyclic, and sigmatropic reactions, were successfully traced back via concerted reaction pathways, and starting materials were computationally obtained with the correct stereochemistry. Furthermore, AFIR was used to predict whether the identified reaction pathway can be expected to occur in good yield relative to other possible reactions of the identified starting material. In order to showcase its practical utility, this state-of-the-art technology was also applied to the retrosynthetic analysis of a natural product with a relatively high number of atoms (52 atoms: endiandric acid C methyl ester), which was first synthesized by Nicolaou in 1982 and provided the corresponding starting polyenes with the correct stereospecificity via three pericyclic reaction cascades (one Diels-Alder reaction as well as 6π and 8π electrocyclic reactions). Moreover, not only systems that obey the Woodward-Hoffmann rules but also systems that violate these rules, such as those recently calculated by Houk, can be retrosynthesized accurately.

A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part II: A ChemTraYzer Study of Chlorinated Dibenzofuran Formation and Decomposition Processes

In our two-paper series, we first present the development of ReaxFF CHOCl parameters using the recently published ParAMS parametrization tool. In this second part, we update the reactive Molecular Dynamics - Quantum Mechanics coupling scheme ChemTraYzer and combine it with our new ReaxFF parameters from Part I to study formation and decomposition processes of chlorinated dibenzofurans. We introduce a self-learning method for recovering failed transition-state searches that improves the overall ChemTraYzer transition-state search success rate by 10 percentage points to a total of 48 %. With ChemTraYzer, we automatically find and quantify more than 500 reactions using transition state theory and DFT. Among the discovered chlorinated dibenzofuran reactions are numerous reactions that are new to the literature. In three case studies, we discuss the set of reactions that are most relevant to the dibenzofuran literature: (i) bimolecular reactions of the chlorinated-dibenzofuran precursors phenoxy radical and 1,3,5-trichlorobenzene, (ii) dibenzofuran chlorination and pyrolysis, and (iii) oxidation of chlorinated dibenzofurans.

Guided discovery of chemical reaction pathways with imposed activation

Computational power and quantum chemical methods have improved immensely since computers were first applied to the study of reactivity, but the de novo prediction of chemical reactions has remained challenging. We show that complex reaction pathways can be efficiently predicted in a guided manner using chemical activation imposed by geometrical constraints of specific reactive modes, which we term imposed activation (IACTA). Our approach is demonstrated on realistic and challenging chemistry, such as a triple cyclization cascade involved in the total synthesis of a natural product, a water-mediated Michael addition, and several oxidative addition reactions of complex drug-like molecules. Notably and in contrast with traditional hand-guided computational chemistry calculations, our method requires minimal human involvement and no prior knowledge of the products or the associated mechanisms. We believe that IACTA will be a transformational tool to screen for chemical reactivity and to study both by-product formation and decomposition pathways in a guided way.

Algorithmic Explorations of Unimolecular and Bimolecular Reaction Spaces

Algorithmic reaction exploration based on transition state searches has already made inroads into many niche applications, but its potential as a general-purpose tool is still largely unrealized. Computational cost and the absence of benchmark problems involving larger molecules remain obstacles to further progress. Here an ultra-low cost exploration algorithm is implemented and used to explore the reactivity of unimolecular and bimolecular reactants, comprising a total of 581 reactions involving 51 distinct reactants. The algorithm discovers all established reaction pathways, where such comparisons are possible, while also revealing a much richer reactivity landscape, including lower barrier reaction pathways and a strong dependence of reaction conformation in the apparent barriers of the reported reactions. The diversity of these benchmarks illustrate that reaction exploration algorithms are approaching general-purpose capability.

An automated method for graph-based chemical space exploration and transition state finding

Algorithms that automatically explore the chemical space have been limited to chemical systems with a low number of atoms due to expensive involved quantum calculations and the large amount of possible reaction pathways. The method described here presents a novel solution to the problem of chemical exploration by generating reaction networks with heuristics based on chemical theory. First, a second version of the reaction network is determined through molecular graph transformations acting upon functional groups of the reacting. Only transformations that break two chemical bonds and form two new ones are considered, leading to a significant performance enhancement compared to previously presented algorithm. Second, energy barriers for this reaction network are estimated through quantum chemical calculations by a growing string method, which can also identify non-octet species missed during the previous step and further define the reaction network. The proposed algorithm has been successfully applied to five different chemical reactions, in all cases identifying the most important reaction pathways.

Graph-Driven Reaction Discovery: Progress, Challenges, and Future Opportunities

Graph-based descriptors, such as bond-order matrices and adjacency matrices, offer a simple and compact way of categorizing molecular structures; furthermore, such descriptors can be readily used to catalog chemical reactions (i.e., bond-making and -breaking). As such, a number of graph-based methodologies have been developed with the goal of automating the process of generating chemical reaction network models describing the possible mechanistic chemistry in a given set of reactant species. Here, we outline the evolution of these graph-based reaction discovery schemes, with particular emphasis on more recent methods incorporating graph-based methods with semiempirical and ab initio electronic structure calculations, minimum-energy path refinements, and transition state searches. Using representative examples from homogeneous catalysis and interstellar chemistry, we highlight how these schemes increasingly act as "virtual reaction vessels" for interrogating mechanistic questions. Finally, we highlight where challenges remain, including issues of chemical accuracy and calculation speeds, as well as the inherent challenge of dealing with the vast size of accessible chemical reaction space.

Enhanced Sampling Aided Design of Molecular Photoswitches

Advances in the evolving field of atomistic simulations promise important insights for the design and fundamental understanding of novel molecular photoswitches. Here, we use state-of-the-art enhanced simulation techniques to unravel the complex, multistep chemistry of donor-acceptor Stenhouse adducts (DASAs). Our reaction discovery workflow consists of enhanced sampling for efficient chemical space exploration, refinement of newly observed pathways with more accurate ab initio electronic structure calculations, and structural modifications to introduce design principles within future generations of DASAs. We showcase our discovery workflow by not only recovering the full photoswitching mechanism of DASA but also predicting a plethora of new plausible thermal pathways and suggesting a way for their experimental validation. Furthermore, we illustrate the tunability of these newly discovered reactions, leading to a potential avenue for controlling DASA dynamics through multiple external stimuli. Overall, these insights could offer alternative routes to increase the efficiency and control of DASA's photoswitching mechanism, providing new elements to design more complex light-responsive materials.

Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods

The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S₀/T₁ minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S₀/S₁ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.

Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks

Fueled by advances in hardware and algorithm design, large-scale automated explorations of chemical reaction space have become possible. Here, we present our approach to an open-source, extensible framework for explorations of chemical reaction mechanisms based on the first-principles of quantum mechanics. It is intended to facilitate reaction network explorations for diverse chemical problems with a wide range of goals such as mechanism elucidation, reaction path optimization, retrosynthetic path validation, reagent design, and microkinetic modeling. The stringent first-principles basis of all algorithms in our framework is key for the general applicability that avoids any restrictions to specific chemical systems. Such an agile framework requires multiple specialized software components of which we present three modules in this work. The key module, Chemoton, drives the exploration of reaction networks. For the exploration itself, we introduce two new algorithms for elementary-step searches that are based on Newton trajectories. The performance of these algorithms is assessed for a variety of reactions characterized by a broad chemical diversity in terms of bonding patterns and chemical elements. Chemoton successfully recovers the vast majority of these. We provide the resulting data, including large numbers of reactions that were not included in our reference set, to be used as a starting point for further explorations and for future reference.

Deep reaction network exploration at a heterogeneous catalytic interface

Characterizing the reaction energies and barriers of reaction networks is central to catalyst development. However, heterogeneous catalytic surfaces pose several unique challenges to automatic reaction network characterization, including large sizes and open-ended reactant sets, that make ad hoc network construction the current state-of-the-art. Here, we show how automated network exploration algorithms can be adapted to the constraints of heterogeneous systems using ethylene oligomerization on silica-supported single-site Ga3+ as a model system. Using only graph-based rules for exploring the network and elementary constraints based on activation energy and size for identifying network terminations, a comprehensive reaction network is generated and validated against standard methods. The algorithm (re)discovers the Ga-alkyl-centered Cossee-Arlman mechanism that is hypothesized to drive major product formation while also predicting several new pathways for producing alkanes and coke precursors. These results demonstrate that automated reaction exploration algorithms are rapidly maturing towards general purpose capability for exploratory catalytic applications.

Automated input structure generation for single-ended reaction path optimizations

The exploration of a reaction network requires highly automated workflows to avoid error-prone and time-consuming manual steps. In this respect, a major bottleneck is the search for transition-state (TS) structures, which frequently fails and, therefore, makes (manual) revision necessary. In this work, we present a technique for obtaining suitable input structures for automated TS searches based on single-ended reaction path optimization algorithms, which makes subsequent TS searches via this method significantly more robust. First, possible input structures are generated based on the spatial alignment of the reactants. The appropriate orientation of reacting groups is achieved via stepwise rotations along selected torsional degrees of freedom. Second, a ranking of the obtained structures is performed according to selected geometric criteria. The main goals are to properly align the reactive atoms, to avoid hindrance within the reaction channel and to resolve steric clashes between the reactants. The developed procedure has been carefully tested on a variety of examples and provides suitable input structures for TS searches within seconds. The method is in daily use in an industrial setting.

Artificial intelligence pathway search to resolve catalytic glycerol hydrogenolysis selectivity

The complex interaction between molecules and catalyst surfaces leads to great difficulties in understanding and predicting the activity and selectivity in heterogeneous catalysis. Here we develop an end-to-end artificial intelligence framework for the activity prediction of heterogeneous catalytic systems (AI-Cat method), which takes simple inputs from names of molecules and metal catalysts and outputs the reaction energy profile from the input molecule to low energy pathway products. The AI-Cat method combines two neural network models, one for predicting reaction patterns and the other for providing the reaction barrier and energy, with a Monte Carlo tree search to resolve the low energy pathways in a reaction network. We then apply AI-Cat to resolve the reaction network of glycerol hydrogenolysis on Cu surfaces, which is a typical selective C-O bond activation system and of key significance for biomass-derived polyol utilization. We show that glycerol hydrogenolysis features a huge reaction network of relevant candidates, containing 420 reaction intermediates and 2467 elementary reactions. Among them, the surface-mediated enol-keto tautomeric resonance is a key step to facilitate the primary C-OH bond breaking and thus selects 1,2-propanediol as the major product on Cu catalysts. 1,3-Propanediol can only be produced under strong acidic conditions and high surface H coverage by following a hydrogenation-dehydration pathway. AI-Cat further discovers six low-energy reaction patterns for C-O bond activation on metals that is of general significance to polyol catalysis. Our results demonstrate that the reaction prediction for complex heterogeneous catalysis is now feasible with AI-based atomic simulation and a Monte Carlo tree search.

Successes and challenges in using machine-learned activation energies in kinetic simulations

The prediction of the thermodynamic and kinetic properties of chemical reactions is increasingly being addressed by machine-learning (ML) methods such as artificial neural networks (ANNs). While a number of recent studies have reported success in predicting chemical reaction activation energies, less attention has focused on how the accuracy of ML predictions filter through to predictions of macroscopic observables. Here, we consider the impact of the uncertainty associated with ML prediction of activation energies on observable properties of chemical reaction networks, as given by microkinetics simulations based on ML-predicted reaction rates. After training an ANN to predict activation energies given standard molecular descriptors for reactants and products alone, we performed microkinetics simulations of three different prototypical reaction networks: formamide decomposition, aldol reactions and decomposition of 3-hydroperoxypropanal. We find that the kinetic modelling predictions can be in excellent agreement with corresponding simulations performed with ab initio calculations, but this is dependent on the inherent energetic landscape of the networks. We use these simulations to suggest some guidelines for when ML-based activation energies can be reliable, and when one should take more care in applications to kinetics modelling.

Fast and automated identification of reactions with low barriers using meta-MD simulations

We test our meta-molecular dynamics (MD) based approach for finding low-barrier (<30 kcal/mol) reactions on uni- and bimolecular reactions extracted from the barrier dataset developed by Grambow, Pattanaik & Green (2020). For unimolecular reactions the meta-MD simulations identify 25 of the 26 products found by Grambow, Pattanaik & Green (2020), while the subsequent semiempirical screening eliminates an additional four reactions due to an overestimation of the reaction energies or estimated barrier heights relative to DFT. In addition, our approach identifies 36 reactions not found by Grambow, Pattanaik & Green (2020), 10 of which are <30 kcal/mol. For bimolecular reactions the meta-MD simulations identify 19 of the 20 reactions found by Grambow, Pattanaik & Green (2020), while the subsequent semiempirical screening eliminates an additional reaction. In addition, we find 34 new low-barrier reactions. For bimolecular reactions we found that it is necessary to “encourage” the reactants to go to previously undiscovered products, by including products found by other MD simulations when computing the biasing potential as well as decreasing the size of the molecular cavity in which the MD occurs, until a reaction is observed. We also show that our methodology can find the correct products for two reactions that are more representative of those encountered in synthetic organic chemistry. The meta-MD hyperparameters used in this study thus appear to be generally applicable to finding low-barrier reactions.

Efficient Reaction Space Exploration with ChemTraYzer-TAD

The development of a reaction model is often a time-consuming process, especially if unknown reactions have to be found and quantified. To alleviate the reaction modeling process, automated procedures for reaction space exploration are highly desired. We present ChemTraYzer-TAD, a new reactive molecular dynamics acceleration technique aimed at efficient reaction space exploration. The new method is based on the basin confinement strategy known from the temperature-accelerated dynamics (TAD) acceleration method. Our method features integrated ChemTraYzer bond-order processing steps for the automatic and on-the-fly determination of the positions of virtual walls in configuration space that confine the system in a potential energy basin. We use the example of 1,3-dioxolane-4-hydroperoxide-2-yl radical oxidation to show that ChemTraYzer-TAD finds more than 100 different parallel reactions for the given set of reactants in less than 2 ns of simulation time. Among the many observed reactions, ChemTraYzer-TAD finds the expected typical low-temperature reactions despite the use of extremely high simulation temperatures up to 5000 K. Our method also finds a new concerted β-scission plus O₂ addition with a lower reaction barrier than the literature-known and so-far dominant β-scission.

CARNOT: a Fragment-Based Direct Molecular Dynamics and Virtual-Reality Simulation Package for Reactive Systems

Traditionally, the study of reaction mechanisms of complex reaction systems such as combustion has been performed on an individual basis by optimizations of transition structure and minimum energy path or by reaction dynamics trajectory calculations for one elementary reaction at a time. It is effective, but time-consuming, whereas important and unexpected processes could have been missed. In this article, we present a direct molecular dynamics (DMD) approach and a virtual-reality simulation program, CARNOT, in which plausible chemical reactions are simulated simultaneously at finite temperature and pressure conditions. A key concept of the present ab initio molecular dynamics method is to partition a large, chemically reactive system into molecular fragments that can be adjusted on the fly of a DMD simulation. The theory represents an extension of the explicit polarization method to reactive events, called ReX-Pol. We propose a highest-and-lowest adapted-spin approximation to define the local spins of individual fragments, rather than treating the entire system by a delocalized wave function. Consequently, the present ab initio DMD can be applied to reactive systems consisting of an arbitrarily varying number of closed and open-shell fragments such as free radicals, zwitterions, and separate ions found in combustion and other reactions. A graph-data structure algorithm was incorporated in CARNOT for the analysis of reaction networks, suitable for reaction mechanism reduction. Employing the PW91 density functional theory and the 6-31+G(d) basis set, the capabilities of the CARNOT program were illustrated by a combustion reaction, consisting of 28 650 atoms, and by reaction network analysis that revealed a range of mechanistic and dynamical events. The method may be useful for applications to other types of complex reactions.

New Tools for Taming Complex Reaction Networks: The Unimolecular Decomposition of Indole Revisited

The level of detail attained in the computational description of reaction mechanisms can be vastly improved through tools for automated chemical space exploration, particularly for systems of small to medium size. Under this approach, the unimolecular decomposition landscape for indole was explored through the automated reaction mechanism discovery program AutoMeKin. Nevertheless, the sheer complexity of the obtained mechanisms might be a hindrance regarding their chemical interpretation. In this spirit, the new Python library amk-tools has been designed to read and manipulate complex reaction networks, greatly simplifying their overall analysis. The package provides interactive dashboards featuring visualizations of the network, the three-dimensional (3D) molecular structures and vibrational normal modes of all chemical species, and the corresponding energy profiles for selected pathways. The combination of the joined mechanism generation and postprocessing workflow with the rich chemistry of indole decomposition enabled us to find new details of the reaction (obtained at the CCSD(T)/aug-cc-pVTZ//M06-2X/MG3S level of theory) that were not reported before: (i) 16 pathways leading to the formation of HCN and NH₃ (via amino radical); (ii) a barrierless reaction between methylene radical and phenyl isocyanide, which might be an operative mechanism under the conditions of the interstellar medium; and (iii) reaction channels leading to both hydrogen cyanide and hydrogen isocyanide, of potential astrochemical interest as the computed HNC/HCN ratios greatly exceed the calculated equilibrium value at very low temperatures. The reported reaction networks can be very valuable to supplement databases of kinetic data, which is of remarkable interest for pyrolysis and astrochemical studies.

Expansive Quantum Mechanical Exploration of Chemical Reaction Paths

Quantum mechanical methods have been well-established for the elucidation of reaction paths of chemical processes and for the explicit dynamics of molecular systems. While they are usually deployed in routine manual calculations on reactions for which some insights are already available (typically from experiment), new algorithms and continuously increasing capabilities of modern computer hardware allow for exploratory open-ended computational campaigns that are unbiased and therefore enable unexpected discoveries. Highly efficient and even automated procedures facilitate systematic approaches toward the exploration of uncharted territory in molecular transformations and dynamics. In this work, we elaborate on such explorative approaches that range from reaction network explorations with (stationary) quantum chemical methods to explorative molecular dynamics and migrant wave packet dynamics. The focus is on recent developments that cover the following strategies. (i) Pruning search options for elementary reaction steps by heuristic rules based on the first-principles of quantum mechanics: Rules are required for reducing the combinatorial explosion of potentially reactive atom pairings, and rooting them in concepts derived from the electronic wave function makes them applicable to any molecular system. (ii) Enforcing reactive events by external biases: Inducing a reaction requires constraints that steer and direct elementary-step searches, which can be formulated in terms of forces, velocities, or supplementary potentials. (iii) Manual steering facilitated by interactive quantum mechanics: As ultrafast quantum chemical methods allow for real-time manual interactions with molecular systems, human-intuition-guided paths can be easily explored with suitable human-machine interfaces. (iv) New approaches for transition-state optimization with continuous curve representations can provide stable schemes to be driven in an automated way by allowing for an efficient tuning of the curve's parameters (instead of a manipulation of a collection of structures along the path), and (v) reactive molecular dynamics and direct wave packet propagation exploit the equations of motion of an underlying mechanical theory (usually, classical Newtonian mechanics or Schrödinger quantum mechanics). Explorative approaches are likely to replace the current state of the art in computational chemistry, because they reduce the human effort to be invested in reaction path elucidations, they are less prone to errors and bias-free, and they cover more extensive regions of the relevant configuration space. As a result, computational investigations that rely on these techniques are more likely to deliver surprising discoveries.

Discovering surface reaction pathways using accelerated molecular dynamics and network analysis tools

We present an automated method that maps surface reaction pathways with no experimental data and with minimal human interventions. In this method, bias potentials promoting surface reactions are applied to enable statistical samplings of the surface reaction within the timescale of ab initio molecular dynamics (MD) simulations, and elementary reactions are extracted automatically using an extension of the method constructed for gas- or liquid-phase reactions. By converting the extracted elementary reaction data to directed graph data, MD trajectories can be efficiently mapped onto reaction pathways using a network analysis tool. To demonstrate the power of the method, it was applied to the steam reforming of methane on the Rh(111) surface and to propane reforming on the Pt(111) and Pt₃Sn(111) surfaces. We discover new energetically favorable pathways for both reactions and reproduce the experimentally-observed materials-dependence of the surface reaction activity and the selectivity for the propane reforming reactions.

We present an automated method that maps surface reaction pathways with no experimental data and with minimal human interventions.

An automatized workflow from molecular dynamic simulation to quantum chemical methods to identify elementary reactions and compute reaction constants

We present an automatized workflow which, starting from molecular dynamics simulations, identifies reaction events, filters them, and prepares them for accurate quantum chemical calculations using, for example, Density Functional Theory (DFT) or Coupled Cluster methods. The capabilities of the automatized workflow are demonstrated by the example of simulations for the combustion of some polycyclic aromatic hydrocarbons (PAHs). It is shown how key elementary reaction candidates are filtered out of a much larger set of redundant reactions and refined further. The molecular species in question are optimized using DFT and reaction energies, barrier heights, and reaction rates are calculated. The setup is general enough to include at this stage configurational sampling, which can be exploited in the future. Using the introduced machinery, we investigate how the observed reaction types depend on the gas atmosphere used in the molecular dynamics simulation. For the re-optimization on the DFT level, we show how the additional information needed to switch from reactive force-field to electronic structure calculations can be filled in and study how well ReaxFF and DFT agree with each other and shine light on the perspective of using more accurate semi-empirical methods in the MD simulation.

Assessing and rationalizing the performance of Hessian update schemes for reaction path Hamiltonian rate calculations

The reaction path Hamiltonian (RPH) can be used to calculate chemical reaction rate constants, going beyond transition-state theory in taking account of recrossing by providing an approximation to the dynamic transmission coefficient. However, the RPH necessitates the calculation of the Hessian matrix at a number of points along the minimum energy path; the associated computational cost stands as a bottleneck in RPH calculations, especially if one is interested in using high-accuracy electronic structure methods. In this work, four different Hessian update schemes (symmetric rank-1, Powell-symmetric Broyden, Bofill, and TS-BFGS updates) are assessed to see whether or not they reliably reproduce calculated transmission coefficients for three different chemical reactions. Based on the reactions investigated, the symmetric rank-1 Hessian update was the least appropriate for RPH construction, giving different transmission coefficients from the standard analytical Hessian approach, as well as inconsistent frequencies and coupling properties. The Bofill scheme, the Powell-symmetric Broyden scheme, and the TS-BFGS scheme were the most reliable Hessian update methods, with transmission coefficients that were in good agreement with those calculated by the standard RPH calculations. The relative accuracy of the different Hessian update schemes is further rationalized by investigating the approximated Coriolis and curvature coupling terms along the reaction-path, providing insight into when these schemes would be expected to work well. Furthermore, the associated computational cost associated with the RPH calculations was substantially reduced by the tested update schemes. Together, these results provide useful rules-of-thumb for using Hessian update schemes in RPH simulations.

Discovering Catalytic Reaction Networks Using Deep Reinforcement Learning from First-Principles

Determining the reaction pathways, which is central to illustrating the working mechanisms of a catalyst, is severely hindered by the high complexity of the reaction and the extreme scarcity of the data. Here, we develop a novel artificial intelligence framework integrating deep reinforcement learning (DRL) techniques with density functional theory simulations to automate the quantitative search and evaluation on the complex catalytic reaction networks from zero knowledge. Our framework quantitatively transforms the first-principles-derived free energy landscape of the chemical reactions to a DRL environment and the corresponding actions. By interacting with this dynamic environment, our model evolves by itself from scratch to a complete reaction path. We demonstrate this framework using the Haber-Bosch process on the most active Fe(111) surface. The new path found by our framework has a lower overall free energy barrier than the previous study based on domain knowledge, demonstrating its outstanding capability in discovering complicated reaction paths. Looking forward, we anticipate that this framework will open the door to exploring the fundamental reaction mechanisms of many catalytic reactions.

AutoMeKin2021: An open-source program for automated reaction discovery

AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin.