1
|
Dékány AÁ, Czakó G. Detailed quasiclassical dynamics of the F - + SiH 3Cl multi-channel reaction. Phys Chem Chem Phys 2024; 26:10008-10020. [PMID: 38482549 DOI: 10.1039/d4cp00048j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/28/2024]
Abstract
We report a detailed quasiclassical trajectory study on the F- + SiH3Cl multi-channel reaction using a full-dimensional ab initio analytical potential energy surface. Reaction probabilities, cross sections, initial attack and scattering angle distributions as well as product relative translational, internal, vibrational, and rotational energy distributions are obtained in the collision energy range of 1-40 kcal mol-1 for the following channels: SiH3F + Cl-, SiH2Cl- + HF, SiH2F- + HCl, SiH2FCl + H-, SiH2 + FHCl-, and SiHFCl- + H2. All the channels are translationally cold indicating indirect mechanisms, except proton transfer (SiH2Cl- + HF), which shows mixed direct-indirect character. The angular distributions vary depending on collision energy and inversion/retention for SiH3F + Cl-. In the case of SiH2Cl- + HF front-side/back-side attack backward-forward/forward scattering preference is found at low/high collision energy. SiH2F- + HCl is formed with isotropic scattering and the preferred angle of attack is similar to the SiH3F + Cl- channel. SiH2FCl + H-/SiH2 + FHCl- favors back-side attack and isotropic/backward scattering, whereas SiHFCl- + H2 does not show significant angular preference.
Collapse
Affiliation(s)
- Attila Á Dékány
- 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
|
de Andrade KN, Peixoto BP, Carneiro JWDM, Fiorot RG. Exploring borderline S N1-S N2 mechanisms: the role of explicit solvation protocols in the DFT investigation of isopropyl chloride. RSC Adv 2024; 14:4692-4701. [PMID: 38318615 PMCID: PMC10841197 DOI: 10.1039/d4ra00066h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2024] [Accepted: 01/30/2024] [Indexed: 02/07/2024] Open
Abstract
Nucleophilic substitution at saturated carbon is a crucial class of organic reactions, playing a pivotal role in various chemical transformations that yield valuable compounds for society. Despite the well-established SN1 and SN2 mechanisms, secondary substrates, particularly in solvolysis reactions, often exhibit a borderline pathway. A molecular-level understanding of these processes is fundamental for developing more efficient chemical transformations. Typically, quantum-chemical simulations of the solvent medium combine explicit and implicit solvation methods. The configuration of explicit molecules can be defined through top-down approaches, such as Monte Carlo (MC) calculations for generating initial configurations, and bottom-up methods that involve user-dependent protocols to add solvent molecules around the substrate. Herein, we investigated the borderline mechanism of the hydrolysis of a secondary substrate, isopropyl chloride (iPrCl), at DFT-M06-2X/aug-cc-pVDZ level, employing explicit and explicit + implicit protocols. Top-down and bottom-up approaches were employed to generate substrate-solvent complexes of varying number (n = 1, 3, 5, 7, 9, and 12) and configurations of H2O molecules. Our findings consistently reveal that regardless of the solvation approach, the hydrolysis of iPrCl follows a loose-SN2-like mechanism with nucleophilic solvent assistance. Increasing the water cluster around the substrate in most cases led to reaction barriers of ΔH‡ ≈ 21 kcal mol-1, with nine water molecules from MC configurations sufficient to describe the reaction. The More O'Ferrall-Jencks plot demonstrates an SN1-like character for all transition state structures, showing a clear merged profile. The fragmentation activation strain analyses indicate that energy barriers are predominantly controlled by solvent-substrate interactions, supported by the leaving group stabilization assessed through CHELPG atomic charges.
Collapse
Affiliation(s)
- Karine Nascimento de Andrade
- Department of Organic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| | - Bárbara Pereira Peixoto
- Department of Organic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| | - José Walkimar de Mesquita Carneiro
- Department of Inorganic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| | - Rodolfo Goetze Fiorot
- Department of Organic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| |
Collapse
|
3
|
Ballay B, Szűcs T, Papp D, Czakó G. Phosphorus-centered ion-molecule reactions: benchmark ab initio characterization of the potential energy surfaces of the X - + PH 2Y [X, Y = F, Cl, Br, I] systems. Phys Chem Chem Phys 2023; 25:28925-28940. [PMID: 37855143 DOI: 10.1039/d3cp03733a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
In the present work we determine the benchmark relative energies and geometries of all the relevant stationary points of the X- + PH2Y [X, Y = F, Cl, Br, I] identity and non-identity reactions using state-of-the-art electronic-structure methods. These phosphorus-centered ion-molecule reactions follow two main reaction routes: bimolecular nucleophilic substitution (SN2), leading to Y- + PH2X, and proton transfer, resulting in HX + PHY- products. The SN2 route can proceed through Walden-inversion, front-side-attack retention, and double-/multiple-inversion pathways. In addition, we also identify the following product channels: H--formation, PH2-- and PH2-formation, 1PH- and 3PH-formation, H2-formation and HY + PHX- formation. The benchmark classical relative energies are obtained by taking into account the core-correlation, scalar relativistic, and post-(T) corrections, which turn out to be necessary to reach subchemical (<1 kcal mol-1) accuracy of the results. Classical relative energies are augmented with zero-point-energy contributions to gain the benchmark adiabatic energies.
Collapse
Affiliation(s)
- Boldizsár Ballay
- 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.
| | - Tímea Szűcs
- 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.
| | - Dóra Papp
- 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
|
4
|
Dékány AÁ, Czakó G. Exploring the versatile reactivity of the F- + SiH3Cl system on a full-dimensional coupled-cluster potential energy surface. J Chem Phys 2023; 158:2895234. [PMID: 37290077 DOI: 10.1063/5.0153083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Accepted: 05/18/2023] [Indexed: 06/10/2023] Open
Abstract
We develop a full-dimensional analytical potential energy surface (PES) for the F- + SiH3Cl reaction using Robosurfer for automatically sampling the configuration space, the robust [CCSD-F12b + BCCD(T) - BCCD]/aug-cc-pVTZ composite level of theory for computing the energy points, and the permutationally invariant polynomial method for fitting. Evolution of the fitting error and the percentage of the unphysical trajectories are monitored as a function of the iteration steps/number of energy points and polynomial order. Quasi-classical trajectory simulations on the new PES reveal rich dynamics resulting in high-probability SN2 (SiH3F + Cl-) and proton-transfer (SiH2Cl- + HF) products as well as several lower-probability channels, such as SiH2F- + HCl, SiH2FCl + H-, SiH2 + FHCl-, SiHFCl- + H2, SiHF + H2 + Cl-, and SiH2 + HF + Cl-. The Walden-inversion and front-side-attack-retention SN2 pathways are found to be competitive, producing nearly racemic products at high collision energies. The detailed atomic-level mechanisms of the various reaction pathways and channels as well as the accuracy of the analytical PES are analyzed along representative trajectories.
Collapse
Affiliation(s)
- Attila Á Dékány
- 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
|
5
|
Christian MS, Nenoff TM, Rimsza JM. Effect of Linker Structure and Functionalization on Secondary Gas Formation in Metal-Organic Frameworks. J Phys Chem A 2023; 127:2881-2888. [PMID: 36947182 DOI: 10.1021/acs.jpca.2c07751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2023]
Abstract
Rare-earth terephthalic acid (BDC)-based metal-organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H2O, NO2, and SO2) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC-gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with -OH, -NH2, and -SH were investigated. The calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.
Collapse
Affiliation(s)
- Matthew S Christian
- Geochemistry Department, Sandia National Laboratories, P.O. Box 5800, Eubank Boulevard SE, Albuquerque, New Mexico 87185, United States
| | - Tina M Nenoff
- Advanced Science & Technology, Sandia National Laboratories, P.O. Box 5800, Eubank Boulevard SE, Albuquerque, New Mexico 87185, United States
| | - Jessica M Rimsza
- Geochemistry Department, Sandia National Laboratories, P.O. Box 5800, Eubank Boulevard SE, Albuquerque, New Mexico 87185, United States
| |
Collapse
|
6
|
Jung Y, Yoo SY, Jin Y, You J, Han S, Yu J, Park Y, Cho SH. Iridium-Catalyzed Chemo-, Diastereo-, and Enantioselective Allyl-Allyl Coupling: Accessing All Four Stereoisomers of (E)-1-Boryl-Substituted 1,5-Dienes by Chirality Pairing. Angew Chem Int Ed Engl 2023; 62:e202218794. [PMID: 36718077 DOI: 10.1002/anie.202218794] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Revised: 01/29/2023] [Accepted: 01/30/2023] [Indexed: 02/01/2023]
Abstract
Here, we report a highly chemo-, diastereo-, and enantioselective allyl-allyl coupling between branched allyl alcohols and α-silyl-substituted allylboronate esters, catalyzed by a chiral iridium complex. The α-silyl-substituted allylboronate esters can be chemoselectively coupled with allyl electrophiles, affording a diverse set of enantioenriched (E)-1-boryl-substituted 1,5-dienes in good yields, with excellent stereoselectivity. By permuting the chiral iridium catalysts and the substrates, we efficiently and selectively obtained all four stereoisomers bearing two consecutive chiral centers. Mechanistic studies via density functional theory calculations revealed the origins of the diastereo- and chemoselectivities, indicating the pivotal roles of the steric interaction, the β-silicon effect, and a rapid desilylation process. Additional synthetic modifications for preparing a variety of enantioenriched compounds containing contiguous chiral centers are also included.
Collapse
Affiliation(s)
- Yongsuk Jung
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 (Republic of, Korea
| | - Seok Yeol Yoo
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141 (Republic of, Korea
| | - Yonghoon Jin
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 (Republic of, Korea
| | - Jaehyun You
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141 (Republic of, Korea
| | - Seungcheol Han
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 (Republic of, Korea
| | - Jeongwoo Yu
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 (Republic of, Korea
| | - Yoonsu Park
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141 (Republic of, Korea
| | - Seung Hwan Cho
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 37673 (Republic of, Korea
| |
Collapse
|
7
|
Zhou B, Zhou Y, Xie D. Accelerated Quantum Mechanics/Molecular Mechanics Simulations via Neural Networks Incorporated with Mechanical Embedding Scheme. J Chem Theory Comput 2023; 19:1157-1169. [PMID: 36724190 DOI: 10.1021/acs.jctc.2c01131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
A powerful tool to study the mechanism of reactions in solutions or enzymes is to perform the ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations. However, the computational cost is too high due to the explicit electronic structure calculations at every time step of the simulation. A neural network (NN) method can accelerate the QM/MM-MD simulations, but it has long been a problem to accurately describe the QM/MM electrostatic coupling by NN in the electrostatic embedding (EE) scheme. In this work, we developed a new method to accelerate QM/MM calculations in the mechanic embedding (ME) scheme. The potentials and partial point charges of QM atoms are first learned in vacuo by the embedded atom neural networks (EANN) approach. MD simulations are then performed on this EANN/MM potential energy surface (PES) to obtain free energy (FE) profiles for reactions, in which the QM/MM electrostatic coupling is treated in the mechanic embedding (ME) scheme. Finally, a weighted thermodynamic perturbation (wTP) corrects the FE profiles in the ME scheme to the EE scheme. For two reactions in water and one in methanol, our simulations reproduced the B3LYP/MM free energy profiles within 0.5 kcal/mol with a speed-up of 30-60-fold. The results show that the strategy of combining EANN potential in the ME scheme with the wTP correction is efficient and reliable for chemical reaction simulations in liquid. Another advantage of our method is that the QM PES is independent of the MM subsystem, so it can be applied to various MM environments as demonstrated by an SN2 reaction studied in water and methanol individually, which used the same EANN PES. The free energy profiles are in excellent accordance with the results obtained from B3LYP/MM-MD simulations. In future, this method will be applied to the reactions of enzymes and their variants.
Collapse
Affiliation(s)
- Boyi Zhou
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yanzi Zhou
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China.,Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
8
|
Zhang F, Zheng W, Yang F, Ma Z, Sun W, Zhao L. Understanding the Reaction Kinetics and Microdynamics between Methylimidazole and Alkyl Thiocyanate for Ionic Liquid Synthesis through Experiments and Theoretical Calculation. Ind Eng Chem Res 2023. [DOI: 10.1021/acs.iecr.2c04323] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/26/2023]
Affiliation(s)
- Fan Zhang
- State Key Laboratory of Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
| | - Weizhong Zheng
- State Key Laboratory of Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
| | - Fan Yang
- State Key Laboratory of Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
| | - Zhihong Ma
- State Key Laboratory of Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
| | - Weizhen Sun
- State Key Laboratory of Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
| | - Ling Zhao
- State Key Laboratory of Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
| |
Collapse
|
9
|
Taylor M, Yu H, Ho J. Predicting Solvent Effects on S N2 Reaction Rates: Comparison of QM/MM, Implicit, and MM Explicit Solvent Models. J Phys Chem B 2022; 126:9047-9058. [PMID: 36300819 DOI: 10.1021/acs.jpcb.2c06000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Solvents are one of the key variables in the optimization of a synthesis yield or properties of a synthesis product. In this paper, contemporary solvent models are applied to predict the rates of SN2 reactions in a range of aqueous and non-aqueous solvents. High-level CCSD(T)/CBS//M06-2X/6-31+G(d) gas phase energies were combined with solvation free energies from SMD, SM12, and ADF-COSMO-RS continuum solvent models, as well as molecular mechanics (MM) explicit solvent models with different atomic charge schemes to predict the rate constants of three SN2 reactions in eight protic and aprotic solvents. It is revealed that the prediction of rate constants in organic solvents is not necessarily less challenging than in water and popular solvent models struggle to predict their rate constants to within 3 log units of experimental values. Among the continuum solvent models, the ADF-COSMO-RS model performed the best in predicting absolute rate contants while the SM12 model was best at predicting relative rate constants with an average accuracy of about 1.5 and 0.8 log units, respectively. The use of computationally more demanding MM explicit solvent models did not translate to improvements in absolute rate constants but was quite effective at predicting relative rate constants due to systematic error cancellation. Free energy barriers obtained from umbrella sampling with explicit solvent QM/MM simulations led to excellent agreement with experimental values, provided that a validated level of theory is used to treat the QM region.
Collapse
Affiliation(s)
- Mackenzie Taylor
- School of Chemistry, The University of New South Wales, Sydney, New South Wales2052, Australia
| | - Haibo Yu
- Molecular Horizons and School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong, New South Wales2522, Australia
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, New South Wales2052, Australia
| |
Collapse
|
10
|
Isamura BK, Lobb KA. A multiscale ONIOM study of the buckminsterfullerene (C 60) Diels-Alder reaction: from model design to reaction path analysis. J Mol Model 2022; 28:327. [PMID: 36138156 DOI: 10.1007/s00894-022-05319-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Accepted: 09/06/2022] [Indexed: 11/29/2022]
Abstract
The hybrid ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) formalism is employed to investigate the Diels-Alder reaction of the buckminsterfullerene C60. Our computations suggest that the ONIOM2(M06-2X/6-31G(d): SVWN/STO-3G) model, enclosing both the diene and the pyracyclene fragment of C60 in the higher-layer, provides a reasonable trade-off between accuracy and computational cost as it comes to predicting reaction energetics. Moreover, the frontier molecular orbital (FMO) theory and activation strain model (ASM) are jointly relied on to rationalize the effect of -OH and -CN substituents on the activation barrier of this reaction. Finally, reaction paths are scrutinized to get insight into the various forces underpinning the process of cycloadduct formation.
Collapse
Affiliation(s)
| | - Kevin Alan Lobb
- Department of Chemistry, Rhodes University, Eastern Cape, Makhanda, 6140, South Africa.,Research Unit in Bioinformatics (RUBi), Rhodes University, Eastern Cape, Makhanda, 6140, South Africa
| |
Collapse
|
11
|
Wu X, Zhang S, Xie J. Investigating the competing E2 and S N2 mechanisms for the microsolvated HO -(H 2O) n=0-4 + CH 3CH 2X (X = Cl, Br, I) reactions. Phys Chem Chem Phys 2022; 24:12993-13005. [PMID: 35582984 DOI: 10.1039/d1cp04010c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We characterized the anti-E2, syn-E2, inv-SN2, and ret-SN2 reaction channels for the reaction of microsolvated HO-(H2O)n anions with CH3CH2X (X = Cl, Br, I), using the CCSD(T)/PP/t//MP2/ECP/d level method, to understand how a solvent influences the competing E2 and SN2 reactions. The calculated sequence of barrier for the four channels is ret-SN2 > syn-E2 > anti-E2 > inv-SN2. The barrier heights increase with incremental hydration as the system transfers from the gas phase to microsolvation, and to bulk solvation (using the PCM implicit solvent model). As the degree of hydration n increases, good correlations have been found between barrier heights and several thermodynamic, geometric and charge parameters, including the reaction enthalpy, proton/ethyl-cation affinity of the hydrated nucleophile, geometric looseness (%L‡) and asymmetry (%AS‡) and charge asymmetry (Δq(X-O)) of the transition structures. Under a molecular orbital scheme, the HOMOs of nucleophiles are stabilized by stepwise hydration, explaining the rise in the barriers. Considering the effect of the leaving group, the barrier heights exhibit linear correlation with the halogen electronegativity and H-acidity of substrate CH3CH2X. In terms of E2/SN2 competition, the barrier difference, , first increases then decreases as the number of explicit water molecules increases, under both microsolvation and bulk solvation conditions, but the inv-SN2 pathway is always favored over the anti-E2 pathway. Energy decomposition analysis attributes the increase of barrier difference to the greater geometric distortion in the anti-E2 transition structure.
Collapse
Affiliation(s)
- Xiangyu Wu
- Key Laboratory of Cluster Science of Ministry of Education, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| | - Shaowen Zhang
- Key Laboratory of Cluster Science of Ministry of Education, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| | - Jing Xie
- Key Laboratory of Cluster Science of Ministry of Education, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| |
Collapse
|
12
|
Bortoli M, Campeggio J, Orian L, Zerbetto M, Polimeno A. Parameter free evaluation of S N2 reaction rates for halide substitution in halomethane. Phys Chem Chem Phys 2022; 24:7474-7480. [PMID: 35274635 DOI: 10.1039/d1cp05220a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We estimate the kinetic constants of a series of archetypal SN2 reactions, i.e., the nucleophilic substitutions of halides in halomethane. A parameter free, multiscale approach recently developed [Campeggio et al., Phys. Chem. Chem. Phys., 2020, 22, 3455] is employed. The protocol relies on quantum mechanical calculations for the description of the energy profile along the intrinsic reaction coordinate, which is then mapped onto a reaction coordinate conveniently built for the reactive process. A Kramers-Klein equation is used to describe the stochastic time evolution of the reaction coordinate and its velocity; friction is parameterized using a hydrodynamic model and Kramers theory is used to derive the rate constant of the reaction. The method is here applied to six SN2 reactions in water at 295.15 K, which differ in the nucleophile and the leaving group. The computed reaction rates are in good agreement with the experimental data and correlate well with the trends observed for the activation energies.
Collapse
Affiliation(s)
- Marco Bortoli
- Institut de Química Computacional i Catàlisi (IQCC) i Departament de Química, Facultat de Ciències, Universitat de Girona, C/M. A. Capmany 69, 17003 Girona, Spain
| | - Jonatan Campeggio
- Consorzio interuniversitario per lo sviluppo dei Sistemi a Grande Interfase, via della Lastruccia 3, 50019 Sesto Fiorentino, Firenze, Italy
| | - Laura Orian
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, Padova, Italy.
| | - Mirco Zerbetto
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, Padova, Italy.
| | - Antonino Polimeno
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, Padova, Italy.
| |
Collapse
|
13
|
Dékány AÁ, Kovács GZ, Czakó G. High-Level Systematic Ab Initio Comparison of Carbon- and Silicon-Centered S N2 Reactions. J Phys Chem A 2021; 125:9645-9657. [PMID: 34709818 PMCID: PMC8591615 DOI: 10.1021/acs.jpca.1c07574] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
We characterize the
stationary points along the Walden inversion,
front-side attack, and double-inversion pathways of the X– + CH3Y and X– + SiH3Y [X,
Y = F, Cl, Br, I] SN2 reactions using chemically accurate
CCSD(T)-F12b/aug-cc-pVnZ [n = D,
T, Q] levels of theory. At the carbon center, Walden inversion dominates
and proceeds via prereaction (X–···H3CY) and postreaction (XCH3···Y–) ion-dipole wells separated by a usually submerged
transition state (X–H3C–Y)−, front-side attack occurs over high barriers, double inversion is
the lowest-energy retention pathway for X = F, and hydrogen- (F–···HCH2Y) and halogen-bonded
(X–···YCH3) complexes
exist in the entrance channel. At the silicon center, Walden inversion
proceeds through a single minimum (X–SiH3–Y)−, the front-side attack is competitive via a usually
submerged transition state separating pre- and postreaction minima
having X–Si–Y angles close to 90°, double inversion
occurs over positive, often high barriers, and hydrogen- and halogen-bonded
complexes are not found. In addition to the SN2 channels
(Y– + CH3X/SiH3X), we report
reaction enthalpies for proton abstraction (HX + CH2Y–/SiH2Y–), hydride substitution
(H– + CH2XY/SiH2XY), XH···Y– complex formation (XH···Y– + 1CH2/1SiH2), and halogen
abstraction (XY + CH3–/SiH3– and XY– + CH3/SiH3).
Collapse
Affiliation(s)
- Attila Á Dékány
- 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
| | - Gyula Z Kovács
- 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
|
14
|
Alkorta I, Elguero J. The SN2 reaction and its relationship with the Walden inversion, the Finkelstein and Menshutkin reactions together with theoretical calculations for the Finkelstein reaction. Struct Chem 2021. [DOI: 10.1007/s11224-021-01805-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
AbstractThis communication gives an overview of the relationships between four reactions that although related were not always perceived as such: SN2, Walden, Finkelstein, and Menshutkin. Binary interactions (SN2 & Walden, SN2 & Menshutkin, SN2 & Finkelstein, Walden & Menshutkin, Walden & Finkelstein, Menshutkin & Finkelstein) were reported. Carbon, silicon, nitrogen, and phosphorus as central atoms and fluorides, chlorides, bromides, and iodides as lateral atoms were considered. Theoretical calculations provide Gibbs free energies that were analyzed with linear models to obtain the halide contributions. The M06-2x DFT computational method and the 6-311++G(d,p) basis set have been used for all atoms except for iodine where the effective core potential def2-TZVP basis set was used. Concerning the central atom pairs, carbon/silicon vs. nitrogen/phosphorus, we reported here for the first time that the effect of valence expansion was known for Si but not for P. Concerning the lateral halogen atoms, some empirical models including the interaction between F and I as entering and leaving groups explain the Gibbs free energies.
Collapse
|
15
|
Jackson R, Zhang W, Pearson J. TSNet: predicting transition state structures with tensor field networks and transfer learning. Chem Sci 2021; 12:10022-10040. [PMID: 34377396 PMCID: PMC8317659 DOI: 10.1039/d1sc01206a] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Accepted: 06/21/2021] [Indexed: 12/14/2022] Open
Abstract
Transition states are among the most important molecular structures in chemistry, critical to a variety of fields such as reaction kinetics, catalyst design, and the study of protein function. However, transition states are very unstable, typically only existing on the order of femtoseconds. The transient nature of these structures makes them incredibly difficult to study, thus chemists often turn to simulation. Unfortunately, computer simulation of transition states is also challenging, as they are first-order saddle points on highly dimensional mathematical surfaces. Locating these points is resource intensive and unreliable, resulting in methods which can take very long to converge. Machine learning, a relatively novel class of algorithm, has led to radical changes in several fields of computation, including computer vision and natural language processing due to its aptitude for highly accurate function approximation. While machine learning has been widely adopted throughout computational chemistry as a lightweight alternative to costly quantum mechanical calculations, little research has been pursued which utilizes machine learning for transition state structure optimization. In this paper TSNet is presented, a new end-to-end Siamese message-passing neural network based on tensor field networks shown to be capable of predicting transition state geometries. Also presented is a small dataset of SN2 reactions which includes transition state structures - the first of its kind built specifically for machine learning. Finally, transfer learning, a low data remedial technique, is explored to understand the viability of pretraining TSNet on widely available chemical data may provide better starting points during training, faster convergence, and lower loss values. Aspects of the new dataset and model shall be discussed in detail, along with motivations and general outlook on the future of machine learning-based transition state prediction.
Collapse
Affiliation(s)
- Riley Jackson
- Department of Chemistry, University of Prince Edward Island Canada
| | - Wenyuan Zhang
- Department of Chemistry, University of Prince Edward Island Canada
| | - Jason Pearson
- Department of Chemistry, University of Prince Edward Island Canada
| |
Collapse
|
16
|
Fugel M, Dittmer A, Kleemiss F, Grabowsky S. On the Role of Hydrogen Bonding in Gas-Phase S N2 Reactions at Silicon. J Phys Chem A 2021; 125:4070-4078. [PMID: 33974418 DOI: 10.1021/acs.jpca.1c00601] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The shape of the potential energy surface (PES) of gas-phase SN2 reactions at silicon is determined by the type of nucleophile, the leaving group, and substituents which remain bonded to silicon. In this study, we present PES scans along the reaction coordinate of six symmetrical SN2 reactions: X- + SiR3X → XSiR3 + X-, where X = Cl or F and R = H, Me, or OMe. While the fluorine systems and the ClSiH3Cl system only give single-well PESs, ClSiMe3Cl and ClSi(OMe)3Cl give triple- and double-well PESs with stable pre- and post-reaction complexes. A complementary bonding analysis (energy decomposition analysis, quantum theory of atoms in molecules, and natural bond orbitals) reveals that the leaving group (X-) is stabilized by hydrogen bonding in the XSiMe3X and XSi(OMe)3X systems. It is shown that this so far neglected stabilizing contribution, along with σ-hole bonding, is responsible for the shapes of the PESs of ClSiMe3Cl and ClSi(OMe)3Cl in the gas phase.
Collapse
Affiliation(s)
- Malte Fugel
- Fachbereich 2-Biologie/ Chemie, Universität Bremen, Leobener Str. 3, 28359 Bremen, Germany
| | - Anneke Dittmer
- Fachbereich 2-Biologie/ Chemie, Universität Bremen, Leobener Str. 3, 28359 Bremen, Germany
| | - Florian Kleemiss
- Fachbereich 2-Biologie/ Chemie, Universität Bremen, Leobener Str. 3, 28359 Bremen, Germany.,Departement für Chemie, Biochemie und Pharmazie, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Simon Grabowsky
- Fachbereich 2-Biologie/ Chemie, Universität Bremen, Leobener Str. 3, 28359 Bremen, Germany.,Departement für Chemie, Biochemie und Pharmazie, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
| |
Collapse
|
17
|
Efficient synthesis of amine-functionalized graphene oxide by ultrasound-assisted reactions and density functional theory mechanistic insight. APPLIED NANOSCIENCE 2021. [DOI: 10.1007/s13204-021-01798-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
18
|
Vermeeren P, Hansen T, Jansen P, Swart M, Hamlin TA, Bickelhaupt FM. A Unified Framework for Understanding Nucleophilicity and Protophilicity in the S N 2/E2 Competition. Chemistry 2020; 26:15538-15548. [PMID: 32866336 PMCID: PMC7756690 DOI: 10.1002/chem.202003831] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2020] [Indexed: 11/11/2022]
Abstract
The concepts of nucleophilicity and protophilicity are fundamental and ubiquitous in chemistry. A case in point is bimolecular nucleophilic substitution (SN 2) and base-induced elimination (E2). A Lewis base acting as a strong nucleophile is needed for SN 2 reactions, whereas a Lewis base acting as a strong protophile (i.e., base) is required for E2 reactions. A complicating factor is, however, the fact that a good nucleophile is often a strong protophile. Nevertheless, a sound, physical model that explains, in a transparent manner, when an electron-rich Lewis base acts as a protophile or a nucleophile, which is not just phenomenological, is currently lacking in the literature. To address this fundamental question, the potential energy surfaces of the SN 2 and E2 reactions of X- +C2 H5 Y model systems with X, Y = F, Cl, Br, I, and At, are explored by using relativistic density functional theory at ZORA-OLYP/TZ2P. These explorations have yielded a consistent overview of reactivity trends over a wide range in reactivity and pathways. Activation strain analyses of these reactions reveal the factors that determine the shape of the potential energy surfaces and hence govern the propensity of the Lewis base to act as a nucleophile or protophile. The concepts of "characteristic distortivity" and "transition state acidity" of a reaction are introduced, which have the potential to enable chemists to better understand and design reactions for synthesis.
Collapse
Affiliation(s)
- Pascal Vermeeren
- Department of Theoretical ChemistryAmsterdam Institute of, Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale, Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - Thomas Hansen
- Department of Theoretical ChemistryAmsterdam Institute of, Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale, Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Leiden Institute of ChemistryLeiden UniversityEinsteinweg 55, 2333CCLeidenThe Netherlands
| | - Paul Jansen
- Department of Theoretical ChemistryAmsterdam Institute of, Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale, Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Laboratory of Physical ChemistryETH ZurichVladimir-Prelog-Weg 28093ZurichSwitzerland
| | - Marcel Swart
- ICREAPg. Lluís Companys 2308010BarcelonaSpain
- IQCC & Dept. QuímicaUniversitat de GironaCampus Montilivi (Ciències)17003GironaSpain
| | - Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Institute of, Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale, Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical ChemistryAmsterdam Institute of, Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale, Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Institute for Molecules and MaterialsRadboud UniversityHeyendaalseweg 1356525AJNijmegenThe Netherlands
| |
Collapse
|
19
|
Hansen T, Vermeeren P, Haim A, van Dorp MJH, Codée JDC, Bickelhaupt FM, Hamlin TA. Regioselectivity of Epoxide Ring-Openings via SN
2 Reactions Under Basic and Acidic Conditions. European J Org Chem 2020. [DOI: 10.1002/ejoc.202000590] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Thomas Hansen
- Department of Theoretical Chemistry; Amsterdam Institute of Molecular and Life Sciences (AIMMS); Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
- Leiden Institute of Chemistry; Leiden University; Einsteinweg 55 2333 CC Leiden The Netherlands
| | - Pascal Vermeeren
- Department of Theoretical Chemistry; Amsterdam Institute of Molecular and Life Sciences (AIMMS); Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Anissa Haim
- Department of Theoretical Chemistry; Amsterdam Institute of Molecular and Life Sciences (AIMMS); Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Maarten J. H. van Dorp
- Department of Theoretical Chemistry; Amsterdam Institute of Molecular and Life Sciences (AIMMS); Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Jeroen D. C. Codée
- Leiden Institute of Chemistry; Leiden University; Einsteinweg 55 2333 CC Leiden The Netherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry; Amsterdam Institute of Molecular and Life Sciences (AIMMS); Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
- Institute for Molecules and Materials (IMM); Radboud University; Heyendaalseweg 135 6525 AJ Nijmegen The Netherlands
| | - Trevor A. Hamlin
- Department of Theoretical Chemistry; Amsterdam Institute of Molecular and Life Sciences (AIMMS); Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| |
Collapse
|
20
|
Bettens T, Alonso M, De Proft F, Hamlin TA, Bickelhaupt FM. Ambident Nucleophilic Substitution: Understanding Non-HSAB Behavior through Activation Strain and Conceptual DFT Analyses. Chemistry 2020; 26:3884-3893. [PMID: 31957943 PMCID: PMC7154642 DOI: 10.1002/chem.202000272] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2020] [Indexed: 01/31/2023]
Abstract
The ability to understand and predict ambident reactivity is key to the rational design of organic syntheses. An approach to understand trends in ambident reactivity is the hard and soft acids and bases (HSAB) principle. The recent controversy over the general validity of this principle prompted us to investigate the competing gas-phase SN 2 reaction channels of archetypal ambident nucleophiles CN- , OCN- , and SCN- with CH3 Cl (SN 2@C) and SiH3 Cl (SN 2@Si), using DFT calculations. Our combined analyses highlight the inability of the HSAB principle to correctly predict the reactivity trends of these simple, model reactions. Instead, we have successfully traced reactivity trends to the canonical orbital-interaction mechanism and the resulting nucleophile-substrate interaction energy. The HOMO-LUMO orbital interactions set the trend in both SN 2@C and SN 2@Si reactions. We provide simple rules for predicting the ambident reactivity of nucleophiles based on our Kohn-Sham molecular orbital analysis.
Collapse
Affiliation(s)
- Tom Bettens
- Eenheid Algemene Chemie (ALGC)Vrije Universiteit BrusselPleinlaan 21050BrusselsBelgium
| | - Mercedes Alonso
- Eenheid Algemene Chemie (ALGC)Vrije Universiteit BrusselPleinlaan 21050BrusselsBelgium
| | - Frank De Proft
- Eenheid Algemene Chemie (ALGC)Vrije Universiteit BrusselPleinlaan 21050BrusselsBelgium
| | - Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Institute for Molecules and Materials (IMM)Radboud University NijmegenHeyendaalseweg 1356525AJNijmegenThe Netherlands
| |
Collapse
|
21
|
Understanding chemical reactivity using the activation strain model. Nat Protoc 2020; 15:649-667. [PMID: 31925400 DOI: 10.1038/s41596-019-0265-0] [Citation(s) in RCA: 160] [Impact Index Per Article: 40.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2019] [Accepted: 10/29/2019] [Indexed: 12/20/2022]
Abstract
Understanding chemical reactivity through the use of state-of-the-art computational techniques enables chemists to both predict reactivity and rationally design novel reactions. This protocol aims to provide chemists with the tools to implement a powerful and robust method for analyzing and understanding any chemical reaction using PyFrag 2019. The approach is based on the so-called activation strain model (ASM) of reactivity, which relates the relative energy of a molecular system to the sum of the energies required to distort the reactants into the geometries required to react plus the strength of their mutual interactions. Other available methods analyze only a stationary point on the potential energy surface, but our methodology analyzes the change in energy along a reaction coordinate. The use of this methodology has been proven to be critical to the understanding of reactions, spanning the realms of the inorganic and organic, as well as the supramolecular and biochemical, fields. This protocol provides step-by-step instructions-starting from the optimization of the stationary points and extending through calculation of the potential energy surface and analysis of the trend-decisive energy terms-that can serve as a guide for carrying out the analysis of any given reaction of interest within hours to days, depending on the size of the molecular system.
Collapse
|
22
|
Savoo N, Laloo JZA, Rhyman L, Ramasami P, Bickelhaupt FM, Poater J. Activation Strain Analyses of Counterion and Solvent Effects on the Ion-Pair S N 2 Reaction of NH 2 - and CH 3 Cl. J Comput Chem 2019; 41:317-327. [PMID: 31713259 DOI: 10.1002/jcc.26104] [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: 08/27/2019] [Revised: 10/21/2019] [Accepted: 10/21/2019] [Indexed: 11/09/2022]
Abstract
We have computationally studied the bimolecular nucleophilic substitution (SN 2) reactions of Mn NH2 (n-1) + CH3 Cl (M+ = Li+ , Na+ , K+ , and MgCl+ ; n = 0, 1) in the gas phase and in tetrahydrofuran solution at OLYP/6-31++G(d,p) using polarizable continuum model implicit solvation. We wish to explore and understand the effect of the metal counterion M+ and of solvation on the reaction profile and the stereochemical preference, that is, backside (SN 2-b) versus frontside attack (SN 2-f). The results were compared to the corresponding ion-pair SN 2 reactions involving F- and OH- nucleophiles. Our analyses with an extended activation strain model of chemical reactivity uncover and explain various trends in SN 2 reactivity along the nucleophiles F- , OH- , and NH 2 - , including solvent and counterion effects. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Nandini Savoo
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius
| | - Jalal Z A Laloo
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius
| | - Lydia Rhyman
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius.,Department of Chemical Sciences, University of Johannesburg, Doornfontein Campus, Johannesburg, 2028, South Africa
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius.,Department of Chemical Sciences, University of Johannesburg, Doornfontein Campus, Johannesburg, 2028, South Africa
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, NL-1081 HV, Amsterdam, The Netherlands.,Institute for Molecules and Materials, Radboud University Nijmegen, NL-6525 AJ, Nijmegen, The Netherlands
| | - Jordi Poater
- Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, 08028, Barcelona, Spain.,ICREA, 08010, Barcelona, Spain
| |
Collapse
|
23
|
Ruff F, Szabó D, Rábai J, Jalsovszky I, Farkas Ö. Mechanism for the reactions of sulfides with hypochlorous acid and
N
‐chlorosulfon‐amides: Formation of solvated chlorosulfonium cation and λ
4
‐sulfane intermediates. J PHYS ORG CHEM 2019. [DOI: 10.1002/poc.4005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Ferenc Ruff
- Department of Organic Chemistry, Institute of ChemistryEötvös Loránd University Budapest Hungary
| | - Dénes Szabó
- Department of Organic Chemistry, Institute of ChemistryEötvös Loránd University Budapest Hungary
| | - József Rábai
- Department of Organic Chemistry, Institute of ChemistryEötvös Loránd University Budapest Hungary
| | - István Jalsovszky
- Department of Organic Chemistry, Institute of ChemistryEötvös Loránd University Budapest Hungary
| | - Ödön Farkas
- Department of Organic Chemistry, Institute of ChemistryEötvös Loránd University Budapest Hungary
| |
Collapse
|
24
|
Fugel M, Ponomarenko MV, Hesse MF, Malaspina LA, Kleemiss F, Sugimoto K, Genoni A, Röschenthaler GV, Grabowsky S. Complementary bonding analysis of the N–Si interaction in pentacoordinated silicon compounds using quantum crystallography. Dalton Trans 2019; 48:16330-16339. [DOI: 10.1039/c9dt02772f] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The unique combination of quantum crystallography and complementary bonding analysis is used to investigate the bonding of pentacoordinated silicon atoms.
Collapse
Affiliation(s)
- Malte Fugel
- University of Bremen
- Department 2: Biology/Chemistry
- Institute of Inorganic Chemistry and Crystallography
- 28359 Bremen
- Germany
| | | | - Maxie F. Hesse
- University of Bremen
- Department 2: Biology/Chemistry
- Institute of Inorganic Chemistry and Crystallography
- 28359 Bremen
- Germany
| | - Lorraine A. Malaspina
- University of Bremen
- Department 2: Biology/Chemistry
- Institute of Inorganic Chemistry and Crystallography
- 28359 Bremen
- Germany
| | - Florian Kleemiss
- University of Bremen
- Department 2: Biology/Chemistry
- Institute of Inorganic Chemistry and Crystallography
- 28359 Bremen
- Germany
| | - Kunihisa Sugimoto
- JASRI/Diffraction & Scattering Division
- Hyogo 679-5198
- Japan
- Institute for Integrated Cell-Material Sciences (iCeMS)
- Kyoto University
| | | | | | - Simon Grabowsky
- University of Bremen
- Department 2: Biology/Chemistry
- Institute of Inorganic Chemistry and Crystallography
- 28359 Bremen
- Germany
| |
Collapse
|
25
|
Cypryk M, Gostyński B, Pokora M. Hydrolysis of trialkoxysilanes catalysed by the fluoride anion. Nucleophilic vs. basic catalysis. NEW J CHEM 2019. [DOI: 10.1039/c9nj03719e] [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/21/2022]
Abstract
A DFT study showed competition between basic and nucleophilic catalysis of alkoxysilane hydrolysis by the fluoride ion.
Collapse
Affiliation(s)
- Marek Cypryk
- Centre of Molecular and Macromolecular Studies
- Polish Academy of Sciences
- Lodz 90-363
- Poland
| | - Bartłomiej Gostyński
- Centre of Molecular and Macromolecular Studies
- Polish Academy of Sciences
- Lodz 90-363
- Poland
| | | |
Collapse
|
26
|
Hamlin TA, Swart M, Bickelhaupt FM. Nucleophilic Substitution (S N 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent. Chemphyschem 2018; 19:1315-1330. [PMID: 29542853 PMCID: PMC6001448 DOI: 10.1002/cphc.201701363] [Citation(s) in RCA: 118] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Indexed: 11/12/2022]
Abstract
The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.
Collapse
Affiliation(s)
- Trevor A. Hamlin
- Department of Theoretical Chemistry andAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
| | - Marcel Swart
- Department of Theoretical Chemistry andAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
- Institut de Química Computacional I Catàlisi and Department de QuímicaUniversitat de Girona17003GironaSpain
- ICREAPg. Lluís Companys 2308010BarcelonaSpain
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry andAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
- Institute of Molecules and Materials (IMM)Radboud UniversityHeyendaalseweg 1356525 AJNijmegenThe Netherlands
| |
Collapse
|
27
|
Hamlin TA, van Beek B, Wolters LP, Bickelhaupt FM. Nucleophilic Substitution in Solution: Activation Strain Analysis of Weak and Strong Solvent Effects. Chemistry 2018; 24:5927-5938. [PMID: 29457865 PMCID: PMC5947303 DOI: 10.1002/chem.201706075] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Indexed: 12/23/2022]
Abstract
We have quantum chemically studied the effect of various polar and apolar solvents on the shape of the potential energy surface (PES) of a diverse collection of archetypal nucleophilic substitution reactions at carbon, silicon, phosphorus, and arsenic by using density functional theory at the OLYP/TZ2P level. In the gas phase, all our model SN 2 reactions have single-well PESs, except for the nucleophilic substitution reaction at carbon (SN 2@C), which has a double-well energy profile. The presence of the solvent can have a significant effect on the shape of the PES and, thus, on the nature of the SN 2 process. Solvation energies, charges on the nucleophile or leaving group, and structural features are compared for the various SN 2 reactions in a spectrum of solvents. We demonstrate how solvation can change the shape of the PES, depending not only on the polarity of the solvent, but also on how the charge is distributed over the interacting molecular moieties during different stages of the reaction. In the case of a nucleophilic substitution at three-coordinate phosphorus, the reaction can be made to proceed through a single-well [no transition state (TS)], bimodal barrier (two TSs), and then through a unimodal transition state (one TS) simply by increasing the polarity of the solvent.
Collapse
Affiliation(s)
- Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
| | - Bas van Beek
- Department of Theoretical ChemistryAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
| | - Lando P. Wolters
- Department of Theoretical ChemistryAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical ChemistryAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
- Institute of Molecules and Materials (IMM)Radboud University NijmegenHeyendaalseweg 1356525 AJNijmegenThe Netherlands
| |
Collapse
|
28
|
Laloo JZA, Rhyman L, Larrañaga O, Ramasami P, Bickelhaupt FM, de Cózar A. Ion-Pair S N 2 Reaction of OH - and CH 3 Cl: Activation Strain Analyses of Counterion and Solvent Effects. Chem Asian J 2018; 13:1138-1147. [PMID: 29437289 DOI: 10.1002/asia.201800082] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2018] [Revised: 02/12/2018] [Indexed: 11/10/2022]
Abstract
We have theoretically studied the non-identity SN 2 reactions of Mn OH(n-1) +CH3 Cl (M+ =Li+ , Na+ , K+ , and MgCl+ ; n=0, 1) in the gas phase and in THF solution at the OLYP/6-31++G(d,p) level using polarizable continuum model (PCM) implicit solvation. We want to explore and understand the effect of the metal counterion M+ and solvation on the reaction profile and the stereoselectivity of these processes. To this end, we have explored the potential energy surfaces of the backside (SN 2-b) and frontside (SN 2-f) pathways. To explain the computed trends, we have carried out analyses with an extended activation strain model (ASM) of chemical reactivity that includes the treatment of solvation effects.
Collapse
Affiliation(s)
- Jalal Z A Laloo
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius
| | - Lydia Rhyman
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius.,Department of Applied Chemistry, University of Johannesburg, Doornfontein Campus, Johannesburg, 2028, South Africa
| | - Olatz Larrañaga
- Departamento de Química Orgánica I, Facultad de Química, Universidad del País Vasco (UPV/EHU) and Donostia International Physics Center (DIPC), P. K. 1072, 20018, San Sebastián-Donostia, Spain
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius.,Department of Applied Chemistry, University of Johannesburg, Doornfontein Campus, Johannesburg, 2028, South Africa
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081, HV, Amsterdam, The Netherlands.,Institute of Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525, AJ, Nijmegen, The Netherlands
| | - Abel de Cózar
- Departamento de Química Orgánica I, Facultad de Química, Universidad del País Vasco (UPV/EHU) and Donostia International Physics Center (DIPC), P. K. 1072, 20018, San Sebastián-Donostia, Spain.,Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081, HV, Amsterdam, The Netherlands.,IKERBASQUE, Basque Foundation for Science, 48013, Bilbao, Spain
| |
Collapse
|
29
|
Proud AJ, Sheppard BJH, Pearson JK. Revealing Electron-Electron Interactions within Lewis Pairs in Chemical Systems. J Am Chem Soc 2018; 140:219-228. [PMID: 29228767 DOI: 10.1021/jacs.7b08935] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The so-called "Lewis pair" is a ubiquitous phenomenon in chemistry and is often used as an intuitive construct to predict and rationalize chemical structure and behavior. Concepts from the very general Valence Shell Electron Pair Repulsion (VSEPR) model to the most esoteric reaction mechanism routinely rely on the notion that electrons tend to exist in pairs and that these pairs can be thought of as being localized to a particular region of space. It is precisely this localization that allows one to intuit how these pairs might behave, generally speaking, so that reasonable predictions may be made regarding molecular structure, intermolecular interactions, property trends, and reaction mechanisms, etc. Of course, it is rather unfortunate that the Lewis model is entirely qualitative and yields no information regarding how any specific electron pair is distributed. Here we demonstrate a novel electronic structure analysis technique that predicts and analyzes precise quantitative details about the relative and absolute distribution of individual electron pairs. This Single Electron Pair Distribution Analysis (SEPDA) reveals quantitative details about the distribution of the well-known Lewis pairs, such as how they are distributed in space and how their relative velocities change in various chemical contexts. We show that these distributions allow one to image the explicitly pairwise electronic behavior of bonds and lone pairs. We further demonstrate how this electronic behavior changes with several conditions to explore the nature of the covalent chemical bond, non-covalent interactions, bond formation, and exotic 3-center-2-electron species. It is shown that indications of the strength of bonded and non-bonded interactions may also be gleaned from such distributions and SEPDA can be used as a tool to differentiate between interaction types. We anticipate that SEPDA will be of broad utility in a wide variety of chemical contexts because it affords a very detailed, visual and intuitive analysis technique that is generally applicable.
Collapse
Affiliation(s)
- Adam Jonathan Proud
- Department of Chemistry, University of Prince Edward Island , Charlottetown, PE C1A 4P3, Canada
| | | | - Jason Kenneth Pearson
- Department of Chemistry, University of Prince Edward Island , Charlottetown, PE C1A 4P3, Canada
| |
Collapse
|
30
|
Kour M, Gupta R, Bansal RK. Experimental and Theoretical Investigation of the Reaction of Amines with Methyl Propiolate. ChemistrySelect 2017. [DOI: 10.1002/slct.201701103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Manjinder Kour
- Department of chemistry the IIS University; Jaipur 302020 India
| | - Raakhi Gupta
- Department of chemistry the IIS University; Jaipur 302020 India
| | - Raj K. Bansal
- Department of chemistry the IIS University; Jaipur 302020 India
| |
Collapse
|
31
|
Hupf E, Olaru M, Raţ CI, Fugel M, Hübschle CB, Lork E, Grabowsky S, Mebs S, Beckmann J. Mapping the Trajectory of Nucleophilic Substitution at Silicon Using aperi-Substituted Acenaphthyl Scaffold. Chemistry 2017; 23:10568-10579. [DOI: 10.1002/chem.201700992] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2017] [Indexed: 11/07/2022]
Affiliation(s)
- Emanuel Hupf
- Institut für Anorganische Chemie und Kristallographie; Fachbereich 2-Biologie/Chemie; Universität Bremen; Leobener Straße NW2 28359 Bremen Germany
- Current address: Department of Chemistry; University of Alberta; Edmonton Alberta T6G 2G2 Canada
| | - Marian Olaru
- Institut für Anorganische Chemie und Kristallographie; Fachbereich 2-Biologie/Chemie; Universität Bremen; Leobener Straße NW2 28359 Bremen Germany
- Centre of Supramolecular Organic and Organometallic Chemistry; Department of Chemistry, Faculty of Chemistry and Chemical Engineering; Babes-Bolyai University; 11 Arany Janos Street 400028 Cluj-Napoca Romania
| | - Ciprian I. Raţ
- Centre of Supramolecular Organic and Organometallic Chemistry; Department of Chemistry, Faculty of Chemistry and Chemical Engineering; Babes-Bolyai University; 11 Arany Janos Street 400028 Cluj-Napoca Romania
| | - Malte Fugel
- Institut für Anorganische Chemie und Kristallographie; Fachbereich 2-Biologie/Chemie; Universität Bremen; Leobener Straße NW2 28359 Bremen Germany
| | | | - Enno Lork
- Institut für Anorganische Chemie und Kristallographie; Fachbereich 2-Biologie/Chemie; Universität Bremen; Leobener Straße NW2 28359 Bremen Germany
| | - Simon Grabowsky
- Institut für Anorganische Chemie und Kristallographie; Fachbereich 2-Biologie/Chemie; Universität Bremen; Leobener Straße NW2 28359 Bremen Germany
| | - Stefan Mebs
- Institut für Experimentalphysik; Freie Universität Berlin; Arnimallee 14 14195 Berlin Germany
| | - Jens Beckmann
- Institut für Anorganische Chemie und Kristallographie; Fachbereich 2-Biologie/Chemie; Universität Bremen; Leobener Straße NW2 28359 Bremen Germany
| |
Collapse
|
32
|
Microsolvation effects on the reactivity of oxy-nucleophiles: the case of gas-phase S N2 reactions of YO -(CH 3OH) n=1,2 towards CH 3Cl. J Mol Model 2017; 23:192. [PMID: 28528446 DOI: 10.1007/s00894-017-3351-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2017] [Accepted: 04/19/2017] [Indexed: 10/19/2022]
Abstract
The modified G4(MP2) method was applied to explore microsolvation effects on the reactivity of four solvated normal oxy-nucleophiles YO-(CH3OH) n=1,2 (Y = CH3, C2H5, FC2H4, ClC2H4), and five α-oxy-nucleophiles YO-(CH3OH) n=1,2 (Y = HO, CH3O, F, Cl, Br), in gas-phase SN2 reactions towards the substrate CH3Cl. Based on a Brønsted-type plot, our calculations reveal that the overall activation barriers of five microsolvated α-oxy-nucleophiles are obviously smaller than the prediction from the correlation line constructed by four normal microsolvated ones to different degrees, and clearly demonstrate the existence of an α-effect in the presence of one or two methanol molecule(s). Moreover, it was found that the α-effect of the mono-methanol microsolvated α-nucleophile is stronger than that of the monohydrated α-nucleophile. However, the α-effect of YO-(CH3OH)2 becomes weaker for Y = HO and CH3O, whereas it becomes stronger for Y = F, Cl, Br than that of YO-(H2O)2, which can be explained by analyses of the activation strain model in the two cases. It was also found that the rationale about the low ionization energy of α-nucleophile inducing the α-effect was not widely significant. Graphical abstract Variation of alpha-effect in the gas-phase SN2 reaction with the microsolvation.
Collapse
|
33
|
Ribaudo G, Bellanda M, Menegazzo I, Wolters LP, Bortoli M, Ferrer-Sueta G, Zagotto G, Orian L. Mechanistic Insight into the Oxidation of Organic Phenylselenides by H 2 O 2. Chemistry 2017; 23:2405-2422. [PMID: 27935210 DOI: 10.1002/chem.201604915] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Indexed: 12/24/2022]
Abstract
The oxidation of organic phenylselenides by H2 O2 is investigated in model compounds, namely, n-butyl phenyl selenide (PhSe(nBu)), bis(phenylselanyl)methane (PhSeMeSePh), diphenyl diselenide (PhSeSePh), and 1,2-bis(phenylselanyl)ethane (PhSeEtSePh). Through a combined experimental (1 H and 77 Se NMR) and computational approach, we characterize the direct oxidation of monoselenide to selenoxide, the stepwise double oxidation of PhSeMeSePh that leads to different diastereomeric diselenoxides, the complete oxidation of the diphenyldiselenide that leads to selenium-selenium bond cleavage, and the subsequent formation of the phenylseleninic product. The oxidation of PhSeEtSePh also results in the formation of phenylseleninic acid along with 1-(vinylseleninyl)benzene, which is derived from a side elimination reaction. The evidence of a direct mechanism, in addition to an autocatalytic mechanism that emerges from kinetic studies, is discussed. By considering our observations of diselenides with chalcogen atoms that are separated by alkyl spacers of different length, a rationale for the advantage of diselenide versus monoselenide catalysts is presented.
Collapse
Affiliation(s)
- Giovanni Ribaudo
- Dipartimento di Scienze del Farmaco, Università degli Studi di Padova, Via Marzolo 5, 35131, Padova, Italy
| | - Massimo Bellanda
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131, Padova, Italy
| | - Ileana Menegazzo
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131, Padova, Italy
| | - Lando P Wolters
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131, Padova, Italy
| | - Marco Bortoli
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131, Padova, Italy
| | - Gerardo Ferrer-Sueta
- Laboratorio de Fisicoquímica Biológica, Instituto de Química Biológica, Facultad de Ciencias, Universidad de la Repúbica, Igua 4225, Montevideo, Uruguay
| | - Giuseppe Zagotto
- Dipartimento di Scienze del Farmaco, Università degli Studi di Padova, Via Marzolo 5, 35131, Padova, Italy
| | - Laura Orian
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131, Padova, Italy
| |
Collapse
|
34
|
Abstract
![]()
Fundamental
principles that determine chemical reactivity and reaction
mechanisms are the very foundation of chemistry and many related fields
of science. Bimolecular nucleophilic substitutions (SN2)
are among the most common and therefore most important reaction types.
In this report, we examine the trends in the SN2 reactions
with respect to increasing electronegativity of the reaction center
by comparing the well-studied backside SN2 Cl– + CH3Cl with similar Cl– substitutions
on the isoelectronic series with the second period elements N, O,
and F in place of C. Relativistic (ZORA) DFT calculations are used
to construct the gas phase reaction potential energy surfaces (PES),
and activation strain analysis, which allows decomposition of the
PES into the geometrical strain and interaction energy, is employed
to analyze the observed trends. We find that SN2@N and
SN2@O have similar PES to the prototypical SN2@C, with the well-defined reaction complex (RC) local minima and
a central barrier, but all stationary points are, respectively, increasingly
stable in energy. The SN2@F, by contrast, exhibits only
a single-well PES with no barrier. Using the activation strain model,
we show that the trends are due to the interaction energy and originate
mainly from the decreasing energy of the empty acceptor orbital (σ*A–Cl) on the reaction center A in the order of C, N,
O, and F. The decreasing steric congestion around the central atom
is also a likely contributor to this trend. Additional decomposition
of the interaction energy using Kohn–Sham molecular orbital
(KS-MO) theory provides further support for this explanation, as well
as suggesting electrostatic energy as the primary reason for the distinct
single-well PES profile for the FCl reaction.
Collapse
Affiliation(s)
- Jan Kubelka
- Department of Chemistry, University of Wyoming , Laramie, Wyoming 82070, United States
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands.,Institute for Molecules and Materials (IMM), Radboud University , Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| |
Collapse
|
35
|
Morera-Boado C, Martínez González M, Miranda-Quintana RA, Suárez M, Martínez-Álvarez R, Martín N, García de la Vega JM. Theoretical Study on the Mechanism of the Thermal Retro-Cycloaddition of Isoxazolinofullerenes. J Phys Chem A 2016; 120:8830-8842. [DOI: 10.1021/acs.jpca.6b07297] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Cercis Morera-Boado
- Laboratorio
de Química, Computacional y Teórica, Facultad de Química, Universidad de la Habana, 10400 La Habana, Cuba
| | - Marco Martínez González
- Laboratorio
de Química, Computacional y Teórica, Facultad de Química, Universidad de la Habana, 10400 La Habana, Cuba
| | - Ramón A. Miranda-Quintana
- Laboratorio
de Química, Computacional y Teórica, Facultad de Química, Universidad de la Habana, 10400 La Habana, Cuba
| | - Margarita Suárez
- Laboratorio
de Síntesis Orgánica, Facultad de Química, Universidad de la Habana, 10400 La Habana, Cuba
| | - Roberto Martínez-Álvarez
- Departamento
de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Nazario Martín
- Departamento
de Química Orgánica I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - José. M. García de la Vega
- Departamento
de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|
36
|
de Paul N. Nziko V, Scheiner S. Effects of Angular Deformation on the Energetics of the S N2 Reaction. European J Org Chem 2016. [DOI: 10.1002/ejoc.201600712] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
| | - Steve Scheiner
- Department of Chemistry and Biochemistry; Utah State University; 84322-0300 Logan UT USA
| |
Collapse
|
37
|
|
38
|
Ehbets J, Lorenzen S, Mahler C, Bertermann R, Berkefeld A, Poater J, Fritz‐Langhals E, Weidner R, Bickelhaupt FM, Tacke R. Synthesis and Hydrolysis of Alkoxy(aminoalkyl)diorganylsilanes of the Formula Type R
2
(RO)Si(CH
2
)
n
NH
2
(R = Alkyl,
n
= 1–3): A Systematic Experimental and Computational Study. Eur J Inorg Chem 2016. [DOI: 10.1002/ejic.201600077] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Julia Ehbets
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg, Germany, http://www‐anorganik.chemie.uni‐wuerzburg.de
| | - Sabine Lorenzen
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg, Germany, http://www‐anorganik.chemie.uni‐wuerzburg.de
| | - Christoph Mahler
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg, Germany, http://www‐anorganik.chemie.uni‐wuerzburg.de
| | - Rüdiger Bertermann
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg, Germany, http://www‐anorganik.chemie.uni‐wuerzburg.de
| | - André Berkefeld
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg, Germany, http://www‐anorganik.chemie.uni‐wuerzburg.de
| | - Jordi Poater
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands, http://www.few.vu.nl/~bickel
- Universitat de Barcelona, Departament de Química Orgànica & Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Catalonia, Spain
| | - Elke Fritz‐Langhals
- Wacker Chemie AG, Consortium für elektrochemische Industrie, Zielstattstrasse 20, 81379 München, Germany
| | - Richard Weidner
- Wacker Chemie AG, Consortium für elektrochemische Industrie, Zielstattstrasse 20, 81379 München, Germany
| | - F. Matthias Bickelhaupt
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands, http://www.few.vu.nl/~bickel
- Radboud University Nijmegen, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Reinhold Tacke
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg, Germany, http://www‐anorganik.chemie.uni‐wuerzburg.de
| |
Collapse
|
39
|
Laloo JZA, Rhyman L, Ramasami P, Bickelhaupt FM, de Cózar A. Ion-Pair SN 2 Substitution: Activation Strain Analyses of Counter-Ion and Solvent Effects. Chemistry 2016; 22:4431-9. [PMID: 26879231 DOI: 10.1002/chem.201504456] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2015] [Indexed: 12/20/2022]
Abstract
The ion-pair SN 2 reactions of model systems MnF(n-1) +CH3Cl(M(+) =Li(+), Na(+), K(+), and MgCl(+); n=0, 1) have been quantum chemically explored by using DFT at the OLYP/6-31++G(d,p) level. The purpose of this study is threefold: 1) to elucidate how the counterion M(+) modifies ion-pair SN 2 reactivity relative to the parent reaction F(-) +CH3Cl; 2) to determine how this influences stereochemical competition between the backside and frontside attacks; and 3) to examine the effect of solvation on these ion-pair SN2 pathways. Trends in reactivity are analyzed and explained by using the activation strain model (ASM) of chemical reactivity. The ASM has been extended to treat reactivity in solution. These findings contribute to a more rational design of tailor-made substitution reactions.
Collapse
Affiliation(s)
- Jalal Z A Laloo
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius
| | - Lydia Rhyman
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837, Mauritius. .,Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, P.O. Box 2457, Riyadh, 11451, Saudi Arabia.
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands. .,Institute of Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg 135, 6525 AJ, Nijmegen, The Netherlands.
| | - Abel de Cózar
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands. .,Departamento de Química Orgánica I, Facultad de Química, Universidad del País Vasco P. K. 1072, 200880, San Sebastián-Donostia, Spain. .,IKERBASQUE, Basque Foundation for Science, 48013, Bilbao, Spain.
| |
Collapse
|
40
|
Bento AP, Solà M, Bickelhaupt FM. E2 and SN2 Reactions of X(-) + CH3CH2X (X = F, Cl); an ab Initio and DFT Benchmark Study. J Chem Theory Comput 2015; 4:929-40. [PMID: 26621234 DOI: 10.1021/ct700318e] [Citation(s) in RCA: 81] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We have computed consistent benchmark potential energy surfaces (PESs) for the anti-E2, syn-E2, and SN2 pathways of X(-) + CH3CH2X with X = F and Cl. This benchmark has been used to evaluate the performance of 31 popular density functionals, covering local-density approximation, generalized gradient approximation (GGA), meta-GGA, and hybrid density-functional theory (DFT). The ab initio benchmark has been obtained by exploring the PESs using a hierarchical series of ab initio methods [up to CCSD(T)] in combination with a hierarchical series of Gaussian-type basis sets (up to aug-cc-pVQZ). Our best CCSD(T) estimates show that the overall barriers for the various pathways increase in the order anti-E2 (X = F) < SN2 (X = F) < SN2 (X = Cl) ∼ syn-E2 (X = F) < anti-E2 (X = Cl) < syn-E2 (X = Cl). Thus, anti-E2 dominates for F(-) + CH3CH2F, and SN2 dominates for Cl(-) + CH3CH2Cl, while syn-E2 is in all cases the least favorable pathway. Best overall agreement with our ab initio benchmark is obtained by representatives from each of the three categories of functionals, GGA, meta-GGA, and hybrid DFT, with mean absolute errors in, for example, central barriers of 4.3 (OPBE), 2.2 (M06-L), and 2.0 kcal/mol (M06), respectively. Importantly, the hybrid functional BHandH and the meta-GGA M06-L yield incorrect trends and qualitative features of the PESs (in particular, an erroneous preference for SN2 over the anti-E2 in the case of F(-) + CH3CH2F) even though they are among the best functionals as measured by their small mean absolute errors of 3.3 and 2.2 kcal/mol in reaction barriers. OLYP and B3LYP have somewhat higher mean absolute errors in central barriers (5.6 and 4.8 kcal/mol, respectively), but the error distribution is somewhat more uniform, and as a consequence, the correct trends are reproduced.
Collapse
Affiliation(s)
- A Patrícia Bento
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands, and Institut de Química Computacional, Universitat de Girona, Campus Montilivi, E-17071 Girona, Catalonia, Spain
| | - Miquel Solà
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands, and Institut de Química Computacional, Universitat de Girona, Campus Montilivi, E-17071 Girona, Catalonia, Spain
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands, and Institut de Química Computacional, Universitat de Girona, Campus Montilivi, E-17071 Girona, Catalonia, Spain
| |
Collapse
|
41
|
Mück FM, Baus JA, Nutz M, Burschka C, Poater J, Bickelhaupt FM, Tacke R. Reactivity of the Donor‐Stabilized Silylenes [
i
PrNC(Ph)N
i
Pr]
2
Si and [
i
PrNC(N
i
Pr
2
)N
i
Pr]
2
Si: Activation of CO
2
and CS
2. Chemistry 2015; 21:16665-72. [DOI: 10.1002/chem.201501788] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2015] [Indexed: 11/11/2022]
Affiliation(s)
- Felix M. Mück
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany), Fax: (+49) 931‐31‐84609
| | - Johannes A. Baus
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany), Fax: (+49) 931‐31‐84609
| | - Marco Nutz
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany), Fax: (+49) 931‐31‐84609
| | - Christian Burschka
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany), Fax: (+49) 931‐31‐84609
| | - Jordi Poater
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands)
| | - F. Matthias Bickelhaupt
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands)
- Radboud University Nijmegen, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen (The Netherlands)
| | - Reinhold Tacke
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany), Fax: (+49) 931‐31‐84609
| |
Collapse
|
42
|
Mück FM, Kloß D, Baus JA, Burschka C, Bertermann R, Poater J, Fonseca Guerra C, Bickelhaupt FM, Tacke R. Stable Four‐Coordinate Guanidinatosilicon(IV) Complexes with SiN
3
El Skeletons (El=S, Se, Te) and SiEl Double Bonds. Chemistry 2015; 21:14011-21. [DOI: 10.1002/chem.201501789] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2015] [Indexed: 11/11/2022]
Affiliation(s)
- Felix M. Mück
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany)
| | - Dorit Kloß
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany)
| | - Johannes A. Baus
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany)
| | - Christian Burschka
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany)
| | - Rüdiger Bertermann
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany)
| | - Jordi Poater
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands)
| | - Célia Fonseca Guerra
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands)
| | - F. Matthias Bickelhaupt
- VU University Amsterdam, Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), De Boelelaan 1083, 1081 HV Amsterdam (The Netherlands)
- Radboud University Nijmegen, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen (The Netherlands)
| | - Reinhold Tacke
- Universität Würzburg, Institut für Anorganische Chemie, Am Hubland, 97074 Würzburg (Germany)
| |
Collapse
|
43
|
Pill MF, Schmidt SW, Beyer MK, Clausen-Schaumann H, Kersch A. A density functional theory model of mechanically activated silyl ester hydrolysis. J Chem Phys 2015; 140:044321. [PMID: 25669537 DOI: 10.1063/1.4862827] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
To elucidate the mechanism of the mechanically activated dissociation of chemical bonds between carboxymethylated amylose (CMA) and silane functionalized silicon dioxide, we have investigated the dissociation kinetics of the bonds connecting CMA to silicon oxide surfaces with density functional calculations including the effects of force, solvent polarizability, and pH. We have determined the activation energies, the pre-exponential factors, and the reaction rate constants of candidate reactions. The weakest bond was found to be the silyl ester bond between the silicon and the alkoxy oxygen atom. Under acidic conditions, spontaneous proton addition occurs close to the silyl ester such that neutral reactions become insignificant. Upon proton addition at the most favored position, the activation energy for bond hydrolysis becomes 31 kJ mol(-1), which agrees very well with experimental observation. Heterolytic bond scission in the protonated molecule has a much higher activation energy. The experimentally observed bi-exponential rupture kinetics can be explained by different side groups attached to the silicon atom of the silyl ester. The fact that different side groups lead to different dissociation kinetics provides an opportunity to deliberately modify and tune the kinetic parameters of mechanically activated bond dissociation of silyl esters.
Collapse
Affiliation(s)
- Michael F Pill
- Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences, Lothstr. 34, 80335 Munich, Germany
| | - Sebastian W Schmidt
- Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences, Lothstr. 34, 80335 Munich, Germany
| | - Martin K Beyer
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Hauke Clausen-Schaumann
- Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences, Lothstr. 34, 80335 Munich, Germany
| | - Alfred Kersch
- Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences, Lothstr. 34, 80335 Munich, Germany
| |
Collapse
|
44
|
Proenza YG, de Souza MAF, Ventura E, do Monte SA, Longo RL. Assessment of density-functionals for describing the X(-) + CH3ONO2 gas-phase reactions with X = F, OH, CH2CN. Phys Chem Chem Phys 2015; 16:26769-78. [PMID: 25372179 DOI: 10.1039/c4cp03674c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The energetics of the ECO2, SN2@C and SN2@N channels of X(-) + CH3ONO2 (X = F, OH, CH2CN) gas-phase reactions were computed using the CCSD(T)/CBS method. This benchmark extends a previous study with X = OH [M. A. F. de Souza et al., J. Am. Chem. Soc., 2012, 134, 19004] and was used to ascertain the accuracy and robustness of nineteen density-functionals for describing these potential energy profiles (PEP) as well as the kinetic product distributions obtained from RRKM calculations. Assessments were based on the mean unsigned error (MUE), the mean signed error (MSE), the #best : #worst (BW) criterion and the statistical confidence interval (CI) for the MSE. In general, double-hybrid (DH) functionals perform better than the range-separated ones, and both are better than the global-hybrid functionals. Based on the MUE and CI criteria the B2GPPLYP, B2PLYP, M08-SO, BMK, ωB97X-D, CAM-B3LYP, M06, M08-HX, ωB97X and B97-K functionals show the best performance in the description of these PEPs. Within this set, the B2GPPLYP functional is the most accurate and robust. The RRKM results indicate that the DHs are the best for describing the selectivities of these reactions. Compared to CCSD(T), the B2PLYP method has a relative error of only ca. 1% for the selectivity and the accuracy to provide the correct conclusion concerning the nonstatistical behavior of these reactions.
Collapse
Affiliation(s)
- Yaicel G Proenza
- Departamento de Química Fundamental, CCEN, Universidade Federal de Pernambuco, Cidade Universitária, Recife, PE 50.740-560, Brazil.
| | | | | | | | | |
Collapse
|
45
|
Liang J, Su Q, Wang Y, Geng Z. SN2 Reaction of IO− + CH3Cl: An Ab Initio and DFT Benchmark Study. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2015. [DOI: 10.1246/bcsj.20140194] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Junxi Liang
- Gansu Key Laboratory of Environmental Friendly Composites and Biomass Utilization, College of Chemical Engineering, Northwest University for Nationalities
| | - Qiong Su
- Gansu Key Laboratory of Environmental Friendly Composites and Biomass Utilization, College of Chemical Engineering, Northwest University for Nationalities
| | - Yanbin Wang
- Gansu Key Laboratory of Environmental Friendly Composites and Biomass Utilization, College of Chemical Engineering, Northwest University for Nationalities
| | - Zhiyuan Geng
- Gansu Key Laboratory of Polymer Materials, College of Chemistry and Chemical Engineering, Key Laboratory of Eco-environment-related Polymer Materials, Ministry of Education, Northwest Normal University
| |
Collapse
|
46
|
Essafi S, Tomasi S, Aggarwal VK, Harvey JN. Homologation of boronic esters with organolithium compounds: a computational assessment of mechanism. J Org Chem 2014; 79:12148-58. [PMID: 25318004 DOI: 10.1021/jo502020e] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Ab initio calculations are reported for the reaction of methyl boronic ester with organolithium reagents with α-leaving groups. The best calculations rely on density functional theory prediction of structures and coupled-cluster theory calculation of accurate potential energies. The results provide strong confirmation of the feasibility of a two-step mechanism with rapid initial formation of a boron-ate complex followed by slower migration of methyl from boron to carbon with loss of the leaving group. The calculated free energy of activation is consistent with observed kinetic behavior, and the calculations provide a framework for exploring substituent and other effects on reactivity. Obtaining reasonable agreement with experiment in this way is not trivial and requires careful treatment of level of theory (density functional theory calculations tend to yield inaccurate results), of conformational complexity, especially for the ate complexes, and of the nature of the microscopic model of reactants and solvent. The methodological challenges and possible pitfalls, many of which are relevant more broadly to computational modeling of organic reaction mechanisms, are discussed in detail.
Collapse
Affiliation(s)
- Stéphanie Essafi
- School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, U.K
| | | | | | | |
Collapse
|
47
|
Shekar S, Brown SN. Mechanism and selectivity of methyl and phenyl migrations in hypervalent silylated iminoquinones. J Org Chem 2014; 79:12047-55. [PMID: 25289545 DOI: 10.1021/jo501888r] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Chlorosilanes R(X)(Y)SiCl (R = Me, Ph; X, Y = Me, Ph, Cl) have been reported to react with Pb(ONO(Q))2 (ONO(Q) = 3,5-di-tert-butyl-1,2-quinone-(3,5-di-tert-butyl-2-oxy-1-phenyl)imine) to give five-coordinate (X)(Y)Si(ON[R]O), in which the R group has migrated from silicon to nitrogen. This migration is intramolecular, as confirmed by the lack of crossover between (CH3)3SiCl and (CD3)3SiCl in their reaction with Pb(ONO(Q))2. Reaction of PhSiMeCl2 takes place with high kinetic stereoselectivity to produce isomer Ph(Cl)Si(ON[Me]O) in which the phenyl is axial in the trigonal bipyramid, which subsequently isomerizes to the thermodynamic isomer with axial chlorine. This indicates that migration takes place preferentially from the stereoisomer of the octahedral intermediate, κ(3)-Ph(CH3)(Cl)Si(ONO(Q)), in which the phenyl and methyl groups are mutually trans, indicating that the observed complete selectivity for methyl over phenyl migration is due to intrinsic differences in migratory aptitude. DFT calculations suggest that migration takes place from this isomer not because it undergoes migration faster than other possible stereoisomers, but because it is formed most rapidly, and migration occurs faster than isomerization.
Collapse
Affiliation(s)
- Sukesh Shekar
- Department of Chemistry and Biochemistry, University of Notre Dame , 251 Nieuwland Science Hall, Notre Dame, Indiana 46556-5670, United States
| | | |
Collapse
|
48
|
Tandem Michael addition of amines to maleic anhydride and 1,3-prototropic shift: experimental and theoretical results. Tetrahedron 2014. [DOI: 10.1016/j.tet.2014.06.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
49
|
Junold K, Nutz M, Baus JA, Burschka C, Fonseca Guerra C, Bickelhaupt FM, Tacke R. The Donor-Stabilized Silylene Bis[N,N′-diisopropylbenzamidinato(−)]silicon(II): Synthesis, Electronic Structure, and Reactivity. Chemistry 2014; 20:9319-29. [DOI: 10.1002/chem.201402483] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2014] [Indexed: 11/09/2022]
|
50
|
Berkefeld A, Guerra CF, Bertermann R, Troegel D, Daiß JO, Stohrer J, Bickelhaupt FM, Tacke R. Silicon α-Effect: A Systematic Experimental and Computational Study of the Hydrolysis of Cα- and Cγ-Functionalized Alkoxytriorganylsilanes of the Formula Type ROSiMe2(CH2)nX (R = Me, Et; n = 1, 3; X = Functional Group). Organometallics 2014. [DOI: 10.1021/om500073m] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- André Berkefeld
- Institut
für Anorganische Chemie, Universität Würzburg, Am Hubland, D-97074 Würzburg, Germany
| | - Célia Fonseca Guerra
- Department
of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling
(ACMM), VU University Amsterdam, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Rüdiger Bertermann
- Institut
für Anorganische Chemie, Universität Würzburg, Am Hubland, D-97074 Würzburg, Germany
| | - Dennis Troegel
- Wacker Chemie AG, Consortium für elektrochemische Industrie, Zielstattstrasse 20, D-81379 München, Germany
| | - Jürgen O. Daiß
- Wacker Chemie AG, Johannes-Hess-Strasse
24, D-84489 Burghausen, Germany
| | - Jürgen Stohrer
- Wacker Chemie AG, Consortium für elektrochemische Industrie, Zielstattstrasse 20, D-81379 München, Germany
| | - F. Matthias Bickelhaupt
- Department
of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling
(ACMM), VU University Amsterdam, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
- Institute
for Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg
135, NL-6525 AJ Nijmegen, The Netherlands
| | - Reinhold Tacke
- Institut
für Anorganische Chemie, Universität Würzburg, Am Hubland, D-97074 Würzburg, Germany
| |
Collapse
|