1
|
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
|
2
|
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
|
3
|
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
|
4
|
García-Andrade X, García Tahoces P, Pérez-Ríos J, Martínez Núñez E. Barrier Height Prediction by Machine Learning Correction of Semiempirical Calculations. J Phys Chem A 2023; 127:2274-2283. [PMID: 36877614 PMCID: PMC10845151 DOI: 10.1021/acs.jpca.2c08340] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 02/19/2023] [Indexed: 03/07/2023]
Abstract
Different machine learning (ML) models are proposed in the present work to predict density functional theory-quality barrier heights (BHs) from semiempirical quantum mechanical (SQM) calculations. The ML models include a multitask deep neural network, gradient-boosted trees by means of the XGBoost interface, and Gaussian process regression. The obtained mean absolute errors are similar to those of previous models considering the same number of data points. The ML corrections proposed in this paper could be useful for rapid screening of the large reaction networks that appear in combustion chemistry or in astrochemistry. Finally, our results show that 70% of the features with the highest impact on model output are bespoke predictors. This custom-made set of predictors could be employed by future Δ-ML models to improve the quantitative prediction of other reaction properties.
Collapse
Affiliation(s)
| | - Pablo García Tahoces
- Department
of Electronics and Computer Science, University
of Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Jesús Pérez-Ríos
- Department
of Physics, Stony Brook University, Stony Brook, New York 11794, United States
- Institute
for Advanced Computational Science, Stony
Brook University, Stony
Brook, New York 11794-3800, United States
| | - Emilio Martínez Núñez
- Department
of Physical Chemistry, University of Santiago
de Compostela, Santiago
de Compostela 15782, Spain
| |
Collapse
|
5
|
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]
|
6
|
A computational and experimental examination of the CID of phosphorylated serine-H +. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023]
|
7
|
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
|
8
|
Unsleber JP, Grimmel SA, Reiher M. Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks. J Chem Theory Comput 2022; 18:5393-5409. [PMID: 35926118 DOI: 10.1021/acs.jctc.2c00193] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Jan P Unsleber
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Stephanie A Grimmel
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
9
|
Garay-Ruiz D, Álvarez-Moreno M, Bo C, Martínez-Núñez E. New Tools for Taming Complex Reaction Networks: The Unimolecular Decomposition of Indole Revisited. ACS PHYSICAL CHEMISTRY AU 2022; 2:225-236. [PMID: 36855573 PMCID: PMC9718323 DOI: 10.1021/acsphyschemau.1c00051] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The level of detail attained in the computational description of reaction mechanisms can be vastly improved through tools for automated chemical space exploration, particularly for systems of small to medium size. Under this approach, the unimolecular decomposition landscape for indole was explored through the automated reaction mechanism discovery program AutoMeKin. Nevertheless, the sheer complexity of the obtained mechanisms might be a hindrance regarding their chemical interpretation. In this spirit, the new Python library amk-tools has been designed to read and manipulate complex reaction networks, greatly simplifying their overall analysis. The package provides interactive dashboards featuring visualizations of the network, the three-dimensional (3D) molecular structures and vibrational normal modes of all chemical species, and the corresponding energy profiles for selected pathways. The combination of the joined mechanism generation and postprocessing workflow with the rich chemistry of indole decomposition enabled us to find new details of the reaction (obtained at the CCSD(T)/aug-cc-pVTZ//M06-2X/MG3S level of theory) that were not reported before: (i) 16 pathways leading to the formation of HCN and NH3 (via amino radical); (ii) a barrierless reaction between methylene radical and phenyl isocyanide, which might be an operative mechanism under the conditions of the interstellar medium; and (iii) reaction channels leading to both hydrogen cyanide and hydrogen isocyanide, of potential astrochemical interest as the computed HNC/HCN ratios greatly exceed the calculated equilibrium value at very low temperatures. The reported reaction networks can be very valuable to supplement databases of kinetic data, which is of remarkable interest for pyrolysis and astrochemical studies.
Collapse
Affiliation(s)
- Diego Garay-Ruiz
- Institute
of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Avinguda Països Catalans,
16, 43007 Tarragona, Spain,Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili (URV), Marcel·lí Domingo s/n, 43007 Tarragona, Spain
| | - Moises Álvarez-Moreno
- Institute
of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Avinguda Països Catalans,
16, 43007 Tarragona, Spain
| | - Carles Bo
- Institute
of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Avinguda Països Catalans,
16, 43007 Tarragona, Spain,Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili (URV), Marcel·lí Domingo s/n, 43007 Tarragona, Spain,
| | - Emilio Martínez-Núñez
- Departmento
de Química Física, Facultade de Química, Universidade de Santiago de Compostela, 15782 Santiago
de Compostela, Spain,
| |
Collapse
|
10
|
Cui Q, Peng J, Xu C, Lan Z. Automatic Approach to Explore the Multireaction Mechanism for Medium-Sized Bimolecular Reactions via Collision Dynamics Simulations and Transition State Searches. J Chem Theory Comput 2022; 18:910-924. [DOI: 10.1021/acs.jctc.1c00795] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qinghai Cui
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
11
|
Steiner M, Reiher M. Autonomous Reaction Network Exploration in Homogeneous and Heterogeneous Catalysis. Top Catal 2022; 65:6-39. [PMID: 35185305 PMCID: PMC8816766 DOI: 10.1007/s11244-021-01543-9] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/17/2021] [Indexed: 12/11/2022]
Abstract
Autonomous computations that rely on automated reaction network elucidation algorithms may pave the way to make computational catalysis on a par with experimental research in the field. Several advantages of this approach are key to catalysis: (i) automation allows one to consider orders of magnitude more structures in a systematic and open-ended fashion than what would be accessible by manual inspection. Eventually, full resolution in terms of structural varieties and conformations as well as with respect to the type and number of potentially important elementary reaction steps (including decomposition reactions that determine turnover numbers) may be achieved. (ii) Fast electronic structure methods with uncertainty quantification warrant high efficiency and reliability in order to not only deliver results quickly, but also to allow for predictive work. (iii) A high degree of autonomy reduces the amount of manual human work, processing errors, and human bias. Although being inherently unbiased, it is still steerable with respect to specific regions of an emerging network and with respect to the addition of new reactant species. This allows for a high fidelity of the formalization of some catalytic process and for surprising in silico discoveries. In this work, we first review the state of the art in computational catalysis to embed autonomous explorations into the general field from which it draws its ingredients. We then elaborate on the specific conceptual issues that arise in the context of autonomous computational procedures, some of which we discuss at an example catalytic system.
Collapse
Affiliation(s)
- Miguel Steiner
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
12
|
Baiardi A, Grimmel SA, Steiner M, Türtscher PL, Unsleber JP, Weymuth T, Reiher M. Expansive Quantum Mechanical Exploration of Chemical Reaction Paths. Acc Chem Res 2022; 55:35-43. [PMID: 34918903 DOI: 10.1021/acs.accounts.1c00472] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Quantum mechanical methods have been well-established for the elucidation of reaction paths of chemical processes and for the explicit dynamics of molecular systems. While they are usually deployed in routine manual calculations on reactions for which some insights are already available (typically from experiment), new algorithms and continuously increasing capabilities of modern computer hardware allow for exploratory open-ended computational campaigns that are unbiased and therefore enable unexpected discoveries. Highly efficient and even automated procedures facilitate systematic approaches toward the exploration of uncharted territory in molecular transformations and dynamics. In this work, we elaborate on such explorative approaches that range from reaction network explorations with (stationary) quantum chemical methods to explorative molecular dynamics and migrant wave packet dynamics. The focus is on recent developments that cover the following strategies. (i) Pruning search options for elementary reaction steps by heuristic rules based on the first-principles of quantum mechanics: Rules are required for reducing the combinatorial explosion of potentially reactive atom pairings, and rooting them in concepts derived from the electronic wave function makes them applicable to any molecular system. (ii) Enforcing reactive events by external biases: Inducing a reaction requires constraints that steer and direct elementary-step searches, which can be formulated in terms of forces, velocities, or supplementary potentials. (iii) Manual steering facilitated by interactive quantum mechanics: As ultrafast quantum chemical methods allow for real-time manual interactions with molecular systems, human-intuition-guided paths can be easily explored with suitable human-machine interfaces. (iv) New approaches for transition-state optimization with continuous curve representations can provide stable schemes to be driven in an automated way by allowing for an efficient tuning of the curve's parameters (instead of a manipulation of a collection of structures along the path), and (v) reactive molecular dynamics and direct wave packet propagation exploit the equations of motion of an underlying mechanical theory (usually, classical Newtonian mechanics or Schrödinger quantum mechanics). Explorative approaches are likely to replace the current state of the art in computational chemistry, because they reduce the human effort to be invested in reaction path elucidations, they are less prone to errors and bias-free, and they cover more extensive regions of the relevant configuration space. As a result, computational investigations that rely on these techniques are more likely to deliver surprising discoveries.
Collapse
Affiliation(s)
- Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A. Grimmel
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L. Türtscher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P. Unsleber
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Thomas Weymuth
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
13
|
Affiliation(s)
- Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie Vladimir-Prelog-Weg 2 8093 Zürich Switzerland
| |
Collapse
|
14
|
Tsutsumi T, Ono Y, Taketsugu T. Visualization of reaction route map and dynamical trajectory in reduced dimension. Chem Commun (Camb) 2021; 57:11734-11750. [PMID: 34642706 DOI: 10.1039/d1cc04667e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In the quantum chemical approach, chemical reaction mechanisms are investigated based on a potential energy surface (PES). Automated reaction path search methods enable us to construct a global reaction route map containing multiple reaction paths corresponding to a series of elementary reaction processes. The on-the-fly molecular dynamics (MD) method provides a classical trajectory exploring the full-dimensional PES based on electronic structure calculations. We have developed two reaction analysis methods, the on-the-fly trajectory mapping method and the reaction space projector (ReSPer) method, by introducing a structural similarity to a pair of geometric structures and revealed dynamic aspects affecting chemical reaction mechanisms. In this review, we will present the details of these analysis methods and discuss the dynamics effects of reaction path curvature and reaction path bifurcation with applications to the CH3OH + OH- collision reaction and the Au5 cluster branching and isomerization reactions.
Collapse
Affiliation(s)
- Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan.
| | - Yuriko Ono
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan. .,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
| |
Collapse
|
15
|
Lucas K, Chen A, Schubmehl M, Kolonko KJ, Barnes GL. Exploring the Effects of Methylation on the CID of Protonated Lysine: A Combined Experimental and Computational Approach. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2021; 32:2675-2684. [PMID: 34677967 DOI: 10.1021/jasms.1c00225] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We report the results of experiments, simulations, and DFT calculations that focus on describing the reaction dynamics observed within the collision-induced dissociation of l-lysine-H+ and its side-chain methylated analogues, Nε-methyl-l-lysine-H+ (Me1-lysine-H+), Nε,Nε-dimethyl-l-lysine-H+ (Me2-lysine-H+), and Nε,Nε,Nε-trimethyl-l-lysine-H+ (Me3-lysine-H+). The major pathways observed in the experimental measurements were m/z 130 and 84, with the former dominant at low collision energies and the latter at intermediate to high collision energies. The m/z 130 peak corresponds to loss of N(CH3)nH3-n, while m/z 84 has the additional loss of H2CO2 likely in the form of H2O + CO. Within the time frame of the direct dynamics simulations, m/z 130 and 101 were the most populous peaks, with the latter identified as an intermediate to m/z 84. The simulations allowed for the determination of several reaction pathways that result in these products. A graph theory analysis enabled the elucidation of the significant structures that compose each peak. Methylation results in the preferential loss of the side-chain amide group and a reduction of cyclic structures within the m/z 84 peak population in simulations.
Collapse
Affiliation(s)
- Kenneth Lucas
- Department of Chemistry and Biochemistry, Siena College, 515 Loudon Road, Loudonville, New York 12211, United States
| | - Amy Chen
- Department of Chemistry and Biochemistry, Siena College, 515 Loudon Road, Loudonville, New York 12211, United States
| | - Megan Schubmehl
- Department of Chemistry and Biochemistry, Siena College, 515 Loudon Road, Loudonville, New York 12211, United States
| | - Kristopher J Kolonko
- Stewart's Advanced Instrumentation and Technology (SAInT) Center, Siena College, 515 Loudon Road, Loudonville, New York 12211, United States
| | - George L Barnes
- Department of Chemistry and Biochemistry, Siena College, 515 Loudon Road, Loudonville, New York 12211, United States
| |
Collapse
|
16
|
Perez-Mellor AF, Spezia R. Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories. J Chem Phys 2021; 155:124103. [PMID: 34598552 DOI: 10.1063/5.0058382] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.
Collapse
Affiliation(s)
- Ariel F Perez-Mellor
- LAMBE UMR8587, Université d'Evry Val d'Essonne, CNRS, CEA, Université Paris-Saclay, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, 91025 Evry, France
| | - Riccardo Spezia
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
17
|
Martínez-Núñez E, Barnes GL, Glowacki DR, Kopec S, Peláez D, Rodríguez A, Rodríguez-Fernández R, Shannon RJ, Stewart JJP, Tahoces PG, Vazquez SA. AutoMeKin2021: An open-source program for automated reaction discovery. J Comput Chem 2021; 42:2036-2048. [PMID: 34387374 DOI: 10.1002/jcc.26734] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 07/16/2021] [Accepted: 07/27/2021] [Indexed: 01/10/2023]
Abstract
AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin.
Collapse
Affiliation(s)
- Emilio Martínez-Núñez
- Department of Physical Chemistry, University of Santiago de Compostela, Santiago de Compostela, Spain
| | - George L Barnes
- Department of Chemistry and Biochemistry, Siena College, Loudonville, New York, USA
| | - David R Glowacki
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | - Sabine Kopec
- Institut de Sciences Moléculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, Orsay, France
| | - Daniel Peláez
- Institut de Sciences Moléculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, Orsay, France
| | - Aurelio Rodríguez
- Galicia Supercomputing Center (CESGA), Santiago de Compostela, Spain
| | | | - Robin J Shannon
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | | | - Pablo G Tahoces
- Department of Electronics and Computer Science, University of Santiago de Compostela, Santiago de Compostela, Spain
| | - Saulo A Vazquez
- Department of Physical Chemistry, University of Santiago de Compostela, Santiago de Compostela, Spain
| |
Collapse
|
18
|
Borges R, Colby SM, Das S, Edison AS, Fiehn O, Kind T, Lee J, Merrill AT, Merz KM, Metz TO, Nunez JR, Tantillo DJ, Wang LP, Wang S, Renslow RS. Quantum Chemistry Calculations for Metabolomics. Chem Rev 2021; 121:5633-5670. [PMID: 33979149 PMCID: PMC8161423 DOI: 10.1021/acs.chemrev.0c00901] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Indexed: 02/07/2023]
Abstract
A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.
Collapse
Affiliation(s)
- Ricardo
M. Borges
- Walter
Mors Institute of Research on Natural Products, Federal University of Rio de Janeiro, Rio de Janeiro 21941-901, Brazil
| | - Sean M. Colby
- Biological
Science Division, Pacific Northwest National
Laboratory, Richland, Washington 99352, United States
| | - Susanta Das
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Arthur S. Edison
- Departments
of Genetics and Biochemistry and Molecular Biology, Complex Carbohydrate
Research Center and Institute of Bioinformatics, University of Georgia, Athens, Georgia 30602, United States
| | - Oliver Fiehn
- West
Coast Metabolomics Center for Compound Identification, UC Davis Genome
Center, University of California, Davis, California 95616, United States
| | - Tobias Kind
- West
Coast Metabolomics Center for Compound Identification, UC Davis Genome
Center, University of California, Davis, California 95616, United States
| | - Jesi Lee
- West
Coast Metabolomics Center for Compound Identification, UC Davis Genome
Center, University of California, Davis, California 95616, United States
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Amy T. Merrill
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Kenneth M. Merz
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Thomas O. Metz
- Biological
Science Division, Pacific Northwest National
Laboratory, Richland, Washington 99352, United States
| | - Jamie R. Nunez
- Biological
Science Division, Pacific Northwest National
Laboratory, Richland, Washington 99352, United States
| | - Dean J. Tantillo
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Lee-Ping Wang
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Shunyang Wang
- West
Coast Metabolomics Center for Compound Identification, UC Davis Genome
Center, University of California, Davis, California 95616, United States
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Ryan S. Renslow
- Biological
Science Division, Pacific Northwest National
Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
19
|
Krettler CA, Thallinger GG. A map of mass spectrometry-based in silico fragmentation prediction and compound identification in metabolomics. Brief Bioinform 2021; 22:6184408. [PMID: 33758925 DOI: 10.1093/bib/bbab073] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2020] [Revised: 01/29/2021] [Accepted: 02/12/2021] [Indexed: 12/27/2022] Open
Abstract
Metabolomics, the comprehensive study of the metabolome, and lipidomics-the large-scale study of pathways and networks of cellular lipids-are major driving forces in enabling personalized medicine. Complicated and error-prone data analysis still remains a bottleneck, however, especially for identifying novel metabolites. Comparing experimental mass spectra to curated databases containing reference spectra has been the gold standard for identification of compounds, but constructing such databases is a costly and time-demanding task. Many software applications try to circumvent this process by utilizing cutting-edge advances in computational methods-including quantum chemistry and machine learning-and simulate mass spectra by performing theoretical, so called in silico fragmentations of compounds. Other solutions concentrate directly on experimental spectra and try to identify structural properties by investigating reoccurring patterns and the relationships between them. The considerable progress made in the field allows recent approaches to provide valuable clues to expedite annotation of experimental mass spectra. This review sheds light on individual strengths and weaknesses of these tools, and attempts to evaluate them-especially in view of lipidomics, when considering complex mixtures found in biological samples as well as mass spectrometer inter-instrument variability.
Collapse
Affiliation(s)
- Christoph A Krettler
- Institute of Biomedical Informatics, Graz University of Technology, Stremayrgasse 16/I, 8010, Graz, Austria.,Omics Center Graz, BioTechMed-Graz, Stiftingtalstrasse 24, 8010, Graz, Austria
| | - Gerhard G Thallinger
- Institute of Biomedical Informatics, Graz University of Technology, Stremayrgasse 16/I, 8010, Graz, Austria.,Omics Center Graz, BioTechMed-Graz, Stiftingtalstrasse 24, 8010, Graz, Austria
| |
Collapse
|
20
|
Tsutsumi T, Ono Y, Arai Z, Taketsugu T. Visualization of the Dynamics Effect: Projection of on-the-Fly Trajectories to the Subspace Spanned by the Static Reaction Path Network. J Chem Theory Comput 2020; 16:4029-4037. [DOI: 10.1021/acs.jctc.0c00018] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Takuro Tsutsumi
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-0810, Japan
| | - Yuriko Ono
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
| | - Zin Arai
- Academy of Emerging Sciences, Chubu University, Kasugai, Aichi 487-8501, Japan
| | - Tetsuya Taketsugu
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-0810, Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
| |
Collapse
|
21
|
Lucas K, Barnes GL. Modeling the Effects of O-Sulfonation on the CID of Serine. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2020; 31:1114-1122. [PMID: 32202776 DOI: 10.1021/jasms.0c00037] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of O-sulfonation on the collision-induced dissociation for serine. Toward this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m/z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m/z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The m/z 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also observed that the RM1 semiempirical method was not able to obtain all of the major peaks seen in experimens for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.
Collapse
Affiliation(s)
- Kenneth Lucas
- Department of Chemistry and Biochemistry Siena College 515 Loudon Road Loudonville, New York 12211, United States
| | - George L Barnes
- Department of Chemistry and Biochemistry Siena College 515 Loudon Road Loudonville, New York 12211, United States
| |
Collapse
|
22
|
Quapp W, Bofill JM. Some Mathematical Reasoning on the Artificial Force Induced Reaction Method. J Comput Chem 2020; 41:629-634. [PMID: 31792984 DOI: 10.1002/jcc.26115] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2019] [Revised: 11/01/2019] [Accepted: 11/04/2019] [Indexed: 11/07/2022]
Abstract
There are works of the Maeda-Morokuma group, which propose the artificial force induced reaction (AFIR) method (Maeda et al., J. Comput. Chem. 2014, 35, 166 and 2018, 39, 233). We study this important method from a theoretical point of view. The understanding of the proposers does not use the barrier breakdown point of the AFIR parameter, which usually is half of the reaction path between the minimum and the transition state which is searched for. Based on a comparison with the theory of Newton trajectories, we could better understand the method. It allows us to follow along some reaction pathways from minimum to saddle point, or vice versa. We discuss some well-known two-dimensional test surfaces where we calculate full AFIR pathways. If one has special AFIR curves at hand, one can also study the behavior of the ansatz. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Wolfgang Quapp
- Leipzig University, Mathematisches Institut, Universität Leipzig, PF 100920, D-04009, Leipzig, Germany
| | - Josep Maria Bofill
- Universitat de Barcelona, Departament de Química Inorgànica i Orgànica, Universitat de Barcelona, and Institut de Química Teòrica i Computacional, Universitat de Barcelona, (IQTCUB), Martí i Franquès, 1, 08028, Barcelona, Spain
| |
Collapse
|
23
|
Lee K, Woo Kim J, Youn Kim W. Efficient Construction of a Chemical Reaction Network Guided By a Monte Carlo Tree Search. CHEMSYSTEMSCHEM 2020. [DOI: 10.1002/syst.201900057] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Kyunghoon Lee
- Department of ChemistryKorea Advanced Institute of Science and Technology (KAIST) 291 Daehak-ro, Yuseong-gu Daejeon 305-701 Korea
| | - Jin Woo Kim
- Department of ChemistryKorea Advanced Institute of Science and Technology (KAIST) 291 Daehak-ro, Yuseong-gu Daejeon 305-701 Korea
| | - Woo Youn Kim
- Department of ChemistryKorea Advanced Institute of Science and Technology (KAIST) 291 Daehak-ro, Yuseong-gu Daejeon 305-701 Korea
- KI for Artificial IntelligenceKorea Advanced Institute of Science and Technology (KAIST) 291 Daehak-ro, Yuseong-gu Daejeon 305-701 Korea
| |
Collapse
|
24
|
Rodríguez-Fernández A, Bonnet L, Crespos C, Larrégaray P, Díez Muiño R. When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111). Phys Chem Chem Phys 2020; 22:22805-22814. [DOI: 10.1039/d0cp02655g] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.
Collapse
Affiliation(s)
| | | | | | | | - Ricardo Díez Muiño
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU)
- 20018 Donostia-SanSebastián
- Spain
- Donostia International Physics Center (DIPC)
- 20018 Donostia-SanSebastián
| |
Collapse
|
25
|
Rodríguez-Fernández A, Bonnet L, Crespos C, Larrégaray P, Díez Muiño R. When Classical Trajectories Get to Quantum Accuracy: The Scattering of H 2 on Pd(111). J Phys Chem Lett 2019; 10:7629-7635. [PMID: 31774684 DOI: 10.1021/acs.jpclett.9b02742] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
When elementary reactive processes occur at such low energies that only a few states of reactants and/or products are available, quantum effects strongly manifest and the standard description of the dynamics within the classical framework fails. We show here, for H2 scattering on Pd(111), that by pseudoquantizing in the spirit of Bohr the relevant final actions of the system, along with adequately treating the diffraction-mediated trapping of the incoming wave, classical simulations achieve an unprecedented agreement with state-of-the-art quantum dynamics calculations.
Collapse
Affiliation(s)
- A Rodríguez-Fernández
- Université de Bordeaux, ISM , UMR 5255, F-33400 Talence , France
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU) , Paseo Manuel de Lardizabal 5 , 20018 Donostia-San Sebastián , Spain
| | - L Bonnet
- Université de Bordeaux, ISM , UMR 5255, F-33400 Talence , France
- CNRS, ISM , UMR 5255, F-33400 Talence , France
| | - C Crespos
- Université de Bordeaux, ISM , UMR 5255, F-33400 Talence , France
- CNRS, ISM , UMR 5255, F-33400 Talence , France
| | - P Larrégaray
- Université de Bordeaux, ISM , UMR 5255, F-33400 Talence , France
- CNRS, ISM , UMR 5255, F-33400 Talence , France
| | - R Díez Muiño
- Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU) , Paseo Manuel de Lardizabal 5 , 20018 Donostia-San Sebastián , Spain
- Donostia International Physics Center (DIPC) , Paseo Manuel de Lardizabal 4 , 20018 Donostia-San Sebastián , Spain
| |
Collapse
|
26
|
De DS, Krummenacher M, Schaefer B, Goedecker S. Finding Reaction Pathways with Optimal Atomic Index Mappings. PHYSICAL REVIEW LETTERS 2019; 123:206102. [PMID: 31809087 DOI: 10.1103/physrevlett.123.206102] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/10/2019] [Indexed: 06/10/2023]
Abstract
Finding complex reaction and transformation pathways involving many intermediate states is, in general, not possible on the density-functional theory level with existing simulation methods, due to the very large number of required energy and force evaluations. For complex reactions, it is not possible to determine which atom in the reactant is mapped onto which atom in the product. Trying out all possible atomic index mappings is not feasible because of the factorial increase in the number of possible mappings. We use a penalty function that is invariant under index permutations to bias the potential energy surface in such a way that it obtains the characteristics of a structure seeker, whose global minimum is the reaction product. By performing a minima-hopping-based global optimization on this biased potential energy surface, we rapidly find intermediate states that lead into the global minimum and allow us to then extract entire reaction pathways. We first demonstrate for a benchmark system, namely, the Lennard-Jones cluster LJ_{38}, that our method finds intermediate states relevant to the lowest energy reaction pathway, and hence we need to consider much fewer intermediate states than previous methods to find the lowest energy reaction pathway. Finally, we apply the method to two real systems, C_{60} and C_{20}H_{20}, and show that the reaction pathways found contain valuable information on how these molecules can be synthesized.
Collapse
Affiliation(s)
- Deb Sankar De
- Department of Physics, Universität Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Marco Krummenacher
- Department of Physics, Universität Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Bastian Schaefer
- Department of Physics, Universität Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Stefan Goedecker
- Department of Physics, Universität Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| |
Collapse
|
27
|
Jara‐Toro RA, Pino GA, Glowacki DR, Shannon RJ, Martínez‐Núñez E. Enhancing Automated Reaction Discovery with Boxed Molecular Dynamics in Energy Space. CHEMSYSTEMSCHEM 2019. [DOI: 10.1002/syst.201900024] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Rafael A. Jara‐Toro
- INIFIQC (CONICET-UNC) Dpto. De Fisicoquímica-Facultad de Ciencias Químicas-Centro Láser de Ciencias MolecularesUniversidad de Córdoba Ciudad Universitaria X50000HUA Córdoba Argentina
| | - Gustavo A. Pino
- INIFIQC (CONICET-UNC) Dpto. De Fisicoquímica-Facultad de Ciencias Químicas-Centro Láser de Ciencias MolecularesUniversidad de Córdoba Ciudad Universitaria X50000HUA Córdoba Argentina
| | - David R. Glowacki
- Centre for Computational Chemistry School of ChemistryUniversity of Bristol Cantock's Close Bristol BS8 1TS UK
| | - Robin J. Shannon
- Centre for Computational Chemistry School of ChemistryUniversity of Bristol Cantock's Close Bristol BS8 1TS UK
| | - Emilio Martínez‐Núñez
- Departmento de Química Física, Facultade de QuímicaUniversidade de Santiago de Compostela 15782 Santiago de Compostela Spain
| |
Collapse
|
28
|
Bonnet L, Larrégaray P, Lara M, Launay JM. Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S( 1D) + X 2, X = H, D. J Phys Chem A 2019; 123:6439-6454. [PMID: 31329443 DOI: 10.1021/acs.jpca.9b04938] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.
Collapse
Affiliation(s)
- L Bonnet
- Université de Bordeaux, ISM , UMR 5255, F-33400 Talence , France.,CNRS , ISM , UMR 5255, F-33400 Talence , France
| | - P Larrégaray
- Université de Bordeaux, ISM , UMR 5255, F-33400 Talence , France.,CNRS , ISM , UMR 5255, F-33400 Talence , France
| | - M Lara
- Departamento de Química Física Aplicada, Facultad de Ciencias , Universidad Autónoma de Madrid , 28049 Madrid , Spain
| | - J-M Launay
- Institut de Physique de Rennes, UMR CNRS 6251 , Université de Rennes I , F-35042 Rennes , France
| |
Collapse
|
29
|
Kim JW, Kim Y, Baek KY, Lee K, Kim WY. Performance of ACE-Reaction on 26 Organic Reactions for Fully Automated Reaction Network Construction and Microkinetic Analysis. J Phys Chem A 2019; 123:4796-4805. [DOI: 10.1021/acs.jpca.9b02161] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jin Woo Kim
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Yeonjoon Kim
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Kyung Yup Baek
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Kyunghoon Lee
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Woo Youn Kim
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
- KAIST Institute for Artificial Intelligence, KAIST 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| |
Collapse
|
30
|
Aguirre NF, Díaz-Tendero S, IdBarkach T, Chabot M, Béroff K, Alcamí M, Martín F. Fully versus constrained statistical fragmentation of carbon clusters and their heteronuclear derivatives. J Chem Phys 2019; 150:144301. [PMID: 30981259 DOI: 10.1063/1.5083864] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Microcanonical Metropolis Monte Carlo (MMMC) method has been shown to describe reasonably well fragmentation of clusters composed of identical atomic species. However, this is not so clear in the case of heteronuclear clusters as some regions of phase space might be inaccessible due to the different mobility of the different atomic species, the existence of large isomerization barriers, or the quite different chemical nature of the possible intermediate species. In this paper, we introduce a constrained statistical model that extends the range of applicability of the MMMC method to such mixed clusters. The method is applied to describe fragmentation of isolated clusters with high, moderate, and no heteronuclear character, namely, CnHm, CnN, and Cn clusters for which experimental fragmentation branching ratios are available in the literature. We show that the constrained statistical model describes fairly well fragmentation of CnHm clusters in contrast with the poor description provided by the fully statistical model. The latter model, however, works pretty well for both Cn and CnN clusters, thus showing that the ultimate reason for this discrepancy is the inability of the MMMC method to selectively explore the whole phase space. This conclusion has driven us to predict the fragmentation patterns of the C4N cluster for which experiments are not yet available.
Collapse
Affiliation(s)
- Néstor F Aguirre
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Sergio Díaz-Tendero
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Tijani IdBarkach
- Institut de Physique Nucléaire d'Orsay, IN2P3-CNRS and Université Paris-Sud, F-91406 Orsay Cedex, France
| | - Marin Chabot
- Institut de Physique Nucléaire d'Orsay, IN2P3-CNRS and Université Paris-Sud, F-91406 Orsay Cedex, France
| | - Karine Béroff
- Institut des Sciences Moléculaires d'Orsay, INP-CNRS and Université Paris-Sud, F-91405 Orsay Cedex, France
| | - Manuel Alcamí
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Fernando Martín
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|