1
|
Laplaza R, Wodrich MD, Corminboeuf C. Overcoming the Pitfalls of Computing Reaction Selectivity from Ensembles of Transition States. J Phys Chem Lett 2024; 15:7363-7370. [PMID: 38990895 DOI: 10.1021/acs.jpclett.4c01657] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
The prediction of reaction selectivity is a challenging task for computational chemistry, not only because many molecules adopt multiple conformations but also due to the exponential relationship between effective activation energies and rate constants. To account for molecular flexibility, an increasing number of methods exist that generate conformational ensembles of transition state (TS) structures. Typically, these TS ensembles are Boltzmann weighted and used to compute selectivity assuming Curtin-Hammett conditions. This strategy, however, can lead to erroneous predictions if the appropriate filtering of the conformer ensembles is not conducted. Here, we demonstrate how any possible selectivity can be obtained by processing the same sets of TS ensembles for a model reaction. To address the burdensome filtering task in a consistent and automated way, we introduce marc, a tool for the modular analysis of representative conformers that aids in avoiding human errors while minimizing the number of reoptimization computations needed to obtain correct reaction selectivity.
Collapse
Affiliation(s)
- Ruben Laplaza
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Matthew D Wodrich
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Clemence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
2
|
Csizi KS, Steiner M, Reiher M. Nanoscale chemical reaction exploration with a quantum magnifying glass. Nat Commun 2024; 15:5320. [PMID: 38909029 PMCID: PMC11193806 DOI: 10.1038/s41467-024-49594-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Accepted: 06/04/2024] [Indexed: 06/24/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Katja-Sophia Csizi
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
- ETH Zurich, NCCR Catalysis, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland.
- ETH Zurich, NCCR Catalysis, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland.
| |
Collapse
|
3
|
Weymuth T, Unsleber JP, Türtscher PL, Steiner M, Sobez JG, Müller CH, Mörchen M, Klasovita V, Grimmel SA, Eckhoff M, Csizi KS, Bosia F, Bensberg M, Reiher M. SCINE-Software for chemical interaction networks. J Chem Phys 2024; 160:222501. [PMID: 38857173 DOI: 10.1063/5.0206974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/09/2024] [Indexed: 06/12/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Thomas Weymuth
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P Unsleber
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L Türtscher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan-Grimo Sobez
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Charlotte H Müller
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Veronika Klasovita
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A Grimmel
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Eckhoff
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Katja-Sophia Csizi
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Francesco Bosia
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Moritz Bensberg
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
4
|
Bensberg M, Reiher M. Uncertainty-Aware First-Principles Exploration of Chemical Reaction Networks. J Phys Chem A 2024; 128:4532-4547. [PMID: 38787736 PMCID: PMC11163430 DOI: 10.1021/acs.jpca.3c08386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2023] [Revised: 05/13/2024] [Accepted: 05/13/2024] [Indexed: 05/26/2024]
Abstract
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.
Collapse
Affiliation(s)
- Moritz Bensberg
- Department of Chemistry and Applied
Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Department of Chemistry and Applied
Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
5
|
Steiner M, Reiher M. A human-machine interface for automatic exploration of chemical reaction networks. Nat Commun 2024; 15:3680. [PMID: 38693117 PMCID: PMC11063077 DOI: 10.1038/s41467-024-47997-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 04/15/2024] [Indexed: 05/03/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
- ETH Zurich, NCCR Catalysis, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland.
- ETH Zurich, NCCR Catalysis, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland.
| |
Collapse
|
6
|
Hutton DJ, Cordes KE, Michel C, Göltl F. Machine Learning-Based Prediction of Activation Energies for Chemical Reactions on Metal Surfaces. J Chem Inf Model 2023; 63:6006-6013. [PMID: 37722106 DOI: 10.1021/acs.jcim.3c00740] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/20/2023]
Abstract
In computational surface catalysis, the calculation of activation energies of chemical reactions is expensive, which, in many cases, limits our ability to understand complex reaction networks. Here, we present a universal, machine learning-based approach for the prediction of activation energies for reactions of C-, O-, and H-containing molecules on transition metal surfaces. We rely on generalized Bronsted-Evans-Polanyi relationships in combination with machine learning-based multiparameter regression techniques to train our model for reactions included in the University of Arizona Reaction database. In our best approach, we find a mean absolute error for activation energies within our test set of 0.14 eV if the reaction energy is known and 0.19 eV if the reaction energy is unknown. We expect that this methodology will often replace the explicit calculation of activation energies within surface catalysis when exploring large reaction networks or screening catalysts for desirable properties in the future.
Collapse
Affiliation(s)
- Daniel J Hutton
- Department of Biosystems Engineering, The University of Arizona, 1177 E. Fourth St., Tucson, Arizona 85719, United States
| | - Kari E Cordes
- Department of Biosystems Engineering, The University of Arizona, 1177 E. Fourth St., Tucson, Arizona 85719, United States
| | - Carine Michel
- ENSL, CNRS, Laboratoire de Chimie UMR 5182, 46 Allée d'Italie, F69364 Lyon, France
| | - Florian Göltl
- Department of Biosystems Engineering, The University of Arizona, 1177 E. Fourth St., Tucson, Arizona 85719, United States
| |
Collapse
|
7
|
Zhao Q, Savoie BM. Deep reaction network exploration of glucose pyrolysis. Proc Natl Acad Sci U S A 2023; 120:e2305884120. [PMID: 37579176 PMCID: PMC10450414 DOI: 10.1073/pnas.2305884120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Accepted: 07/03/2023] [Indexed: 08/16/2023] Open
Abstract
Resolving the reaction networks associated with biomass pyrolysis is central to understanding product selectivity and aiding catalyst design to produce more valuable products. However, even the pyrolysis network of relatively simple [Formula: see text]-D-glucose remains unresolved due to its significant complexity in terms of the depth of the network and the number of major products. Here, a transition-state-guided reaction exploration has been performed that provides complete pathways to most significant experimental pyrolysis products of [Formula: see text]-D-glucose. The resulting reaction network involves over 31,000 reactions and transition states computed at the semiempirical quantum chemistry level and approximately 7,000 kinetically relevant reactions and transition states characterized with density function theory, comprising the largest reaction network reported for biomass pyrolysis. The exploration was conducted using graph-based rules to explore the reactivities of intermediates and an adaption of the Dijkstra algorithm to identify kinetically relevant intermediates. This simple exploration policy surprisingly (re)identified pathways to most major experimental pyrolysis products, many intermediates proposed by previous computational studies, and also identified new low-barrier reaction mechanisms that resolve outstanding discrepancies between reaction pathways and yields in isotope labeling experiments. This network also provides explanatory pathways for the high yield of hydroxymethylfurfural and the reaction pathway that contributes most to the formation of hydroxyacetaldehyde during glucose pyrolysis. Due to the limited domain knowledge required to generate this network, this approach should also be transferable to other complex reaction network prediction problems in biomass pyrolysis.
Collapse
Affiliation(s)
- Qiyuan Zhao
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN47906
| | - Brett M. Savoie
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN47906
| |
Collapse
|
8
|
Toniato A, Unsleber JP, Vaucher AC, Weymuth T, Probst D, Laino T, Reiher M. Quantum chemical data generation as fill-in for reliability enhancement of machine-learning reaction and retrosynthesis planning. DIGITAL DISCOVERY 2023; 2:663-673. [PMID: 37312681 PMCID: PMC10259370 DOI: 10.1039/d3dd00006k] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/21/2023] [Accepted: 03/09/2023] [Indexed: 06/15/2023]
Abstract
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.
Collapse
Affiliation(s)
- Alessandra Toniato
- Laboratory of Physical Chemistry, ETH Zurich Vladimir-Prelog-Weg 2 8093 Zurich Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), ETH Zurich Vladimir-Prelog-Weg 1-5/10 8093 Zurich Switzerland
- IBM Research Europe 8803 Rüschlikon Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), IBM Research 8803 Rüschlikon Switzerland
| | - Jan P Unsleber
- Laboratory of Physical Chemistry, ETH Zurich Vladimir-Prelog-Weg 2 8093 Zurich Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), ETH Zurich Vladimir-Prelog-Weg 1-5/10 8093 Zurich Switzerland
| | - Alain C Vaucher
- IBM Research Europe 8803 Rüschlikon Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), IBM Research 8803 Rüschlikon Switzerland
| | - Thomas Weymuth
- Laboratory of Physical Chemistry, ETH Zurich Vladimir-Prelog-Weg 2 8093 Zurich Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), ETH Zurich Vladimir-Prelog-Weg 1-5/10 8093 Zurich Switzerland
| | - Daniel Probst
- IBM Research Europe 8803 Rüschlikon Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), IBM Research 8803 Rüschlikon Switzerland
| | - Teodoro Laino
- IBM Research Europe 8803 Rüschlikon Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), IBM Research 8803 Rüschlikon Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich Vladimir-Prelog-Weg 2 8093 Zurich Switzerland
- National Center for Competence in Research-Catalysis (NCCR Catalysis), ETH Zurich Vladimir-Prelog-Weg 1-5/10 8093 Zurich Switzerland
| |
Collapse
|
9
|
Kraka E, Antonio JJ, Freindorf M. Reaction mechanism - explored with the unified reaction valley approach. Chem Commun (Camb) 2023; 59:7151-7165. [PMID: 37233449 DOI: 10.1039/d3cc01576a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
One of the ultimate goals of chemistry is to understand and manipulate chemical reactions, which implies the ability to monitor the reaction and its underlying mechanism at an atomic scale. In this article, we introduce the Unified Reaction Valley Approach (URVA) as a tool for elucidating reaction mechanisms, complementing existing computational procedures. URVA combines the concept of the potential energy surface with vibrational spectroscopy and describes a chemical reaction via the reaction path and the surrounding reaction valley traced out by the reacting species on the potential energy surface on their way from the entrance to the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting species is registered by a change in the normal vibrational modes spanning the reaction valley and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction, with curvature minima reflecting minimal change and curvature maxima indicating the location of important chemical events such as bond breaking/formation, charge polarization and transfer, rehybridization, etc. A decomposition of the path curvature into internal coordinate components or other coordinates of relevance for the reaction under consideration, provides comprehensive insight into the origin of the chemical changes taking place. After giving an overview of current experimental and computational efforts to gain insight into the mechanism of a chemical reaction and presenting the theoretical background of URVA, we illustrate how URVA works for three diverse processes, (i) [1,3] hydrogen transfer reactions; (ii) α-keto-amino inhibitor for SARS-CoV-2 Mpro; (iii) Rh-catalyzed cyanation. We hope that this article will inspire our computational colleagues to add URVA to their repertoire and will serve as an incubator for new reaction mechanisms to be studied in collaboration with our experimental experts in the field.
Collapse
Affiliation(s)
- Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA.
| | - Juliana J Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA.
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA.
| |
Collapse
|
10
|
Unsleber JP. Accelerating Reaction Network Explorations with Automated Reaction Template Extraction and Application. J Chem Inf Model 2023. [PMID: 37216641 DOI: 10.1021/acs.jcim.3c00102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jan P Unsleber
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
11
|
Pasquini M, Stenta M. LinChemIn: SynGraph-a data model and a toolkit to analyze and compare synthetic routes. J Cheminform 2023; 15:41. [PMID: 37005691 PMCID: PMC10067316 DOI: 10.1186/s13321-023-00714-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Accepted: 03/20/2023] [Indexed: 04/04/2023] Open
Abstract
BACKGROUND The increasing amount of chemical reaction data makes traditional ways to navigate its corpus less effective, while the demand for novel approaches and instruments is rising. Recent data science and machine learning techniques support the development of new ways to extract value from the available reaction data. On the one side, Computer-Aided Synthesis Planning tools can predict synthetic routes in a model-driven approach; on the other side, experimental routes can be extracted from the Network of Organic Chemistry, in which reaction data are linked in a network. In this context, the need to combine, compare and analyze synthetic routes generated by different sources arises naturally. RESULTS Here we present LinChemIn, a python toolkit that allows chemoinformatics operations on synthetic routes and reaction networks. Wrapping some third-party packages for handling graph arithmetic and chemoinformatics and implementing new data models and functionalities, LinChemIn allows the interconversion between data formats and data models and enables route-level analysis and operations, including route comparison and descriptors calculation. Object-Oriented Design principles inspire the software architecture, and the modules are structured to maximize code reusability and support code testing and refactoring. The code structure should facilitate external contributions, thus encouraging open and collaborative software development. CONCLUSIONS The current version of LinChemIn allows users to combine synthetic routes generated from various tools and analyze them, and constitutes an open and extensible framework capable of incorporating contributions from the community and fostering scientific discussion. Our roadmap envisages the development of sophisticated metrics for routes evaluation, a multi-parameter scoring system, and the implementation of an entire "ecosystem" of functionalities operating on synthetic routes. LinChemIn is freely available at https://github.com/syngenta/linchemin.
Collapse
Affiliation(s)
- Marta Pasquini
- Syngenta Crop Protection AG, Schaffhauserstrasse, 4332, Stein, AG, Switzerland.
| | - Marco Stenta
- Syngenta Crop Protection AG, Schaffhauserstrasse, 4332, Stein, AG, Switzerland
| |
Collapse
|
12
|
Zhao Q, Vaddadi SM, Woulfe M, Ogunfowora LA, Garimella SS, Isayev O, Savoie BM. Comprehensive exploration of graphically defined reaction spaces. Sci Data 2023; 10:145. [PMID: 36935430 PMCID: PMC10025260 DOI: 10.1038/s41597-023-02043-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Accepted: 02/27/2023] [Indexed: 03/21/2023] Open
Abstract
Existing reaction transition state (TS) databases are comparatively small and lack chemical diversity. Here, this data gap has been addressed using the concept of a graphically-defined model reaction to comprehensively characterize a reaction space associated with C, H, O, and N containing molecules with up to 10 heavy (non-hydrogen) atoms. The resulting dataset is composed of 176,992 organic reactions possessing at least one validated TS, activation energy, heat of reaction, reactant and product geometries, frequencies, and atom-mapping. For 33,032 reactions, more than one TS was discovered by conformational sampling, allowing conformational errors in TS prediction to be assessed. Data is supplied at the GFN2-xTB and B3LYP-D3/TZVP levels of theory. A subset of reactions were recalculated at the CCSD(T)-F12/cc-pVDZ-F12 and ωB97X-D2/def2-TZVP levels to establish relative errors. The resulting collection of reactions and properties are called the Reaction Graph Depth 1 (RGD1) dataset. RGD1 represents the largest and most chemically diverse TS dataset published to date and should find immediate use in developing novel machine learning models for predicting reaction properties.
Collapse
Affiliation(s)
- Qiyuan Zhao
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Sai Mahit Vaddadi
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Michael Woulfe
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Lawal A Ogunfowora
- Department of Chemistry, Purdue University, West Lafayette, IN, 47906, USA
| | - Sanjay S Garimella
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Brett M Savoie
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA.
| |
Collapse
|
13
|
Concentration‐Flux‐Steered Mechanism Exploration with an Organocatalysis Application. Isr J Chem 2023. [DOI: 10.1002/ijch.202200123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023]
|
14
|
Unsleber JP, Liu H, Talirz L, Weymuth T, Mörchen M, Grofe A, Wecker D, Stein CJ, Panyala A, Peng B, Kowalski K, Troyer M, Reiher M. High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation. J Chem Phys 2023; 158:084803. [PMID: 36859110 DOI: 10.1063/5.0136526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Jan P Unsleber
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Hongbin Liu
- Microsoft Quantum, Redmond, Washington 98052, USA
| | | | - Thomas Weymuth
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Adam Grofe
- Microsoft Quantum, Redmond, Washington 98052, USA
| | - Dave Wecker
- Microsoft Quantum, Redmond, Washington 98052, USA
| | - Christopher J Stein
- Department of Chemistry, TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, D-85748 Garching, Germany
| | - Ajay Panyala
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Bo Peng
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Karol Kowalski
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | | | - Markus Reiher
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
15
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
16
|
Abstract
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.
Collapse
Affiliation(s)
- Maike Vahl
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany.
| | - Jonny Proppe
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany.
| |
Collapse
|
17
|
Türtscher PL, Reiher M. Pathfinder─Navigating and Analyzing Chemical Reaction Networks with an Efficient Graph-Based Approach. J Chem Inf Model 2023; 63:147-160. [PMID: 36515968 PMCID: PMC9832502 DOI: 10.1021/acs.jcim.2c01136] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
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.
Collapse
|
18
|
Wen M, Spotte-Smith EWC, Blau SM, McDermott MJ, Krishnapriyan AS, Persson KA. Chemical reaction networks and opportunities for machine learning. NATURE COMPUTATIONAL SCIENCE 2023; 3:12-24. [PMID: 38177958 DOI: 10.1038/s43588-022-00369-z] [Citation(s) in RCA: 18] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Accepted: 11/08/2022] [Indexed: 01/06/2024]
Abstract
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.
Collapse
Affiliation(s)
- Mingjian Wen
- Chemical and Biomolecular Engineering, University of Houston, Houston, TX, USA
- Energy Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Evan Walter Clark Spotte-Smith
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Materials Science and Engineering, University of California, Berkeley, Berkeley, CA, USA
| | - Samuel M Blau
- Energy Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Matthew J McDermott
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Materials Science and Engineering, University of California, Berkeley, Berkeley, CA, USA
| | - Aditi S Krishnapriyan
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, CA, USA
- Electrical Engineering and Computer Science, University of California, Berkeley, Berkeley, CA, USA
| | - Kristin A Persson
- Materials Science and Engineering, University of California, Berkeley, Berkeley, CA, USA.
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, CA, USA.
| |
Collapse
|
19
|
Hashemi A, Bougueroua S, Gaigeot MP, Pidko EA. ReNeGate: A Reaction Network Graph-Theoretical Tool for Automated Mechanistic Studies in Computational Homogeneous Catalysis. J Chem Theory Comput 2022; 18:7470-7482. [PMID: 36321652 DOI: 10.1021/acs.jctc.2c00404] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
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.
Collapse
Affiliation(s)
- Ali Hashemi
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, Delft 2629 HZ, The Netherlands
| | - Sana Bougueroua
- Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement (LAMBE) UMR8587, Universite Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, Evry-Courcouronnes 91025, France
| | - Marie-Pierre Gaigeot
- Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement (LAMBE) UMR8587, Universite Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, Evry-Courcouronnes 91025, France
| | - Evgeny A Pidko
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, Delft 2629 HZ, The Netherlands
| |
Collapse
|
20
|
Xu T, Wang Y, Xiong Z, Wang Y, Zhou Y, Li X. A Rising 2D Star: Novel MBenes with Excellent Performance in Energy Conversion and Storage. NANO-MICRO LETTERS 2022; 15:6. [PMID: 36472760 PMCID: PMC9727130 DOI: 10.1007/s40820-022-00976-5] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Accepted: 10/26/2022] [Indexed: 06/17/2023]
Abstract
As a flourishing member of the two-dimensional (2D) nanomaterial family, MXenes have shown great potential in various research areas. In recent years, the continued growth of interest in MXene derivatives, 2D transition metal borides (MBenes), has contributed to the emergence of this 2D material as a latecomer. Due to the excellent electrical conductivity, mechanical properties and electrical properties, thus MBenes attract more researchers' interest. Extensive experimental and theoretical studies have shown that they have exciting energy conversion and electrochemical storage potential. However, a comprehensive and systematic review of MBenes applications has not been available so far. For this reason, we present a comprehensive summary of recent advances in MBenes research. We started by summarizing the latest fabrication routes and excellent properties of MBenes. The focus will then turn to their exciting potential for energy storage and conversion. Finally, a brief summary of the challenges and opportunities for MBenes in future practical applications is presented.
Collapse
Affiliation(s)
- Tianjie Xu
- Hubei Province Key Laboratory of Science in Metallurgical Process, Wuhan University of Science and Technology, Wuhan, 430081, People's Republic of China
| | - Yuhua Wang
- Hubei Province Key Laboratory of Science in Metallurgical Process, Wuhan University of Science and Technology, Wuhan, 430081, People's Republic of China.
| | - Zuzhao Xiong
- Hubei Province Key Laboratory of Science in Metallurgical Process, Wuhan University of Science and Technology, Wuhan, 430081, People's Republic of China
| | - Yitong Wang
- Hubei Province Key Laboratory of Science in Metallurgical Process, Wuhan University of Science and Technology, Wuhan, 430081, People's Republic of China
| | - Yujin Zhou
- Hubei Province Key Laboratory of Science in Metallurgical Process, Wuhan University of Science and Technology, Wuhan, 430081, People's Republic of China
| | - Xifei Li
- Institute of Advanced Electrochemical Energy and School of Materials Science and Engineering, Xi'an University of Technology, Xi'an, 710048, People's Republic of China.
- Center for International Cooperation On Designer Low-Carbon and Environmental Materials (CDLCEM), Zhengzhou University, Zhengzhou, 450001, Henan, People's Republic of China.
| |
Collapse
|
21
|
Zhao Q, Xu Y, Greeley J, Savoie BM. Deep reaction network exploration at a heterogeneous catalytic interface. Nat Commun 2022; 13:4860. [PMID: 35982057 PMCID: PMC9388529 DOI: 10.1038/s41467-022-32514-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Accepted: 08/03/2022] [Indexed: 11/09/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Qiyuan Zhao
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Yinan Xu
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Jeffrey Greeley
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA.
| | - Brett M Savoie
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA.
| |
Collapse
|
22
|
Cohen M, Vlachos DG. Modified Energy Span Analysis of Catalytic Parallel Pathways and Selectivity. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c01991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Maximilian Cohen
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
| | - Dionisios G. Vlachos
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
| |
Collapse
|
23
|
Duan C, Nandy A, Adamji H, Roman-Leshkov Y, Kulik HJ. Machine Learning Models Predict Calculation Outcomes with the Transferability Necessary for Computational Catalysis. J Chem Theory Comput 2022; 18:4282-4292. [PMID: 35737587 DOI: 10.1021/acs.jctc.2c00331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Virtual high-throughput screening (VHTS) and machine learning (ML) have greatly accelerated the design of single-site transition-metal catalysts. VHTS of catalysts, however, is often accompanied with a high calculation failure rate and wasted computational resources due to the difficulty of simultaneously converging all mechanistically relevant reactive intermediates to expected geometries and electronic states. We demonstrate a dynamic classifier approach, i.e., a convolutional neural network that monitors geometry optimizations on the fly, and exploit its good performance and transferability in identifying geometry optimization failures for catalyst design. We show that the dynamic classifier performs well on all reactive intermediates in the representative catalytic cycle of the radical rebound mechanism for the conversion of methane to methanol despite being trained on only one reactive intermediate. The dynamic classifier also generalizes to chemically distinct intermediates and metal centers absent from the training data without loss of accuracy or model confidence. We rationalize this superior model transferability as arising from the use of electronic structure and geometric information generated on-the-fly from density functional theory calculations and the convolutional layer in the dynamic classifier. When used in combination with uncertainty quantification, the dynamic classifier saves more than half of the computational resources that would have been wasted on unsuccessful calculations for all reactive intermediates being considered.
Collapse
Affiliation(s)
- Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Husain Adamji
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Yuriy Roman-Leshkov
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
24
|
Eisenstein O. From the Felkin‐Anh Rule to the Grignard Reaction: an Almost Circular 50 Year Adventure in the World of Molecular Structures and Reaction Mechanisms with Computational Chemistry**. Isr J Chem 2022. [DOI: 10.1002/ijch.202100138] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Odile Eisenstein
- ICGM, Univ. Montpellier, CNRS, ENSCM, Montpellier, 34095 France Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences University of Oslo Oslo 0315 Norway
| |
Collapse
|
25
|
Trunschke A. Prospects and challenges for autonomous catalyst discovery viewed from an experimental perspective. Catal Sci Technol 2022. [DOI: 10.1039/d2cy00275b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Autonomous catalysis research requires elaborate integration of operando experiments into automated workflows. Suitable experimental data for analysis by artificial intelligence can be measured more readily according to standard operating procedures.
Collapse
Affiliation(s)
- Annette Trunschke
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Department of Inorganic Chemistry, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
26
|
Laplaza R, Sobez JG, Wodrich MD, Reiher M, Corminboeuf C. The (not so) simple prediction of enantioselectivity – a pipeline for high-fidelity computations. Chem Sci 2022; 13:6858-6864. [PMID: 35774159 PMCID: PMC9200111 DOI: 10.1039/d2sc01714h] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 05/17/2022] [Indexed: 11/21/2022] Open
Abstract
The computation of reaction selectivity represents an appealing complementary route to experimental studies and a powerful means to refine catalyst design strategies. Accurately establishing the selectivity of reactions facilitated by molecular catalysts, however, remains a challenging task for computational chemistry. The small free energy differences that lead to large variations in the enantiomeric ratio (er) represent particularly tricky quantities to predict with sufficient accuracy to be helpful for prioritizing experiments. Further complicating this problem is the fact that standard approaches typically consider only one or a handful of conformers identified through human intuition as pars pro toto of the conformational space. Obviously, this assumption can potentially lead to dramatic failures should key energetic low-lying structures be missed. Here, we introduce a multi-level computational pipeline leveraging the graph-based Molassembler library to construct an ensemble of molecular catalysts. The manipulation and interpretation of molecules as graphs provides a powerful and direct route to tailored functionalization and conformer generation that facilitates high-throughput mechanistic investigations of chemical reactions. The capabilities of this approach are validated by examining a Rh(iii) catalyzed asymmetric C–H activation reaction and assessing the limitations associated with the underlying ligand design model. Specifically, the presence of remarkably flexible chiral Cp ligands, which induce the experimentally observed high level of selectivity, present a rich configurational landscape where multiple unexpected conformations contribute to the reported enantiomeric ratios (er). Using Molassembler, we show that considering about 20 transition state conformations per catalysts, which are generated with little human intervention and are not tied to “back-of-the-envelope” models, accurately reproduces experimental er values with limited computational expense. The computation of reaction selectivity represents an appealing complementary route to experimental studies and a powerful mean to refine catalyst design strategies.![]()
Collapse
Affiliation(s)
- Rubén Laplaza
- Laboratory for Computational Molecular Design (LCMD), Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Jan-Grimo Sobez
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), ETH Zürich, Vladimir-Prelog-Weg 1-5/10, 8093 Zürich, Switzerland
| | - Matthew D. Wodrich
- Laboratory for Computational Molecular Design (LCMD), Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), ETH Zürich, Vladimir-Prelog-Weg 1-5/10, 8093 Zürich, Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design (LCMD), Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- National Center for Competence in Research-Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| |
Collapse
|