1
|
Tasi DA, Czakó G. Benchmark ab initio characterization of the complex potential energy surfaces of the HOO - + CH 3Y [Y = F, Cl, Br, I] reactions. Phys Chem Chem Phys 2024; 26:16048-16059. [PMID: 38779842 DOI: 10.1039/d4cp01071j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The α-effect is a well-known phenomenon in organic chemistry, and is related to the enhanced reactivity of nucleophiles involving one or more lone-pair electrons adjacent to the nucleophilic center. The gas-phase bimolecular nucleophilic substitution (SN2) reactions of α-nucleophile HOO- with methyl halides have been thoroughly investigated experimentally and theoretically; however, these investigations have mainly focused on identifying and characterizing the α-effect of HOO-. Here, we perform the first comprehensive high-level ab initio mapping for the HOO- + CH3Y [Y = F, Cl, Br and I] reactions utilizing the modern explicitly-correlated CCSD(T)-F12b method with the aug-cc-pVnZ [n = 2-4] basis sets. The present ab initio characterization considers five distinct product channels of SN2: (CH3OOH + Y-), proton abstraction (CH2Y- + H2O2), peroxide ion substitution (CH3OO- + HY), SN2-induced elimination (CH2O + HY + HO-) and SN2-induced rearrangement (CH2(OH)O- + HY). Moreover, besides the traditional back-side attack Walden inversion, the pathways of front-side attack, double inversion and halogen-bond complex formation have also been explored for SN2. With regard to the Walden inversion of HOO- + CH3Cl, the previously unaddressed discrepancies concerning the geometry of the corresponding transition state are clarified. For the HOO- + CH3F reaction, the recently identified SN2-induced elimination is found to be more exothermic than the SN2 channel, submerged by ∼36 kcal mol-1. The accuracy of our high-level ab initio calculations performed in the present study is validated by the fact that our new benchmark 0 K reaction enthalpies show excellent agreement with the experimental data in nearly all cases.
Collapse
Affiliation(s)
- Domonkos A Tasi
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
2
|
Teramoto H, Saito T, Aoki M, Murayama B, Kobayashi M, Nakamura T, Taketsugu T. Reproducing the Reaction Route Map on the Shape Space from Its Quotient by the Complete Nuclear Permutation-Inversion Group. J Chem Theory Comput 2023; 19:5886-5896. [PMID: 37642714 DOI: 10.1021/acs.jctc.3c00500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
This study develops an algorithm to reproduce reaction route maps (RRMs) in the shape space from the outputs of potential search algorithms. To demonstrate the algorithm, global reaction route mapping is utilized as a potential search algorithm, but the proposed algorithm should work with other potential search algorithms in principle. The proposed algorithm does not require any encoding of the molecular configurations and is thus applicable to complicated realistic molecules for which efficient encoding is not readily available. We show that subgraphs of an RRM mapped to each other by the action of the symmetry group are isomorphic and also provide an algorithm to compute the set of feasible transformations in the sense of Longuet-Higgins. We demonstrate the proposed algorithm in toy models and in more realistic molecules. Finally, we remark on absolute rate theory from our perspective.
Collapse
Affiliation(s)
- Hiroshi Teramoto
- Faculty of Engineering Science, Kansai University, Suita 564-8680, Japan
| | - Takuya Saito
- Department of Mathematics, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
- Department of Economics and Statistics, University of Turin, 10124 Turin, Italy
| | - Masamitsu Aoki
- Department of Mathematics, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
| | - Burai Murayama
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-0810, Japan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
- WPI-ICReDD, Hokkaido University, Sapporo 001-0021, Japan
| | - Takenobu Nakamura
- National Institute of Advanced Industrial Science and Technology, Tsukuba 305-8568, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
- WPI-ICReDD, Hokkaido University, Sapporo 001-0021, Japan
| |
Collapse
|
3
|
Tasi DA, Michaelsen T, Wester R, Czakó G. Quasi-classical trajectory study of the OH - + CH 3I reaction: theory meets experiment. Phys Chem Chem Phys 2023; 25:4005-4014. [PMID: 36649119 DOI: 10.1039/d2cp05553h] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Regarding OH- + CH3I, several studies have focused on the dynamics of the reaction. Here, high-level quasi-classical trajectory simulations are carried out at four different collision energies on our recently developed potential energy surface. In all, more than half a million trajectories are performed, and for the first time, the detailed quasi-classical trajectory results are compared with the reanalysed crossed-beam ion imaging experiments. Concerning the previously reported direct dynamics study of OH- + CH3I, a better agreement can be obtained between the revised experiment and our novel theoretical results. Furthermore, in the present work, the benchmark geometries, frequencies and relative energies of the stationary points are also determined for the OH- + CH3I proton-abstraction channel along with the earlier characterized SN2 channel.
Collapse
Affiliation(s)
- Domonkos A Tasi
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Tim Michaelsen
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Gábor Czakó
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
4
|
Tsutsumi T, Ono Y, Taketsugu T. Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene. J Chem Theory Comput 2022; 18:7483-7495. [PMID: 36351076 DOI: 10.1021/acs.jctc.2c00560] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.
Collapse
Affiliation(s)
- Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo060-0810, Japan.,L-Station, Creative Research Institution (CRI), Hokkaido University, Sapporo060-0812, Japan
| | - Yuriko Ono
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo001-0021, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo001-0021, Japan
| |
Collapse
|
5
|
Abstract
Differences in entropies of competing transition states can direct kinetic selectivity. Understanding and modeling such entropy differences at the molecular level is complicated by the fact that entropy is statistical in nature; i.e., it depends on multiple vibrational states of transition structures, the existence of multiple dynamically accessible pathways past these transition structures, and contributions from multiple transition structures differing in conformation/configuration. The difficulties associated with modeling each of these contributors are discussed here, along with possible solutions, all with an eye toward the development of portable qualitative models of use to experimentalists aiming to design reactions that make use of entropy to control kinetic selectivity.
Collapse
Affiliation(s)
- Dean J Tantillo
- Department of Chemistry, University of California-Davis, 1 Shields Ave, Davis, California 95616, United States
| |
Collapse
|