1
|
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
|
2
|
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
|
3
|
Lovitt CF, Capra NE, Lastowski RJ, Girolami GS. Steric and Electronic Analyses of Ligand Effects on the Stability of σ-Methane Coordination Complexes: A DFT Study. Organometallics 2022. [DOI: 10.1021/acs.organomet.2c00434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Charity Flener Lovitt
- School of Science, Technology, Engineering & Mathematics, University of Washington Bothell, 18115 Campus Way NE, Bothell, Washington 98011, United States
- School of Chemical Sciences, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Nicolas E. Capra
- School of Chemical Sciences, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - R. Joseph Lastowski
- School of Chemical Sciences, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Gregory S. Girolami
- School of Chemical Sciences, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
Hansen T, Sun X, Dalla Tiezza M, van Zeist WJ, Poater J, Hamlin TA, Bickelhaupt FM. C(spn)-X (n = 1-3) Bond Activation by Palladium. Chemistry 2021; 28:e202103953. [PMID: 34958486 PMCID: PMC9306469 DOI: 10.1002/chem.202103953] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Indexed: 11/09/2022]
Abstract
We have studied the palladium-mediated activation of C( sp n )-X bonds (n = 1-3 and X = H, CH 3 , Cl) in archetypal model substrates H 3 C-CH 2 -X, H 2 C=CH-X and HC≡C-X by catalysts PdL n with L n = no ligand, Cl - , and (PH 3 ) 2 , using relativistic density functional theory at ZORA-BLYP/TZ2P. The oxidative addition barrier decreases along this series, even though the strength of the bonds increases going from C( sp 3 )-X, to C( sp 2 )-X, to C( sp )-X. Activation strain and matching energy decomposition analyses reveal that the decreased oxidative addition barrier going from sp 3 to sp 2 to sp , originates from a reduction in the destabilizing steric (Pauli) repulsion between catalyst and substrate. This is the direct consequence of the decreasing coordination number of the carbon atom in C( sp n )-X, which goes from four, to three, to two along this series. The associated net stabilization of the catalyst-substrate interaction dominates the trend in strain energy which indeed becomes more destabilizing along this same series as the bond becomes stronger from C( sp 3 )-X to C( sp )-X.
Collapse
Affiliation(s)
- Thomas Hansen
- Vrije Universiteit Amsterdam, Theoretical Chemistry, NETHERLANDS
| | - Xiaobo Sun
- Vrije Universiteit Amsterdam, Theoretical Chemistry, NETHERLANDS
| | | | | | - Jordi Poater
- University of Barcelona: Universitat de Barcelona, Inorganic and organic chemistry, SPAIN
| | - Trevor A Hamlin
- Vrije Universiteit Amsterdam, Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, De Boelelaan 1083, 1081 HV, Amsterdam, NETHERLANDS
| | | |
Collapse
|
7
|
de Azevedo Santos L, Ramalho TC, Hamlin TA, Bickelhaupt FM. Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study. J Comput Chem 2021; 42:688-698. [PMID: 33543482 PMCID: PMC7986859 DOI: 10.1002/jcc.26489] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2020] [Revised: 01/15/2021] [Accepted: 01/18/2021] [Indexed: 12/27/2022]
Abstract
We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.
Collapse
Affiliation(s)
- Lucas de Azevedo Santos
- Department of Theoretical Chemistry, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamAmsterdamNetherlands
- Department of Chemistry, Institute of Natural SciencesFederal University of LavrasLavrasBrazil
| | - Teodorico C. Ramalho
- Department of Chemistry, Institute of Natural SciencesFederal University of LavrasLavrasBrazil
- Center for Basic and Applied ResearchUniversity Hradec KraloveHradec KraloveCzech Republic
| | - Trevor A. Hamlin
- Department of Theoretical Chemistry, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamAmsterdamNetherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamAmsterdamNetherlands
- Institute for Molecules and MaterialsRadboud University NijmegenNijmegenNetherlands
| |
Collapse
|
8
|
Roglans A, Pla-Quintana A, Solà M. Mechanistic Studies of Transition-Metal-Catalyzed [2 + 2 + 2] Cycloaddition Reactions. Chem Rev 2020; 121:1894-1979. [DOI: 10.1021/acs.chemrev.0c00062] [Citation(s) in RCA: 66] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Anna Roglans
- Institut de Quı́mica Computacional i Catàlisi (IQCC) and Departament de Quı́mica, Universitat de Girona, C/Maria Aurèlia Capmany, 69, E-17003, Girona, Catalonia, Spain
| | - Anna Pla-Quintana
- Institut de Quı́mica Computacional i Catàlisi (IQCC) and Departament de Quı́mica, Universitat de Girona, C/Maria Aurèlia Capmany, 69, E-17003, Girona, Catalonia, Spain
| | - Miquel Solà
- Institut de Quı́mica Computacional i Catàlisi (IQCC) and Departament de Quı́mica, Universitat de Girona, C/Maria Aurèlia Capmany, 69, E-17003, Girona, Catalonia, Spain
| |
Collapse
|
9
|
Zhang L, Yamazaki K, Leitch JA, Manzano R, Atkinson VAM, Hamlin TA, Dixon DJ. Dual catalytic enantioselective desymmetrization of allene-tethered cyclohexanones. Chem Sci 2020; 11:7444-7450. [PMID: 34123026 PMCID: PMC8159440 DOI: 10.1039/d0sc02878a] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2020] [Accepted: 06/23/2020] [Indexed: 11/30/2022] Open
Abstract
The construction of enantioenriched azabicyclo[3.3.1]nonan-6-one heterocycles via an enantioselective desymmetrization of allene-linked cyclohexanones, enabled through a dual catalytic system, that provides synchronous activation of the cyclohexanone with a chiral prolinamide and the allene with a copper(i) co-catalyst to deliver the stereodefined bicyclic core, is described. Successful application to oxygen analogues was also achieved, thereby providing a new enantioselective synthetic entry to architecturally complex bicyclic ethereal frameworks. The mechanistic pathway and the origin of enantio- and diastereoselectivities has been uncovered using density functional theory (DFT) calculations.
Collapse
Affiliation(s)
- Lin Zhang
- Department of Chemistry, Chemistry Research Laboratory, University of Oxford 12 Mansfield Road Oxford UK
| | - Ken Yamazaki
- Department of Chemistry, Chemistry Research Laboratory, University of Oxford 12 Mansfield Road Oxford UK
| | - Jamie A Leitch
- Department of Chemistry, Chemistry Research Laboratory, University of Oxford 12 Mansfield Road Oxford UK
| | - Ruben Manzano
- Department of Chemistry, Chemistry Research Laboratory, University of Oxford 12 Mansfield Road Oxford UK
| | - Victoria A M Atkinson
- Department of Chemistry, Chemistry Research Laboratory, University of Oxford 12 Mansfield Road Oxford UK
| | - 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
| | - Darren J Dixon
- Department of Chemistry, Chemistry Research Laboratory, University of Oxford 12 Mansfield Road Oxford UK
| |
Collapse
|
10
|
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
|
11
|
Li J, Meng X, Hu C, Su Z. Cooperative Catalysis of Chiral Guanidine and Rh 2(OAc) 4 in Asymmetric O-H Insertion of Carboxylic Acid: A Theoretical Investigation. J Org Chem 2019; 84:15020-15031. [PMID: 31675228 DOI: 10.1021/acs.joc.8b02850] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We reported a mechanistic study on asymmetric O-H insertion reaction of α-diazoester with carboxylic acid using Rh2(OAc)4/chiral guanidine-amide as the cocatalyst by density functional theory [B3LYP-D3(BJ)/def2-TZVP//B3LYP-D3(BJ)/[6-31G**, SDD] (SMD, Et2O)]. The catalytic reaction included two stages: (i) formation of Rh-carbene species, subsequently by the construction of C-O bond forming enol and (ii) chiral guanidinium salt-assisted H-transfer to the enol. In cooperative catalysis, Rh2(OAc)4 helped to form an enol intermediate via high-reactivity Rh-carbene species, while the in situ-formed guanidium carboxylate acted as a chiral proton shuttle to construct a hydrogen bonding net for the stereo-determinant protonation. The repulsions between the phenyl group of the enol intermediate and the cyclohexyl as well as the ortho-substituted isopropyl group of chiral guanidine played important roles in controlling stereoselectivity. A disadvantageous steric arrangement in si-face attack weakened the stabilizing electrostatic and orbital interaction of reacting species in the H-transfer step, enhancing the pathway to form a predominant product with R-configuration in the two competing pathways. A model was proposed to explain the asymmetric induction of chiral guanidine-amide in protonation.
Collapse
Affiliation(s)
- Jing Li
- Key Laboratory of Green Chemistry and Technology, Ministry of Education, College of Chemistry , Sichuan University , Chengdu , Sichuan 610064 , P.R. China
| | - Xiangxiang Meng
- Key Laboratory of Green Chemistry and Technology, Ministry of Education, College of Chemistry , Sichuan University , Chengdu , Sichuan 610064 , P.R. China
| | - Changwei Hu
- Key Laboratory of Green Chemistry and Technology, Ministry of Education, College of Chemistry , Sichuan University , Chengdu , Sichuan 610064 , P.R. China
| | - Zhishan Su
- Key Laboratory of Green Chemistry and Technology, Ministry of Education, College of Chemistry , Sichuan University , Chengdu , Sichuan 610064 , P.R. China
| |
Collapse
|
12
|
Carlsen R, Jenkins JR, Huang TCJ, Pugh SL, Ess DH. Paddle Ball Dynamics during Conversion of a Rh–Methyl Hydride Complex to a Rh–Methane σ-Complex through Reductive Coupling. Organometallics 2019. [DOI: 10.1021/acs.organomet.8b00936] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ryan Carlsen
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States
| | - Jordan R. Jenkins
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States
| | - Tsung-Chiang Johnny Huang
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States
| | - Samuel L. Pugh
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States
| | - Daniel H. Ess
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, United States
| |
Collapse
|
13
|
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
|
14
|
Nitsch J, Wolters LP, Fonseca Guerra C, Bickelhaupt FM, Steffen A. Enhanced π-Back-Donation as a Way to Higher Coordination Numbers in d10[M(NHC)n] Complexes: A DFT Study. Chemistry 2016; 23:614-622. [DOI: 10.1002/chem.201603861] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2016] [Indexed: 11/07/2022]
Affiliation(s)
- Jörn Nitsch
- Institut für Anorganische Chemie; Universität Würzburg; Am Hubland 97074 Würzburg Germany
| | - Lando P. Wolters
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM); Vrije Universiteit Amsterdam; De Boelelaan 1083 1081 HV Amsterdam The Netherlands
| | - Célia Fonseca Guerra
- 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
| | - Andreas Steffen
- Institut für Anorganische Chemie; Universität Würzburg; Am Hubland 97074 Würzburg Germany
| |
Collapse
|
15
|
Ma PP, Wang YC, Wang WX, Deng ZP, Niu GP, Wang XL, Li S, Zhang YW. Theoretical investigation on activation of CH and CC bonds of 2-butyne by gas-phase Nb atom. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.03.033] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
16
|
Wang J, Liu L, Wilson AK. Oxidative Cleavage of the β-O-4 Linkage of Lignin by Transition Metals: Catalytic Properties and the Performance of Density Functionals. J Phys Chem A 2016; 120:737-46. [DOI: 10.1021/acs.jpca.5b08854] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiaqi Wang
- Department of Chemistry and
Center for Advanced Scientific Computing and Modeling, University of North Texas, Denton, Texas 76203-5017, United States
| | - Lily Liu
- Department of Chemistry and
Center for Advanced Scientific Computing and Modeling, University of North Texas, Denton, Texas 76203-5017, United States
| | - Angela K. Wilson
- Department of Chemistry and
Center for Advanced Scientific Computing and Modeling, University of North Texas, Denton, Texas 76203-5017, United States
| |
Collapse
|
17
|
de Jong GT, Bickelhaupt FM. Catalytic Carbon-Halogen Bond Activation: Trends in Reactivity, Selectivity, and Solvation. J Chem Theory Comput 2015; 3:514-29. [PMID: 26637032 DOI: 10.1021/ct600342j] [Citation(s) in RCA: 95] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have theoretically studied the oxidative addition of all halomethanes CH3X (with X = F, Cl, Br, I, At) to Pd and PdCl(-), using both nonrelativistic and zeroth-order-regular-approximation-relativistic density functional theory at BLYP/QZ4P. Our study covers the gas phase as well as the condensed phase (water), where solvent effects are described with the conductor-like screening model. The activation of the C*-X bond may proceed via two stereochemically different pathways: (i) direct oxidative insertion (OxIn) which goes with retention of the configuration at C* and (ii) an alternative SN2 pathway which goes with inversion of the configuration at C*. In the gas phase, for Pd, the OxIn pathway has the lowest reaction barrier for all CH3X's. Anion assistance, that is, going from Pd to PdCl(-), changes the preference for all CH3X's from OxIn to the SN2 pathway. Gas-phase reaction barriers for both pathways to C-X activation generally decrease as X descends in group 17. Two striking solvent effects are (i) the shift in reactivity of Pd + CH3X from OxIn to SN2 in the case of the smaller halogens, F and Cl, and (ii) the shift in reactivity of PdCl(-) + CH3X in the opposite direction, that is, from SN2 to OxIn, in the case of the heavier halogens, I and At. We use the activation strain model to arrive at a qualitative understanding of how the competition between OxIn and SN2 pathways is determined by the halogen atom in the activated C-X bond, by anion assistance, and by solvation.
Collapse
Affiliation(s)
- G Theodoor de Jong
- Afdeling Theoretische Chemie, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - F Matthias Bickelhaupt
- Afdeling Theoretische Chemie, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
18
|
Sienkiewicz N, Strzelec K, Pospiech P, Cypryk M, Szmechtyk T. New palladium catalyst immobilized on epoxy resin: synthesis, characterization and catalytic activity. Appl Organomet Chem 2015. [DOI: 10.1002/aoc.3389] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Natalia Sienkiewicz
- Institute of Polymer and Dye Technology, Faculty of Chemistry; Lodz University of Technology; Stefanowskiego 12/16 90-924 Lodz Poland
| | - Krzysztof Strzelec
- Institute of Polymer and Dye Technology, Faculty of Chemistry; Lodz University of Technology; Stefanowskiego 12/16 90-924 Lodz Poland
| | - Piotr Pospiech
- Center of Molecular and Macromolecular Studies; Polish Academy of Sciences; Sienkiewicza 112 90-363 Lodz Poland
| | - Marek Cypryk
- Center of Molecular and Macromolecular Studies; Polish Academy of Sciences; Sienkiewicza 112 90-363 Lodz Poland
| | - Tomasz Szmechtyk
- Institute of Polymer and Dye Technology, Faculty of Chemistry; Lodz University of Technology; Stefanowskiego 12/16 90-924 Lodz Poland
| |
Collapse
|
19
|
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
|
20
|
Rahmani L, Ketabi S. Solvation of alanine and histidine functionalized carbon nanotubes in aqueous media: A Monte Carlo simulation study. J Mol Liq 2015. [DOI: 10.1016/j.molliq.2015.04.021] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
|
21
|
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
|
22
|
Performance of recent density functionals to discriminate between olefin and nitrogen binding to palladium. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1533-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
23
|
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]
|
24
|
Grüber R, Fleurat-Lessard P. Density Functional Study of Indole Formation by an Intramolecular Heck Reaction. Organometallics 2014. [DOI: 10.1021/om5001309] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Raymond Grüber
- Université de Lyon, Institut de Chimie de Lyon,
École Normale Supérieure
de Lyon, CNRS, 46 allée d’Italie, 69364 Lyon Cedex 07, France
| | - Paul Fleurat-Lessard
- Université de Lyon, Institut de Chimie de Lyon,
École Normale Supérieure
de Lyon, CNRS, 46 allée d’Italie, 69364 Lyon Cedex 07, France
| |
Collapse
|
25
|
Orian L, Swart M, Bickelhaupt FM. Indenyl effect due to metal slippage? Computational exploration of rhodium-catalyzed acetylene [2+2+2] cyclotrimerization. Chemphyschem 2014; 15:219-28. [PMID: 24273097 DOI: 10.1002/cphc.201300934] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2013] [Indexed: 11/07/2022]
Abstract
The mechanism of CpRh (Cp=cyclopentadienyl) and IndRh (Ind=indenyl)-catalyzed acetylene [2+2+2] cyclotrimerization has been revisited aiming at finding an explanation for the better performance of the latter catalyst found experimentally. The hypothesis that an ancillary ligand of the precatalyst remains bonded to the metal center throughout the whole catalytic cycle, based on the experimental evidence that the nature of this ligand can exert some control in cocyclotrimerization of different alkynes, is considered. Strong hapticity variations occur in both the CpRh- and IndRh-catalyzed processes. As the Ind ligand undergoes a more facile slippage than Cp, the energy profile is far smoother in the IndRh-catalyzed cyclotrimerization. This difference in the energetics of the process translates into an enhanced activity of the IndRh catalyst, in nice agreement with experiment.
Collapse
Affiliation(s)
- Laura Orian
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35129 Padova (Italy).
| | | | | |
Collapse
|
26
|
Orian L, Wolters LP, Bickelhaupt FM. In silico design of heteroaromatic half-sandwich RhI catalysts for acetylene [2+2+2] cyclotrimerization: evidence of a reverse indenyl effect. Chemistry 2013; 19:13337-47. [PMID: 24038672 DOI: 10.1002/chem.201301990] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2013] [Indexed: 11/09/2022]
Abstract
A mechanistic density functional theory study of acetylene [2+2+2] cyclotrimerization to benzene catalyzed by Rh(I) half metallocenes is presented. The catalyst fragment contains a heteroaromatic ligand, that is, the 1,2-azaborolyl (Ab) or the 3a,7a-azaborindenyl (Abi) anions, which are isostructural and isoelectronic to the hydrocarbon cyclopentadienyl (Cp) and indenyl (Ind) anions, respectively, but differ from the last ones on having two adjacent carbon atoms replaced with a boron and a nitrogen atom. The better performance of either the classic hydrocarbon or the heteroaromatic catalysts is found to depend on the different mechanistic paths that can be envisioned for the process. The present analyses uncover and explain general structure-reactivity relationships that may serve as rational design principles. In particular, we provide evidence of a reverse indenyl effect.
Collapse
Affiliation(s)
- Laura Orian
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35129 Padova (Italy).
| | | | | |
Collapse
|
27
|
Munz D, Meyer D, Strassner T. Methane CH Activation by Palladium Complexes with Chelating Bis(NHC) Ligands: A DFT Study. Organometallics 2013. [DOI: 10.1021/om400232u] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Dominik Munz
- Physikalische Organische
Chemie, TU Dresden, Bergstrasse
66, 01069 Dresden, Germany
| | - Dirk Meyer
- Physikalische Organische
Chemie, TU Dresden, Bergstrasse
66, 01069 Dresden, Germany
| | - Thomas Strassner
- Physikalische Organische
Chemie, TU Dresden, Bergstrasse
66, 01069 Dresden, Germany
| |
Collapse
|
28
|
Nizovtsev AS. Activation of C-H bond in methane by Pd atom from the bonding evolution theory perspective. J Comput Chem 2013; 34:1917-24. [DOI: 10.1002/jcc.23345] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2013] [Revised: 05/08/2013] [Accepted: 05/10/2013] [Indexed: 12/31/2022]
Affiliation(s)
- Anton S. Nizovtsev
- Nikolaev Institute of Inorganic Chemistry; Siberian Branch of the Russian Academy of Sciences; Academician Lavrentiev Avenue 3; Novosibirsk; 630090; Russian Federation
| |
Collapse
|
29
|
Steinmetz M, Grimme S. Benchmark study of the performance of density functional theory for bond activations with (ni,pd)-based transition-metal catalysts. ChemistryOpen 2013; 2:115-24. [PMID: 24551548 PMCID: PMC3703816 DOI: 10.1002/open.201300012] [Citation(s) in RCA: 125] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2013] [Indexed: 11/25/2022] Open
Abstract
The performance of 23 density functionals, including one LDA, four GGAs, three meta-GGAs, three hybrid GGAs, eight hybrid meta-GGAs, and ten double-hybrid functionals, was investigated for the computation of activation energies of various covalent main-group single bonds by four catalysts: Pd, PdCl−, PdCl2, and Ni (all in the singlet state). A reactant complex, the barrier, and reaction energy were considered, leading to 164 energy data points for statistical analysis. Extended Gaussian AO basis sets were used in all calculations. The best functional for the complete benchmark set relative to estimated CCSD(T)/CBS reference data is PBE0-D3, with an MAD value of 1.1 kcal mol−1 followed by PW6B95-D3, the double hybrid PWPB95-D3, and B3LYP-D3 (1.9 kcal mol−1 each). The other tested hybrid meta-GGAs perform less well (M06-HF: 7.0 kcal mol−1; M06-2X: 6.3 kcal mol−1; M06: 4.9 kcal mol−1) for the investigated reactions. In the Ni case, some double hybrids show larger errors due to partial breakdown of the perturbative treatment for the correlation energy in cases with difficult electronic structures (partial multi-reference character). Only double hybrids either with very low amounts of perturbative correlation (e.g., PBE0-DH) or that use the opposite-spin correlation component only (e.g., PWPB95) seem to be more robust. We also investigated the effect of the D3 dispersion correction. While the barriers are not affected by this correction, significant and mostly positive results were observed for reaction energies. Furthermore, six very recently proposed double-hybrid functionals were analyzed regarding the influence of the amount of Fock exchange as well as the type of perturbative correlation treatment. According to these results, double hybrids with <50–60 % of exact exchange and ∼30 % perturbative correlation perform best.
Collapse
Affiliation(s)
- Marc Steinmetz
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn Beringstr. 4, 53115 Bonn (Germany)
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn Beringstr. 4, 53115 Bonn (Germany)
| |
Collapse
|
30
|
Abstract
We have investigated the molecular geometries of a series of dicoordinated d(10)-transition-metal complexes ML2 (M=Co(-), Rh(-), Ir(-), Ni, Pd, Pt, Cu(+), Ag(+), Au(+); L=NH3, PH3, CO) using relativistic density functional theory (DFT) at ZORA-BLYP/TZ2P. Not all complexes have the expected linear ligand-metal-ligand (L-M-L) angle: this angle varies from 180° to 128.6° as a function of the metal as well as the ligands. Our main objective is to present a detailed explanation why ML2 complexes can become bent. To this end, we have analyzed the bonding mechanism in ML2 as a function of the L-M-L angle using quantitative Kohn-Sham molecular orbital (MO) theory in combination with an energy decomposition analysis (EDA) scheme. The origin of bent L-M-L structures is π backdonation. In situations of strong π backdonation, smaller angles increase the overlap of the ligand's acceptor orbital with a higher-energy donor orbital on the metal-ligand fragment, and therefore favor π backdonation, resulting in additional stabilization. The angle of the complexes thus depends on the balance between this additional stabilization and increased steric repulsion that occurs as the complexes are bent.
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) ; Institute for Molecules and Materials, Radboud University Nijmegen Heyendaalseweg 135, 6525 AJ Nijmegen (The Netherlands)
| |
Collapse
|
31
|
Seth M, Ziegler T, Steinmetz M, Grimme S. Modeling Transition Metal Reactions with Range-Separated Functionals. J Chem Theory Comput 2013; 9:2286-99. [DOI: 10.1021/ct301112m] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Michael Seth
- Department of Chemistry, University of Calgary, University Drive 2500, Calgary,
AB T2N-1N4, Canada
| | - Tom Ziegler
- Department of Chemistry, University of Calgary, University Drive 2500, Calgary,
AB T2N-1N4, Canada
| | - Marc Steinmetz
- Mulliken
Center for Theoretical Chemistry, Institute for Physical and Theoretical
Chemistry, Universität Bonn, Beringstrasse
4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken
Center for Theoretical Chemistry, Institute for Physical and Theoretical
Chemistry, Universität Bonn, Beringstrasse
4, D-53115 Bonn, Germany
| |
Collapse
|
32
|
Zarzycki B, Zell T, Schmidt D, Radius U. P4Activation at Ni0: Selective Formation of an NHC-Stabilized, Dinuclear Nickel Complex [Ni2(iPr2Im)4(μ,η2:2-P2)]. Eur J Inorg Chem 2013. [DOI: 10.1002/ejic.201201368] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
33
|
Zarzycki B, Bickelhaupt FM, Radius U. Symmetrical P4 cleavage at cobalt half sandwich complexes [(η5-C5H5)Co(L)] (L = CO, NHC) – a computational case study on the mechanism of symmetrical P4 degradation to P2 ligands. Dalton Trans 2013; 42:7468-81. [DOI: 10.1039/c3dt50267h] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
Sanhueza IA, Wagner AM, Sanford MS, Schoenebeck F. On the role of anionic ligands in the site-selectivity of oxidative C–H functionalization reactions of arenes. Chem Sci 2013. [DOI: 10.1039/c3sc00017f] [Citation(s) in RCA: 81] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
35
|
Yang HQ, Hu CW, Gao C, Yang MY, Li FM, Li CQ, Li XY. Theoretical study on the gas-phase reaction mechanism between palladium monoxide and methane. J Comput Chem 2011; 32:3440-55. [PMID: 21919016 DOI: 10.1002/jcc.21926] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2011] [Revised: 07/16/2011] [Accepted: 07/27/2011] [Indexed: 11/06/2022]
Abstract
The gas-phase reaction mechanism between palladium monoxide and methane has been theoretically investigated on the singlet and triplet state potential energy surfaces (PESs) at the CCSD(T)/AVTZ//B3LYP/6-311+G(2d, 2p), SDD level. The major reaction channel leads to the products PdCH(2) + H(2)O, whereas the minor channel results in the products Pd + CH(3)OH, CH(2)OPd + H(2), and PdOH + CH(3). The minimum energy reaction pathway for the formation of main products (PdCH(2) + H(2)O), involving one spin inversion, prefers to start at the triplet state PES and afterward proceed along the singlet state PES, where both CH(3)PdOH and CH(3)Pd(O)H are the critical intermediates. Furthermore, the rate-determining step is RS-CH(3) PdOH → RS-2-TS1cb → RS-CH(2)Pd(H)OH with the rate constant of k = 1.48 × 10(12) exp(-93,930/RT). For the first C-H bond cleavage, both the activation strain ΔE(≠)(strain) and the stabilizing interaction ΔE(≠)(int) affect the activation energy ΔE(≠), with ΔE(≠)(int) in favor of the direct oxidative insertion. On the other hand, in the PdCH(2) + H(2) O reaction, the main products are Pd + CH(3)OH, and CH(3)PdOH is the energetically preferred intermediate. In the CH(2)OPd + H(2) reaction, the main products are Pd + CH(3)OH with the energetically preferred intermediate H(2)PdOCH(2). In the Pd + CH(3)OH reaction, the main products are CH(2)OPd + H(2), and H(2)PdOCH(2) is the energetically predominant intermediate. The intermediates, PdCH(2), H(2) PdCO, and t-HPdCHO are energetically preferred in the PdC + H(2), PdCO + H(2), and H(2)Pd + CO reactions, respectively. Besides, PdO toward methane activation exhibits higher reaction efficiency than the atom Pd and its first-row congener NiO.
Collapse
Affiliation(s)
- Hua-Qing Yang
- College of Chemical Engineering, Sichuan University, Chengdu, Sichuan 610065, People's Republic of China
| | | | | | | | | | | | | |
Collapse
|
36
|
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
|
37
|
|
38
|
Zhang X, Schwarz H. Bonding in cationic MOH n + (M = K − La, Hf − Rn; n = 0–2): DFT performances and periodic trends. Theor Chem Acc 2010. [DOI: 10.1007/s00214-010-0861-0] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
39
|
Flener-Lovitt C, Woon DE, Dunning TH, Girolami GS. A DFT and ab Initio Benchmarking Study of Metal−Alkane Interactions and the Activation of Carbon−Hydrogen Bonds. J Phys Chem A 2009; 114:1843-51. [DOI: 10.1021/jp9058033] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Charity Flener-Lovitt
- School of Chemical Sciences, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801
| | - David E. Woon
- School of Chemical Sciences, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801
| | - Thom H. Dunning
- School of Chemical Sciences, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801
| | - Gregory S. Girolami
- School of Chemical Sciences, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801
| |
Collapse
|
40
|
Osuna SÃ, Swart M, Baerends EJ, Bickelhaupt FM, Solà M. Homolytic versus Heterolytic Dissociation of Alkalimetal Halides: The Effect of Microsolvation. Chemphyschem 2009; 10:2955-65. [DOI: 10.1002/cphc.200900480] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
41
|
Hicks JD, Hyde AM, Cuezva AM, Buchwald SL. Pd-catalyzed N-arylation of secondary acyclic amides: catalyst development, scope, and computational study. J Am Chem Soc 2009; 131:16720-34. [PMID: 19886610 PMCID: PMC2805443 DOI: 10.1021/ja9044357] [Citation(s) in RCA: 191] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
We report the efficient N-arylation of acyclic secondary amides and related nucleophiles with aryl nonaflates, triflates, and chlorides. This method allows for easy variation of the aromatic component in tertiary aryl amides. A new biaryl phosphine with P-bound 3,5-(bis)trifluoromethylphenyl groups was found to be uniquely effective for this amidation. The critical aspects of the ligand were explored through synthetic, mechanistic, and computational studies. Systematic variation of the ligand revealed the importance of (1) a methoxy group on the aromatic carbon of the "top ring" ortho to the phosphorus and (2) two highly electron-withdrawing P-bound 3,5-(bis)trifluoromethylphenyl groups. Computational studies suggest the electron-deficient nature of the ligand is important in facilitating amide binding to the LPd(II)(Ph)(X) intermediate.
Collapse
Affiliation(s)
- Jacqueline D. Hicks
- Department of Chemistry, Massachusetts Institute of Technology, 77
Massachusetts Avenue, Cambridge, MA 02139
| | | | | | - Stephen L. Buchwald
- Department of Chemistry, Massachusetts Institute of Technology, 77
Massachusetts Avenue, Cambridge, MA 02139
| |
Collapse
|
42
|
Abstract
We have computationally explored C–X bond activation by the group-11 transition-metal cations Cu+, Ag+, and Au+, and, for comparison, Pd, using relativistic density functional theory (DFT) at ZORA-BLYP/TZ2P. Oxidative insertion of the second-row transition-metal species Ag+ and Pd leads, for a given bond, to the highest overall reaction barriers. On the other hand, if we compare the different bonds oxidative insertion into the C–F bond is associated with (one of the) highest overall barriers whereas insertion into the C–Cl bond leads to the lowest overall barrier for any transition metal. The main trends in reactivity are rationalized using the activation strain model of chemical reactivity, which is an extension of the fragment approach to reaction profiles. In this model, the shape of the reaction profile ΔE(ζ) and the height of the overall reaction barrier ΔE≠ = ΔE(ζ=ζTS) are interpreted in terms of the strain energy ΔEstrain(ζ) associated with deforming the reactants along the reaction coordinate ζ plus the interaction energy ΔEint(ζ) between these deformed reactants: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ).
Collapse
Affiliation(s)
- G. Theodoor de Jong
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - 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
| |
Collapse
|
43
|
|
44
|
Cypryk M, Pospiech P, Strzelec K, Sobczak JW. Soluble Alkylthiopolysiloxane-Supported Palladium Catalysts for the Heck Reaction. PHOSPHORUS SULFUR 2009. [DOI: 10.1080/10426500902947930] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Marek Cypryk
- a Center of Molecular and Macromolecular Studies , Polish Academy of Sciences , Łódź, Poland
| | - Piotr Pospiech
- a Center of Molecular and Macromolecular Studies , Polish Academy of Sciences , Łódź, Poland
| | | | - Janusz W. Sobczak
- c Institute of Physical Chemistry, Polish Academy of Sciences , Warsaw, Poland
| |
Collapse
|
45
|
Effect of Basis Set Superposition Error on the MP2 Relative Energies of Gold Cluster Au6. B KOREAN CHEM SOC 2009. [DOI: 10.5012/bkcs.2009.30.4.794] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
46
|
Bento AP, Bickelhaupt FM. Frontside versus Backside S(N)2 substitution at group 14 atoms: origin of reaction barriers and reasons for their absence. Chem Asian J 2008; 3:1783-92. [PMID: 18712744 DOI: 10.1002/asia.200800065] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We have theoretically studied the gas-phase nucleophilic substitution at group-14 atoms (S(N)2@A) in the model reactions of Cl(-)+AH(3)Cl (A=C, Si, Ge, Sn, and Pb) using relativistic density functional theory (DFT) at ZORA-OLYP/TZ2P. Firstly, we wish to explore and understand how the reaction coordinate zeta, and potential energy surfaces (PES) along zeta, vary as the center of nucleophilic attack changes from carbon to the heavier group-14 atoms. Secondly, a comparison between the more common backside reaction (S(N)2-b) and the frontside pathway (S(N)2-f) is performed. The S(N)2-b reaction is found to have a central barrier for A=C, but none for the other group-14 atoms, A=Si-Pb. Relativistic effects destabilize reactant complexes and transition species by up to 10 kcal mol(-1) (for S(N)2-f@Pb), but they do not change relative heights of barriers. We also address the nature of the transformation in the frontside S(N)2-f reactions in terms of turnstile rotation versus Berry-pseudorotation mechanism.
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, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
47
|
Vicent C, Feliz M, Llusar R. Intrinsic Gas-Phase Reactivity toward Methanol of Trinuclear Tungsten W3S4 Complexes Bearing W−X (X = Br, OH) Groups. J Phys Chem A 2008; 112:12550-8. [DOI: 10.1021/jp804263q] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Cristian Vicent
- Serveis Centrals d′Instrumentació Cientifica and Departament de Química Física i Analítica, Universitat Jaume I, Av. Sos Baynat s/n, 12071 Castelló, Spain
| | - Marta Feliz
- Serveis Centrals d′Instrumentació Cientifica and Departament de Química Física i Analítica, Universitat Jaume I, Av. Sos Baynat s/n, 12071 Castelló, Spain
| | - Rosa Llusar
- Serveis Centrals d′Instrumentació Cientifica and Departament de Química Física i Analítica, Universitat Jaume I, Av. Sos Baynat s/n, 12071 Castelló, Spain
| |
Collapse
|
48
|
Orian L, Zeist WJV, Bickelhaupt FM. Linkage Isomerism of Nitriles in Rhodium Half-Sandwich Metallacycles. Organometallics 2008. [DOI: 10.1021/om8004614] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Laura Orian
- Dip. Scienze Chimiche Università degli Studi di Padova, Via Marzolo 1, 35129 Padova, Italy, and Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Willem-Jan van Zeist
- Dip. Scienze Chimiche Università degli Studi di Padova, Via Marzolo 1, 35129 Padova, Italy, and Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - F. Matthias Bickelhaupt
- Dip. Scienze Chimiche Università degli Studi di Padova, Via Marzolo 1, 35129 Padova, Italy, and Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
49
|
Tsai HHG, Chung MW, Chou YK, Hou DR. Interplay of Hydrogenation and Dehydrogenation in Isoindoline and Indoline Isomers: A Density Functional Theory Study. J Phys Chem A 2008; 112:5278-85. [DOI: 10.1021/jp7121246] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hui-Hsu Gavin Tsai
- Department of Chemistry, National Central University, Jhong-Li City Tao-Yuan County 32001, Taiwan
| | - Ming-Weng Chung
- Department of Chemistry, National Central University, Jhong-Li City Tao-Yuan County 32001, Taiwan
| | - Yi-Kang Chou
- Department of Chemistry, National Central University, Jhong-Li City Tao-Yuan County 32001, Taiwan
| | - Duen-Ren Hou
- Department of Chemistry, National Central University, Jhong-Li City Tao-Yuan County 32001, Taiwan
| |
Collapse
|
50
|
Guo Z, Ke Z, Phillips DL, Zhao C. Intrinsic Reaction Coordinate Analysis of the Activation of CH4 by Molybdenum Atoms: A Density Functional Theory Study of the Crossing Seams of the Potential Energy Surfaces. Organometallics 2007. [DOI: 10.1021/om7007452] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zhen Guo
- Department of Chemistry, The University of Hong Kong, Pokfulam Road, Hong Kong, P. R. China, and School of Chemistry and Chemical Engineering, Sun Yat-Sen University, Guangzhou 510275, P.R. China
| | - Zhuofeng Ke
- Department of Chemistry, The University of Hong Kong, Pokfulam Road, Hong Kong, P. R. China, and School of Chemistry and Chemical Engineering, Sun Yat-Sen University, Guangzhou 510275, P.R. China
| | - David Lee Phillips
- Department of Chemistry, The University of Hong Kong, Pokfulam Road, Hong Kong, P. R. China, and School of Chemistry and Chemical Engineering, Sun Yat-Sen University, Guangzhou 510275, P.R. China
| | - Cunyuan Zhao
- Department of Chemistry, The University of Hong Kong, Pokfulam Road, Hong Kong, P. R. China, and School of Chemistry and Chemical Engineering, Sun Yat-Sen University, Guangzhou 510275, P.R. China
| |
Collapse
|