1
|
Martí C, Devereux C, Najm HN, Zádor J. Evaluation of Rate Coefficients in the Gas Phase Using Machine-Learned Potentials. J Phys Chem A 2024. [PMID: 38427974 DOI: 10.1021/acs.jpca.3c07872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2024]
Abstract
We assess the capability of machine-learned potentials to compute rate coefficients by training a neural network (NN) model and applying it to describe the chemical landscape on the C5H5 potential energy surface, which is relevant to molecular weight growth in combustion and interstellar media. We coupled the resulting NN with an automated kinetics workflow code, KinBot, to perform all necessary calculations to compute the rate coefficients. The NN is benchmarked exhaustively by evaluating its performance at the various stages of the kinetics calculations: from the electronic energy through the computation of zero point energy, barrier heights, entropic contributions, the portion of the PES explored, and finally the overall rate coefficients as formulated by transition state theory.
Collapse
Affiliation(s)
- Carles Martí
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| | - Christian Devereux
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| | - Habib N Najm
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| | - Judit Zádor
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| |
Collapse
|
2
|
Zádor J, Martí C, Van de Vijver R, Johansen SL, Yang Y, Michelsen HA, Najm HN. Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces. J Phys Chem A 2023; 127:565-588. [PMID: 36607817 DOI: 10.1021/acs.jpca.2c06558] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. We bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.
Collapse
Affiliation(s)
- Judit Zádor
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | - Carles Martí
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | | | - Sommer L Johansen
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | - Yoona Yang
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| | - Hope A Michelsen
- Department of Mechanical Engineering, University of Colorado Boulder, Boulder80309, Colorado, United States
| | - Habib N Najm
- Combustion Research Facility, Sandia National Laboratories, Livermore94550, California, United States
| |
Collapse
|
3
|
Kuntz D, Wilson AK. Machine learning, artificial intelligence, and chemistry: how smart algorithms are reshaping simulation and the laboratory. PURE APPL CHEM 2022. [DOI: 10.1515/pac-2022-0202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
Machine learning and artificial intelligence are increasingly gaining in prominence through image analysis, language processing, and automation, to name a few applications. Machine learning is also making profound changes in chemistry. From revisiting decades-old analytical techniques for the purpose of creating better calibration curves, to assisting and accelerating traditional in silico simulations, to automating entire scientific workflows, to being used as an approach to deduce underlying physics of unexplained chemical phenomena, machine learning and artificial intelligence are reshaping chemistry, accelerating scientific discovery, and yielding new insights. This review provides an overview of machine learning and artificial intelligence from a chemist’s perspective and focuses on a number of examples of the use of these approaches in computational chemistry and in the laboratory.
Collapse
Affiliation(s)
- David Kuntz
- Department of Chemistry , University of North Texas , Denton , TX 76201 , USA
| | - Angela K. Wilson
- Department of Chemistry , Michigan State University , East Lansing , MI 48824 , USA
| |
Collapse
|
4
|
Moberg DR, Jasper AW. Permutationally Invariant Polynomial Expansions with Unrestricted Complexity. J Chem Theory Comput 2021; 17:5440-5455. [PMID: 34469127 DOI: 10.1021/acs.jctc.1c00352] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. The parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.
Collapse
Affiliation(s)
- Daniel R Moberg
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Ahren W Jasper
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
5
|
Lu J, Zhang H, Yu J, Shan D, Qi J, Chen J, Song H, Yang M. Predicting Rate Constants of Hydroxyl Radical Reactions with Alkanes Using Machine Learning. J Chem Inf Model 2021; 61:4259-4265. [PMID: 34378385 DOI: 10.1021/acs.jcim.1c00809] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The hydrogen abstraction reactions of the hydroxyl radical with alkanes play an important role in combustion chemistry and atmospheric chemistry. However, site-specific reaction constants are difficult to obtain experimentally and theoretically. Recently, machine learning has proved its ability to predict chemical properties. In this work, a machine learning approach is developed to predict the temperature-dependent site-specific rate constants of the title reactions. Multilayered neural network (NN) models are developed by training the site-specific rate constants of 11 reactions, and several schemes are designed to improve the prediction accuracy. The results show that the proposed NN models are robust in predicting the site-specific and overall rate constants.
Collapse
Affiliation(s)
- Junhui Lu
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Huimin Zhang
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Jinhui Yu
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Dezun Shan
- School of Physical Science and Technology, Wuhan University, Wuhan 430072, China
| | - Ji Qi
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Jiawen Chen
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Hongwei Song
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Minghui Yang
- Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.,University of Chinese Academy of Sciences, Beijing 100049, China.,Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan 430071, China
| |
Collapse
|
6
|
Ferguson AL, Hachmann J, Miller TF, Pfaendtner J. The Journal of Physical Chemistry A/ B/ C Virtual Special Issue on Machine Learning in Physical Chemistry. J Phys Chem A 2021; 124:9113-9118. [PMID: 33147969 DOI: 10.1021/acs.jpca.0c09205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
7
|
Ferguson AL, Hachmann J, Miller TF, Pfaendtner J. The Journal of Physical Chemistry A/ B/ C Virtual Special Issue on Machine Learning in Physical Chemistry. J Phys Chem B 2021; 124:9767-9772. [PMID: 33147970 DOI: 10.1021/acs.jpcb.0c09206] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
8
|
Manzhos S, Carrington T. Neural Network Potential Energy Surfaces for Small Molecules and Reactions. Chem Rev 2020; 121:10187-10217. [PMID: 33021368 DOI: 10.1021/acs.chemrev.0c00665] [Citation(s) in RCA: 110] [Impact Index Per Article: 27.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We review progress in neural network (NN)-based methods for the construction of interatomic potentials from discrete samples (such as ab initio energies) for applications in classical and quantum dynamics including reaction dynamics and computational spectroscopy. The main focus is on methods for building molecular potential energy surfaces (PES) in internal coordinates that explicitly include all many-body contributions, even though some of the methods we review limit the degree of coupling, due either to a desire to limit computational cost or to limited data. Explicit and direct treatment of all many-body contributions is only practical for sufficiently small molecules, which are therefore our primary focus. This includes small molecules on surfaces. We consider direct, single NN PES fitting as well as more complex methods that impose structure (such as a multibody representation) on the PES function, either through the architecture of one NN or by using multiple NNs. We show how NNs are effective in building representations with low-dimensional functions including dimensionality reduction. We consider NN-based approaches to build PESs in the sums-of-product form important for quantum dynamics, ways to treat symmetry, and issues related to sampling data distributions and the relation between PES errors and errors in observables. We highlight combinations of NNs with other ideas such as permutationally invariant polynomials or sums of environment-dependent atomic contributions, which have recently emerged as powerful tools for building highly accurate PESs for relatively large molecular and reactive systems.
Collapse
Affiliation(s)
- Sergei Manzhos
- Centre Énergie Matériaux Télécommunications, Institut National de la Recherche Scientifique, 1650, Boulevard Lionel-Boulet, Varennes, Québec City, Québec J3X 1S2, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston Ontario K7L 3N6, Canada
| |
Collapse
|
9
|
Quiton SJ, Mitra U, Mallikarjun Sharada S. A matrix completion algorithm to recover modes orthogonal to the minimum energy path in chemical reactions. J Chem Phys 2020; 153:054122. [PMID: 32770897 DOI: 10.1063/5.0018326] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Structured statistical methods are promising for recovering or completing information from noisy and incomplete data with high fidelity. In particular, matrix completion exploits underlying structural properties such as rank or sparsity. Our objective is to employ matrix completion to reduce computational effort associated with the calculation of multiple quantum chemical Hessians, which are necessary for identification of temperature-dependent free energy maxima under canonical variational transition state theory (VTST). We demonstrate proof-of-principle of an algebraic variety-based matrix completion method for recovering missing elements in a matrix of transverse Hessian eigenvalues constituting the minimum energy path (MEP) of a reaction. The algorithm, named harmonic variety-based matrix completion (HVMC), utilizes the fact that the points lying on the MEP of a reaction step constitute an algebraic variety in the reaction path Hamiltonian representation. We demonstrate that, with as low as 30% random sampling of matrix elements for the largest system in our test set (46 atoms), the complete matrix of eigenvalues can be recovered. We further establish algorithm performance for VTST rate calculations by quantifying errors in zero-point energies and vibrational free energies. Motivated by this success, we outline next steps toward developing a practical HVMC algorithm, which utilizes a gradient-based sampling protocol for low-cost VTST rate computations.
Collapse
Affiliation(s)
- Stephen Jon Quiton
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, USA
| | - Urbashi Mitra
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|