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.

Reaction Networks Resemble Low-Dimensional Regular Lattices

The computational exploration, manipulation, and design of complex chemical reactions face fundamental challenges related to the high-dimensional nature of potential energy surfaces (PESs) that govern reactivity. Accurately modeling complex reactions is crucial for understanding the chemical processes involved in, for example, organocatalysis, autocatalytic cycles, and one-pot molecular assembly. Our prior research demonstrated that discretizing PESs using heuristics based on bond breaking and bond formation produces a reaction network representation with a low-dimensional structure (metric space). We now find that these stoichiometry-preserving reaction networks possess additional, though approximate, structure and resemble low-dimensional regular lattices with a small amount of random edge rewiring. The heuristics-based discretization thus generates a nonlinear dimensionality reduction by a factor of 10 with an a posteriori error measure (probability of random rewiring). The structure becomes evident through a comparative analysis of CHNO reaction networks of varying stoichiometries against a panel of size-matched generative network models, taking into account their local, metric, and global properties. The generative models include random networks (Erdős-Rényi and bipartite random networks), regular lattices (periodic and nonperiodic), and network models with a tunable level of "randomness" (Watts-Strogatz graphs and regular lattices with random rewiring). The CHNO networks are simultaneously closely matched in all these properties by 3-4-dimensional regular lattices with 10% or less of edges randomly rewired. The effective dimensionality reduction is found to be independent of the system size, stoichiometry, and ruleset, suggesting that search and sampling algorithms for PESs of complex chemical reactions can be effectively leveraged.

Nanoscale chemical reaction exploration with a quantum magnifying glass

Nanoscopic systems exhibit diverse molecular substructures by which they facilitate specific functions. Theoretical models of them, which aim at describing, understanding, and predicting these capabilities, are difficult to build. Viable quantum-classical hybrid models come with specific challenges regarding atomistic structure construction and quantum region selection. Moreover, if their dynamics are mapped onto a state-to-state mechanism such as a chemical reaction network, its exhaustive exploration will be impossible due to the combinatorial explosion of the reaction space. Here, we introduce a "quantum magnifying glass" that allows one to interactively manipulate nanoscale structures at the quantum level. The quantum magnifying glass seamlessly combines autonomous model parametrization, ultra-fast quantum mechanical calculations, and automated reaction exploration. It represents an approach to investigate complex reaction sequences in a physically consistent manner with unprecedented effortlessness in real time. We demonstrate these features for reactions in bio-macromolecules and metal-organic frameworks, diverse systems that highlight general applicability.

SCINE-Software for chemical interaction networks

The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.

Identification of Reaction Network Hypotheses for Complex Feedstocks from Spectroscopic Measurements with Minimal Human Intervention

In this work, we detail an automated reaction network hypothesis generation protocol for processes involving complex feedstocks where information about the species and reactions involved is unknown. Our methodology is process agnostic and can be utilized in any reactive process with spectroscopic measurements that provide information on the evolution of the components in the mixture. We decompose the mixture spectra to obtain spectroscopic signatures of the individual components and use a 1-D convolutional neural network to automatically identify functional groups indicated by them. We employ atom-atom mapping to automatically recover reaction rules that are applied on candidate molecules identified from chemistry databases through fingerprint similarity. The method is tested on synthetic data and on spectroscopic measurements of lab-scale batch hydrothermal liquefaction (HTL) of biomass to determine the accuracy of prediction across datasets of varying complexities. Our methodology is able to identify reaction network hypotheses containing reaction networks close to the ground truth in the case of synthetic data, and we are also able to recover candidate molecules and reaction networks close to the ones reported in the previous literature studies for biomass pyrolysis.

Uncertainty-Aware First-Principles Exploration of Chemical Reaction Networks

Exploring large chemical reaction networks with automated exploration approaches and accurate quantum chemical methods can require prohibitively large computational resources. Here, we present an automated exploration approach that focuses on the kinetically relevant part of the reaction network by interweaving (i) large-scale exploration of chemical reactions, (ii) identification of kinetically relevant parts of the reaction network through microkinetic modeling, (iii) quantification and propagation of uncertainties, and (iv) reaction network refinement. Such an uncertainty-aware exploration of kinetically relevant parts of a reaction network with automated accuracy improvement has not been demonstrated before in a fully quantum mechanical approach. Uncertainties are identified by local or global sensitivity analysis. The network is refined in a rolling fashion during the exploration. Moreover, the uncertainties are considered during kinetically steering of a rolling reaction network exploration. We demonstrate our approach for Eschenmoser-Claisen rearrangement reactions. The sensitivity analysis identifies that only a small number of reactions and compounds are essential for describing the kinetics reliably, resulting in efficient explorations without sacrificing accuracy and without requiring prior knowledge about the chemistry unfolding.

Benchmarking machine-readable vectors of chemical reactions on computed activation barriers

In recent years, there has been a surge of interest in predicting computed activation barriers, to enable the acceleration of the automated exploration of reaction networks. Consequently, various predictive approaches have emerged, ranging from graph-based models to methods based on the three-dimensional structure of reactants and products. In tandem, many representations have been developed to predict experimental targets, which may hold promise for barrier prediction as well. Here, we bring together all of these efforts and benchmark various methods (Morgan fingerprints, the DRFP, the CGR representation-based Chemprop, SLATMd, B2Rl2, EquiReact and language model BERT + RXNFP) for the prediction of computed activation barriers on three diverse datasets.

A human-machine interface for automatic exploration of chemical reaction networks

Autonomous reaction network exploration algorithms offer a systematic approach to explore mechanisms of complex chemical processes. However, the resulting reaction networks are so vast that an exploration of all potentially accessible intermediates is computationally too demanding. This renders brute-force explorations unfeasible, while explorations with completely pre-defined intermediates or hard-wired chemical constraints, such as element-specific coordination numbers, are not flexible enough for complex chemical systems. Here, we introduce a STEERING WHEEL to guide an otherwise unbiased automated exploration. The STEERING WHEEL algorithm is intuitive, generally applicable, and enables one to focus on specific regions of an emerging network. It also allows for guiding automated data generation in the context of mechanism exploration, catalyst design, and other chemical optimization challenges. The algorithm is demonstrated for reaction mechanism elucidation of transition metal catalysts. We highlight how to explore catalytic cycles in a systematic and reproducible way. The exploration objectives are fully adjustable, allowing one to harness the STEERING WHEEL for both structure-specific (accurate) calculations as well as for broad high-throughput screening of possible reaction intermediates.

Interoperable workflows by exchanging grid-based data between quantum-chemical program packages

Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability.

Physical Chemistry Models for Chemical Research in the XXth and XXIst Centuries

Thermodynamic hypotheses and models are the touchstone for chemical results, but the actual models based on time-invariance, which have performed efficiently in the development of chemistry, are nowadays invalid for the interpretation of the behavior of complex systems exhibiting nonlinear kinetics and with matter and energy exchange flows with the surroundings. Such fields of research will necessarily foment and drive the use of thermodynamic models based on the description of irreversibility at the macroscopic level, instead of the current models which are strongly anchored in microreversibility.

CREST-A program for the exploration of low-energy molecular chemical space

Conformer-rotamer sampling tool (CREST) is an open-source program for the efficient and automated exploration of molecular chemical space. Originally developed in Pracht et al. [Phys. Chem. Chem. Phys. 22, 7169 (2020)] as an automated driver for calculations at the extended tight-binding level (xTB), it offers a variety of molecular- and metadynamics simulations, geometry optimization, and molecular structure analysis capabilities. Implemented algorithms include automated procedures for conformational sampling, explicit solvation studies, the calculation of absolute molecular entropy, and the identification of molecular protonation and deprotonation sites. Calculations are set up to run concurrently, providing efficient single-node parallelization. CREST is designed to require minimal user input and comes with an implementation of the GFNn-xTB Hamiltonians and the GFN-FF force-field. Furthermore, interfaces to any quantum chemistry and force-field software can easily be created. In this article, we present recent developments in the CREST code and show a selection of applications for the most important features of the program. An important novelty is the refactored calculation backend, which provides significant speed-up for sampling of small or medium-sized drug molecules and allows for more sophisticated setups, for example, quantum mechanics/molecular mechanics and minimum energy crossing point calculations.

Network topology mapping of chemical compounds space

We define bipartite and monopartite relational networks of chemical elements and compounds using two different datasets of inorganic chemical and material compounds, as well as study their topology. We discover that the connectivity between elements and compounds is distributed exponentially for materials, and with a fat tail for chemicals. Compounds networks show similar distribution of degrees, and feature a highly-connected club due to oxygen . Chemical compounds networks appear more modular than material ones, while the communities detected reveal different dominant elements specific to the topology. We successfully reproduce the connectivity of the empirical chemicals and materials networks by using a family of fitness models, where the fitness values are derived from the abundances of the elements in the aggregate compound data. Our results pave the way towards a relational network-based understanding of the inherent complexity of the vast chemical knowledge atlas, and our methodology can be applied to other systems with the ingredient-composite structure.

Exploring Chemical Space Using Ab Initio Hyperreactor Dynamics

In recent years, first-principles exploration of chemical reaction space has provided valuable insights into intricate reaction networks. Here, we introduce ab initio hyperreactor dynamics, which enables rapid screening of the accessible chemical space from a given set of initial molecular species, predicting new synthetic routes that can potentially guide subsequent experimental studies. For this purpose, different hyperdynamics derived bias potentials are applied along with pressure-inducing spherical confinement of the molecular system in ab initio molecular dynamics simulations to efficiently enhance reactivity under mild conditions. To showcase the advantages and flexibility of the hyperreactor approach, we present a systematic study of the method's parameters on a HCN toy model and apply it to a recently introduced experimental model for the prebiotic formation of glycinal and acetamide in interstellar ices, which yields results in line with experimental findings. In addition, we show how the developed framework enables the study of complicated transitions like the first step of a nonenzymatic DNA nucleoside synthesis in an aqueous environment, where the molecular fragmentation problem of earlier nanoreactor approaches is avoided.

Scientific Deep Machine Learning Concepts for the Prediction of Concentration Profiles and Chemical Reaction Kinetics: Consideration of Reaction Conditions

Emerging concepts from scientific deep machine learning such as physics-informed neural networks (PINNs) enable a data-driven approach for the study of complex kinetic problems. We present an extended framework that combines the advantages of PINNs with the detailed consideration of experimental parameter variations for the simulation and prediction of chemical reaction kinetics. The approach is based on truncated Taylor series expansions for the underlying fundamental equations, whereby the external variations can be interpreted as perturbations of the kinetic parameters. Accordingly, our method allows for an efficient consideration of experimental parameter settings and their influence on the concentration profiles and reaction kinetics. A particular advantage of our approach, in addition to the consideration of univariate and multivariate parameter variations, is the robust model-based exploration of the parameter space to determine optimal reaction conditions in combination with advanced reaction insights. The benefits of this concept are demonstrated for higher-order chemical reactions including catalytic and oscillatory systems in combination with small amounts of training data. All predicted values show a high level of accuracy, demonstrating the broad applicability and flexibility of our approach.

A query-response causal analysis of reaction events in biochemical reaction networks

The stochastic kinetics of biochemical reaction networks is described by a chemical master equation (CME) and the underlying laws of mass action. Assuming network-free simulations of the rule-based models of biochemical reaction networks (BRNs), this paper departs from the usual analysis of network dynamics as the time-dependent distributions of chemical species counts, and instead considers statistically evaluating the sequences of reaction events generated from the stochastic simulations. The reaction event-time series can be used for reaction clustering, identifying rare events, and recognizing the periods of increased or steady-state activity. However, the main aim of this paper is to device an effective method for identifying causally and anti-causally related sub-sequences of reaction events using their empirical probabilities. This allows discovering some of the causal dynamics of BRNs as well as uncovering their short-term deterministic behaviors. In particular, it is proposed that the reaction sub-sequences that are conditionally nearly certain or nearly uncertain can be considered as being causally related. Moreover, since the time-ordering of reaction events is locally irrelevant, the reaction sub-sequences can be transformed into the reaction sets or multi-sets. The distance metrics can be then used to define the equivalences among the reaction events. The proposed method for identifying the causally related reaction sub-sequences has been implemented as a computationally efficient query-response mechanism. The method was evaluated for five models of genetic networks in seven defined numerical experiments. The models were simulated in BioNetGen using the open-source network-free simulator NFsim. This simulator had to be modified first to allow recording the traces of reaction events, and it is available in the Github repository, ploskot/nfsim_1.20. The generated event time-series were analyzed with Python and Matlab scripts. The whole process of data generation, analysis and visualization has been nearly fully automated using shell scripts. This demonstrates the opportunities for substantially increasing the research productivity by creating automated data generation and processing pipelines.

Diffusion-based generative AI for exploring transition states from 2D molecular graphs

The exploration of transition state (TS) geometries is crucial for elucidating chemical reaction mechanisms and modeling their kinetics. Recently, machine learning (ML) models have shown remarkable performance for prediction of TS geometries. However, they require 3D conformations of reactants and products often with their appropriate orientations as input, which demands substantial efforts and computational cost. Here, we propose a generative approach based on the stochastic diffusion method, namely TSDiff, for prediction of TS geometries just from 2D molecular graphs. TSDiff outperforms the existing ML models with 3D geometries in terms of both accuracy and efficiency. Moreover, it enables to sample various TS conformations, because it learns the distribution of TS geometries for diverse reactions in training. Thus, TSDiff finds more favorable reaction pathways with lower barrier heights than those in the reference database. These results demonstrate that TSDiff shows promising potential for an efficient and reliable TS exploration.

Accurate transition state generation with an object-aware equivariant elementary reaction diffusion model

Transition state search is key in chemistry for elucidating reaction mechanisms and exploring reaction networks. The search for accurate 3D transition state structures, however, requires numerous computationally intensive quantum chemistry calculations due to the complexity of potential energy surfaces. Here we developed an object-aware SE(3) equivariant diffusion model that satisfies all physical symmetries and constraints for generating sets of structures-reactant, transition state and product-in an elementary reaction. Provided reactant and product, this model generates a transition state structure in seconds instead of hours, which is typically required when performing quantum-chemistry-based optimizations. The generated transition state structures achieve a median of 0.08 Å root mean square deviation compared to the true transition state. With a confidence scoring model for uncertainty quantification, we approach an accuracy required for reaction barrier estimation (2.6 kcal mol-1) by only performing quantum chemistry-based optimizations on 14% of the most challenging reactions. We envision usefulness for our approach in constructing large reaction networks with unknown mechanisms.

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.

Δ2 machine learning for reaction property prediction

The emergence of Δ-learning models, whereby machine learning (ML) is used to predict a correction to a low-level energy calculation, provides a versatile route to accelerate high-level energy evaluations at a given geometry. However, Δ-learning models are inapplicable to reaction properties like heats of reaction and activation energies that require both a high-level geometry and energy evaluation. Here, a Δ2-learning model is introduced that can predict high-level activation energies based on low-level critical-point geometries. The Δ2 model uses an atom-wise featurization typical of contemporary ML interatomic potentials (MLIPs) and is trained on a dataset of ∼167 000 reactions, using the GFN2-xTB energy and critical-point geometry as a low-level input and the B3LYP-D3/TZVP energy calculated at the B3LYP-D3/TZVP critical point as a high-level target. The excellent performance of the Δ2 model on unseen reactions demonstrates the surprising ease with which the model implicitly learns the geometric deviations between the low-level and high-level geometries that condition the activation energy prediction. The transferability of the Δ2 model is validated on several external testing sets where it shows near chemical accuracy, illustrating the benefits of combining ML models with readily available physical-based information from semi-empirical quantum chemistry calculations. Fine-tuning of the Δ2 model on a small number of Gaussian-4 calculations produced a 35% accuracy improvement over DFT activation energy predictions while retaining xTB-level cost. The Δ2 model approach proves to be an efficient strategy for accelerating chemical reaction characterization with minimal sacrifice in prediction accuracy.

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.

Quantum chemical data generation as fill-in for reliability enhancement of machine-learning reaction and retrosynthesis planning

Data-driven synthesis planning has seen remarkable successes in recent years by virtue of modern approaches of artificial intelligence that efficiently exploit vast databases with experimental data on chemical reactions. However, this success story is intimately connected to the availability of existing experimental data. It may well occur in retrosynthetic and synthesis design tasks that predictions in individual steps of a reaction cascade are affected by large uncertainties. In such cases, it will, in general, not be easily possible to provide missing data from autonomously conducted experiments on demand. However, first-principles calculations can, in principle, provide missing data to enhance the confidence of an individual prediction or for model retraining. Here, we demonstrate the feasibility of such an ansatz and examine resource requirements for conducting autonomous first-principles calculations on demand.

Accelerating Reaction Network Explorations with Automated Reaction Template Extraction and Application

Autonomously exploring chemical reaction networks with first-principles methods can generate vast data. Especially autonomous explorations without tight constraints risk getting trapped in regions of reaction networks that are not of interest. In many cases, these regions of the networks are only exited once fully searched. Consequently, the required human time for analysis and computer time for data generation can make these investigations unfeasible. Here, we show how simple reaction templates can facilitate the transfer of chemical knowledge from expert input or existing data into new explorations. This process significantly accelerates reaction network explorations and improves cost-effectiveness. We discuss the definition of the reaction templates and their generation based on molecular graphs. The resulting simple filtering mechanism for autonomous reaction network investigations is exemplified with a polymerization reaction.

Combustion, Chemistry, and Carbon Neutrality

Combustion is a reactive oxidation process that releases energy bound in chemical compounds used as fuels─energy that is needed for power generation, transportation, heating, and industrial purposes. Because of greenhouse gas and local pollutant emissions associated with fossil fuels, combustion science and applications are challenged to abandon conventional pathways and to adapt toward the demand of future carbon neutrality. For the design of efficient, low-emission processes, understanding the details of the relevant chemical transformations is essential. Comprehensive knowledge gained from decades of fossil-fuel combustion research includes general principles for establishing and validating reaction mechanisms and process models, relying on both theory and experiments with a suite of analytic monitoring and sensing techniques. Such knowledge can be advantageously applied and extended to configure, analyze, and control new systems using different, nonfossil, potentially zero-carbon fuels. Understanding the impact of combustion and its links with chemistry needs some background. The introduction therefore combines information on exemplary cultural and technological achievements using combustion and on nature and effects of combustion emissions. Subsequently, the methodology of combustion chemistry research is described. A major part is devoted to fuels, followed by a discussion of selected combustion applications, illustrating the chemical information needed for the future.

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

Machine Learning in Unmanned Systems for Chemical Synthesis

Chemical synthesis is state-of-the-art, and, therefore, it is generally based on chemical intuition or experience of researchers. The upgraded paradigm that incorporates automation technology and machine learning (ML) algorithms has recently been merged into almost every subdiscipline of chemical science, from material discovery to catalyst/reaction design to synthetic route planning, which often takes the form of unmanned systems. The ML algorithms and their application scenarios in unmanned systems for chemical synthesis were presented. The prospects for strengthening the connection between reaction pathway exploration and the existing automatic reaction platform and solutions for improving autonomation through information extraction, robots, computer vision, and intelligent scheduling were proposed.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

The computational road to reactivity scales

Reactivity scales are useful research tools for chemists, both experimental and computational. However, to determine the reactivity of a single molecule, multiple measurements need to be carried out, which is a time-consuming and resource-intensive task. In this Tutorial Review, we present alternative approaches for the efficient generation of quantitative structure-reactivity relationships that are based on quantum chemistry, supervised learning, and uncertainty quantification. First published in 2002, we observe a tendency for these relationships to become not only more predictive but also more interpretable over time.

Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces

Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. We bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.

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.

Benzoquinoline Chemical Space: A Helpful Approach in Antibacterial and Anticancer Drug Design

Benzoquinolines are used in many drug design projects as starting molecules subject to derivatization. This computational study aims to characterize e benzoquinone drug space to ease future drug design processes based on these molecules. The drug space is composed of all benzoquinones, which are active on topoisomerase II and ATP synthase. Topological, chemical, and bioactivity spaces are explored using computational methodologies based on virtual screening and scaffold hopping and molecular docking, respectively. Topological space is a geometrical space in which the elements composing it can be defined as a set of neighbors (which satisfy a particular axiom). In such space, a chemical space can be defined as the property space spanned by all possible molecules and chemical compounds adhering to a given set of construction principles and boundary conditions. In this chemical space, the potentially pharmacologically active molecules form the bioactivity space. Results show a poly-morphological chemical space that suggests distinct characteristics. The chemical space is correlated with properties such as steric energy, the number of hydrogen bonds, the presence of halogen atoms, and membrane permeability-related properties. Lastly, novel chemical compounds (such as oxadiazole methybenzamide and floro methylcyclohexane diene) with drug-like potential, active on TOPO II and ATP synthase have been identified.

Chemically Driven Multimodal Locomotion of Active, Flexible Sheets

The inhibitor-promoter feedback loop is a vital component in regulatory pathways that controls functionality in living systems. In this loop, the production of chemical A at one site promotes the production of chemical B at another site, but B inhibits the production of A. In solution, differences in the volumes of the reactants and products of this reaction can generate buoyancy-driven fluid flows, which will deform neighboring soft material. To probe the intrinsic interrelationship among chemistry, hydrodynamics, and fluid-structure interactions, we model a bio-inspired system where a flexible sheet immersed in solution encompasses two spatially separated catalytic patches, which drive the A-B inhibitor-promotor reaction. The convective rolls of fluid generated above the patches can circulate inward or outward depending on the chemical environment. Within the regime displaying chemical oscillations, the dynamic fluid-structure interactions morph the shape of the sheet to periodically "fly", "crawl", or "swim" along the bottom of the confining chamber, revealing an intimate coupling between form and function in this system. The oscillations in the sheet's motion in turn affect the chemical oscillations in the solution. In the regime with non-oscillatory chemistry, the induced flow still morphs the shape of the sheet, but now, the fluid simply translates the sheet along the length of the chamber. The findings reveal the potential for enzymatic reactions in the body to generate hydrodynamic behavior that modifies the shape of neighboring soft tissue, which in turn modifies both the fluid dynamics and the enzymatic reaction. The findings indicate that this non-linear dynamic behavior can be playing a critical role in the functioning of regulatory pathways in living systems.

Pathfinder─Navigating and Analyzing Chemical Reaction Networks with an Efficient Graph-Based Approach

While the field of first-principles explorations into chemical reaction space has been continuously growing, the development of strategies for analyzing resulting chemical reaction networks (CRNs) is lagging behind. A CRN consists of compounds linked by reactions. Analyzing how these compounds are transformed into one another based on kinetic modeling is a nontrivial task. Here, we present the graph-optimization-driven algorithm and program Pathfinder to allow for such an analysis of a CRN. The CRN for this work has been obtained with our open-source Chemoton reaction network exploration software. Chemoton probes reactive combinations of compounds for elementary steps and sorts them into reactions. By encoding these reactions of the CRN as a graph consisting of compound and reaction vertices and adding information about activation barriers as well as required reagents to the edges of the graph yields a complete graph-theoretical representation of the CRN. Since the probabilities of the formation of compounds depend on the starting conditions, the consumption of any compound during a reaction must be accounted for to reflect the availability of reagents. To account for this, we introduce compound costs to reflect compound availability. Simultaneously, the determined compound costs rank the compounds in the CRN in terms of their probability to be formed. This ranking then allows us to probe easily accessible compounds in the CRN first for further explorations into yet unexplored terrain. We first illustrate the working principle on an abstract small CRN. Afterward, Pathfinder is demonstrated in the example of the disproportionation of iodine with water and the comproportionation of iodic acid and hydrogen iodide. Both processes are analyzed within the same CRN, which we construct with our autonomous first-principles CRN exploration software Chemoton [Unsleber, J. P.; J. Chem. Theory Comput. 2022, 18, 5393-5409] guided by Pathfinder.

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.

Chemical reaction networks and opportunities for machine learning

Chemical reaction networks (CRNs), defined by sets of species and possible reactions between them, are widely used to interrogate chemical systems. To capture increasingly complex phenomena, CRNs can be leveraged alongside data-driven methods and machine learning (ML). In this Perspective, we assess the diverse strategies available for CRN construction and analysis in pursuit of a wide range of scientific goals, discuss ML techniques currently being applied to CRNs and outline future CRN-ML approaches, presenting scientific and technical challenges to overcome.

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.

ReNeGate: A Reaction Network Graph-Theoretical Tool for Automated Mechanistic Studies in Computational Homogeneous Catalysis

Exploration of the chemical reaction space of chemical transformations in multicomponent mixtures is one of the main challenges in contemporary computational chemistry. To remove expert bias from mechanistic studies and to discover new chemistries, an automated graph-theoretical methodology is proposed, which puts forward a network formalism of homogeneous catalysis reactions and utilizes a network analysis tool for mechanistic studies. The method can be used for analyzing trajectories with single and multiple catalytic species and can provide unique conformers of catalysts including multinuclear catalyst clusters along with other catalytic mixture components. The presented three-step approach has the integrated ability to handle multicomponent catalytic systems of arbitrary complexity (mixtures of reactants, catalyst precursors, ligands, additives, and solvents). It is not limited to predefined chemical rules, does not require prealignment of reaction mixture components consistent with a reaction coordinate, and is not agnostic to the chemical nature of transformations. Conformer exploration, reactive event identification, and reaction network analysis are the main steps taken for identifying the pathways in catalytic systems given the starting precatalytic reaction mixture as the input. Such a methodology allows us to efficiently explore catalytic systems in realistic conditions for either previously observed or completely unknown reactive events in the context of a network representing different intermediates. Our workflow for the catalytic reaction space exploration exclusively focuses on the identification of thermodynamically feasible conversion channels, representative of the (secondary) catalyst deactivation or inhibition paths, which are usually most difficult to anticipate based solely on expert chemical knowledge. Thus, the expert bias is sought to be removed at all steps, and the chemical intuition is limited to the choice of the thermodynamic constraint imposed by the applicable experimental conditions in terms of threshold energy values for allowed transformations. The capabilities of the proposed methodology have been tested by exploring the reactivity of Mn complexes relevant for catalytic hydrogenation chemistry to verify previously postulated activation mechanisms and unravel unexpected reaction channels relevant to rare deactivation events.

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.

Fully Automated Generation of Prebiotically Relevant Reaction Networks from Optimized Nanoreactor Simulations

The nanoreactor approach first introduced by the group of Martı́nez [Wang et al. Nat. Chem. 2014, 6, 1044-1048] has recently attracted much attention because of its ability to accelerate the discovery of reaction pathways. Here, we provide a comprehensive study of various simulation parameters and present an alternative implementation for the reactivity-enhancing spherical constraint function, as well as for the detection of reaction events. In this context, a fully automated postsimulation evaluation procedure based on RDKit and NetworkX analysis is introduced. The chemical and physical robustness of the procedure is examined by investigating the reactivity of selected homogeneous systems. The optimized procedure is applied at the GFN2-xTB level of theory to a system composed of HCN molecules and argon atoms, acting as a buffer, yielding prebiotically plausible primary and secondary precursors for the synthesis of RNA. Furthermore, the formose reaction network is explored leading to numerous sugar precursors. The discovered compounds reflect experimental findings; however, new synthetic routes and a large collection of exotic, highly reactive molecules are observed, highlighting the predictive power of the nanoreactor approach for unraveling the reactive manifold.

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.

Network Analysis Reveals Spatial Clustering and Annotation of Complex Chemical Spaces: Application to Astrochemistry

How are molecules linked to each other in complex systems? In a proof-of-concept study, we have developed the method mol2net (https://zenodo.org/record/7025094) to generate and analyze the molecular network of complex astrochemical data (from high-resolution Orbitrap MS¹ analysis of H₂O:CH₃OH:NH₃ interstellar ice analogs) in a data-driven and unsupervised manner, without any prior knowledge about chemical reactions. The molecular network is clustered according to the initial NH₃ content and unlocked HCN, NH₃, and H₂O as spatially resolved key transformations. In comparison with the PubChem database, four subsets were annotated: (i) saturated C-backbone molecules without N, (ii) saturated N-backbone molecules, (iii) unsaturated C-backbone molecules without N, and (iv) unsaturated N-backbone molecules. These findings were validated with previous results (e.g., identifying the two major graph components as previously described N-poor and N-rich molecular groups) but with additional information about subclustering, key transformations, and molecular structures, and thus, the structural characterization of large complex organic molecules in interstellar ice analogs has been significantly refined.

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.

O-Acetyl Migration within the Sialic Acid Side Chain: A Mechanistic Study Using the Ab Initio Nanoreactor

Many disease-causing viruses target sialic acids on the surface of host cells. Some viruses bind preferentially to sialic acids with O-acetyl modification at the hydroxyl group of C7, C8, or C9 on the glycerol-like side chain. Studies of proteins binding to sialosides containing O-acetylated sialic acids are crucial in understanding the related diseases but experimentally difficult due to the lability of the ester group. We recently showed that O-acetyl migration among hydroxyl groups of C7, C8, and C9 in sialic acids occurs in all directions in a pH-dependent manner. In the current study, we elucidate a full mechanistic pathway for the migration of O-acetyl among C7, C8, and C9. We used an ab initio nanoreactor to explore potential reaction pathways and density functional theory, pK_a calculations, and umbrella sampling to investigate elementary steps of interest. We found that when a base is present, migration is easy in any direction and involves three key steps: deprotonation of the hydroxyl group, cyclization between the two carbons, and the migration of the O-acetyl group. This dynamic equilibrium may play a defensive role against pathogens that evolve to gain entry to the cell by binding selectively to one acetylation state.

Computational Prediction of Speciation Diagrams and Nucleation Mechanisms: Molecular Vanadium, Niobium, and Tantalum Oxide Nanoclusters in Solution

Understanding the aqueous speciation of molecular metal-oxo-clusters plays a key role in different fields such as catalysis, electrochemistry, nuclear waste recycling, and biochemistry. To describe the speciation accurately, it is essential to elucidate the underlying self-assembly processes. Herein, we apply a computational method to predict the speciation and formation mechanisms of polyoxovanadates, -niobates, and -tantalates. While polyoxovanadates have been widely studied, polyoxoniobates and -tantalates lack the same level of understanding. First, we propose a pentavanadate cluster ([V₅O₁₄]^3-) as a key intermediate for the formation of the decavanadate. Our computed phase speciation diagram is in particularly good agreement with the experiments. Second, we report the formation constants of the heptaniobate, [Nb₇O₂₂]^9-, decaniobate, [Nb₁₀O₂₈]^6-, and tetracosaniobate [H₉Nb₂₄O₇₂]^15-. Additionally, we compute the speciation and phase diagram of niobium, which so far was restricted to Lindqvist derivates. Finally, we predict the formation constant of the decatantalate ([Ta₁₀O₂₆]^6-) in water, even though it had only been synthesized in toluene. Furthermore, we also calculate the corresponding speciation and phase diagrams for polyoxotantalates. Overall, we show that our method can be successfully applied to different families of molecular metal oxides without any need for readjustments; therefore, it can be regarded as a trustworthy tool for exploring polyoxometalates' chemistry.

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.

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.

The chemistry of AlF and CaF production in buffer gas sources

In this work, we explore the role of chemical reactions on the properties of buffer gas cooled molecular beams. In particular, we focus on scenarios relevant to the formation of AlF and CaF via chemical reactions between the Ca and Al atoms ablated from a solid target in an atmosphere of a fluorine-containing gas, in this case, SF6 and NF3. Reactions are studied following an ab initio molecular dynamics approach, and the results are rationalized following a tree-shaped reaction model based on Bayesian inference. We find that NF3 reacts more efficiently with hot metal atoms to form monofluoride molecules than SF6. In addition, when using NF3, the reaction products have lower kinetic energy, requiring fewer collisions to thermalize with the cryogenic helium. Furthermore, we find that the reaction probability for AlF formation is much higher than for CaF across a broad range of temperatures.