Automatic Approach to Explore the Multireaction Mechanism for Medium-Sized Bimolecular Reactions via Collision Dynamics Simulations and Transition State Searches

Reactant Discovery with an Ab Initio Nanoreactor: Exploration of Astrophysical N-Heterocycle Precursors and Formation Pathways

The incorporation of nitrogen atoms into cyclic compounds is essential for terrestrial life; nitrogen-containing (N-)heterocycles make up DNA and RNA nucleobases, several amino acids, B vitamins, porphyrins, and other components of biomolecules. The discovery of these molecules on meteorites with non-terrestrial isotopic abundances supports the hypothesis of exogenous delivery of prebiotic material to early Earth; however, there has been no detection of these species in interstellar environments, indicating that there is a need for greater knowledge of their astrochemical formation and destruction pathways. Here, we present results of simulations of gas-phase pyrrole and pyridine formation from an ab initio nanoreactor, a first-principles molecular dynamics simulation method that accelerates reaction discovery by applying non-equilibrium forces that are agnostic to individual reaction coordinates. Using the nanoreactor in a retrosynthetic mode, starting with the N-heterocycle of interest and a radical leaving group, then considering the discovered reaction pathways in reverse, a rich landscape of N-heterocycle-forming reactivity can be found. Several of these reaction pathways, when mapped to their corresponding minimum energy paths, correspond to novel barrierless formation pathways for pyridine and pyrrole, starting from both detected and hypothesized astrochemical precursors. This study demonstrates how first-principles reaction discovery can build mechanistic knowledge in astrochemical environments as well as in early Earth models such as Titan's atmosphere where N-heterocycles have been tentatively detected.

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.

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.

Discover, Sample, and Refine: Exploring Chemistry with Enhanced Sampling Techniques

Over the last few decades, enhanced sampling methods have been continuously improved. Here, we exploit this progress and propose a modular workflow for blind reaction discovery and determination of reaction paths. In a three-step strategy, at first we use a collective variable derived from spectral graph theory in conjunction with the explore variant of the on-the-fly probability enhanced sampling method to drive reaction discovery runs. Once different chemical products are determined, we construct an ad-hoc neural network-based collective variable to improve sampling, and finally we refine the results using the free energy perturbation theory and a more accurate Hamiltonian. We apply this strategy to both intramolecular and intermolecular reactions. Our workflow requires minimal user input and extends the power of ab initio molecular dynamics to explore and characterize the reaction space.