1
|
Hwang W, Austin SL, Blondel A, Boittier ED, Boresch S, Buck M, Buckner J, Caflisch A, Chang HT, Cheng X, Choi YK, Chu JW, Crowley MF, Cui Q, Damjanovic A, Deng Y, Devereux M, Ding X, Feig MF, Gao J, Glowacki DR, Gonzales JE, Hamaneh MB, Harder ED, Hayes RL, Huang J, Huang Y, Hudson PS, Im W, Islam SM, Jiang W, Jones MR, Käser S, Kearns FL, Kern NR, Klauda JB, Lazaridis T, Lee J, Lemkul JA, Liu X, Luo Y, MacKerell AD, Major DT, Meuwly M, Nam K, Nilsson L, Ovchinnikov V, Paci E, Park S, Pastor RW, Pittman AR, Post CB, Prasad S, Pu J, Qi Y, Rathinavelan T, Roe DR, Roux B, Rowley CN, Shen J, Simmonett AC, Sodt AJ, Töpfer K, Upadhyay M, van der Vaart A, Vazquez-Salazar LI, Venable RM, Warrensford LC, Woodcock HL, Wu Y, Brooks CL, Brooks BR, Karplus M. CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed. J Phys Chem B 2024; 128:9976-10042. [PMID: 39303207 DOI: 10.1021/acs.jpcb.4c04100] [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/22/2024]
Abstract
Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.
Collapse
Affiliation(s)
- Wonmuk Hwang
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas 77843, United States
- Department of Materials Science and Engineering, Texas A&M University, College Station, Texas 77843, United States
- Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843, United States
- Center for AI and Natural Sciences, Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Steven L Austin
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Arnaud Blondel
- Institut Pasteur, Université Paris Cité, CNRS UMR3825, Structural Bioinformatics Unit, 28 rue du Dr. Roux F-75015 Paris, France
| | - Eric D Boittier
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Wahringerstrasse 17, 1090 Vienna, Austria
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, Cleveland, Ohio 44106, United States
| | - Joshua Buckner
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zürich, CH-8057 Zürich, Switzerland
| | - Hao-Ting Chang
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, ROC
| | - Xi Cheng
- Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Yeol Kyo Choi
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jhih-Wei Chu
- Institute of Bioinformatics and Systems Biology, Department of Biological Science and Technology, Institute of Molecular Medicine and Bioengineering, and Center for Intelligent Drug Systems and Smart Bio-devices (IDS2B), National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, ROC
| | - Michael F Crowley
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| | - Ana Damjanovic
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Yuqing Deng
- Shanghai R&D Center, DP Technology, Ltd., Shanghai 201210, China
| | - Mike Devereux
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Xinqiang Ding
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Michael F Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Jiali Gao
- School of Chemical Biology & Biotechnology, Peking University Shenzhen Graduate School, Shenzhen, Guangdong 518055, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 518055, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - David R Glowacki
- CiTIUS Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, 15705 Santiago de Compostela, Spain
| | - James E Gonzales
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas 77843, United States
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Mehdi Bagerhi Hamaneh
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, Cleveland, Ohio 44106, United States
| | | | - Ryan L Hayes
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, Irvine, California 92697, United States
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, United States
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
| | - Yandong Huang
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Phillip S Hudson
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
- Medicine Design, Pfizer Inc., Cambridge, Massachusetts 02139, United States
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Shahidul M Islam
- Department of Chemistry, Delaware State University, Dover, Delaware 19901, United States
| | - Wei Jiang
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Michael R Jones
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Fiona L Kearns
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Nathan R Kern
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, Institute for Physical Science and Technology, Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Themis Lazaridis
- Department of Chemistry, City College of New York, New York, New York 10031, United States
| | - Jinhyuk Lee
- Disease Target Structure Research Center, Korea Research Institute of Bioscience and Biotechnology, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology, Daejeon 34141, Republic of Korea
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| | - Xiaorong Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yun Luo
- Department of Biotechnology and Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California 91766, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Dan T Major
- Department of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden
| | - Victor Ovchinnikov
- Harvard University, Department of Chemistry and Chemical Biology, Cambridge, Massachusetts 02138, United States
| | - Emanuele Paci
- Dipartimento di Fisica e Astronomia, Universitá di Bologna, Bologna 40127, Italy
| | - Soohyung Park
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Amanda R Pittman
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Carol Beth Post
- Borch Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, Indiana 47907, United States
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Yifei Qi
- School of Pharmacy, Fudan University, Shanghai 201203, China
| | | | - Daniel R Roe
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Benoit Roux
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | | | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Alexander J Sodt
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Arjan van der Vaart
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | | | - Richard M Venable
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Luke C Warrensford
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Yujin Wu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Charles L Brooks
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Martin Karplus
- Harvard University, Department of Chemistry and Chemical Biology, Cambridge, Massachusetts 02138, United States
- Laboratoire de Chimie Biophysique, ISIS, Université de Strasbourg, 67000 Strasbourg, France
| |
Collapse
|
2
|
Shannon R, Blitz MA, Seakins PW. Solving the OH + Glyoxal Problem: A Complete Theoretical Description of Post-Transition-State Energy Deposition in Activated Systems. J Phys Chem A 2024; 128:1501-1510. [PMID: 38377581 PMCID: PMC10910583 DOI: 10.1021/acs.jpca.3c07823] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 01/25/2024] [Accepted: 01/31/2024] [Indexed: 02/22/2024]
Abstract
Activated chemistry in coupled reaction systems has broadened our understanding of the chemical kinetics. In the case of intermediates formed in gas phase abstraction reactions (e.g., OH + HC(O)C(O)H (glyoxal) →HC(O)CO + H2O), it is particularly crucial to understand how the reaction energy is partitioned between product species as this determines the propensity for a given product to undergo "prompt" dissociation (e.g., HC(O)CO → HCO + CO) before the excess reaction energy is removed. An example of such an activated system is the OH + glyoxal + O2 coupled reaction system. In this work, we develop a molecular dynamics pipeline, which, combined with a master equation analysis, accurately models previous experimental measurements. This new work resolves previous complexities and discrepancies from earlier master equation modeling for this reaction system. The detailed molecular dynamics approach employed here is a powerful new tool for modeling challenging activated reaction systems.
Collapse
Affiliation(s)
- Robin Shannon
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, U.K.
| | - Mark A. Blitz
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, U.K.
- National
Centre for Atmospheric Science, University
of Leeds, Leeds LS2 9JT, U.K.
| | - Paul W. Seakins
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, U.K.
| |
Collapse
|
3
|
Meuwly M. Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─ Quo Vadis?. J Phys Chem B 2022; 126:2155-2167. [PMID: 35286087 DOI: 10.1021/acs.jpcb.2c00212] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
4
|
Zhang X, Lefebvre PL, Harvey JN. Effect of solvent motions on the dynamics of the Diels-Alder reaction. Phys Chem Chem Phys 2021; 24:1120-1130. [PMID: 34928279 DOI: 10.1039/d1cp05272a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
How solvent motions affect the dynamics of chemical reactions in which the solute undergoes a substantial shape change is a fundamental but elusive issue. This work utilizes reactive simulation and Grote-Hynes theory to explore the effect of solvent motions on the dynamics of the Diels-Alder reaction (in the reverse direction, this reaction involves very substantial solute expansion) in aprotic solvents. The results reveal that the solvent environment is not sufficiently constraining to influence transition state passage dynamics, with the calculated transmission coefficients being close to unity. Even when solvent motions are suppressed or artificially slowed down, the solvent only affects the reaction dynamics in the transition state region to a very small extent. The only notable effect of solvent occurs far from the transition state region and corresponds to caging of the reactants within the reactant well.
Collapse
Affiliation(s)
- Xiaoyong Zhang
- Theoretical and Computational Chemistry, Department of Chemistry, KU Leuven Celestijnenlaan 200F, 3001, Leuven, Belgium.
| | - Pierre-Louis Lefebvre
- Theoretical and Computational Chemistry, Department of Chemistry, KU Leuven Celestijnenlaan 200F, 3001, Leuven, Belgium. .,Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| | - Jeremy N Harvey
- Theoretical and Computational Chemistry, Department of Chemistry, KU Leuven Celestijnenlaan 200F, 3001, Leuven, Belgium.
| |
Collapse
|
5
|
Zhang X, Vázquez SA, Harvey JN. Vibrational Energy Relaxation of Deuterium Fluoride in d-Dichloromethane: Insights from Different Potentials. J Chem Theory Comput 2021; 17:1277-1289. [PMID: 33550803 DOI: 10.1021/acs.jctc.0c01059] [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/28/2022]
Abstract
Vibrationally excited deuterium fluoride (DF) formed by fluorine atom reaction with a solvent was found (Science, 2015, 347, 530) to relax rapidly (less than 10 ps) in acetonitrile-d3 (CD3CN) and dichloromethane-d2 (CD2Cl2). However, insights into how CD2Cl2 facilitates this energy relaxation have so far been lacking, given the weak interaction between DF and a single CD2Cl2. In this work, we report the results of reactive simulations with a two-state reactive empirical valence bond (EVB) potential to study the energy deposited into nascent DF after transition-state passage and of nonequilibrium molecular dynamics simulations using multiple different potential energy functions to model the relaxation dynamics. For these second simulations, we used the standard Merck molecular force field (MMFF) potential, an MMFF-based covalent-ionic empirical valence bond (EVB) potential (EVBCI), a newly developed potential [referred to as MMFF(rDF)] which extends upon the MMFF potential by making the DF/CD2Cl2 interaction depend on the value of the D-F bond stretching coordinate and by taking the anisotropic charge distribution of the solvent molecules into account, the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) potential, and the quantum mechanics/molecular mechanics (QM/MM) potential. The relaxation is revealed to be highly sensitive to the potential used. Neither standard MMFF nor EVBCI reproduces the experimentally observed rapid relaxation dynamics, and they also fail to provide a good description of the interaction potential between DF and CD2Cl2 as calculated using CCSD(T)-F12. This is attributed to the use of a point-charge model for the solute and to failing to model the anisotropic electrostatic properties of CD2Cl2. The MMFF(rDF), AMOEBA, and QM/MM potentials all reproduce the CCSD(T)-F12 two-body DF---CD2Cl2 interaction potential rather well but only with the QM/MM approach is fast vibrational relaxation obtained (lifetimes of ∼288, ∼186, and ∼8 ps, respectively), which we attribute to differences in the solute-solvent local structure. With QM/MM, a unique "many-body" interaction pattern in which DF is in close contact with two solvent Cl atoms and more than three solvent D atoms is found, but this structure is not seen with other potentials. The QM/MM dynamics also display enhanced solute-solvent interactions with vibrationally excited DF that induce a DF band redshift and hence a resonant overlap with solvent C-D modes, which facilitate the intermolecular energy transfer. Our work also suggests that potentials used to model energy relaxation need to capture the fine structure of solute-solvent interactions and not just the two-body part.
Collapse
Affiliation(s)
- Xiaoyong Zhang
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven B-3001, Belgium
| | - Saulo A Vázquez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Jeremy N Harvey
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven B-3001, Belgium
| |
Collapse
|
6
|
Amabilino S, Bratholm LA, Bennie SJ, O’Connor MB, Glowacki DR. Training atomic neural networks using fragment-based data generated in virtual reality. J Chem Phys 2020; 153:154105. [DOI: 10.1063/5.0015950] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Affiliation(s)
- Silvia Amabilino
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Intangible Realities Laboratory, University of Bristol, Bristol BS8 1UB, United Kingdom
| | - Lars A. Bratholm
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Intangible Realities Laboratory, University of Bristol, Bristol BS8 1UB, United Kingdom
| | - Simon J. Bennie
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Intangible Realities Laboratory, University of Bristol, Bristol BS8 1UB, United Kingdom
| | - Michael B. O’Connor
- Intangible Realities Laboratory, University of Bristol, Bristol BS8 1UB, United Kingdom
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
| | - David R. Glowacki
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Intangible Realities Laboratory, University of Bristol, Bristol BS8 1UB, United Kingdom
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
| |
Collapse
|
7
|
Scivetti I, Sen K, Elena AM, Todorov I. Reactive Molecular Dynamics at Constant Pressure via Nonreactive Force Fields: Extending the Empirical Valence Bond Method to the Isothermal-Isobaric Ensemble. J Phys Chem A 2020; 124:7585-7597. [PMID: 32820921 DOI: 10.1021/acs.jpca.0c05461] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Empirical Valence Bond (EVB) method offers a suitable framework to obtain reactive potentials through the coupling of nonreactive force fields. In this formalism, most of the implemented coupling terms are built using functional forms that depend on spatial coordinates, while parameters are fitted against reference data to model the change of chemistry between the participating nonreactive states. In this work, we demonstrate that the use of such coupling terms precludes the computation of the stress tensor for condensed phase systems and prevents the possibility to carry out EVB molecular dynamics in the isothermal-isobaric (NPT) ensemble. Alternatively, we make use of coupling terms that depend on the energy gaps, defined as the energy differences between the participating nonreactive force fields, and derive a general expression for the EVB stress tensor suitable for computation. Implementation of this new methodology is tested for a model of a single reactive malonaldehyde solvated in nonreactive water. Mass densities and probability distributions for the values of the energy gaps computed in the NPT ensemble reveal a negligible role of the reactive potential in the limit of low concentrated solutions, thus corroborating for the first time the validity of approximations based on the canonical NVT ensemble, customarily adopted for EVB simulations. The presented formalism also aims to contribute to future implementations and extensions of the EVB method to research the limit of highly concentrated solutions.
Collapse
Affiliation(s)
- Ivan Scivetti
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K.,Department of Chemistry, University of Liverpool, Liverpool L69 3BX, U.K
| | - Kakali Sen
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Alin M Elena
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Ilian Todorov
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| |
Collapse
|
8
|
Koner D, Salehi SM, Mondal P, Meuwly M. Non-conventional force fields for applications in spectroscopy and chemical
reaction dynamics. J Chem Phys 2020; 153:010901. [DOI: 10.1063/5.0009628] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Seyedeh Maryam Salehi
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Padmabati Mondal
- Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra
Pradesh, India
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel,
Switzerland and Department of Chemistry, Brown University, Providence, Rhode Island, USA
| |
Collapse
|
9
|
Wang CH, Masunov AE, Allison TC, Chang S, Lim C, Jin Y, Vasu SS. Molecular Dynamics of Combustion Reactions in Supercritical Carbon Dioxide. 6. Computational Kinetics of Reactions between Hydrogen Atom and Oxygen Molecule H + O 2 ⇌ HO + O and H + O 2 ⇌ HO 2. J Phys Chem A 2019; 123:10772-10781. [PMID: 31820644 DOI: 10.1021/acs.jpca.9b08789] [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/30/2022]
Abstract
Reactions of the hydrogen atom and the oxygen molecule are among the most important ones in the hydrogen and hydrocarbon oxidation mechanisms, including combustion in a supercritical CO2 (sCO2) environment, known as oxy-combustion or the Allam cycle. Development of these energy technologies requires understanding of chemical kinetics of H + O2 ⇌ HO + O and H + O2 ⇌ HO2 in high pressures and concentrations of CO2. Here, we combine quantum treatment of the reaction system by the transition state theory with classical molecular dynamics simulation and the multistate empirical valence bonding method to treat environmental effects. Potential of mean force in the sCO2 solvent at various temperatures 1000-2000 K and pressures 100-400 atm was obtained. The reaction rate for H + O2 ⇌ HO + O was found to be pressure-independent and described by the extended Arrhenius equation 4.23 × 10-7 T-0.73 exp(-21 855.2 cal/mol/RT) cm3/molecule/s, while the reaction rate H + O2 ⇌ HO2 is pressure-dependent and can be expressed as 5.22 × 10-2 T-2.86 exp(-7247.4 cal/mol/RT) cm3/molecule/s at 300 atm.
Collapse
Affiliation(s)
- Chun-Hung Wang
- NanoScience Technology Center , University of Central Florida , 12424 Research Parkway , Orlando , Florida 32826 , United States
| | - Artëm E Masunov
- NanoScience Technology Center , University of Central Florida , 12424 Research Parkway , Orlando , Florida 32826 , United States.,School of Modeling, Simulation, and Training , University of Central Florida , 3100 Technology Parkway , Orlando , Florida 32816 , United States.,Department of Chemistry , University of Central Florida , 4111 Libra Drive , Orlando , Florida 32816 , United States.,South Ural State University , Lenin pr. 76 , Chelyabinsk 454080 , Russia.,National Research Nuclear University MEPhI , Kashirskoye shosse 31 , Moscow 115409 , Russia
| | - Timothy C Allison
- Southwest Research Institute , San Antonio , Texas 78238 , United States
| | - Sungho Chang
- KEPCO Research Institute , Daejeon 34050 , Korea
| | - Chansun Lim
- Hanwha Power Systems , Seongnam , Gyeonggi 13488 , Korea
| | - Yuin Jin
- Hanwha Power Systems , Seongnam , Gyeonggi 13488 , Korea
| | - Subith S Vasu
- Center for Advanced Turbomachinery and Energy Research (CATER), Mechanical and Aerospace Engineering , University of Central Florida , Orlando , Florida 32816 , United States
| |
Collapse
|
10
|
Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk. J Phys Chem B 2019; 123:9846-9861. [PMID: 31647873 DOI: 10.1021/acs.jpcb.9b03258] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H2O)100H+, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.
Collapse
|
11
|
O'Connor MB, Bennie SJ, Deeks HM, Jamieson-Binnie A, Jones AJ, Shannon RJ, Walters R, Mitchell TJ, Mulholland AJ, Glowacki DR. Interactive molecular dynamics in virtual reality from quantum chemistry to drug binding: An open-source multi-person framework. J Chem Phys 2019; 150:220901. [PMID: 31202243 DOI: 10.1063/1.5092590] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
As molecular scientists have made progress in their ability to engineer nanoscale molecular structure, we face new challenges in our ability to engineer molecular dynamics (MD) and flexibility. Dynamics at the molecular scale differs from the familiar mechanics of everyday objects because it involves a complicated, highly correlated, and three-dimensional many-body dynamical choreography which is often nonintuitive even for highly trained researchers. We recently described how interactive molecular dynamics in virtual reality (iMD-VR) can help to meet this challenge, enabling researchers to manipulate real-time MD simulations of flexible structures in 3D. In this article, we outline various efforts to extend immersive technologies to the molecular sciences, and we introduce "Narupa," a flexible, open-source, multiperson iMD-VR software framework which enables groups of researchers to simultaneously cohabit real-time simulation environments to interactively visualize and manipulate the dynamics of molecular structures with atomic-level precision. We outline several application domains where iMD-VR is facilitating research, communication, and creative approaches within the molecular sciences, including training machines to learn potential energy functions, biomolecular conformational sampling, protein-ligand binding, reaction discovery using "on-the-fly" quantum chemistry, and transport dynamics in materials. We touch on iMD-VR's various cognitive and perceptual affordances and outline how these provide research insight for molecular systems. By synergistically combining human spatial reasoning and design insight with computational automation, technologies such as iMD-VR have the potential to improve our ability to understand, engineer, and communicate microscopic dynamical behavior, offering the potential to usher in a new paradigm for engineering molecules and nano-architectures.
Collapse
Affiliation(s)
- Michael B O'Connor
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Simon J Bennie
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Helen M Deeks
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Alexander Jamieson-Binnie
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Alex J Jones
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Robin J Shannon
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Rebecca Walters
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Thomas J Mitchell
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - David R Glowacki
- Intangible Realities Laboratory, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
12
|
Wang CH, Panteleev SV, Masunov AE, Allison TC, Chang S, Lim C, Jin Y, Vasu SS. Molecular Dynamics of Combustion Reactions in Supercritical Carbon Dioxide. Part 5: Computational Study of Ethane Dissociation and Recombination Reactions C 2H 6 ⇌ CH 3 + CH 3. J Phys Chem A 2019; 123:4776-4784. [PMID: 31034229 DOI: 10.1021/acs.jpca.9b02302] [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/30/2022]
Abstract
Fossil fuel oxy-combustion is an emerging technology where the habitual nitrogen diluent is replaced by high-pressure supercritical CO2 (sCO2), which increases the efficiency of energy conversion. In this study, the chemical kinetics of the combustion reaction C2H6 ⇌ CH3 + CH3 in the sCO2 environment is predicted at 30-1000 atm and 1000-2000 K. We adopt a multiscale approach, where the reactive complex is treated quantum mechanically in rigid rotor/harmonic oscillator approximation, while environment effects at different densities are taken into account by the potential of mean force, produced with classical molecular dynamics (MD). Here, we used boxed MD, where enhanced sampling of infrequent events of barrier crossing is accomplished without application of the bias potential. The multistate empirical valence bond model is applied to describe free radical formation accurately at the cost of the classical force field. Predicted rates at low densities agree well with the literature data. Rate constants at 300 atm are 2.41 × 1014 T-0.20 exp(-77.03 kcal/mol/ RT) 1/s for ethane dissociation and 8.44 × 10-19 T1.42 exp(19.89 kcal/mol/ RT) cm3/molecule/s for methyl-methyl recombination.
Collapse
Affiliation(s)
- Chun-Hung Wang
- NanoScienece Technology Center , University of Central Florida , 12424 Research Parkway , Orlando , Florida 32826 , United States
| | - Sergey V Panteleev
- N. I. Lobachevsky State University of Nizhny Novgorod , Gagarin Av. 23 , Nizhny Novgorod 603950 , Russia
| | - Artëm E Masunov
- NanoScienece Technology Center , University of Central Florida , 12424 Research Parkway , Orlando , Florida 32826 , United States.,Department of Chemistry , University of Central Florida , 4111 Libra Dr. , Orlando , Florida 32816 , United States.,South Ural State University , Lenin pr. 76 , Chelyabinsk 454080 , Russia.,National Research Nuclear University MEPhI , Kashirskoye shosse 31 , Moscow 115409 , Russia
| | - Timothy C Allison
- Southwest Research Institute , San Antonio , Texas 78238 , United States
| | - Sungho Chang
- KEPCO Research Institute , Daejeon 34050 , Korea
| | - Chansun Lim
- Hanwha Power Systems , Seongnam , Gyeonggi 13488 , Korea
| | - Yuin Jin
- Hanwha Power Systems , Seongnam , Gyeonggi 13488 , Korea
| | - Subith S Vasu
- Center for Advanced Turbomachinery and Energy Research (CATER), Mechanical and Aerospace Engineering , University of Central Florida , Orlando , Florida 32816 , United States
| |
Collapse
|
13
|
Shannon RJ, Hornung B, Tew DP, Glowacki DR. Anharmonic Molecular Mechanics: Ab Initio Based Morse Parametrizations for the Popular MM3 Force Field. J Phys Chem A 2019; 123:2991-2999. [PMID: 30793911 DOI: 10.1021/acs.jpca.8b12006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.
Collapse
Affiliation(s)
- R J Shannon
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K.,Department of Mechanical Engineering , Stanford University , 452 Escondido Mall , Stanford , California 94305 , United States
| | - B Hornung
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K
| | - D P Tew
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K
| | - D R Glowacki
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K.,Department of Computer Science , University of Bristol , Bristol BS8 1UB , U.K
| |
Collapse
|
14
|
Amabilino S, Bratholm LA, Bennie SJ, Vaucher AC, Reiher M, Glowacki DR. Training Neural Nets To Learn Reactive Potential Energy Surfaces Using Interactive Quantum Chemistry in Virtual Reality. J Phys Chem A 2019; 123:4486-4499. [DOI: 10.1021/acs.jpca.9b01006] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Silvia Amabilino
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Lars A. Bratholm
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Simon J. Bennie
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Alain C. Vaucher
- Laboratory of Physical Chemistry, ETH Zurich, Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Zurich, Switzerland
| | | |
Collapse
|
15
|
Molecular dynamics of combustion reactions in supercritical carbon dioxide. Part 4: boxed MD study of formyl radical dissociation and recombination. J Mol Model 2019; 25:35. [PMID: 30631947 DOI: 10.1007/s00894-018-3912-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Accepted: 12/18/2018] [Indexed: 10/27/2022]
Abstract
Fossil fuel oxy-combustion is an emergent technology where habitual nitrogen diluent is replaced by high pressure (supercritical) carbon dioxide. The supercritical state of CO2 increases the efficiency of the energy conversion and the absence of nitrogen from the reaction mixture reduces pollution by NOx. However, the effects of a supercritical environment on elementary reactions kinetics are not well understood at present. We used boxed molecular dynamics simulations at the QM/MM theory level to predict the kinetics of dissociation/recombination reaction HCO• + [M] ↔ H• + CO + [M], an important elementary step in many combustion processes. A wide range of temperatures (400-1600 K) and pressures (0.3-1000 atm) were studied. Potentials of mean force were plotted and used to predict activation free energies and rate constants. Based on the data obtained, extended Arrhenius equation parameters were fitted and tabulated. The apparent activation energy for the recombination reaction becomes negative above 30 atm. As the temperature increased, the pressure effect on the rate constant decreased. While at 400 K the pressure increase from 0.3 atm to 300 atm accelerated the dissociation reaction by a factor of 250, at 1600 K the same pressure increase accelerated this reaction by a factor of 100. Graphical abstract Formyl radical surrounded by carbon dioxide molecules.
Collapse
|
16
|
Zhang X, Harvey JN. EVB and polarizable MM study of energy relaxation in fluorine–acetonitrile reactions. Phys Chem Chem Phys 2019; 21:14331-14340. [DOI: 10.1039/c8cp06686h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Many-body effects can impact on rates of energy transfer from a ‘hot’ DF solute to acetonitrile solvent.
Collapse
Affiliation(s)
- Xiaoyong Zhang
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry
- KU Leuven
- B-3001 Leuven
- Belgium
| | - Jeremy N. Harvey
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry
- KU Leuven
- B-3001 Leuven
- Belgium
| |
Collapse
|
17
|
Classical trajectory studies of collisional energy transfer. ACTA ACUST UNITED AC 2019. [DOI: 10.1016/b978-0-444-64207-3.00003-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
18
|
Kulkarni Y, Kamerlin SCL. Computational physical organic chemistry using the empirical valence bond approach. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.07.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
19
|
Garcia-Meseguer R, Carpenter BK. Re-Evaluating the Transition State for Reactions in Solution. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800841] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Barry K. Carpenter
- School of Chemistry; Cardiff University; CF10 3AT Cardiff United Kingdom
| |
Collapse
|
20
|
Meuwly M. Reactive molecular dynamics: From small molecules to proteins. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1386] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Markus Meuwly
- Department of Chemistry University of Basel Basel Switzerland
- Department of Chemistry Brown University Providence Rhode Island
| |
Collapse
|
21
|
Athokpam B, Ramesh SG. Alkyl hydrogen atom abstraction reactions of the CN radical with ethanol. J Chem Phys 2018; 148:134503. [PMID: 29626852 DOI: 10.1063/1.5021634] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a study of the abstraction of alkyl hydrogen atoms from the β and α positions of ethanol by the CN radical in solution using the Empirical Valence Bond (EVB) method. We have built separate 2 × 2 EVB models for the Hβ and Hα reactions, where the atom transfer is parameterized using ab initio calculations. The intra- and intermolecular potentials of the reactant and product molecules were modelled with the General AMBER Force Field, with some modifications. We have carried out the dynamics in water and chloroform, which are solvents of contrasting polarity. We have computed the potential of mean force for both abstractions in each of the solvents. They are found to have a small and early barrier along the reaction coordinate with a large energy release. Analyzing the solvent structure around the reaction system, we have found two solvents to have little effect on either reaction. Simulating the dynamics from the transition state, we also study the fate of the energies in the HCN vibrational modes. The HCN molecule is born vibrationally hot in the CH stretch in both reactions and additionally in the HCN bends for the Hα abstraction reaction. In the early stage of the dynamics, we find that the CN stretch mode gains energy at the expense of the energy in CH stretch mode.
Collapse
Affiliation(s)
- Bijyalaxmi Athokpam
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| | - Sai G Ramesh
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
22
|
Essafi S, Harvey JN. Rates of Molecular Vibrational Energy Transfer in Organic Solutions. J Phys Chem A 2018; 122:3535-3540. [DOI: 10.1021/acs.jpca.7b12563] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stéphanie Essafi
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Jeremy N. Harvey
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| |
Collapse
|
23
|
Abstract
The dynamics of chemical reactions in liquid solutions are now amenable to direct study using ultrafast laser spectroscopy techniques and advances in computer simulation methods. The surrounding solvent affects the chemical reaction dynamics in numerous ways, which include: (i) formation of complexes between reactants and solvent molecules; (ii) modifications to transition state energies and structures relative to the reactants and products; (iii) coupling between the motions of the reacting molecules and the solvent modes, and exchange of energy; (iv) solvent caging of reactants and products; and (v) structural changes to the solvation shells in response to the changing chemical identity of the solutes, on timescales which may be slower than the reactive events. This article reviews progress in the study of bimolecular chemical reaction dynamics in solution, concentrating on reactions which occur on ground electronic states. It illustrates this progress with reference to recent experimental and computational studies, and considers how the various ways in which a solvent affects the chemical reaction dynamics can be unravelled. Implications are considered for research in fields such as mechanistic synthetic chemistry.
Collapse
Affiliation(s)
- Andrew J Orr-Ewing
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, UK.
| |
Collapse
|
24
|
Schmid MH, Das AK, Landis CR, Meuwly M. Multi-State VALBOND for Atomistic Simulations of Hypervalent Molecules, Metal Complexes, and Reactions. J Chem Theory Comput 2018; 14:3565-3578. [DOI: 10.1021/acs.jctc.7b01210] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Maurus H. Schmid
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Akshaya Kumar Das
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Clark R. Landis
- Department of Chemistry, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
25
|
Panteleev SV, Masunov AE, Vasu SS. Molecular Dynamics Study of Combustion Reactions in a Supercritical Environment. Part 2: Boxed MD Study of CO + OH → CO2 + H Reaction Kinetics. J Phys Chem A 2018; 122:897-908. [DOI: 10.1021/acs.jpca.7b09774] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sergey V. Panteleev
- NanoScienece
Technology Center, University of Central Florida, 12424 Research
Parkway, Suite 400, Orlando, Florida 32826, United States
- N. I. Lobachevsky State University of Nizhny Novgorod, Gagarin Av. 23, Nizhny Novgorod 603950, Russia
| | - Artëm E. Masunov
- NanoScienece
Technology Center, University of Central Florida, 12424 Research
Parkway, Suite 400, Orlando, Florida 32826, United States
- Department
of Chemistry, and Department of Physics, University of Central Florida, 4111 Libra Drive, Orlando, Florida 32816, United States
- South Ural State University, Lenin pr. 76, Chelyabinsk 454080, Russia
- National Research Nuclear University MEPhI, Kashirskoye shosse 31, Moscow, 115409, Russia
| | - Subith S. Vasu
- Center for
Advanced Turbomachinery and Energy Research (CATER), Mechanical and
Aerospace Engineering, University of Central Florida, Orlando, Florida 32816, United States
| |
Collapse
|
26
|
Shannon R, Glowacki DR. A Simple “Boxed Molecular Kinetics” Approach To Accelerate Rare Events in the Stochastic Kinetic Master Equation. J Phys Chem A 2018; 122:1531-1541. [DOI: 10.1021/acs.jpca.7b12521] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Robin Shannon
- Mechanical Engineering, Stanford University, Stanford, California 94305, United States
- School of Chemistry, University of Bristol, Bristol, BS8 1TS, United Kingdom
| | - David R. Glowacki
- Mechanical Engineering, Stanford University, Stanford, California 94305, United States
- School of Chemistry, University of Bristol, Bristol, BS8 1TS, United Kingdom
- Department of Computer Science, University of Bristol, Bristol, BS8 1UB, United Kingdom
| |
Collapse
|
27
|
Purg M, Kamerlin SCL. Empirical Valence Bond Simulations of Organophosphate Hydrolysis: Theory and Practice. Methods Enzymol 2018; 607:3-51. [DOI: 10.1016/bs.mie.2018.06.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
28
|
Pratihar S, Ma X, Xie J, Scott R, Gao E, Ruscic B, Aquino AJA, Setser DW, Hase WL. Post-transition state dynamics and product energy partitioning following thermal excitation of the F⋯HCH2CN transition state: Disagreement with experiment. J Chem Phys 2017; 147:144301. [DOI: 10.1063/1.4985894] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Subha Pratihar
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Xinyou Ma
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Jing Xie
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Rebecca Scott
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Eric Gao
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA and Computation Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Adelia J. A. Aquino
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
- School of Pharmaceutical Science and Technology, Tianjin University, Tianjin 300072, People’s Republic of China
- Institute for Soil Research University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - Donald W. Setser
- Institute for Soil Research University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - William L. Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| |
Collapse
|
29
|
Spezia R, Martínez-Nuñez E, Vazquez S, Hase WL. Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:20170035. [PMID: 28320909 PMCID: PMC5360905 DOI: 10.1098/rsta.2017.0035] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 01/24/2017] [Indexed: 06/06/2023]
Abstract
In this Introduction, we show the basic problems of non-statistical and non-equilibrium phenomena related to the papers collected in this themed issue. Over the past few years, significant advances in both computing power and development of theories have allowed the study of larger systems, increasing the time length of simulations and improving the quality of potential energy surfaces. In particular, the possibility of using quantum chemistry to calculate energies and forces 'on the fly' has paved the way to directly study chemical reactions. This has provided a valuable tool to explore molecular mechanisms at given temperatures and energies and to see whether these reactive trajectories follow statistical laws and/or minimum energy pathways. This themed issue collects different aspects of the problem and gives an overview of recent works and developments in different contexts, from the gas phase to the condensed phase to excited states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.
Collapse
Affiliation(s)
- Riccardo Spezia
- Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, CEA CNRS Université Paris Saclay, 91025 Evry, France
- LAMBE, Université d'Evry, 91025 Evry, France
| | - Emilio Martínez-Nuñez
- Departamento de Química Física and Centro Singular de Investigación en Química, Biológica y Materiales Moleculares (CIQUS), Campus Vida, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | - Saulo Vazquez
- Departamento de Química Física and Centro Singular de Investigación en Química, Biológica y Materiales Moleculares (CIQUS), Campus Vida, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | - William L Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX 79409, USA
| |
Collapse
|
30
|
Glowacki DR, Rodgers WJ, Shannon R, Robertson SH, Harvey JN. Reaction and relaxation at surface hotspots: using molecular dynamics and the energy-grained master equation to describe diamond etching. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:rsta.2016.0206. [PMID: 28320908 PMCID: PMC5360904 DOI: 10.1098/rsta.2016.0206] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Accepted: 01/11/2017] [Indexed: 06/06/2023]
Abstract
The extent to which vibrational energy transfer dynamics can impact reaction outcomes beyond the gas phase remains an active research question. Molecular dynamics (MD) simulations are the method of choice for investigating such questions; however, they can be extremely expensive, and therefore it is worth developing cheaper models that are capable of furnishing reasonable results. This paper has two primary aims. First, we investigate the competition between energy relaxation and reaction at 'hotspots' that form on the surface of diamond during the chemical vapour deposition process. To explore this, we developed an efficient reactive potential energy surface by fitting an empirical valence bond model to higher-level ab initio electronic structure theory. We then ran 160 000 NVE trajectories on a large slab of diamond, and the results are in reasonable agreement with experiment: they suggest that energy dissipation from surface hotspots is complete within a few hundred femtoseconds, but that a small fraction of CH3 does in fact undergo dissociation prior to the onset of thermal equilibrium. Second, we developed and tested a general procedure to formulate and solve the energy-grained master equation (EGME) for surface chemistry problems. The procedure we outline splits the diamond slab into system and bath components, and then evaluates microcanonical transition-state theory rate coefficients in the configuration space of the system atoms. Energy transfer from the system to the bath is estimated using linear response theory from a single long MD trajectory, and used to parametrize an energy transfer function which can be input into the EGME. Despite the number of approximations involved, the surface EGME results are in reasonable agreement with the NVE MD simulations, but considerably cheaper. The results are encouraging, because they offer a computationally tractable strategy for investigating non-equilibrium reaction dynamics at surfaces for a broader range of systems.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.
Collapse
Affiliation(s)
- David R Glowacki
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, UK
- Department of Mechanical Engineering, Stanford University, 452 Escondido Mall, Stanford, CA 94305, USA
| | - W J Rodgers
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Robin Shannon
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- Department of Mechanical Engineering, Stanford University, 452 Escondido Mall, Stanford, CA 94305, USA
| | - Struan H Robertson
- Dassault Systémes BIOVIA, 334 Cambridge Science Park, Cambridge CB4 0WN, UK
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Heverlee, Belgium
| |
Collapse
|
31
|
Shin JY, Shaloski MA, Crim FF, Case AS. First Evidence of Vibrationally Driven Bimolecular Reactions in Solution: Reactions of Br Atoms with Dimethylsulfoxide and Methanol. J Phys Chem B 2017; 121:2486-2494. [DOI: 10.1021/acs.jpcb.7b00035] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jae Yoon Shin
- Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - Michael A. Shaloski
- Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - F. Fleming Crim
- Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - Amanda S. Case
- Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| |
Collapse
|
32
|
Sisto A, Stross C, van der Kamp MW, O’Connor M, McIntosh-Smith S, Johnson GT, Hohenstein EG, Manby FR, Glowacki DR, Martinez TJ. Atomistic non-adiabatic dynamics of the LH2 complex with a GPU-accelerated ab initio exciton model. Phys Chem Chem Phys 2017; 19:14924-14936. [DOI: 10.1039/c7cp00492c] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present GPU-accelerated ab initio molecular dynamics simulations of nonadiabatic dynamics in the LH2 complex in full atomistic detail.
Collapse
Affiliation(s)
- Aaron Sisto
- PULSE Institute and Department of Chemistry
- Stanford University
- Stanford
- USA
- SLAC National Accelerator Laboratory
| | - Clem Stross
- School of Chemistry
- University of Bristol
- Bristol
- UK
| | | | - Michael O’Connor
- School of Chemistry
- University of Bristol
- Bristol
- UK
- Department of Computer Science
| | | | - Graham T. Johnson
- California Institute for Quantitative Biosciences (QB3)
- University of California
- San Francisco
- USA
- Department of Bioengineering and Therapeutic Sciences
| | | | | | - David R. Glowacki
- School of Chemistry
- University of Bristol
- Bristol
- UK
- Department of Computer Science
| | - Todd J. Martinez
- PULSE Institute and Department of Chemistry
- Stanford University
- Stanford
- USA
- SLAC National Accelerator Laboratory
| |
Collapse
|
33
|
Carpenter BK, Harvey JN, Orr-Ewing AJ. The Study of Reactive Intermediates in Condensed Phases. J Am Chem Soc 2016; 138:4695-705. [DOI: 10.1021/jacs.6b01761] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Barry K. Carpenter
- School
of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, U.K
| | - Jeremy N. Harvey
- Department
of Chemistry, KU Leuven, Celestijnen Laan 200F, B-3001 Heverlee, Belgium
| | - Andrew J. Orr-Ewing
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| |
Collapse
|
34
|
Preston TJ, Hornung B, Pandit S, Harvey JN, Orr-Ewing AJ. Dynamical Effects and Product Distributions in Simulated CN + Methane Reactions. J Phys Chem A 2016; 120:4672-82. [DOI: 10.1021/acs.jpca.5b09487] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Thomas J. Preston
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| | - Balázs Hornung
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| | - Shubhrangshu Pandit
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| | - Jeremy N. Harvey
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven (Heverlee), Belgium
| | - Andrew J. Orr-Ewing
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| |
Collapse
|
35
|
O'Connor M, Paci E, McIntosh-Smith S, Glowacki DR. Adaptive free energy sampling in multidimensional collective variable space using boxed molecular dynamics. Faraday Discuss 2016; 195:395-419. [DOI: 10.1039/c6fd00138f] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The past decade has seen the development of a new class of rare event methods in which molecular configuration space is divided into a set of boundaries/interfaces, and then short trajectories are run between boundaries. For all these methods, an important concern is how to generate boundaries. In this paper, we outline an algorithm for adaptively generating boundaries along a free energy surface in multi-dimensional collective variable (CV) space, building on the boxed molecular dynamics (BXD) rare event algorithm. BXD is a simple technique for accelerating the simulation of rare events and free energy sampling which has proven useful for calculating kinetics and free energy profiles in reactive and non-reactive molecular dynamics (MD) simulations across a range of systems, in both NVT and NVE ensembles. Two key developments outlined in this paper make it possible to automate BXD, and to adaptively map free energy and kinetics in complex systems. First, we have generalized BXD to multidimensional CV space. Using strategies from rigid-body dynamics, we have derived a simple and general velocity-reflection procedure that conserves energy for arbitrary collective variable definitions in multiple dimensions, and show that it is straightforward to apply BXD to sampling in multidimensional CV space so long as the Cartesian gradients ∇CV are available. Second, we have modified BXD to undertake on-the-fly statistical analysis during a trajectory, harnessing the information content latent in the dynamics to automatically determine boundary locations. Such automation not only makes BXD considerably easier to use; it also guarantees optimal boundaries, speeding up convergence. We have tested the multidimensional adaptive BXD procedure by calculating the potential of mean force for a chemical reaction recently investigated using both experimental and computational approaches – i.e., F + CD3CN → DF + D2CN in both the gas phase and a strongly coupled explicit CD3CN solvent. The results obtained using multidimensional adaptive BXD agree well with previously published experimental and computational results, providing good evidence for its reliability.
Collapse
Affiliation(s)
- Mike O'Connor
- School of Chemistry
- University of Bristol
- Bristol BS8 1TS, UK
- Department of Computer Science
- University of Bristol
| | - Emanuele Paci
- Astbury Centre for Structural Molecular Biology
- University of Leeds
- Leeds, UK
| | | | - David R. Glowacki
- School of Chemistry
- University of Bristol
- Bristol BS8 1TS, UK
- Department of Computer Science
- University of Bristol
| |
Collapse
|
36
|
Koyama D, Coulter P, Grubb MP, Greetham GM, Clark IP, Orr-Ewing AJ. Reaction Dynamics of CN Radicals in Acetonitrile Solutions. J Phys Chem A 2015; 119:12924-34. [DOI: 10.1021/acs.jpca.5b10720] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daisuke Koyama
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| | - Philip Coulter
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| | - Michael P. Grubb
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| | - Gregory M. Greetham
- Central
Laser Facility, Research Complex at Harwell, Science and Technology
Facilities Council, Rutherford Appleton Laboratory, Harwell Oxford, Didcot, Oxfordshire OX11 0QX, U.K
| | - Ian P. Clark
- Central
Laser Facility, Research Complex at Harwell, Science and Technology
Facilities Council, Rutherford Appleton Laboratory, Harwell Oxford, Didcot, Oxfordshire OX11 0QX, U.K
| | - Andrew J. Orr-Ewing
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| |
Collapse
|
37
|
Hornung B, Preston TJ, Pandit S, Harvey JN, Orr-Ewing AJ. Computational Study of Competition between Direct Abstraction and Addition-Elimination in the Reaction of Cl Atoms with Propene. J Phys Chem A 2015; 119:9452-64. [PMID: 26288318 DOI: 10.1021/acs.jpca.5b07052] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quasi-classical trajectory calculations on a newly constructed and full-dimensionality potential energy surface (PES) examine the dynamics of the reaction of Cl atoms with propene. The PES is an empirical valence bond (EVB) fit to high-level ab initio energies and incorporates deep potential energy wells for the 1-chloropropyl and 2-chloropropyl radicals, a direct H atom abstraction route to HCl + allyl radical (CH2CHCH2(•)) products (Δ(r)H(298K)(⊖) = −63.1 kJ mol(-1)), and a pathway connecting these regions. In total, 94 000 successful reactive trajectories were used to compute distributions of angular scattering and HCl vibrational and rotational level populations. These measures of the reaction dynamics agree satisfactorily with available experimental data. The dominant reaction pathway is direct abstraction of a hydrogen atom from the methyl group of propene occurring in under 500 fs. Less than 10% of trajectories follow an addition–elimination route via the two isomeric chloropropyl radicals. Large amplitude motions of the Cl about the propene molecular framework couple the addition intermediates to the direct abstraction pathway. The EVB method provides a good description of the complicated PES for the Cl + propene reaction despite fitting to a limited number of ab initio points, with the further advantage that dynamics specific to certain mechanisms can be studied in isolation by switching off coupling terms in the EVB matrix connecting different regions of the PES.
Collapse
Affiliation(s)
- Balázs Hornung
- School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Thomas J Preston
- School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Shubhrangshu Pandit
- School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven , Celestijnenlaan 200F, B-3001 Leuven (Heverlee), Belgium
| | - Andrew J Orr-Ewing
- School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|