1
|
Durán R, Barrales-Martínez C, Santana-Romo F, Rodríguez DF, Zacconi FC, Herrera B. Substitution Effects in Aryl Halides and Amides into the Reaction Mechanism of Ullmann-Type Coupling Reactions. Molecules 2024; 29:1770. [PMID: 38675590 PMCID: PMC11051942 DOI: 10.3390/molecules29081770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Revised: 04/03/2024] [Accepted: 04/10/2024] [Indexed: 04/28/2024] Open
Abstract
In this article, we present a comprehensive computational investigation into the reaction mechanism of N-arylation of substituted aryl halides through Ullmann-type coupling reactions. Our computational findings, obtained through DFT ωB97X-D/6-311G(d,p) and ωB97X-D/LanL2DZ calculations, reveal a direct relation between the previously reported experimental reaction yields and the activation energy of haloarene activation, which constitutes the rate-limiting step in the overall coupling process. A detailed analysis of the reaction mechanism employing the Activation Strain Model indicates that the strain in the substituted iodoanilines is the primary contributor to the energy barrier, representing an average of 80% of the total strain energy. Additional analysis based on conceptual Density Functional Theory (DFT) suggests that the nucleophilicity of the nitrogen in the lactam is directly linked to the activation energies. These results provide valuable insights into the factors influencing energetic barriers and, consequently, reaction yields. These insights enable the rational modification of reactants to optimize the N-arylation process.
Collapse
Affiliation(s)
- Rocío Durán
- Instituto de Investigación Interdisciplinaria (I3), Vicerrectoría Académica, Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Lircay, Talca 3460000, Chile; (R.D.); (C.B.-M.)
| | - César Barrales-Martínez
- Instituto de Investigación Interdisciplinaria (I3), Vicerrectoría Académica, Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Lircay, Talca 3460000, Chile; (R.D.); (C.B.-M.)
| | - Fabián Santana-Romo
- Escuela de Química, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile; (F.S.-R.); (D.F.R.)
- Departamento de Ciencias Exactas, Universidad de las Fuerzas Armadas ESPE, Sangolquí 171103, Ecuador
| | - Diego F. Rodríguez
- Escuela de Química, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile; (F.S.-R.); (D.F.R.)
| | - Flavia C. Zacconi
- Escuela de Química, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile; (F.S.-R.); (D.F.R.)
- Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago 7820436, Chile
| | - Barbara Herrera
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile
| |
Collapse
|
2
|
Sun X, Hansen T, Poater J, Hamlin TA, Bickelhaupt FM. Rational design of iron catalysts for C-X bond activation. J Comput Chem 2023; 44:495-505. [PMID: 35137432 PMCID: PMC10078697 DOI: 10.1002/jcc.26818] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Revised: 01/16/2022] [Accepted: 01/20/2022] [Indexed: 01/10/2023]
Abstract
We have quantum chemically studied the iron-mediated CX bond activation (X = H, Cl, CH3 ) by d8 -FeL4 complexes using relativistic density functional theory at ZORA-OPBE/TZ2P. We find that by either modulating the electronic effects of a generic iron-catalyst by a set of ligands, that is, CO, BF, PH3 , BN(CH3 )2 , or by manipulating structural effects through the introduction of bidentate ligands, that is, PH2 (CH2 )n PH2 with n = 6-1, one can significantly decrease the reaction barrier for the CX bond activation. The combination of both tuning handles causes a decrease of the CH activation barrier from 10.4 to 4.6 kcal mol-1 . Our activation strain and Kohn-Sham molecular orbital analyses reveal that the electronic tuning works via optimizing the catalyst-substrate interaction by introducing a strong second backdonation interaction (i.e., "ligand-assisted" interaction), while the mechanism for structural tuning is mainly caused by the reduction of the required activation strain because of the pre-distortion of the catalyst. In all, we present design principles for iron-based catalysts that mimic the favorable behavior of their well-known palladium analogs in the bond-activation step of cross-coupling reactions.
Collapse
Affiliation(s)
- Xiaobo Sun
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Barcelona, Spain
| | - Thomas Hansen
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Barcelona, Spain.,Leiden Institute of Chemistry, Leiden University, Leiden, The Netherlands
| | - Jordi Poater
- Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Barcelona, Spain.,ICREA, Barcelona, Spain
| | - Trevor A Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Friedrich Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Institute for Molecules and Materials (IMM), Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
3
|
Moloto BP, Vermeeren P, Tiezza MD, Bouwens T, Esterhuysen C, Hamlin TA, Bickelhaupt FM. Palladium-catalyzed activation of H nA–AH n bonds (AH n = CH 3, NH 2, OH, F). PURE APPL CHEM 2023. [DOI: 10.1515/pac-2022-1004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Abstract
We have quantum chemically studied activation of HnA–AHn bonds (AHn = CH3, NH2, OH, F) by PdLn catalysts with Ln = no ligand, PH3, (PH3)2, using relativistic density functional theory at ZORA-BLYP/TZ2P. The activation energy associated with the oxidative addition step decreases from H3C–CH3 to H2N–NH2 to HO–OH to F–F, where the activation of the F–F bond is barrierless. Activation strain and Kohn–Sham molecular orbital analyses reveal that the enhanced reactivity along this series of substrates originates from a combination of (i) reduced activation strain due to a weaker HnA–AHn bond; (ii) decreased Pauli repulsion as a result of a difference in steric shielding of the HnA–AHn bond; and (iii) enhanced backbonding interaction between the occupied 4d atomic orbitals of the palladium catalyst and σ* acceptor orbital of the substrate.
Collapse
Affiliation(s)
- Bryan Phuti Moloto
- Department of Theoretical Chemistry , Amsterdam Institute of Molecular and Life Sciences (AIMMS), and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam , The Netherlands , URL:
- Department of Chemistry and Polymer Science , Stellenbosch University , Private Bag X1 , Matieland , Stellenbosch , 7602 , South Africa
| | - Pascal Vermeeren
- Department of Theoretical Chemistry , Amsterdam Institute of Molecular and Life Sciences (AIMMS), and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam , The Netherlands , URL:
| | - Marco Dalla Tiezza
- Department of Theoretical Chemistry , Amsterdam Institute of Molecular and Life Sciences (AIMMS), and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam , The Netherlands , URL:
| | - Tessel Bouwens
- Department of Theoretical Chemistry , Amsterdam Institute of Molecular and Life Sciences (AIMMS), and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam , The Netherlands , URL:
| | - Catharine Esterhuysen
- Department of Chemistry and Polymer Science , Stellenbosch University , Private Bag X1 , Matieland , Stellenbosch , 7602 , South Africa
| | - Trevor A. Hamlin
- Department of Theoretical Chemistry , Amsterdam Institute of Molecular and Life Sciences (AIMMS), and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam , The Netherlands , URL:
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry , Amsterdam Institute of Molecular and Life Sciences (AIMMS), and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam , De Boelelaan 1083, 1081 HV Amsterdam , The Netherlands , URL:
- Institute for Molecules and Materials (IMM), Radboud University , Heyendaalseweg 135, 6525 AJ Nijmegen , The Netherlands
- Department of Chemical Sciences , University of Johannesburg , Auckland Park , Johannesburg 2006 , South Africa
| |
Collapse
|
4
|
Yu S, Tiekink EH, Vermeeren P, Bickelhaupt FM, Hamlin TA. How Bases Catalyze Diels-Alder Reactions. Chemistry 2023; 29:e202203121. [PMID: 36330879 PMCID: PMC10108159 DOI: 10.1002/chem.202203121] [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: 10/06/2022] [Revised: 11/03/2022] [Accepted: 11/04/2022] [Indexed: 11/06/2022]
Abstract
We have quantum chemically studied the base-catalyzed Diels-Alder (DA) reaction between 3-hydroxy-2-pyrone and N-methylmaleimide using dispersion-corrected density functional theory. The uncatalyzed reaction is slow and is preceded by the extrusion of CO2 via a retro-DA reaction. Base catalysis, for example, by triethylamine, lowers the reaction barrier up to 10 kcal mol-1 , causing the reaction to proceed smoothly at low temperature, which quenches the expulsion of CO2 , yielding efficient access to polyoxygenated natural compounds. Our activation strain analyses reveal that the base accelerates the DA reaction via two distinct electronic mechanisms: i) by the HOMO-raising effect, which enhances the normal electron demand orbital interaction; and ii) by donating charge into 3-hydroxy-2-pyrone which accumulates in its reactive region and promotes strongly stabilizing secondary electrostatic interactions with N-methylmaleimide.
Collapse
Affiliation(s)
- Song Yu
- 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, (TheNetherlands
| | - Eveline H Tiekink
- 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, (TheNetherlands
| | - 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, (TheNetherlands
| | - 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, (TheNetherlands.,Institute for Molecules and Materials (IMM), Radboud University, Heyendaalseweg 135, 6525 AJ, Nijmegen (The, Netherlands.,Department of Chemical Sciences, University of Johannesburg, Auckland Park, Johannesburg, 2006, South Africa
| | - 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, (TheNetherlands
| |
Collapse
|
5
|
Yang MC, Zhang ZF, Su MD. Understanding the Reactivity of Combination Reactions of Intramolecular Geminal Group 13 Element/Phosphorus and Gallium/Group 15 Element Frustrated Lewis Pairs with CS 2. Inorg Chem 2022; 61:12959-12976. [PMID: 35930703 DOI: 10.1021/acs.inorgchem.2c00259] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The reactions of CS2 captured by intramolecular geminal G13/P-based (G13 = group 13 elements) and Ga/G15-based (G15 = group 15 elements) frustrated Lewis pairs have been theoretically examined by using density functional theory (DFT) computations. With regard to the nine FLP-related compounds, our DFT calculated results reveal that only Al/P-Rea and Ga/P-Rea can kinetically and thermodynamically precede the energetically feasible combination reactions with CS2 to form the five-membered heterocyclic adducts. Our activation strain model analyses on the nine aforementioned model molecules indicate that the atomic radius of the Lewis acceptor (G13) and the Lewis donor (G15) plays a role in controlling their barrier heights to obtain good orbital overlaps among G13/P-Rea, Ga/G15-Rea, and CS2. Our theoretical observations based on the energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) approach strongly indicate that the donor-acceptor bonding (i.e., singlet-singlet bonding) rather than the electron-sharing bonding (i.e., triplet-triplet bonding) plays a central role in determining the bonding conditions of the transition states, G13/P-TS and Ga/G15-TS. In addition, the theoretical evidence obtained by the frontier molecular orbital theory and EDA-NOCV analyses reveals that the best description for the bonding natures of the combination reactions of intramolecular geminal G13/P-Rea and Ga/G15-Rea with CS2 is the lone pair(G15) → p-π*(C) interaction rather than the p-π*(G13) ← p-π(S) interaction. Moreover, our present DFT computations concerning the calculated structures and corresponding relative energetics of the stationary points connected with the aforementioned sophisticated approaches are in accordance with the Hammond postulate.
Collapse
Affiliation(s)
- Ming-Chung Yang
- Department of Applied Chemistry, National Chiayi University, Chiayi 60004, Taiwan
| | - Zheng-Feng Zhang
- Department of Applied Chemistry, National Chiayi University, Chiayi 60004, Taiwan
| | - Ming-Der Su
- Department of Applied Chemistry, National Chiayi University, Chiayi 60004, Taiwan.,Department of Medicinal and Applied Chemistry, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
| |
Collapse
|
6
|
Hansen T, Sun X, Dalla Tiezza M, van Zeist W, van Stralen JNP, Geerke DP, Wolters LP, Poater J, Hamlin TA, Bickelhaupt FM. C−X Bond Activation by Palladium: Steric Shielding versus Steric Attraction. Chemistry 2022; 28:e202201093. [PMID: 35420229 PMCID: PMC9401605 DOI: 10.1002/chem.202201093] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2022] [Indexed: 11/07/2022]
Abstract
The C−X bond activation (X = H, C) of a series of substituted C(n°)−H and C(n°)−C(m°) bonds with C(n°) and C(m°) = H3C− (methyl, 0°), CH3H2C− (primary, 1°), (CH3)2HC− (secondary, 2°), (CH3)3C− (tertiary, 3°) by palladium were investigated using relativistic dispersion‐corrected density functional theory at ZORA‐BLYP‐D3(BJ)/TZ2P. The effect of the stepwise introduction of substituents was pinpointed at the C−X bond on the bond activation process. The C(n°)−X bonds become substantially weaker going from C(0°)−X, to C(1°)−X, to C(2°)−X, to C(3°)−X because of the increasing steric repulsion between the C(n°)‐ and X‐group. Interestingly, this often does not lead to a lower barrier for the C(n°)−X bond activation. The C−H activation barrier, for example, decreases from C(0°)−X, to C(1°)−X, to C(2°)−X and then increases again for the very crowded C(3°)−X bond. For the more congested C−C bond, in contrast, the activation barrier always increases as the degree of substitution is increased. Our activation strain and matching energy decomposition analyses reveal that these differences in C−H and C−C bond activation can be traced back to the opposing interplay between steric repulsion across the C−X bond versus that between the catalyst and substrate.
Collapse
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
- Departament de Química Inorgànicai i Orgànica & IQTCUB Universitat de Barcelona Martí i Franquès 1-11 08028 Barcelona Spain
| | - Xiaobo Sun
- 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
- Departament de Química Inorgànicai i Orgànica & IQTCUB Universitat de Barcelona Martí i Franquès 1-11 08028 Barcelona Spain
| | - Marco Dalla Tiezza
- 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
| | - Willem‐Jan van Zeist
- 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
| | - Joost N. P. van Stralen
- 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
| | - Daan P. Geerke
- 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
| | - Lando P. Wolters
- 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
| | - Jordi Poater
- Departament de Química Inorgànicai i Orgànica & IQTCUB Universitat de Barcelona Martí i Franquès 1-11 08028 Barcelona Spain
- ICREA Passeig Lluís Companys 23 08010 Barcelona Spain
| | - 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
| | - 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
| |
Collapse
|
7
|
Abstract
We have quantum chemically explored the competition between the SN2 and SN2' pathways for X- + H2C═CHCH2Y (X, Y = F, Cl, Br, I) using a combined relativistic density functional theory and coupled-cluster theory approach. Bimolecular nucleophilic substitution reactions at allylic systems, i.e., Cγ═Cβ-Cα-Y, bearing a leaving-group at the α-position, proceed either via a direct attack at the α-carbon (SN2) or via an attack at the γ-carbon, involving a concerted allylic rearrangement (SN2'), in both cases leading to the expulsion of the leaving-group. Herein, we provide a physically sound model to rationalize under which circumstances a nucleophile will follow either the aliphatic SN2 or allylic SN2' pathway. Our activation strain analyses expose the underlying physical factors that steer the SN2/SN2' competition and, again, demonstrate that the concepts of a reaction's "characteristic distortivity" and "transition state acidity" provide explanations and design tools for understanding and predicting reactivity trends in organic synthesis.
Collapse
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.,Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, 08028 Barcelona, Spain
| | - 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
| | - Lea de Jong
- 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
| | - 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
|
8
|
Moloto BP, Vermeeren P, Dalla Tiezza M, Esterhuysen C, Bickelhaupt FM, Hamlin TA. Palladium‐Catalyzed Activation of Carbon–Halogen Bonds: Electrostatics‐Controlled Reactivity. European J Org Chem 2022. [DOI: 10.1002/ejoc.202200722] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | | | | | | | | | - Trevor A. Hamlin
- Vrije Universiteit Amsterdam Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling De Boelelaan 1083 1081 HV Amsterdam NETHERLANDS
| |
Collapse
|
9
|
Sorbelli D, Belanzoni P, Belpassi L, Lee JW, Ciancaleoni G. An ETS-NOCV-based computational strategies for the characterization of concerted transition states involving CO 2. J Comput Chem 2022; 43:717-727. [PMID: 35194805 PMCID: PMC9303928 DOI: 10.1002/jcc.26829] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2021] [Revised: 01/25/2022] [Accepted: 02/07/2022] [Indexed: 11/18/2022]
Abstract
Due to the presence of both a slightly acidic carbon and a slightly basic oxygen, carbon dioxide is often involved in concerted transition states (TSs) with two (or more) different molecular events interlaced in the same step. The possibility of isolating and quantitatively evaluating each molecular event would be important to characterize and understand the reaction mechanism in depth. This could be done, in principle, by measuring the relevant distances in the optimized TS, but often distances are not accurate enough, especially in the presence of many simultaneous processes. Here, we have applied the Extended Transition State‐Natural Orbital for Chemical Valence‐method (ETS‐NOCV), also in combination with the Activation Strain Model (ASM) and Energy Decomposition Analysis (EDA), to separate and quantify these molecular events at the TS of both organometallic and organic reactions. For the former, we chose the decomposition of formic acid to CO2 by an iridium catalyst, and for the latter, a CO2‐mediated transamidation and its chemical variations (hydro‐ and aminolysis of an ester) as case studies. We demonstrate that the one‐to‐one mapping between the “molecular events” and the ETS‐NOCV components is maintained along the entire lowest energy path connecting reactants and products around the TS, thus enabling a detailed picture on the relative importance of each interacting component. The methodology proposed here provides valuable insights into the effect of different chemical substituents on the reaction mechanism and promises to be generally applicable for any concerted TSs.
Collapse
Affiliation(s)
- Diego Sorbelli
- Department of Chemistry, Biology and Biotechnology, University of Perugia, Perugia, I-06123, Italy
| | - Paola Belanzoni
- Department of Chemistry, Biology and Biotechnology, University of Perugia, Perugia, I-06123, Italy.,CNR Institute of Chemical Science and Technologies "Giulio Natta" (CNR-SCITEC), c/o Department of Chemistry, Biology and Biotechnology, University of Perugia, Perugia, I-06123, Italy
| | - Leonardo Belpassi
- CNR Institute of Chemical Science and Technologies "Giulio Natta" (CNR-SCITEC), c/o Department of Chemistry, Biology and Biotechnology, University of Perugia, Perugia, I-06123, Italy
| | - Ji-Woong Lee
- Department of Chemistry, University of Copenhagen, Copenhagen, Ø 2100, Denmark.,Nanoscience Center, University of Copenhagen, Copenhagen, Ø 2100, Denmark
| | - Gianluca Ciancaleoni
- Department of Chemistry and Industrial Chemistry, University of Pisa, Pisa, I-56124, Italy.,CIRCC, Bari, Italy
| |
Collapse
|
10
|
Wang T, Yu Z, Liu T. Theoretical study on the highly diastereoselective palladium-catalyzed cascade carbocyclization of enallene. MOLECULAR CATALYSIS 2021. [DOI: 10.1016/j.mcat.2021.111994] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
11
|
Hansen T, Vermeeren P, Bickelhaupt FM, Hamlin TA. Origin of the α‐Effect in S
N
2 Reactions. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202106053] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/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
| | - 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 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
|
12
|
Hansen T, Vermeeren P, Bickelhaupt FM, Hamlin TA. Origin of the α-Effect in S N 2 Reactions. Angew Chem Int Ed Engl 2021; 60:20840-20848. [PMID: 34087047 PMCID: PMC8518820 DOI: 10.1002/anie.202106053] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 05/27/2021] [Indexed: 11/07/2022]
Abstract
The α-effect is a term used to explain the dramatically enhanced reactivity of α-nucleophiles (R-Y-X:- ) compared to their parent normal nucleophile (R-X:- ) by deviating from the classical Brønsted-type reactivity-basicity relationship. The exact origin of this effect is, however, still heavily under debate. In this work, we have quantum chemically analyzed the α-effect of a set of anionic nucleophiles, including O-, N- and S-based normal and α-nucleophiles, participating in an SN 2 reaction with ethyl chloride using relativistic density functional theory at ZORA-OLYP/QZ4P. Our activation strain and Kohn-Sham molecular orbital analyses identified two criteria an α-nucleophile needs to fulfill in order to show α-effect: (i) a small HOMO lobe on the nucleophilic center, pointing towards the substrate, to reduce the repulsive occupied-occupied orbital overlap and hence (steric) Pauli repulsion with the substrate; and (ii) a sufficiently high energy HOMO to overcome the loss of favorable HOMO-LUMO orbital overlap with the substrate, as a consequence of the first criterion, by reducing the HOMO-LUMO orbital energy gap. If one of these two criteria is not fulfilled, one can expect no α-effect or inverse α-effect.
Collapse
Affiliation(s)
- Thomas Hansen
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
- Leiden Institute of ChemistryLeiden UniversityEinsteinweg 552333 CCLeidenThe Netherlands
| | - Pascal Vermeeren
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe 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 10831081 HVAmsterdamThe Netherlands
- Institute for Molecules and MaterialsRadboud UniversityHeyendaalseweg 1356525 AJNijmegenThe Netherlands
| | - Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
| |
Collapse
|
13
|
Yu S, Bickelhaupt FM, Hamlin TA. Switch From Pauli-Lowering to LUMO-Lowering Catalysis in Brønsted Acid-Catalyzed Aza-Diels-Alder Reactions. ChemistryOpen 2021; 10:784-789. [PMID: 34351072 PMCID: PMC8340067 DOI: 10.1002/open.202100172] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Revised: 07/20/2021] [Indexed: 12/16/2022] Open
Abstract
Brønsted acid-catalyzed inverse-electron demand (IED) aza-Diels-Alder reactions between 2-aza-dienes and ethylene were studied using quantum chemical calculations. The computed activation energy systematically decreases as the basic sites of the diene progressively become protonated. Our activation strain and Kohn-Sham molecular orbital analyses traced the origin of this enhanced reactivity to i) "Pauli-lowering catalysis" for mono-protonated 2-aza-dienes due to the induction of an asynchronous, but still concerted, reaction pathway that reduces the Pauli repulsion between the reactants; and ii) "LUMO-lowering catalysis" for multi-protonated 2-aza-dienes due to their highly stabilized LUMO(s) and more concerted synchronous reaction path that facilitates more efficient orbital overlaps in IED interactions. In all, we illustrate how the novel concept of "Pauli-lowering catalysis" can be overruled by the traditional concept of "LUMO-lowering catalysis" when the degree of LUMO stabilization is extreme as in the case of multi-protonated 2-aza-dienes.
Collapse
Affiliation(s)
- Song Yu
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands) and
| | - F. Matthias Bickelhaupt
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands) and
- Institute for Molecules and Materials (IMM)Radboud UniversityHeyendaalseweg 1356525 AJNijmegenThe Netherlands
| | - Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands) and
| |
Collapse
|
14
|
Vermeeren P, Doppert MT, Bickelhaupt FM, Hamlin TA. How metallylenes activate small molecules. Chem Sci 2021; 12:4526-4535. [PMID: 34163718 PMCID: PMC8179460 DOI: 10.1039/d0sc05987k] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Accepted: 02/02/2021] [Indexed: 11/29/2022] Open
Abstract
We have studied the activation of dihydrogen by metallylenes using relativistic density functional theory (DFT). Our detailed activation strain and Kohn-Sham molecular orbital analyses have quantified the physical factors behind the decreased reactivity of the metallylene on going down Group 14, from carbenes to stannylenes. Along this series, the reactivity decreases due to a worsening of the back-donation interaction between the filled lone-pair orbital of the metallylene and the σ*-orbital of H2, which, therefore, reduces the metallylene-substrate interaction and increases the reaction barrier. As the metallylene ligand is varied from nitrogen to phosphorus to arsenic a significant rate enhancement is observed for the activation of H2 due to (i) a reduced steric (Pauli) repulsion between the metallylene and the substrate; and (ii) less activation strain, as the metallylene becomes increasingly more predistorted. Using a rationally designed metallylene with an optimal Group 14 atom and ligand combination, we show that a number of small molecules (i.e. HCN, CO2, H2, NH3) may also be readily activated. For the first time, we show the ability of our H2 activated designer metallylenes to hydrogenate unsaturated hydrocarbons. The results presented herein will serve as a guide for the rational design of metallylenes toward the activation of small molecules and subsequent reactions.
Collapse
Affiliation(s)
- 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
| | - Michael T Doppert
- 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
| | - 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
|
15
|
Fukaya K, Saito A, Nakajima N, Urabe D. Computational Analysis of the Selective Formation of the C4α-C8' Bond in the Intermolecular Coupling of Catechin Derivatives. J Org Chem 2020; 85:5010-5018. [PMID: 32149508 DOI: 10.1021/acs.joc.0c00261] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Procyanidin B3 is a natural flavonoid composed of two catechins connected via a C4α-C8' bond. The couplings of catechin derivatives, promoted by Lewis acids, have been widely applied to the syntheses of procyanidin B3 and related flavonoids because the reactions construct the C4α-C8' bond in a highly stereo- and regioselective manner. However, the structural complexity of the catechin derivatives has complicated the exploration of a detailed mechanism for this selectivity. Here, we report the results of a computational study to provide plausible origins for the selective C4α-C8' bond formation of catechin derivatives 1 and 2 by using models 5 and 7. Although a systematic search did not provide SN2-like transition states, we successfully identified transition states TS-A, TS-B, and TS-C for the SN1-type C4α-C8', C4β-C8', and C4α-C6' bond formations, respectively, from a total of 233 transition states to justify the stereo- and regioselectivity of the experimental results. The analysis of these structures by NCIPLOT mapping and the distortion/interaction strain model suggests that the eclipsed interaction at the forming C-C bond between the electrophile and the nucleophile destabilizes TS-B, while the strain of the electrophile destabilizes TS-C. Consequently, the C4α-C8' bond is formed via the lowest energy transition state TS-A.
Collapse
Affiliation(s)
- Keisuke Fukaya
- Biotechnology Research Center and Department of Biotechnology, Toyama Prefectural University, 5180 Kurokawa, Imizu, Toyama 939-0398, Japan
| | - Akiko Saito
- Graduate School of Engineering, Osaka Electro-Communication University, 18-8 Hatsu-cho, Neyagawa, Osaka 572-8530, Japan
| | - Noriyuki Nakajima
- Biotechnology Research Center and Department of Pharmaceutical Engineering, Toyama Prefectural University, 5180 Kurokawa, Imizu, Toyama 939-0398, Japan
| | - Daisuke Urabe
- Biotechnology Research Center and Department of Biotechnology, Toyama Prefectural University, 5180 Kurokawa, Imizu, Toyama 939-0398, Japan
| |
Collapse
|
16
|
Durán R, Herrera B. Theoretical Study of the Mechanism of Catalytic Enanteoselective N-H and O-H Insertion Reactions. J Phys Chem A 2020; 124:2-11. [PMID: 31809051 DOI: 10.1021/acs.jpca.9b07274] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Theoretical density functional theory (DFT) calculations were carried out to study bond insertion reactions using a copper(I)-Box-carbenoid as a bond activator. In order to understand the reaction mechanism where N-H and O-H bonds actively participate, the reaction force (RF) and activation strain model (ASM) were used. Results indicate that the first step of the reaction is barrierless for both bond insertions (N-H and O-H), and the second step of the insertion reaction in the phenol (O-H bond) is favored kinetically and thermodynamically with regard to the aniline substrate (N-H bond). The enantioselectivity is driven by the ligand of the catalyst by steric repulsion, favoring the formation of the R isomer. The analysis of the reaction force and ASM exhibited that the higher energy barrier in aniline is mainly due to a higher W2 contribution together with repulsive interactions, which hinders the insertion process.
Collapse
Affiliation(s)
- Rocío Durán
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia , Pontificia Universidad Católica de Chile , Av. Vicuña Mackenna, 4860 Macul, Santiago , Chile
| | - Barbara Herrera
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia , Pontificia Universidad Católica de Chile , Av. Vicuña Mackenna, 4860 Macul, Santiago , Chile
| |
Collapse
|
17
|
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
|
18
|
Sun X, Soini TM, Poater J, Hamlin TA, Bickelhaupt FM. PyFrag 2019-Automating the exploration and analysis of reaction mechanisms. J Comput Chem 2019; 40:2227-2233. [PMID: 31165500 PMCID: PMC6771738 DOI: 10.1002/jcc.25871] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2019] [Revised: 05/16/2019] [Accepted: 05/16/2019] [Indexed: 02/06/2023]
Abstract
We present a substantial update to the PyFrag 2008 program, which was originally designed to perform a fragment-based activation strain analysis along a provided potential energy surface. The original PyFrag 2008 workflow facilitated the characterization of reaction mechanisms in terms of the intrinsic properties, such as strain and interaction, of the reactants. The new PyFrag 2019 program has automated and reduced the time-consuming and laborious task of setting up, running, analyzing, and visualizing computational data from reaction mechanism studies to a single job. PyFrag 2019 resolves three main challenges associated with the automated computational exploration of reaction mechanisms: it (1) computes the reaction path by carrying out multiple parallel calculations using initial coordinates provided by the user; (2) monitors the entire workflow process; and (3) tabulates and visualizes the final data in a clear way. The activation strain and canonical energy decomposition results that are generated relate the characteristics of the reaction profile in terms of intrinsic properties (strain, interaction, orbital overlaps, orbital energies, populations) of the reactant species. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Xiaobo Sun
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale ModelingVrije Universiteit AmsterdamDe Boelelaan 1083, 1081 HVAmsterdamNetherlands
| | - Thomas M. Soini
- Software for Chemistry & Materials B.V.De Boelelaan 1083, 1081 HVAmsterdamNetherlands
| | - Jordi Poater
- ICREA, Pg. Lluís Companys 23, 08010 Barcelona, Spain and Departament de Química Inorgànica i Orgànica & IQTCUBUniversitat de Barcelona08028BarcelonaCataloniaSpain
| | - Trevor A. Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale ModelingVrije Universiteit AmsterdamDe Boelelaan 1083, 1081 HVAmsterdamNetherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale ModelingVrije Universiteit AmsterdamDe Boelelaan 1083, 1081 HVAmsterdamNetherlands
- Institute for Molecules and MaterialsRadboud UniversityHeyendaalseweg 135, 6525 AJNijmegenNetherlands
| |
Collapse
|
19
|
Narsaria AK, Hamlin TA, Lammertsma K, Bickelhaupt FM. Dual Activation of Aromatic Diels-Alder Reactions. Chemistry 2019; 25:9902-9912. [PMID: 31111976 PMCID: PMC6771859 DOI: 10.1002/chem.201901617] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Revised: 05/20/2019] [Indexed: 11/20/2022]
Abstract
The unusually fast Diels-Alder reactions of [5]cyclophanes were analyzed by DFT at the BLYP-D3(BJ)/TZ2P level of theory. The computations were guided by an integrated activation-strain and Kohn-Sham molecular orbital analysis. It is revealed why both [5]metacyclophane and [5]paracyclophane exhibit a significant rate enhancement compared to their planar benzene analogue. The activation strain analyses revealed that the enhanced reactivity originates from 1) predistortion of the aromatic core resulting in a reduced activation strain of the aromatic diene, and/or 2) enhanced interaction with the dienophile through a distortion-controlled lowering of the HOMO-LUMO gap within the diene. Both of these physical mechanisms and thus the rate of Diels-Alder cycloaddition can be tuned through different modes of geometrical distortion (meta versus para bridging) and by heteroatom substitution in the aromatic ring. Judicious choice of the bridge and heteroatom in the aromatic core enables effective tuning of the aromatic Diels-Alder reactivity to achieve activation barriers as low as 2 kcal mol-1 , which is an impressive 35 kcal mol-1 lower than that of benzene.
Collapse
Affiliation(s)
- Ayush K. Narsaria
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - Trevor A. Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - Koop Lammertsma
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Department of ChemistryUniversity of JohannesburgAuckland ParkJohannesburg2006South Africa
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Institute for Molecules and Materials (IMM)Radboud UniversityHeyendaalseweg 1356525AJNijmegenThe Netherlands
| |
Collapse
|
20
|
Hamlin TA, Levandowski BJ, Narsaria AK, Houk KN, Bickelhaupt FM. Structural Distortion of Cycloalkynes Influences Cycloaddition Rates both by Strain and Interaction Energies. Chemistry 2019; 25:6342-6348. [PMID: 30779472 PMCID: PMC6519225 DOI: 10.1002/chem.201900295] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2019] [Revised: 02/16/2019] [Indexed: 12/27/2022]
Abstract
The reactivities of 2‐butyne, cycloheptyne, cyclooctyne, and cyclononyne in the 1,3‐dipolar cycloaddition reaction with methyl azide were evaluated through DFT calculations at the M06‐2X/6‐311++G(d)//M06‐2X/6‐31+G(d) level of theory. Computed activation free energies for the cycloadditions of cycloalkynes are 16.5–22.0 kcal mol−1 lower in energy than that of the acyclic 2‐butyne. The strained or predistorted nature of cycloalkynes is often solely used to rationalize this significant rate enhancement. Our distortion/interaction–activation strain analysis has been revealed that the degree of geometrical predistortion of the cycloalkyne ground‐state geometries acts to enhance reactivity compared with that of acyclic alkynes through three distinct mechanisms, not only due to (i) a reduced strain or distortion energy, but also to (ii) a smaller HOMO–LUMO gap, and (iii) an enhanced orbital overlap, which both contribute to more stabilizing orbital interactions.
Collapse
Affiliation(s)
- Trevor A Hamlin
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands
| | - Brian J Levandowski
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA, 90095, USA
| | - Ayush K Narsaria
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands
| | - Kendall N Houk
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA, 90095, USA
| | - 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 for Molecules and Materials (IMM), Radboud University of Nijmegen, Heyendaalseweg 135, 6525 AJ, Nijmegen, The Netherlands
| |
Collapse
|
21
|
Yu S, de Bruijn HM, Svatunek D, Hamlin TA, Bickelhaupt FM. Factors Controlling the Diels-Alder Reactivity of Hetero-1,3-Butadienes. ChemistryOpen 2018; 7:995-1004. [PMID: 30524925 PMCID: PMC6276106 DOI: 10.1002/open.201800193] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Indexed: 12/29/2022] Open
Abstract
We have quantum chemically explored the Diels-Alder reactivities of a systematic series of hetero-1,3-butadienes with ethylene by using density functional theory at the BP86/TZ2P level. Activation strain analyses provided physical insight into the factors controlling the relative cycloaddition reactivity of aza- and oxa-1,3-butadienes. We find that dienes with a terminal heteroatom, such as 2-propen-1-imine (NCCC) or acrolein (OCCC), are less reactive than the archetypal 1,3-butadiene (CCCC), primarily owing to weaker orbital interactions between the more electronegative heteroatoms with ethylene. Thus, the addition of a second heteroatom at the other terminal position (NCCN and OCCO) further reduces the reactivity. However, the introduction of a nitrogen atom in the backbone (CNCC) leads to enhanced reactivity, owing to less Pauli repulsion resulting from polarization of the diene HOMO in CNCC towards the nitrogen atom and away from the terminal carbon atom. The Diels-Alder reactions of ethenyl-diazene (NNCC) and 1,3-diaza-butadiene (NCNC), which contain heteroatoms at both the terminal and backbone positions, are much more reactive due to less activation strain compared to CCCC.
Collapse
Affiliation(s)
- Song Yu
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM) Vrije Universiteit Amsterdam De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Hans M de Bruijn
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM) Vrije Universiteit Amsterdam De Boelelaan 1083 1081 HV Amsterdam The Netherlands
- Leiden Institute of Chemistry, Gorlaeus Laboratories Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Dennis Svatunek
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM) Vrije Universiteit Amsterdam De Boelelaan 1083 1081 HV Amsterdam The Netherlands
- Institut für Angewandte Synthesechemie Technische Universität Wien (TU Wien) Getreidemarkt 9 1060 Vienna Austria
| | - Trevor A Hamlin
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM) Vrije Universiteit Amsterdam De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - 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 for Molecules and Materials (IMM) Radboud University Heyendaalseweg 135 6525 AJ Nijmegen The Netherlands
| |
Collapse
|
22
|
Andryushchenko VA, Chekmarev SF. Modeling of Multicolor Single-Molecule Förster Resonance Energy-Transfer Experiments on Protein Folding. J Phys Chem B 2018; 122:10678-10685. [PMID: 30383961 DOI: 10.1021/acs.jpcb.8b07737] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Using a coarse-grained, Cα-model of BBL protein, a multicolor single-molecule Förster resonance energy transfer (FRET) experiment is modeled. Three fluorophores are introduced, which, for simplicity, are associated with Cα beads. Two fluorophores are placed at the ends of protein chain and the third one at the middle of the chain. The free-energy surfaces (FESs) depending on the interfluorophore distances and on the FRET efficiencies corresponding to these distances have been constructed and compared with the FESs depending on the conventional collective variables, such as the fraction of native contacts and radius of gyration. It has been found that multicolor experiments can successfully resolve all essential BBL states that are revealed by the conventional FESs. The resolution of these states with the FRET-efficiency histogram is found to be successful if the energy transfer is measured between the fluorophores at the BBL ends. We also show that, although the present model construct of BBL is very simple, it captures some characteristic features of the single-molecule FRET experiments, such as the pattern of the FRET-efficiency histograms and their evolution with the denaturant concentration.
Collapse
Affiliation(s)
- Vladimir A Andryushchenko
- Institute of Thermophysics , SB RAS , 630090 Novosibirsk , Russia.,Department of Physics , Novosibirsk State University , 630090 Novosibirsk , Russia
| | - Sergei F Chekmarev
- Institute of Thermophysics , SB RAS , 630090 Novosibirsk , Russia.,Department of Physics , Novosibirsk State University , 630090 Novosibirsk , Russia
| |
Collapse
|
23
|
Hamlin TA, Svatunek D, Yu S, Ridder L, Infante I, Visscher L, Bickelhaupt FM. Elucidating the Trends in Reactivity of Aza-1,3-Dipolar Cycloadditions. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800572] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Trevor A. Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Dennis Svatunek
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
- Institute of Applied Synthetic Chemistry; Technische Universität Wien (TU Wien); Getreidemarkt 9 1060 Vienna Austria
| | - Song Yu
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Lars Ridder
- Netherlands eScience Center; Science Park 140 1098 XG Amsterdam The Netherlands
| | - Ivan Infante
- Institute for Molecules and Materials (IMM); Radboud University; Heyendaalseweg 135 6525 AJ Nijmegen The Netherlands
| | - Lucas Visscher
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - 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
|
24
|
Arylic C-X Bond Activation by Palladium Catalysts: Activation Strain Analyses of Reactivity Trends. Sci Rep 2018; 8:10729. [PMID: 30013049 PMCID: PMC6048108 DOI: 10.1038/s41598-018-28998-3] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Accepted: 07/03/2018] [Indexed: 11/09/2022] Open
Abstract
We have quantum chemically explored arylic carbon-substituent bond activation via oxidative insertion of a palladium catalyst in C6H5X + PdLn model systems (X = H, Cl, CH3; Ln = no ligand, PH3, (PH3)2, PH2C2H4PH2) using relativistic density functional theory at ZORA-BLYP/TZ2P. Besides exploring reactivity trends and comparing them to aliphatic C-X activation, we aim at uncovering the physical factors behind the activity and selectivity. Our results show that barriers for arylic C-X activation are lower than those for the corresponding aliphatic C-X bonds. However, trends along bonds or upon variation of ligands are similar. Thus, bond activation barriers increase along C-Cl < C-H < C-C and along Pd < Pd(PH3) or Pd(PH2C2H4PH2) < Pd(PH3)2. Activation strain analyses in conjunction with quantitative molecular orbital theory trace these trends to the rigidity and bonding capability of the various C-X bonds, model catalysts, and ligands.
Collapse
|
25
|
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
|
26
|
Merzoud L, Saal A, Moussaoui R, Ouamerali O, Morell C, Chermette H. Fluorine substituent effect on the stereochemistry of catalyzed and non-catalyzed Diels–Alder reactions. The case of R-butenone with cyclopentadiene: a computational assessment of the mechanism. Phys Chem Chem Phys 2018; 20:16102-16116. [DOI: 10.1039/c8cp00985f] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
A semiquantitative prediction of stereoselectivity due to substitutions of dienophile is obtained.
Collapse
Affiliation(s)
- Lynda Merzoud
- Laboratory of Computational and Theoretical Chemistry and Photonics
- USTHB University
- Algeria
- Département de Chimie
- UMMTO University of Tizi–Ouzou
| | - Amar Saal
- Laboratory of Computational and Theoretical Chemistry and Photonics
- USTHB University
- Algeria
- Département de Chimie
- UMMTO University of Tizi–Ouzou
| | - Ramdane Moussaoui
- Laboratoire de Chimie Appliquée et de Génie Chimique
- Université M. Mammeri
- Tizi Ouzou
- Algeria
| | - Ourida Ouamerali
- Laboratory of Computational and Theoretical Chemistry and Photonics
- USTHB University
- Algeria
| | - Christophe Morell
- Université de Lyon
- Université Claude Bernard Lyon 1
- ENS-Lyon
- Institut des Sciences Analytiques
- 69622 Villeurbanne Cedex
| | - Henry Chermette
- Université de Lyon
- Université Claude Bernard Lyon 1
- ENS-Lyon
- Institut des Sciences Analytiques
- 69622 Villeurbanne Cedex
| |
Collapse
|
27
|
Bickelhaupt FM, Houk KN. Das Distortion/Interaction‐Activation‐Strain‐Modell zur Analyse von Reaktionsgeschwindigkeiten. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201701486] [Citation(s) in RCA: 176] [Impact Index Per Article: 25.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- F. Matthias Bickelhaupt
- Department of Theoretical Chemistry und Amsterdam Center for Multiscale Modeling (ACMM) Vrije Universiteit Amsterdam De Boelelaan 1083 1081 HV Amsterdam Niederlande
- Institute of Molecules and Materials (IMM) Radboud University Heyendaalseweg 135 6525 AJ Nijmegen Niederlande
| | - Kendall N. Houk
- Department of Chemistry and Biochemistry und Department of Chemical and Biomolecular Engineering University of California Los Angeles CA 90095-1569 USA
| |
Collapse
|
28
|
Bickelhaupt FM, Houk KN. Analyzing Reaction Rates with the Distortion/Interaction-Activation Strain Model. Angew Chem Int Ed Engl 2017; 56:10070-10086. [PMID: 28447369 PMCID: PMC5601271 DOI: 10.1002/anie.201701486] [Citation(s) in RCA: 948] [Impact Index Per Article: 135.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2017] [Revised: 04/10/2017] [Indexed: 12/21/2022]
Abstract
The activation strain or distortion/interaction model is a tool to analyze activation barriers that determine reaction rates. For bimolecular reactions, the activation energies are the sum of the energies to distort the reactants into geometries they have in transition states plus the interaction energies between the two distorted molecules. The energy required to distort the molecules is called the activation strain or distortion energy. This energy is the principal contributor to the activation barrier. The transition state occurs when this activation strain is overcome by the stabilizing interaction energy. Following the changes in these energies along the reaction coordinate gives insights into the factors controlling reactivity. This model has been applied to reactions of all types in both organic and inorganic chemistry, including substitutions and eliminations, cycloadditions, and several types of organometallic reactions.
Collapse
Affiliation(s)
- 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 of Molecules and Materials (IMM), Radboud University, Heyendaalseweg 135, 6525 AJ, Nijmegen, The Netherlands
| | - Kendall N Houk
- Department of Chemistry and Biochemistry and Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, CA, 90095-1569, USA
| |
Collapse
|
29
|
Jing FQ, Cao JW, Liu XJ, Hu YF, Ma HT, Bian WS. Theoretical Study on Mechanism and Kinetics of Reaction of O(3P) with Propane. CHINESE J CHEM PHYS 2016. [DOI: 10.1063/1674-0068/29/cjcp1603042] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
30
|
Wolters LP, Koekkoek R, Bickelhaupt FM. Role of Steric Attraction and Bite-Angle Flexibility in Metal-Mediated C–H Bond Activation. ACS Catal 2015. [DOI: 10.1021/acscatal.5b01354] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Lando P. Wolters
- Department
of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Rick Koekkoek
- Department
of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - F. Matthias Bickelhaupt
- Department
of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University, 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
| |
Collapse
|
31
|
Wolters LP, Bickelhaupt FM. Selective C-H and C-C Bond Activation: Electronic Regimes as a Tool for Designing d(10) MLn Catalysts. Chem Asian J 2015. [PMID: 26218844 DOI: 10.1002/asia.201500368] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We wish to understand how a transition-metal catalyst can be rationally designed so as to selectively activate one particular bond in a substrate, herein, C-H and C-C bonds in ethane. To this end, we quantum chemically analyzed the activity and selectivity of a large series of model catalysts towards ethane and, for comparison, methane, by using the activation strain model and quantitative molecular orbital theory. The model catalysts comprise d(10) MLn complexes with coordination numbers n=0, 1, and 2; metal centers M=Co(-), Rh(-), Ir(-), Ni, Pd, Pt, Cu(+), Ag(+), and Au(+); and ligands L=NH3, PH3, and CO. Our analyses reveal that rather subtle electronic differences between bonds can be exploited to induce a lower barrier for activating one or the other, depending, among other factors, on the catalysts electronic regime (i.e., s-regime versus d-regime catalysts). Interestingly, the concepts and design principles emerging from this work can also be applied to the more challenging problem of differentiating between activation of the C-H bonds in ethane versus those in methane.
Collapse
Affiliation(s)
- Lando P Wolters
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands. .,Radboud University Nijmegen, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ, Nijmegen, The Netherlands.
| |
Collapse
|
32
|
Wolters LP, Bickelhaupt FM. The activation strain model and molecular orbital theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2015; 5:324-343. [PMID: 26753009 PMCID: PMC4696410 DOI: 10.1002/wcms.1221] [Citation(s) in RCA: 239] [Impact Index Per Article: 26.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/16/2015] [Revised: 04/05/2015] [Accepted: 04/07/2015] [Indexed: 12/23/2022]
Abstract
The activation strain model is a powerful tool for understanding reactivity, or inertness, of molecular species. This is done by relating the relative energy of a molecular complex along the reaction energy profile to the structural rigidity of the reactants and the strength of their mutual interactions: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). We provide a detailed discussion of the model, and elaborate on its strong connection with molecular orbital theory. Using these approaches, a causal relationship is revealed between the properties of the reactants and their reactivity, e.g., reaction barriers and plausible reaction mechanisms. This methodology may reveal intriguing parallels between completely different types of chemical transformations. Thus, the activation strain model constitutes a unifying framework that furthers the development of cross-disciplinary concepts throughout various fields of chemistry. We illustrate the activation strain model in action with selected examples from literature. These examples demonstrate how the methodology is applied to different research questions, how results are interpreted, and how insights into one chemical phenomenon can lead to an improved understanding of another, seemingly completely different chemical process. WIREs Comput Mol Sci 2015, 5:324-343. doi: 10.1002/wcms.1221.
Collapse
Affiliation(s)
- Lando P Wolters
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), VU University AmsterdamAmsterdam, The Netherlands; Dipartimento di Scienze Chimiche, Università degli Studi di PadovaPadova, Italy
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), VU University AmsterdamAmsterdam, The Netherlands; Institute of Molecules and Materials (IMM), Radboud University NijmegenNijmegen, The Netherlands
| |
Collapse
|
33
|
Wolters LP, van Zeist WJ, Bickelhaupt FM. New Concepts for Designing d10-M(L)nCatalysts: d Regime, s Regime and Intrinsic Bite-Angle Flexibility. Chemistry 2014; 20:11370-81. [DOI: 10.1002/chem.201403237] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2014] [Indexed: 11/09/2022]
|
34
|
Driving and retarding forces in a chemical reaction. J Mol Model 2014; 20:2351. [PMID: 25037489 DOI: 10.1007/s00894-014-2351-0] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2013] [Accepted: 06/15/2014] [Indexed: 10/25/2022]
Abstract
The reaction force F(ξ) is the negative gradient of the potential energy of a chemical process along the intrinsic reaction coordinate ξ. We extend the rigorous concept of F(ξ) to the "activation strain model" of Bickelhaupt et al., to formulate the "strain" force Fstr(ξ) that retards a reaction and the "interaction" force Fint(ξ) that drives it. These are investigated for a group of Diels-Alder cycloadditions. The results fully support the interpretation of the minimum of F(ξ) as defining the beginning of the transition from deformed reactants to eventual products.
Collapse
|
35
|
van Zeist WJ, Bickelhaupt FM. Steric nature of the bite angle. A closer and a broader look. Dalton Trans 2011; 40:3028-38. [PMID: 21331411 DOI: 10.1039/c0dt01550d] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The bite angle (ligand-metal-ligand angle) is known to greatly influence the activity of catalytically active transition-metal complexes towards bond activation. Here, we have computationally explored how and why the bite angle has such effects in a wide range of prototypical C-X bonds and palladium complexes, using relativistic density functional theory at ZORA-BLYP/TZ2P. Our model reactions cover the substrates H(3)C-X (with X = H, CH(3), Cl) and, among others, the model catalysts, Pd[PH(2)(CH(2))(n)PH(2)] (with n = 2-6) and Pd[PR(2)(CH(2))(n)PR(2)] (n = 2-4 and R = Me, Ph, t-Bu, Cl), Pd(PH(3))X(-) (X = Cl, Br, I), as well as palladium complexes of chelating and non-chelating N-heterocyclic carbenes. The purpose is to elaborate on an earlier finding that bite-angle effects have a predominantly (although not exclusively) steric nature: a smaller bite angle makes more room for coordinating a substrate by bending away the ligands. Indeed, the present results further consolidate this steric picture by revealing its occurrence in this broader range of model reactions and by identifying and quantifying the exact working mechanism through activation strain analyses.
Collapse
Affiliation(s)
- Willem-Jan van Zeist
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands
| | | |
Collapse
|
36
|
van Zeist WJ, Bickelhaupt FM. The activation strain model of chemical reactivity. Org Biomol Chem 2010; 8:3118-27. [PMID: 20490400 DOI: 10.1039/b926828f] [Citation(s) in RCA: 538] [Impact Index Per Article: 38.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Herein, we provide an account of the activation strain model of chemical reactivity and its recent applications. In this model, the potential energy surface DeltaE(zeta) along the reaction coordinate zeta is decomposed into the strain DeltaE(strain)(zeta) of the increasingly deformed reactants plus the interaction DeltaE(int)(zeta) between these deformed reactants, i.e., DeltaE(zeta) = DeltaE(strain)(zeta) + DeltaE(int)(zeta). The purpose of this fragment-based approach is to arrive at a qualitative understanding, based on accurate calculations, of the trends in activation barriers and transition-state geometries (e.g., early or late along the reaction coordinate) in terms of the reactants' properties. The usage of the activation strain model is illustrated by a number of concrete applications, by us and others, in the fields of catalysis and organic chemistry.
Collapse
Affiliation(s)
- Willem-Jan van Zeist
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, VU University, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands
| | | |
Collapse
|
37
|
|