How the Lewis Base F– Catalyzes the 1,3-Dipolar Cycloaddition between Carbon Dioxide and Nitrilimines
- Dennis SvatunekDennis SvatunekDepartment of Theoretical Chemistry, Amsterdam Institute of Molecular and Life Sciences (AIMSS), Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The NetherlandsInstitute of Applied Synthetic Chemistry, TU Wien (Vienna University of Technology), A-1060, Vienna, AustriaDepartment of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, Los Angeles, United StatesMore by Dennis Svatunek
- ,
- Thomas HansenThomas HansenDepartment of Theoretical Chemistry, Amsterdam Institute of Molecular and Life Sciences (AIMSS), Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The NetherlandsLeiden Institute of Chemistry, Leiden University, Einsteinweg 55, 2333 CC Leiden, The NetherlandsMore by Thomas Hansen
- ,
- Kendall N. HoukKendall N. HoukDepartment of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, Los Angeles, United StatesMore by Kendall N. Houk
- , and
- Trevor A. Hamlin*Trevor A. Hamlin*E-mail: [email protected]Department of Theoretical Chemistry, Amsterdam Institute of Molecular and Life Sciences (AIMSS), Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The NetherlandsMore by Trevor A. Hamlin
Abstract

The mechanism of the Lewis base F– catalyzed 1,3-dipolar cycloaddition between CO2 and nitrilimines is interrogated using DFT calculations. F– activates the nitrilimine, not CO2 as proposed in the literature, and imparts a significant rate enhancement for the cycloaddition. The origin of this catalysis is in the strength of the primary orbital interactions between the reactants. The Lewis base activated nitrilimine–F– has high-lying filled FMOs. The smaller FMO-LUMO gap promotes a rapid nucleophilic attack and overall cycloaddition with CO2.
This publication is licensed under
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
Scheme 1

Scheme 2

Figure 1

Figure 1. (a) Reaction profiles (ΔGtoluene in kcal mol–1) for the three reaction pathways of the studied 1,3-dipolar cycloaddition between carbon dioxide and nitrilimine 1, including (i) uncatalyzed cycloaddition (black), (ii) formation of CO2F– followed by stepwise cycloaddition (blue), and (iii) formation of 2 followed by stepwise cycloaddition (red), computed at SMD(toluene)-M06-2X-D3/def2-TZVP; TS = transition state and INT = intermediate. (b) Key transition state structures with key bond lengths (in Å) for the three reaction pathways.
Figure 2

Figure 2. (a) Distortion/interaction–activation strain model analysis; and (b) energy decomposition analysis of the cycloaddition reaction of 1 (black) and 2 (red) with CO2 (transition states indicated with a dot). (13) (c) Frontier molecular orbital diagram of the most important FMO-LUMOCO2 orbital interaction with the calculated energy gaps, orbital overlaps, and the S2/Δε terms, at consistent geometries with a N···C bond length of 2.14 Å. Computed at SMD(toluene)-M06-2X-D3/def2-TZVP using autoDIAS (14) for (a) and M06-2X-D3/TZ2P//SMD(toluene)-M06-2X-D3/def2-TZVP using PyFrag (15) for (b,c).
Scheme 3

Methods
Computational Details
Distortion/Interaction–Activation Strain and Energy Decomposition Analysis


Voronoi Deformation Density

Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.joc.0c02963.
Additional computational results; and Cartesian coordinates, energies, and number of imaginary frequencies of all stationary points (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.
Acknowledgments
D.S. is grateful to the Christiana Hörbiger Award (TU Wien) and the “Hochschuljubiläumsstiftung der Stadt Wien” (H-331849/2018) for financial support. K.N.H. is grateful to the National Science Foundation of the US for financial support of this research (CHE-1764328). Quantum chemical calculations were performed on the Vienna Scientific Cluster (Austria), the Hoffman2 Cluster (UCLA Institute for Digital Research and Education), and the Cartesius supercomputer (SURFsara, Amsterdam). The authors also acknowledge Pascal Vermeeren of the Vrije Universiteit Amsterdam for fruitful discussions.
References
This article references 25 other publications.
- 1(a) Milani, B.; Licini, G.; Clot, E.; Albrecht, M. Small Molecule Activation. Dalton Trans. 2016, 45, 14419– 14420, DOI: 10.1039/C6DT90140A[Crossref], [PubMed], [CAS], Google Scholar.1ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFWisb3F&md5=806efb4c52eec36651c19c47687080faSmall molecule activationMilani, Barbara; Licini, Giulia; Clot, Eric; Albrecht, MartinDalton Transactions (2016), 45 (37), 14419-14420CODEN: DTARAF; ISSN:1477-9226. (Royal Society of Chemistry)There is no expanded citation for this reference.(b) Zhang, W.; Lai, W.; Cao, R. Energy-Related Small Molecule Activation Reactions: Oxygen Reduction and Hydrogen and Oxygen Evolution Reactions Catalyzed by Porphyrin- and Corrole-Based Systems. Chem. Rev. 2017, 117, 3717– 3797, DOI: 10.1021/acs.chemrev.6b00299[ACS Full Text.
], [CAS], Google Scholar
1bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xhsl2gtbnO&md5=04523ae018db1238aece4e24e7012ef3Energy-Related Small Molecule Activation Reactions: Oxygen Reduction and Hydrogen and Oxygen Evolution Reactions Catalyzed by Porphyrin- and Corrole-Based SystemsZhang, Wei; Lai, Wenzhen; Cao, RuiChemical Reviews (Washington, DC, United States) (2017), 117 (4), 3717-3797CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)Globally increasing energy demands and environmental concerns related to the use of fossil fuels have stimulated extensive research to identify new energy systems and economies that are sustainable, clean, low cost, and environmentally benign. Hydrogen generation from solar-driven water splitting is a promising strategy to store solar energy in chem. bonds. The subsequent combustion of hydrogen in fuel cells produces elec. energy, and the only exhaust is water. These two reactions compose an ideal process to provide clean and sustainable energy. In such a process, a hydrogen evolution reaction (HER), an oxygen evolution reaction (OER) during water splitting, and an oxygen redn. reaction (ORR) as a fuel cell cathodic reaction are key steps that affect the efficiency of the overall energy conversion. Catalysts play key roles in this process by improving the kinetics of these reactions. Porphyrin-based and corrole-based systems are versatile and can efficiently catalyze the ORR, OER, and HER. Because of the significance of energy-related small mol. activation, this review covers recent progress in hydrogen evolution, oxygen evolution, and oxygen redn. reactions catalyzed by porphyrins and corroles.(c) Gómez, J. E.; Kleij, A. W. Catalytic Nonreductive Valorization of Carbon Dioxide into Fine Chemicals. In Adv. Organomet. Chem.; Pérez, P. J., Ed.; Academic Press: 2019; Vol. 71, pp 175– 226.(d) Chen, L.; Chen, G.; Leung, C.-F.; Cometto, C.; Robert, M.; Lau, T.-C. Molecular Quaterpyridine-Based Metal Complexes for Small Molecule Activation: Water Splitting and CO2 Reduction. Chem. Soc. Rev. 2020, 49, 7271– 7283, DOI: 10.1039/D0CS00927J[Crossref], [PubMed], [CAS], Google Scholar1dhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhvVKltrjF&md5=e98173e996d0191dcc4a7396e2ea48c5Molecular quaterpyridine-based metal complexes for small molecule activation: water splitting and CO2 reductionChen, Lingjing; Chen, Gui; Leung, Chi-Fai; Cometto, Claudio; Robert, Marc; Lau, Tai-ChuChemical Society Reviews (2020), 49 (20), 7271-7283CODEN: CSRVBR; ISSN:0306-0012. (Royal Society of Chemistry)A review. Artificial photosynthesis is considered as one of the most promising strategies for solar-to-fuel conversion through sunlight-driven water splitting and CO2 redn. This tutorial describes recent developments in the use of metal quaterpyridine complexes as electrocatalyts and photocatalysts for artificial photosynthesis. - 2(a) Liu, Q.; Wu, L.; Jackstell, R.; Beller, M. Using Carbon Dioxide as a Building Block in Organic Synthesis. Nat. Commun. 2015, 6, 5933, DOI: 10.1038/ncomms6933[Crossref], [PubMed], [CAS], Google Scholar.2ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhtF2lurjM&md5=e3ba6b5612253c0c6ebec3d937f137f8Using carbon dioxide as a building block in organic synthesisLiu, Qiang; Wu, Lipeng; Jackstell, Ralf; Beller, MatthiasNature Communications (2015), 6 (), 5933pp.CODEN: NCAOBW; ISSN:2041-1723. (Nature Publishing Group)A review. The most recent advances made in the area of CO2 valorization-turning CO2 into a useful chem. feedstock-under mild conditions has been reviewed. A special focus is given on the reaction modes for the CO2 activation and its application as C1 building block in org. synthesis. The following subjects will be addressed in this review : (1) novel transformations using carbon dioxide (briefly summarized); (2) different reaction modes for CO2 activation (main focus of this review); and (3) potential new applications of CO2 valorization.(b) Vogt, C.; Monai, M.; Sterk, E. B.; Palle, J.; Melcherts, A. E. M.; Zijlstra, B.; Groeneveld, E.; Berben, P. H.; Boereboom, J. M.; Hensen, E. J. M.; Meirer, F.; Filot, I. A. W.; Weckhuysen, B. M. Understanding Carbon Dioxide Activation and Carbon-Carbon Coupling over Nickel. Nat. Commun. 2019, 10, 5330, DOI: 10.1038/s41467-019-12858-3[Crossref], [PubMed], [CAS], Google Scholar.2bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3Mfjt1ajtg%253D%253D&md5=4fee60b8521239e02af56b1978b41b28Understanding carbon dioxide activation and carbon-carbon coupling over nickelVogt Charlotte; Monai Matteo; Sterk Ellen B; Palle Jonas; Melcherts Angela E M; Boereboom Jelle M; Meirer Florian; Weckhuysen Bert M; Zijlstra Bart; Hensen Emiel J M; Filot Ivo A W; Groeneveld Esther; Berben Peter HNature communications (2019), 10 (1), 5330 ISSN:.Carbon dioxide is a desired feedstock for platform molecules, such as carbon monoxide or higher hydrocarbons, from which we will be able to make many different useful, value-added chemicals. Its catalytic hydrogenation over abundant metals requires the amalgamation of theoretical knowledge with materials design. Here we leverage a theoretical understanding of structure sensitivity, along with a library of different supports, to tune the selectivity of methanation in the Power-to-Gas concept over nickel. For example, we show that carbon dioxide hydrogenation over nickel can and does form propane, and that activity and selectivity can be tuned by supporting different nickel particle sizes on various oxides. This theoretical and experimental toolbox is not only useful for the highly selective production of methane, but also provides new insights for carbon dioxide activation and subsequent carbon-carbon coupling towards value-added products thereby reducing the deleterious effects of this environmentally harmful molecule.(c) Schilling, W.; Das, S. Transition Metal-Free Synthesis of Carbamates Using CO2 as the Carbon Source. ChemSusChem 2020, 13, 6246– 6258[PubMed], [CAS], Google Scholar2chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXitFKjtrrJ&md5=a6465c2b7cba8c074554ac90f2ce1172Transition Metal-Free Synthesis of Carbamates Using CO2 as the Carbon SourceSchilling, Waldemar; Das, ShoubhikChemSusChem (2020), 13 (23), 6246-6258CODEN: CHEMIZ; ISSN:1864-5631. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. Utilization of carbon dioxide as a C1 synthon is highly attractive for the synthesis of valuable chems. However, activation of CO2 is highly challenging, owing to its thermodn. stability and kinetic inertness. With this in mind, several strategies have been developed for the generation of carbon-heteroatom bonds. Among these, formation of C-N bonds is highly attractive, esp., when carbamates can be synthesized directly from CO2. This Minireview focuses on transition metal-free approaches for the fixation of CO2 to generate carbamates for the prodn. of fine chems. and pharmaceuticals. Within the past decade, transition metal-free approaches have gained increasing attention, but traditional reviews have rarely focused on these approaches. Direct comparisons between such methods have been even more scarce. This Minireview seeks to address this discrepancy.
- 3Guo, C.-X.; Zhang, W.-Z.; Zhang, N.; Lu, X.-B. 1,3-Dipolar Cycloaddition of Nitrile Imine with Carbon Dioxide: Access to 1,3,4-Oxadiazole-2(3H)-ones. J. Org. Chem. 2017, 82, 7637, DOI: 10.1021/acs.joc.7b00963[ACS Full Text
], [CAS], Google Scholar
3https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVGgt73L&md5=0d77612903d20b3a9087cefd6e7096831,3-Dipolar Cycloaddition of Nitrile Imine with Carbon Dioxide: Access to 1,3,4-Oxadiazole-2(3H)-onesGuo, Chun-Xiao; Zhang, Wen-Zhen; Zhang, Ning; Lu, Xiao-BingJournal of Organic Chemistry (2017), 82 (14), 7637-7642CODEN: JOCEAH; ISSN:0022-3263. (American Chemical Society)Efficient synthesis of 1,3,4-oxadiazole-2(3H)-one was achieved by CsF/18-crown-6 mediated 1,3-dipolar cycloaddn. of nitrile imine and 2.0 MPa of CO2. CsF/18-crown-6 played a key role in enhancing the reactivity of CO2 as a 1,3-dipolarophile. The practical utility of this transition-metal-free approach to 1,3,4-oxadiazole-2(3H)-one is highlighted by the convenient synthesis of a com. herbicide Oxadiazon and a MAO B inhibitor. - 4(a) Zhang, X.; Gross, U.; Seppelt, K. Fluorocarbonate, [FCO2]−: Preparation and Structure. Angew. Chem., Int. Ed. Engl. 1995, 34, 1858– 1860, DOI: 10.1002/anie.199518581[Crossref], [CAS], Google Scholar.4ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXotF2ku78%253D&md5=1cb058cd84cf7b7967dca42bc0b48a22Fluorocarbonate, [FCO2]-: preparation and structureZhang, Xiongzhi; Gross, Udo; Seppelt, KonradAngewandte Chemie, International Edition in English (1995), 34 (17), 1858-60CODEN: ACIEAY; ISSN:0570-0833. (VCH)CO2 reacts with F- salts of Me4N+ or Me3(Me3CCH2)N+ or 1,1,3,3,5,5-hexamethylpiperidinium (pip+) to give [FCO2]- salts. Elemental anal. of pip+[FCO2] and NMR of Me4N[FCO2] confirm the identity of the anion. For the latter, multinuclear NMR, IR and Raman data are given. Ab initio calcns. for enthalpy of reaction of F- with CO, N2O, or H2C:C:CH2 indicate the fluoride affinity of CO2 may be the approx. lower limit for a reaction.(b) Arnold, D. W.; Bradforth, S. E.; Kim, E. H.; Neumark, D. M. Study of halogen-carbon dioxide clusters and the fluoroformyloxyl radical by photodetachment of X–(CO2) (X= I, Cl, Br) and FCO2–. J. Chem. Phys. 1995, 102, 3493, DOI: 10.1063/1.468575[Crossref], [CAS], Google Scholar4bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXjvFGqtbY%253D&md5=2d1eeb7b405437e16aa326c4d58c8ed9Study of halogen-carbon dioxide clusters and the fluoroformyloxyl radical by photodetachment of X-(CO2) (X = I,Cl,Br) and FCO-2Arnold, Don W.; Bradforth, Stephen E.; Kim, Eun H.; Neumark, Daniel M.Journal of Chemical Physics (1995), 102 (9), 3493-509CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Photoelectron spectra have been measured for the anions X-(CO2), with X = I, Br, Cl, and F. The vibrationally resolved spectra show that I-(CO2), Br-(CO2), and Cl-(CO2) are primarily electrostatically bound clusters, although the charge-quadrupole interaction is strong enough to distort the CO2 mol. by as much as 10° [in Cl-(CO2)]. Ab initio calcns. and electrostatic models are used to describe the geometry and bonding of these clusters. The photoelectron spectrum of FCO2- is qual. different and shows transitions to both the ~X 2B2 ground and the ~A2A2 first excited electronic states of the covalently bound FCO2 radical. The previously unobserved ~A 2A2 state is measured to lie 0.579 eV above the ground state. Vibrational frequencies are assigned with the assistance of ab initio calcns. The FCO2 heat of formation is detd. to be ΔfH0298(FCO2)= -85.2±2.8 kcal/mol. While both FCO2- and FCO2 are more strongly bound than the other halide-CO2 clusters, the C-F bonds are very weak relative to C-F bounds found in other halocarbon compds.
- 5(a) Ukai, K.; Aoki, M.; Takaya, J.; Iwasawa, N. Rhodium(I)-Catalyzed Carboxylation of Aryl- and Alkenylboronic Esters with CO2. J. Am. Chem. Soc. 2006, 128, 8706, DOI: 10.1021/ja061232m[ACS Full Text.
], [CAS], Google Scholar
5ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XmtVSjsbo%253D&md5=0c880ab4e3f61a8cf0f5e6ddab46953cRhodium(I)-Catalyzed Carboxylation of Aryl- and Alkenylboronic Esters with CO2Ukai, Kazutoshi; Aoki, Masao; Takaya, Jun; Iwasawa, NobuharuJournal of the American Chemical Society (2006), 128 (27), 8706-8707CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)When the esters of arylboronic acids with 2,2-dimethyl-1,3-propanol were treated with a catalytic amt. of [Rh(OH)(cod)]2 in the presence of 1,3-bis(diphenylphosphino)propane and CsF in dioxane at 60 °C under carbon dioxide atm., the benzoic acid derivs. were obtained in good yields. Reactions of alkenylboronic esters also proceeded under similar conditions to give α,β-unsatd. carboxylic acids. As these boronic esters are now easily available through coupling or direct borylation reactions, this method would be a useful method for the prepn. of various functionalized aryl- and alkenyl carboxylic acids. For example, the rhodium-catalyzed reaction of a boronic acid ester, i.e., 5,5-dimethyl-2-phenyl-1,3,2-dioxaborinane, with carbon dioxide gave benzoic acid.(b) Li, S.; Ma, S. Quadri-Synergetic Effect for Highly Effective Carbon Dioxide Fixation and Its Application to Indoloquinolinone. Adv. Synth. Catal. 2012, 354, 2387– 2394, DOI: 10.1002/adsc.201200469[Crossref], [CAS], Google Scholar.5bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38Xht1Glu7fE&md5=b15476302f7b12a0de48cb35cefd67a2Quadri-Synergetic Effect for Highly Effective Carbon Dioxide Fixation and Its Application to IndoloquinolinoneLi, Suhua; Ma, ShengmingAdvanced Synthesis & Catalysis (2012), 354 (13), 2387-2394, S2387/1-S2387/65CODEN: ASCAF7; ISSN:1615-4150. (Wiley-VCH Verlag GmbH & Co. KGaA)The first copper-catalyzed cyclic anti-nucleometallation-carboxylation of 2-alkynylanilines with carbon dioxide in the presence of dimethylzinc (ZnMe2) and cesium fluoride (CsF) for the effective synthesis of indolyl-3-carboxylic acids and indolodihydropyran-2-one is described. Through a mechanistic study, it is unveiled that the metal ions Cu2+, Zn2+, Cs+ and F- are working together for this CO2-based highly efficient carboxylation.(c) Takaya, J.; Tadami, S.; Ukai, K.; Iwasawa, N. Copper (I)-catalyzed carboxylation of aryl-and alkenylboronic esters. Org. Lett. 2008, 10, 2697– 2700, DOI: 10.1021/ol800829q[ACS Full Text], [CAS], Google Scholar
5chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXmsVKhtbg%253D&md5=80cfe34f459e75c4618d72e8245dac2eCopper(I)-Catalyzed Carboxylation of Aryl- and Alkenylboronic EstersTakaya, Jun; Tadami, Satoshi; Ukai, Kazutoshi; Iwasawa, NobuharuOrganic Letters (2008), 10 (13), 2697-2700CODEN: ORLEF7; ISSN:1523-7060. (American Chemical Society)The copper(I)-catalyzed carboxylation reaction of aryl- and alkenylboronic esters proceeded smoothly under CO2 to give the corresponding carboxylic acids in good yield. This reaction showed wide generality with higher functional group tolerance compared to the corresponding Rh(I)-catalyzed reaction. - 6Denmark, S. E.; Beutner, G. L. Lewis Base Catalysis in Organic Synthesis. Angew. Chem., Int. Ed. 2008, 47, 1560– 1638, DOI: 10.1002/anie.200604943[Crossref], [CAS], Google Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXis1Krtr8%253D&md5=234cf65fa06fdaaabd490440bd2d7ce4Lewis base catalysis in organic synthesisDenmark, Scott E.; Beutner, Gregory L.Angewandte Chemie, International Edition (2008), 47 (9), 1560-1638CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. The legacy of Gilbert Newton Lewis (1875-1946) pervades the lexicon of chem. bonding and reactivity; the power of his concept of donor-acceptor bonding is evident in the eponymous foundations of electron-pair acceptors (Lewis acids) and donors (Lewis bases). His discovery ushered in the use of Lewis acids as reagents and catalysts for org. reactions. However, in recent years, the recognition that Lewis bases can also serve in this capacity has grown enormously. Most importantly, it has become increasingly apparent that the behavior of Lewis bases as agents for promoting chem. reactions is not merely as an electronic complement of the cognate Lewis acids: in fact Lewis bases are capable of enhancing both the electrophilic and nucleophilic character of mols. to which they are bound. This diversity of behavior leads to a remarkable versatility for the catalysis of reactions by Lewis bases.
- 7Murillo, F.; Barroso, J.; de los Santos, M. G.; Ávila, G.; Pan, S.; Fernández-Herrera, M. A.; Merino, G. Revisiting the Formation Mechanism of 1,3,4-Oxadiazole-2(3H)-ones from Hydrazonyl Chloride and Carbon Dioxide. J. Org. Chem. 2018, 83, 13045– 13050, DOI: 10.1021/acs.joc.8b01676[ACS Full Text
], [CAS], Google Scholar
7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhvVKgtb7L&md5=a2af3f4295d724218efd82995fcb7c29Revisiting the Formation Mechanism of 1,3,4-Oxadiazole-2(3H)-ones from Hydrazonyl Chloride and Carbon DioxideMurillo, Fernando; Barroso, Jorge; de los Santos, Maria G.; Avila, Gustavo; Pan, Sudip; Fernandez-Herrera, Maria A.; Merino, GabrielJournal of Organic Chemistry (2018), 83 (21), 13045-13050CODEN: JOCEAH; ISSN:0022-3263. (American Chemical Society)The reaction mechanism for the synthesis of 1,3,4-oxadiazole-2(3H)-ones from hydrazonyl chloride and CO2 in the presence of CsF/18-crown-6 and toluene, is revisited via d. functional theory computations. Although this reaction was earlier classified as a 1,3-dipolar cycloaddn., we found some competing pathways involved therein. The mechanisms including the (F-CO2)- anion and the nitrile imine intermediate are some options. The dimerization of nitrile imine is another competing mechanism in this reaction. Our results show that the most favorable mechanism proceeds via a stepwise pathway without involving any nitrile imine intermediate or the (F-CO2)- anion. The F- anion, resulting from the formation of a complex between 18-crown-6 and Cs+ cation, deprotonates the nitrile imine precursor easily, which acts then as a nucleophilic anion, enhancing the reactivity of CO2 toward it. The mechanism for the reaction with COS, an isoelectronic analog of CO2, is also explored. - 8(a) Wang, X. S.; Lee, Y. J.; Liu, W. R. The nitrilimine-alkene cycloaddition is an ultra rapid click reaction. Chem. Commun. 2014, 50, 3176– 3179, DOI: 10.1039/C3CC48682F[Crossref], [PubMed], [CAS], Google Scholar.8ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXjtFSnsb8%253D&md5=110e7d55cf4cdb059796fad9c551e83dThe nitrilimine-alkene cycloaddition is an ultra rapid click reactionWang, Xiaoshan Shayna; Lee, Yan-Jiun; Liu, Wenshe R.Chemical Communications (Cambridge, United Kingdom) (2014), 50 (24), 3176-3179CODEN: CHCOFS; ISSN:1359-7345. (Royal Society of Chemistry)The transient formation of nitrilimine in aq. conditions is greatly influenced by pH and chloride. In basic conditions (pH 10) with no chloride, a diarylnitrilimine precursor readily ionizes to form diarylnitrilimine that reacts almost instantly with an acrylamide-contg. protein and fluorescently labels it.(b) Molteni, G.; Ponti, A. The Nitrilimine-Alkene Cycloaddition Regioselectivity Rationalized by Density Functional Theory Reactivity Indices. Molecules 2017, 22, 202, DOI: 10.3390/molecules22020202[Crossref], [CAS], Google Scholar.8bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXovVyrsbY%253D&md5=6911064fc37031e100a76f6c48fa66e6The nitrilimine-alkene cycloaddition regioselectivity rationalized by density functional theory reactivity indicesMolteni, Giorgio; Ponti, AlessandroMolecules (2017), 22 (2), 202/1-202/12CODEN: MOLEFW; ISSN:1420-3049. (MDPI AG)Conventional frontier MO theory is not able to satisfactorily explain the regioselectivity outcome of the nitrilimine-alkene cycloaddn. We considered that conceptual d. functional theory (DFT) could be an effective theor. framework to rationalize the regioselectivity of the title reaction. Several nitrilimine-alkene cycloaddns. were analyzed, for which we could find regioselectivity data in the literature. We computed DFT reactivity indexes at the B3LYP/6-311G(2d,p)//B3LYP/6-31G(d,p) and employed the grand potential stabilization criterion to calc. the preferred regioisomer. Exptl. and calcd. regioselectivity agree in the vast majority of cases. It was concluded that predominance of a single regioisomer can be obtained by maximizing (i) the chem. p.d. between nitrilimine and alkene and (ii) the local softness difference between the reactive at. sites within each reactant. Such maximization can be achieved by carefully selecting the substituents on both reactants.(c) Shawali, A. S. Chemoselectivity in 1,3-dipolar cycloaddition reactions of nitrilimines with multifunctionalized dipolarophiles. Curr. Org. Chem. 2014, 18, 598– 614, DOI: 10.2174/1385272819666140201002900[Crossref], [CAS], Google Scholar8chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXmsVWmtrg%253D&md5=ee2c4aa64224c841e816496fa766445eChemoselectivity in 1,3-Dipolar Cycloaddition Reactions of Nitrilimines with Multifunctionalized DipolarophilesShawali, Ahmad S.Current Organic Chemistry (2014), 18 (5), 598-614CODEN: CORCFE; ISSN:1385-2728. (Bentham Science Publishers Ltd.)A review. This review presents a survey of literature reports dealing with both site and peri-selectivities in cycloaddn. reactions of nitrilimines with multifunctionalized dipolarophiles. The literature results covered in this review, during the period from 1970 to mid 2013, demonstrated that chemo-selectivity plays an important role in synthetic design.
- 9(a) Bickelhaupt, F. M.; Houk, K. N. Analyzing Reaction Rates with the Distortion/Interaction-Activation Strain Model. Angew. Chem., Int. Ed. 2017, 56, 10070– 10086, DOI: 10.1002/anie.201701486[Crossref], [CAS], Google Scholar.9ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtFOlt7nL&md5=6e86c7cfaf7b09e58c0a50a589601578Analyzing Reaction Rates with the Distortion/Interaction-Activation Strain ModelBickelhaupt, F. Matthias; Houk, Kendall N.Angewandte Chemie, International Edition (2017), 56 (34), 10070-10086CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)The activation strain or distortion/interaction model is a tool to analyze activation barriers that det. reaction rates. For bimol. 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 mols. The energy required to distort the mols. 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 org. and inorg. chem., including substitutions and eliminations, cycloaddns., and several types of organometallic reactions.Bickelhaupt, F. M.; Houk, K. N. Angew. Chem. 2017, 129, 10204– 10221, DOI: 10.1002/ange.201701486 .(b) Vermeeren, P.; van der Lubbe, S. C. C.; Fonseca Guerra, C.; Bickelhaupt, F. M.; Hamlin, T. A. Understanding Chemical Reactivity Using the Activation Strain Model. Nat. Protoc. 2020, 15, 649– 667, DOI: 10.1038/s41596-019-0265-0[Crossref], [PubMed], [CAS], Google Scholar.9bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXntFOisg%253D%253D&md5=419807df5cfe254a7efa44da869ac4d7Understanding chemical reactivity using the activation strain modelVermeeren, Pascal; van der Lubbe, Stephanie C. C.; Fonseca Guerra, Celia; Bickelhaupt, F. Matthias; Hamlin, Trevor A.Nature Protocols (2020), 15 (2), 649-667CODEN: NPARDW; ISSN:1750-2799. (Nature Research)Understanding chem. 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 chem. 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 mol. 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 methodol. analyzes the change in energy along a reaction coordinate. The use of this methodol. has been proven to be crit. to the understanding of reactions, spanning the realms of the inorg. and org., as well as the supramol. and biochem., fields. This protocol provides step-by-step instructions-starting from the optimization of the stationary points and extending through calcn. of the potential energy surface and anal. of the trend-decisive energy terms-that can serve as a guide for carrying out the anal. of any given reaction of interest within hours to days, depending on the size of the mol. system.(c) Fernández, I.; Bickelhaupt, F. M. The activation strain model and molecular orbital theory: understanding and designing chemical reactions. Chem. Soc. Rev. 2014, 43, 4953– 4967, DOI: 10.1039/C4CS00055B[Crossref], [PubMed], [CAS], Google Scholar9chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhtValsrfN&md5=13f0b1e56b57a8605b030db8c0d5340aThe activation strain model and molecular orbital theory: understanding and designing chemical reactionsFernandez, Israel; Bickelhaupt, F. MatthiasChemical Society Reviews (2014), 43 (14), 4953-4967CODEN: CSRVBR; ISSN:0306-0012. (Royal Society of Chemistry)A review. In this Tutorial Review, we make the point that a true understanding of trends in reactivity (as opposed to measuring or simply computing them) requires a causal reactivity model. To this end, we present and discuss the Activation Strain Model (ASM). The ASM establishes the desired causal relationship between reaction barriers, on one hand, and the properties of reactants and characteristics of reaction mechanisms, on the other hand. In the ASM, the potential energy surface ΔE(ζ) along the reaction coordinate ζ is decompd. into the strain ΔEstrain(ζ) of the reactants that become increasingly deformed as the reaction proceeds, plus the interaction ΔEint(ζ) between these deformed reactants, i.e., ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). The ASM can be used in conjunction with any quantum chem. program. An anal. of the method and its application to problems in org. and organometallic chem. illustrate the power of the ASM as a unifying concept and a tool for rational design of reactants and catalysts.
- 10Bickelhaupt, F. M.; Baerends, E. J. Kohn-Sham Density Functional Theory: Predicting and Understanding Chemistry. In Reviews in Computational Chemistry; Lipkowitz, K. B., Boyd, D. B., Eds.; Wiley-VCH: New York, 2000; Vol. 15, pp 1– 86.
- 11(a) van Meer, R.; Gritsenko, O. V.; Baerends, E. J. Physical Meaning of Virtual Kohn-Sham Orbitals and Orbital Energies: An Ideal Basis for the Description of Molecular Excitations. J. Chem. Theory Comput. 2014, 10, 4432– 4441, DOI: 10.1021/ct500727c[ACS Full Text
], [CAS], Google Scholar
11ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhsFWjsrjL&md5=7142e461c840d366cba04ca734f4278dPhysical Meaning of Virtual Kohn-Sham Orbitals and Orbital Energies: An Ideal Basis for the Description of Molecular Excitationsvan Meer, R.; Gritsenko, O. V.; Baerends, E. J.Journal of Chemical Theory and Computation (2014), 10 (10), 4432-4441CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)In recent years, several benchmark studies on the performance of large sets of functionals in time-dependent d. functional theory (TDDFT) calcns. of excitation energies have been performed. The tested functionals do not approx. exact Kohn-Sham orbitals and orbital energies closely. We highlight the advantages of (close to) exact Kohn-Sham orbitals and orbital energies for a simple description, very often as just a single orbital-to-orbital transition, of mol. excitations. Benchmark calcns. are performed for the statistical av. of orbital potentials (SAOP) functional for the potential [J. Chem. Phys. 2000, 112, 1344; 2001, 114, 652], which approximates the true Kohn-Sham potential much better than LDA, GGA, mGGA, and hybrid potentials do. An accurate Kohn-Sham potential not only performs satisfactorily for calcd. vertical excitation energies of both valence and Rydberg transitions, but also exhibits appealing properties of the KS orbitals including occupied orbital energies close to ionization energies, virtual-occupied orbital energy gaps very close to excitation energies, realistic shapes of virtual orbitals, leading to straightforward interpretation of most excitations as single orbital transitions. We stress that such advantages are completely lost in time-dependent Hartree-Fock and partly in hybrid approaches. Many excitations and excitation energies calcd. with local d., generalized gradient, and hybrid functionals are spurious. There is, with an accurate KS, or even the LDA or GGA potentials, nothing problematic about the "band gap" in mols.: the HOMO-LUMO gap is close to the first excitation energy (the optical gap). - 12(a) Yu, S.; Vermeeren, P.; van Dommelen, K.; Bickelhaupt, F. M.; Hamlin, T. A. Understanding the 1,3-Dipolar Cycloadditions of Allenes. Chem. - Eur. J. 2020, 26, 11529– 11539, DOI: 10.1002/chem.202000857[Crossref], [PubMed], [CAS], Google Scholar.12ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhsFGhs73K&md5=d5177b3c6ce4715505e99e8b9d846dffUnderstanding the 1,3-Dipolar Cycloadditions of AllenesYu, Song; Vermeeren, Pascal; van Dommelen, Kevin; Bickelhaupt, F. Matthias; Hamlin, Trevor A.Chemistry - A European Journal (2020), 26 (50), 11529-11539CODEN: CEUJED; ISSN:0947-6539. (Wiley-VCH Verlag GmbH & Co. KGaA)We have quantum chem. studied the reactivity, site-, and regioselectivity of the 1,3-dipolar cycloaddn. between Me azide and various allenes, including the archetypal allene propadiene, heteroallenes, and cyclic allenes, by using d. functional theory (DFT). The 1,3-dipolar cycloaddn. reactivity of linear (hetero)allenes decreases as the no. of heteroatoms in the allene increases, and formation of the 1,5-adduct is, in all cases, favored over the 1,4-adduct. Both effects find their origin in the strength of the primary orbital interactions. The cycloaddn. reactivity of cyclic allenes was also investigated, and the increased predistortion of allenes, that results upon cyclization, leads to systematically lower activation barriers not due to the expected variations in the strain energy, but instead from the differences in the interaction energy. The geometric predistortion of cyclic allenes enhances the reactivity compared to linear allenes through a unique mechanism that involves a smaller HOMO-LUMO gap, which manifests as more stabilizing orbital interactions.(b) Hamlin, T. A.; Levandowski, B. J.; Narsaria, A. K.; Houk, K. N.; Bickelhaupt, F. M. Bickelhaupt. Structural Distortion of Cycloalkynes Influences Cycloaddition Rates by both Strain and Interaction Energies. Chem. - Eur. J. 2019, 25, 6342– 6348, DOI: 10.1002/chem.201900295[Crossref], [PubMed], [CAS], Google Scholar.12bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXmtFGks78%253D&md5=45443ff96577c491f65e9ac1e983cfd3Structural Distortion of Cycloalkynes Influences Cycloaddition Rates both by Strain and Interaction EnergiesHamlin, Trevor A.; Levandowski, Brian J.; Narsaria, Ayush K.; Houk, Kendall N.; Bickelhaupt, F. MatthiasChemistry - A European Journal (2019), 25 (25), 6342-6348CODEN: CEUJED; ISSN:0947-6539. (Wiley-VCH Verlag GmbH & Co. KGaA)The reactivities of 2-butyne, cycloheptyne, cyclooctyne, and cyclononyne in the 1,3-dipolar cycloaddn. reaction with Me azide were evaluated through DFT calcns. at the M06-2X/6-311++G(d)//M06-2X/6-31+G(d) level of theory. Computed activation free energies for the cycloaddns. 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 anal. 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.(c) Hamlin, T. A.; Svatunek, D.; Yu, S.; Ridder, L.; Infante, I.; Visscher, L.; Bickelhaupt, F. M. Elucidating the Trends in Reactivity of Aza-1,3-Dipolar Cycloadditions. Eur. J. Org. Chem. 2019, 2019, 378– 386, DOI: 10.1002/ejoc.201800572[Crossref], [CAS], Google Scholar12chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsFKqtrjI&md5=f3243b5ec83ba4fe8cfa4711292358eaElucidating the Trends in Reactivity of Aza-1,3-Dipolar CycloadditionsHamlin, Trevor A.; Svatunek, Dennis; Yu, Song; Ridder, Lars; Infante, Ivan; Visscher, Lucas; Bickelhaupt, F. MatthiasEuropean Journal of Organic Chemistry (2019), 2019 (2-3), 378-386CODEN: EJOCFK; ISSN:1099-0690. (Wiley-VCH Verlag GmbH & Co. KGaA)This report describes a d. functional theory investigation into the reactivities of a series of aza-1,3-dipoles with ethylene at the BP86/TZ2P level. A benchmark study was carried out using QMflows, a newly developed program for automated workflows of quantum chem. calcns. In total, 24 1,3-dipolar cycloaddn. (1,3-DCA) reactions were benchmarked using the highly accurate G3B3 method as a ref. We screened a no. of exchange and correlation functionals, including PBE, OLYP, BP86, BLYP, both with and without explicit dispersion corrections, to assess their accuracies and to det. which of these computationally efficient functionals performed the best for calcg. the energetics for cycloaddn. reactions. The BP86/TZ2P method produced the smallest errors for the activation and reaction enthalpies. Then, to understand the factors controlling the reactivity in these reactions, seven archetypal aza-1,3-dipolar cycloaddns. were investigated using the activation strain model and energy decompn. anal. Our investigations highlight the fact that differences in activation barrier for these 1,3-DCA reactions do not arise from differences in strain energy of the dipole, as previously proposed. Instead, relative reactivities originate from differences in interaction energy. Anal. of the 1,3-dipole-dipolarophile interactions reveals the reactivity trends primarily result from differences in the extent of the primary orbital interactions.
- 13
Differences in the dispersion curves were neglectable and were omitted for clarity.
There is no corresponding record for this reference. - 14Svatunek, D.; Houk, K. N. autoDIAS: a python tool for an automated distortion/interaction activation strain analysis. J. Comput. Chem. 2019, 40, 2509, DOI: 10.1002/jcc.26023[Crossref], [PubMed], [CAS], Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtlejtrfI&md5=52541195d19bc0c9afa1c1da7a34cfc7autoDIAS: a python tool for an automated distortion/interaction activation strain analysisSvatunek, Dennis; Houk, Kendall N.Journal of Computational Chemistry (2019), 40 (28), 2509-2515CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)The distortion/interaction activation strain (DIAS) anal. is a powerful tool for the investigation of energy barriers. However, setup and data anal. of such a calcn. can be cumbersome and requires lengthy intervention of the user. We present autoDIAS, a python tool for the automated setup, performance, and data extn. of the DIAS anal., including automated detection of fragments and relevant geometric parameters. © 2019 Wiley Periodicals, Inc.
- 15Sun, X.; Soini, T. M.; Poater, J.; Hamlin, T. A.; Bickelhaupt, F. M. PyFrag 2019—Automating the exploration and analysis of reaction mechanisms. J. Comput. Chem. 2019, 40, 2227, DOI: 10.1002/jcc.25871[Crossref], [PubMed], [CAS], Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtFWlsrrN&md5=732dcd3dcee95362fd072fd71191236fPyFrag 2019-Automating the exploration and analysis of reaction mechanismsSun, Xiaobo; Soini, Thomas M.; Poater, Jordi; Hamlin, Trevor A.; Bickelhaupt, F. MatthiasJournal of Computational Chemistry (2019), 40 (25), 2227-2233CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)We present a substantial update to the PyFrag 2008 program, which was originally designed to perform a fragment-based activation strain anal. 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 assocd. with the automated computational exploration of reaction mechanisms: it (1) computes the reaction path by carrying out multiple parallel calcns. 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 decompn. 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.
- 16(a) Pracht, P.; Bohle, F.; Grimme, S. Automated exploration of the low-energy chemical space with fast quantum chemical methods. Phys. Chem. Chem. Phys. 2020, 22, 7169– 7192, DOI: 10.1039/C9CP06869D[Crossref], [PubMed], [CAS], Google Scholar.16ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXjsVSmtr0%253D&md5=2eb2432e4672d91f5fe1acb28ad2c1a0Automated exploration of the low-energy chemical space with fast quantum chemical methodsPracht, Philipp; Bohle, Fabian; Grimme, StefanPhysical Chemistry Chemical Physics (2020), 22 (14), 7169-7192CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We propose and discuss an efficient scheme for the in silico sampling for parts of the mol. chem. space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodn. ensembles at a quantum chem. level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated min. energy structures normally form the basis of further, mostly DFT computational work, such as the calcn. of spectra or macroscopic properties. By using basic quantum chem. methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chem. methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized mols.) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chem. and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.(b) Grimme, S.; Bannwarth, C.; Dohm, S.; Hansen, A.; Pisarek, J.; Pracht, P.; Neese, F. Fully Automated Quantum-Chemistry-Based Computation of Spin-Spin-Coupled Nuclear Magnetic Resonance Spectra. Angew. Chem., Int. Ed. 2017, 56, 14763– 14769, DOI: 10.1002/anie.201708266[Crossref], [CAS], Google Scholar16bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhs1emu7jI&md5=d73a66caa9934ea2e7ec995376b39a3eFully Automated Quantum-Chemistry-Based Computation of Spin-Spin-Coupled Nuclear Magnetic Resonance SpectraGrimme, Stefan; Bannwarth, Christoph; Dohm, Sebastian; Hansen, Andreas; Pisarek, Jana; Pracht, Philipp; Seibert, Jakob; Neese, FrankAngewandte Chemie, International Edition (2017), 56 (46), 14763-14769CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)The authors present a composite procedure for the quantum-chem. computation of spin-spin-coupled 1H NMR spectra for general, flexible mols. in soln. that is based on four main steps, namely conformer/rotamer ensemble (CRE) generation by the fast tight-binding method GFN-xTB and a newly developed search algorithm, computation of the relative free energies and NMR parameters, and solving the spin Hamiltonian. In this way the NMR-specific nuclear permutation problem is solved, and the correct spin symmetries were obtained. Energies, shielding consts., and spin-spin couplings are computed at state-of-the-art DFT levels with continuum solvation. A few (in)org. and transition-metal complexes are presented, and very good, unprecedented agreement between the theor. and exptl. spectra was achieved. The approach is routinely applicable to systems with up to 100-150 atoms and may open new avenues for the detailed (conformational) structure elucidation of, for example, natural products or drug mols.
- 17Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M. Gaussian 09, rev. D.01; Gaussian Inc.: Wallingford, CT, 2009.Google ScholarThere is no corresponding record for this reference.
- 18Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215, DOI: 10.1007/s00214-007-0310-x[Crossref], [CAS], Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXltFyltbY%253D&md5=c31d6f319d7c7a45aa9b716220e4a422The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionalsZhao, Yan; Truhlar, Donald G.Theoretical Chemistry Accounts (2008), 120 (1-3), 215-241CODEN: TCACFW; ISSN:1432-881X. (Springer GmbH)We present two new hybrid meta exchange-correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amt. of nonlocal exchange (2X), and it is parametrized only for nonmetals. The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree-Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree-Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochem., four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for mol. excitation energies. We also illustrate the performance of these 17 methods for three databases contg. 40 bond lengths and for databases contg. 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochem., kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chem. and for noncovalent interactions.
- 19Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297– 3305, DOI: 10.1039/b508541a[Crossref], [PubMed], [CAS], Google Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpsFWgu7o%253D&md5=a820fb6055c993b50c405ba0fc62b194Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracyWeigend, Florian; Ahlrichs, ReinhartPhysical Chemistry Chemical Physics (2005), 7 (18), 3297-3305CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 mols. representing (nearly) all elements-except lanthanides-in their common oxidn. states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, d. functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
- 20Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Generalized born solvation model SM12. J. Phys. Chem. B 2009, 113, 6378, DOI: 10.1021/jp810292n[ACS Full Text
], [CAS], Google Scholar
20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXksV2is74%253D&md5=54931a64c70d28445ee53876a8b1a4b9Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface TensionsMarenich, Aleksandr V.; Cramer, Christopher J.; Truhlar, Donald G.Journal of Physical Chemistry B (2009), 113 (18), 6378-6396CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)We present a new continuum solvation model based on the quantum mech. charge d. of a solute mol. interacting with a continuum description of the solvent. The model is called SMD, where the "D" stands for "d." to denote that the full solute electron d. is used without defining partial at. charges. "Continuum" denotes that the solvent is not represented explicitly but rather as a dielec. medium with surface tension at the solute-solvent boundary. SMD is a universal solvation model, where "universal" denotes its applicability to any charged or uncharged solute in any solvent or liq. medium for which a few key descriptors are known (in particular, dielec. const., refractive index, bulk surface tension, and acidity and basicity parameters). The model separates the observable solvation free energy into two main components. The first component is the bulk electrostatic contribution arising from a self-consistent reaction field treatment that involves the soln. of the nonhomogeneous Poisson equation for electrostatics in terms of the integral-equation-formalism polarizable continuum model (IEF-PCM). The cavities for the bulk electrostatic calcn. are defined by superpositions of nuclear-centered spheres. The second component is called the cavity-dispersion-solvent-structure term and is the contribution arising from short-range interactions between the solute and solvent mols. in the first solvation shell. This contribution is a sum of terms that are proportional (with geometry-dependent proportionality consts. called at. surface tensions) to the solvent-accessible surface areas of the individual atoms of the solute. The SMD model has been parametrized with a training set of 2821 solvation data including 112 aq. ionic solvation free energies, 220 solvation free energies for 166 ions in acetonitrile, methanol, and DMSO, 2346 solvation free energies for 318 neutral solutes in 91 solvents (90 nonaq. org. solvents and water), and 143 transfer free energies for 93 neutral solutes between water and 15 org. solvents. The elements present in the solutes are H, C, N, O, F, Si, P, S, Cl, and Br. The SMD model employs a single set of parameters (intrinsic at. Coulomb radii and at. surface tension coeffs.) optimized over six electronic structure methods: M05-2X/MIDI!6D, M05-2X/6-31G*, M05-2X/6-31+G**, M05-2X/cc-pVTZ, B3LYP/6-31G*, and HF/6-31G*. Although the SMD model has been parametrized using the IEF-PCM protocol for bulk electrostatics, it may also be employed with other algorithms for solving the nonhomogeneous Poisson equation for continuum solvation calcns. in which the solute is represented by its electron d. in real space. This includes, for example, the conductor-like screening algorithm. With the 6-31G* basis set, the SMD model achieves mean unsigned errors of 0.6-1.0 kcal/mol in the solvation free energies of tested neutrals and mean unsigned errors of 4 kcal/mol on av. for ions with either Gaussian03 or GAMESS. - 21Goerigk, L.; Hansen, A.; Bauer, C.; Ehrlich, S.; Najibi, A.; Grimme, S. A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions. Phys. Chem. Chem. Phys. 2017, 19, 32184, DOI: 10.1039/C7CP04913G[Crossref], [PubMed], [CAS], Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhslajtLnF&md5=f9393c9e3907336c4da053743797f8dfA look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactionsGoerigk, Lars; Hansen, Andreas; Bauer, Christoph; Ehrlich, Stephan; Najibi, Asim; Grimme, StefanPhysical Chemistry Chemical Physics (2017), 19 (48), 32184-32215CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We present the GMTKN55 benchmark database for general main group thermochem., kinetics and noncovalent interactions. Compared to its popular predecessor GMTKN30, it allows assessment across a larger variety of chem. problems - with 13 new benchmark sets being presented for the first time - and it also provides ref. values of significantly higher quality for most sets. GMTKN55 comprises 1505 relative energies based on 2462 single-point calcns. and it is accessible to the user community via a dedicated website. Herein, we demonstrate the importance of better ref. values, and we re-emphasize the need for London-dispersion corrections in d. functional theory (DFT) treatments of thermochem. problems, including Minnesota methods. We assessed 217 variations of dispersion-cor. and -uncorrected d. functional approxns., and carried out a detailed anal. of 83 of them to identify robust and reliable approaches. Double-hybrid functionals are the most reliable approaches for thermochem. and noncovalent interactions, and they should be used whenever tech. feasible. These are, in particular, DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), and B2GPPLYP-D3(BJ). The best hybrids are ωB97X-V, M052X-D3(0), and ωB97X-D3, but we also recommend PW6B95-D3(BJ) as the best conventional global hybrid. At the meta-generalized-gradient (meta-GGA) level, the SCAN-D3(BJ) method can be recommended. Other meta-GGAs are outperformed by the GGA functionals revPBE-D3(BJ), B97-D3(BJ), and OLYP-D3(BJ). We note that many popular methods, such as B3LYP, are not part of our recommendations. In fact, with our results we hope to inspire a change in the user community's perception of common DFT methods. We also encourage method developers to use GMTKN55 for cross-validation studies of new methodologies.
- 22Ribeiro, R. F.; Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Use of solution-phase vibrational frequencies in continuum models for the free energy of solvation. J. Phys. Chem. B 2011, 115, 14556, DOI: 10.1021/jp205508z[ACS Full Text
], [CAS], Google Scholar
22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhsFSjtr3O&md5=3a164fbab7255d92e1099064e7f72261Use of Solution-Phase Vibrational Frequencies in Continuum Models for the Free Energy of SolvationRibeiro, Raphael F.; Marenich, Aleksandr V.; Cramer, Christopher J.; Truhlar, Donald G.Journal of Physical Chemistry B (2011), 115 (49), 14556-14562CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)We find that vibrational contributions to a solute's free energy are in general insensitive to whether the solute vibrational frequencies are computed in the gas phase or in soln. In most cases, the difference is smaller than the intrinsic error in solvation free energies assocd. with the continuum approxn. to solvation modeling, although care must be taken to avoid spurious results assocd. with limitations in the quantum-mech. harmonic-oscillator approxn. for very low-frequency mol. vibrations. We compute solute vibrational partition functions in aq. and carbon tetrachloride soln. and compare them to gas-phase mol. partition functions computed with the same level of theory and the same quasiharmonic approxn. for the diverse and extensive set of mols. and ions included in the training set of the SMD continuum solvation model, and we find mean unsigned differences in vibrational contributions to the solute free energy of only about 0.2 kcal/mol. On the basis of these results and a review of the theory, we conclude, in contrast to previous work, that using partition functions computed for mols. optimized in soln. is a correct and useful approach for averaging over solute degrees of freedom when computing free energies of solutes in soln., and it is moreover recommended for cases where liq. and gas-phase solute structures differ appreciably or when stationary points present in liq. soln. do not exist in the gas phase, for which we provide some examples. When gas-phase and soln.-phase geometries and frequencies are similar, the use of gas-phase geometries and frequencies is a useful approxn. - 23te Velde, G.; Bickelhaupt, F. M.; Baerends, E. J.; Fonseca Guerra, C.; van Gisbergen, S. J. A.; Snijders, J. G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931, DOI: 10.1002/jcc.1056[Crossref], [CAS], Google Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXjtlGntrw%253D&md5=314e7e942de9b28e664afc5adb2f574fChemistry with ADFTe Velde, G.; Bickelhaupt, F. M.; Baerends, E. J.; Fonseca Guerra, C.; Van Gisbergen, S. J. A.; Snijders, J. G.; Ziegler, T.Journal of Computational Chemistry (2001), 22 (9), 931-967CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)A review with 241 refs. We present the theor. and tech. foundations of the Amsterdam D. Functional (ADF) program with a survey of the characteristics of the code (numerical integration, d. fitting for the Coulomb potential, and STO basis functions). Recent developments enhance the efficiency of ADF (e.g., parallelization, near order-N scaling, QM/MM) and its functionality (e.g., NMR chem. shifts, COSMO solvent effects, ZORA relativistic method, excitation energies, frequency-dependent (hyper)polarizabilities, at. VDD charges). In the Applications section we discuss the phys. model of the electronic structure and the chem. bond, i.e., the Kohn-Sham MO (MO) theory, and illustrate the power of the Kohn-Sham MO model in conjunction with the ADF-typical fragment approach to quant. understand and predict chem. phenomena. We review the "Activation-strain TS interaction" (ATS) model of chem. reactivity as a conceptual framework for understanding how activation barriers of various types of (competing) reaction mechanisms arise and how they may be controlled, for example, in org. chem. or homogeneous catalysis. Finally, we include a brief discussion of exemplary applications in the field of biochem. (structure and bonding of DNA) and of time-dependent d. functional theory (TDDFT) to indicate how this development further reinforces the ADF tools for the anal. of chem. phenomena.
- 24Legault, C. Y. CYLview, 1.0b; Université de Sherbrooke: Canada, Sherbrooke, QC, 2009, www.cylview.org.Google ScholarThere is no corresponding record for this reference.
- 25Fonseca Guerra, C.; Handgraaf, J. W.; Baerends, E. J.; Bickelhaupt, F. M. Voronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysis?. J. Comput. Chem. 2004, 25, 189– 210, DOI: 10.1002/jcc.10351[Crossref], [PubMed], [CAS], Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD3srmvVOktw%253D%253D&md5=eb0605ed105ca7e7f2064a26620fa1cdVoronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysisFonseca Guerra Celia; Handgraaf Jan-Willem; Baerends Evert Jan; Bickelhaupt F MatthiasJournal of computational chemistry (2004), 25 (2), 189-210 ISSN:0192-8651.We present the Voronoi Deformation Density (VDD) method for computing atomic charges. The VDD method does not explicitly use the basis functions but calculates the amount of electronic density that flows to or from a certain atom due to bond formation by spatial integration of the deformation density over the atomic Voronoi cell. We compare our method to the well-known Mulliken, Hirshfeld, Bader, and Weinhold [Natural Population Analysis (NPA)] charges for a variety of biological, organic, and inorganic molecules. The Mulliken charges are (again) shown to be useless due to heavy basis set dependency, and the Bader charges (and often also the NPA charges) are not realistic, yielding too extreme values that suggest much ionic character even in the case of covalent bonds. The Hirshfeld and VDD charges, which prove to be numerically very similar, are to be recommended because they yield chemically meaningful charges. We stress the need to use spatial integration over an atomic domain to get rid of basis set dependency, and the need to integrate the deformation density in order to obtain a realistic picture of the charge rearrangement upon bonding. An asset of the VDD charges is the transparency of the approach owing to the simple geometric partitioning of space. The deformation density based charges prove to conform to chemical experience.
Cited By
This article is cited by 15 publications.
- Ze-Jie Lv, Philipp D. Engel, Lukas Alig, Sandipan Maji, Max C. Holthausen, Sven Schneider. Stabilizing Doubly Deprotonated Diazomethane: Isolable Complexes with CN22– and CN2– Radical Ligands. Journal of the American Chemical Society 2022, 144 (48) , 21872-21877. https://doi.org/10.1021/jacs.2c10743
- Jennifer P. Hernández-Mancera, Natalia Rojas-Valencia, Francisco Núñez-Zarur. Rationalizing the Substituent Effects in Diels–Alder Reactions of Triazolinediones with Anthracene. The Journal of Physical Chemistry A 2022, 126 (38) , 6657-6667. https://doi.org/10.1021/acs.jpca.2c04970
- Jorge Escorihuela, Lawrence M. Wolf. Computational Study on the Co-Mediated Intramolecular Pauson–Khand Reaction of Fluorinated and Chiral N-Tethered 1,7-Enynes. Organometallics 2022, 41 (18) , 2525-2534. https://doi.org/10.1021/acs.organomet.2c00227
- Thomas Hansen, Jasper C. Roozee, F. Matthias Bickelhaupt, Trevor A. Hamlin. How Solvation Influences the SN2 versus E2 Competition. The Journal of Organic Chemistry 2022, 87 (3) , 1805-1813. https://doi.org/10.1021/acs.joc.1c02354
- Ángel Martín Pendás, Evelio Francisco, Dimas Suárez, Aurora Costales, Natalia Díaz, Julen Munárriz, Tomás Rocha-Rinza, José Manuel Guevara-Vela. Atoms in molecules in real space: a fertile field for chemical bonding. Physical Chemistry Chemical Physics 2023, 25 (15) , 10231-10262. https://doi.org/10.1039/D2CP05540F
- Xiaobo Sun, Thomas Hansen, Jordi Poater, Trevor A. Hamlin, Friedrich Matthias Bickelhaupt. Rational design of iron catalysts for C – X bond activation. Journal of Computational Chemistry 2023, 44 (4) , 495-505. https://doi.org/10.1002/jcc.26818
- Asma Hassan, Muhammad Ismail, Ali H. Reshak, Zeshan Zada, Abdul Ahad Khan, Fazal Ur Rehman M, Muhammad Arif, Khadija Siraj, Shafqat Zada, G. Murtaza, Muhammad M. Ramli. Effect of heteroatoms on structural, electronic and spectroscopic properties of polyfuran, polythiophene and polypyrrole: A hybrid DFT approach. Journal of Molecular Structure 2023, 1274 , 134484. https://doi.org/10.1016/j.molstruc.2022.134484
- Song Yu, Eveline H. Tiekink, Pascal Vermeeren, F. Matthias Bickelhaupt, Trevor A. Hamlin. How Bases Catalyze Diels‐Alder Reactions. Chemistry – A European Journal 2023, 29 (7) https://doi.org/10.1002/chem.202203121
- Javier E. Alfonso‐Ramos, Ruben Van Lommel, David Hernández‐Castillo, Frank De Proft, Roy González‐Alemán, Erik V. Van der Eycken, Gerardo M. Ojeda‐Carralero. Origins of the Reactivity in 1,3‐Dipolar Cycloadditions of Acyl Isocyanide Ylides. European Journal of Organic Chemistry 2023, 26 (3) https://doi.org/10.1002/ejoc.202201028
- Thomas Hansen, Alba Nin‐Hill, Jeroen D. C. Codée, Trevor A. Hamlin, Carme Rovira. Rational Tuning of the Reactivity of Three‐Membered Heterocycle Ring Openings via S N 2 Reactions. Chemistry – A European Journal 2022, 28 (60) https://doi.org/10.1002/chem.202201649
- Pablo García-Aznar, Jorge Escorihuela. Computational insights into the inverse electron-demand Diels–Alder reaction of norbornenes with 1,2,4,5-tetrazines: norbornene substituents’ effects on the reaction rate. Organic & Biomolecular Chemistry 2022, 20 (32) , 6400-6412. https://doi.org/10.1039/D2OB01121B
- Thomas Hansen, Xiaobo Sun, Marco Dalla Tiezza, Willem‐Jan van Zeist, Joost N. P. van Stralen, Daan P. Geerke, Lando P. Wolters, Jordi Poater, Trevor A. Hamlin, F. Matthias Bickelhaupt. C−X Bond Activation by Palladium: Steric Shielding versus Steric Attraction. Chemistry – A European Journal 2022, 28 (44) https://doi.org/10.1002/chem.202201093
- Pascal Vermeeren, Trevor A. Hamlin, F. Matthias Bickelhaupt. How Ionization Catalyzes Diels–Alder Reactions. Chemistry – A European Journal 2022, 28 (40) https://doi.org/10.1002/chem.202201620
- Thomas Hansen, Xiaobo Sun, Marco Dalla Tiezza, Willem‐Jan Zeist, Jordi Poater, Trevor A. Hamlin, F. M. Bickelhaupt. C( sp n )−X (n=1–3) Bond Activation by Palladium. Chemistry – A European Journal 2022, 28 (26) https://doi.org/10.1002/chem.202103953
- Pascal Vermeeren, Trevor A. Hamlin, F. Matthias Bickelhaupt. Chemical reactivity from an activation strain perspective. Chemical Communications 2021, 57 (48) , 5880-5896. https://doi.org/10.1039/D1CC02042K
Abstract
Scheme 1
Scheme 1. Proposed Mechanism of the Fluoride-Catalyzed 1,3-Dipolar Cycloaddition between CO2 and Nitrilimines by Lu and Co-workers, in Which the Lewis Base F– Together with CO2 Forms the Activated Dipolarophile (i.e., CO2F–); (3) R1 and R2 = Aryl or AlkylScheme 2
Scheme 2. Possible Reaction Pathways for the 1,3-Dipolar Cycloaddition between Carbon Dioxide and Nitrilimine 1 in the Presence of the Lewis Base F–Figure 1
Figure 1. (a) Reaction profiles (ΔGtoluene in kcal mol–1) for the three reaction pathways of the studied 1,3-dipolar cycloaddition between carbon dioxide and nitrilimine 1, including (i) uncatalyzed cycloaddition (black), (ii) formation of CO2F– followed by stepwise cycloaddition (blue), and (iii) formation of 2 followed by stepwise cycloaddition (red), computed at SMD(toluene)-M06-2X-D3/def2-TZVP; TS = transition state and INT = intermediate. (b) Key transition state structures with key bond lengths (in Å) for the three reaction pathways.
Figure 2
Figure 2. (a) Distortion/interaction–activation strain model analysis; and (b) energy decomposition analysis of the cycloaddition reaction of 1 (black) and 2 (red) with CO2 (transition states indicated with a dot). (13) (c) Frontier molecular orbital diagram of the most important FMO-LUMOCO2 orbital interaction with the calculated energy gaps, orbital overlaps, and the S2/Δε terms, at consistent geometries with a N···C bond length of 2.14 Å. Computed at SMD(toluene)-M06-2X-D3/def2-TZVP using autoDIAS (14) for (a) and M06-2X-D3/TZ2P//SMD(toluene)-M06-2X-D3/def2-TZVP using PyFrag (15) for (b,c).
Scheme 3
Scheme 3. Novel Mechanism Emerging from Our Study for the Lewis Base F– Catalyzed 1,3-Dipolar Cycloaddition of CO2 to Nitrilimines, Where F– Activates the Dipole, Instead of the DipolarophileReferences
ARTICLE SECTIONSThis article references 25 other publications.
- 1(a) Milani, B.; Licini, G.; Clot, E.; Albrecht, M. Small Molecule Activation. Dalton Trans. 2016, 45, 14419– 14420, DOI: 10.1039/C6DT90140A[Crossref], [PubMed], [CAS], Google Scholar.1ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFWisb3F&md5=806efb4c52eec36651c19c47687080faSmall molecule activationMilani, Barbara; Licini, Giulia; Clot, Eric; Albrecht, MartinDalton Transactions (2016), 45 (37), 14419-14420CODEN: DTARAF; ISSN:1477-9226. (Royal Society of Chemistry)There is no expanded citation for this reference.(b) Zhang, W.; Lai, W.; Cao, R. Energy-Related Small Molecule Activation Reactions: Oxygen Reduction and Hydrogen and Oxygen Evolution Reactions Catalyzed by Porphyrin- and Corrole-Based Systems. Chem. Rev. 2017, 117, 3717– 3797, DOI: 10.1021/acs.chemrev.6b00299[ACS Full Text.
], [CAS], Google Scholar
1bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xhsl2gtbnO&md5=04523ae018db1238aece4e24e7012ef3Energy-Related Small Molecule Activation Reactions: Oxygen Reduction and Hydrogen and Oxygen Evolution Reactions Catalyzed by Porphyrin- and Corrole-Based SystemsZhang, Wei; Lai, Wenzhen; Cao, RuiChemical Reviews (Washington, DC, United States) (2017), 117 (4), 3717-3797CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)Globally increasing energy demands and environmental concerns related to the use of fossil fuels have stimulated extensive research to identify new energy systems and economies that are sustainable, clean, low cost, and environmentally benign. Hydrogen generation from solar-driven water splitting is a promising strategy to store solar energy in chem. bonds. The subsequent combustion of hydrogen in fuel cells produces elec. energy, and the only exhaust is water. These two reactions compose an ideal process to provide clean and sustainable energy. In such a process, a hydrogen evolution reaction (HER), an oxygen evolution reaction (OER) during water splitting, and an oxygen redn. reaction (ORR) as a fuel cell cathodic reaction are key steps that affect the efficiency of the overall energy conversion. Catalysts play key roles in this process by improving the kinetics of these reactions. Porphyrin-based and corrole-based systems are versatile and can efficiently catalyze the ORR, OER, and HER. Because of the significance of energy-related small mol. activation, this review covers recent progress in hydrogen evolution, oxygen evolution, and oxygen redn. reactions catalyzed by porphyrins and corroles.(c) Gómez, J. E.; Kleij, A. W. Catalytic Nonreductive Valorization of Carbon Dioxide into Fine Chemicals. In Adv. Organomet. Chem.; Pérez, P. J., Ed.; Academic Press: 2019; Vol. 71, pp 175– 226.(d) Chen, L.; Chen, G.; Leung, C.-F.; Cometto, C.; Robert, M.; Lau, T.-C. Molecular Quaterpyridine-Based Metal Complexes for Small Molecule Activation: Water Splitting and CO2 Reduction. Chem. Soc. Rev. 2020, 49, 7271– 7283, DOI: 10.1039/D0CS00927J[Crossref], [PubMed], [CAS], Google Scholar1dhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhvVKltrjF&md5=e98173e996d0191dcc4a7396e2ea48c5Molecular quaterpyridine-based metal complexes for small molecule activation: water splitting and CO2 reductionChen, Lingjing; Chen, Gui; Leung, Chi-Fai; Cometto, Claudio; Robert, Marc; Lau, Tai-ChuChemical Society Reviews (2020), 49 (20), 7271-7283CODEN: CSRVBR; ISSN:0306-0012. (Royal Society of Chemistry)A review. Artificial photosynthesis is considered as one of the most promising strategies for solar-to-fuel conversion through sunlight-driven water splitting and CO2 redn. This tutorial describes recent developments in the use of metal quaterpyridine complexes as electrocatalyts and photocatalysts for artificial photosynthesis. - 2(a) Liu, Q.; Wu, L.; Jackstell, R.; Beller, M. Using Carbon Dioxide as a Building Block in Organic Synthesis. Nat. Commun. 2015, 6, 5933, DOI: 10.1038/ncomms6933[Crossref], [PubMed], [CAS], Google Scholar.2ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhtF2lurjM&md5=e3ba6b5612253c0c6ebec3d937f137f8Using carbon dioxide as a building block in organic synthesisLiu, Qiang; Wu, Lipeng; Jackstell, Ralf; Beller, MatthiasNature Communications (2015), 6 (), 5933pp.CODEN: NCAOBW; ISSN:2041-1723. (Nature Publishing Group)A review. The most recent advances made in the area of CO2 valorization-turning CO2 into a useful chem. feedstock-under mild conditions has been reviewed. A special focus is given on the reaction modes for the CO2 activation and its application as C1 building block in org. synthesis. The following subjects will be addressed in this review : (1) novel transformations using carbon dioxide (briefly summarized); (2) different reaction modes for CO2 activation (main focus of this review); and (3) potential new applications of CO2 valorization.(b) Vogt, C.; Monai, M.; Sterk, E. B.; Palle, J.; Melcherts, A. E. M.; Zijlstra, B.; Groeneveld, E.; Berben, P. H.; Boereboom, J. M.; Hensen, E. J. M.; Meirer, F.; Filot, I. A. W.; Weckhuysen, B. M. Understanding Carbon Dioxide Activation and Carbon-Carbon Coupling over Nickel. Nat. Commun. 2019, 10, 5330, DOI: 10.1038/s41467-019-12858-3[Crossref], [PubMed], [CAS], Google Scholar.2bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3Mfjt1ajtg%253D%253D&md5=4fee60b8521239e02af56b1978b41b28Understanding carbon dioxide activation and carbon-carbon coupling over nickelVogt Charlotte; Monai Matteo; Sterk Ellen B; Palle Jonas; Melcherts Angela E M; Boereboom Jelle M; Meirer Florian; Weckhuysen Bert M; Zijlstra Bart; Hensen Emiel J M; Filot Ivo A W; Groeneveld Esther; Berben Peter HNature communications (2019), 10 (1), 5330 ISSN:.Carbon dioxide is a desired feedstock for platform molecules, such as carbon monoxide or higher hydrocarbons, from which we will be able to make many different useful, value-added chemicals. Its catalytic hydrogenation over abundant metals requires the amalgamation of theoretical knowledge with materials design. Here we leverage a theoretical understanding of structure sensitivity, along with a library of different supports, to tune the selectivity of methanation in the Power-to-Gas concept over nickel. For example, we show that carbon dioxide hydrogenation over nickel can and does form propane, and that activity and selectivity can be tuned by supporting different nickel particle sizes on various oxides. This theoretical and experimental toolbox is not only useful for the highly selective production of methane, but also provides new insights for carbon dioxide activation and subsequent carbon-carbon coupling towards value-added products thereby reducing the deleterious effects of this environmentally harmful molecule.(c) Schilling, W.; Das, S. Transition Metal-Free Synthesis of Carbamates Using CO2 as the Carbon Source. ChemSusChem 2020, 13, 6246– 6258[PubMed], [CAS], Google Scholar2chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXitFKjtrrJ&md5=a6465c2b7cba8c074554ac90f2ce1172Transition Metal-Free Synthesis of Carbamates Using CO2 as the Carbon SourceSchilling, Waldemar; Das, ShoubhikChemSusChem (2020), 13 (23), 6246-6258CODEN: CHEMIZ; ISSN:1864-5631. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. Utilization of carbon dioxide as a C1 synthon is highly attractive for the synthesis of valuable chems. However, activation of CO2 is highly challenging, owing to its thermodn. stability and kinetic inertness. With this in mind, several strategies have been developed for the generation of carbon-heteroatom bonds. Among these, formation of C-N bonds is highly attractive, esp., when carbamates can be synthesized directly from CO2. This Minireview focuses on transition metal-free approaches for the fixation of CO2 to generate carbamates for the prodn. of fine chems. and pharmaceuticals. Within the past decade, transition metal-free approaches have gained increasing attention, but traditional reviews have rarely focused on these approaches. Direct comparisons between such methods have been even more scarce. This Minireview seeks to address this discrepancy.
- 3Guo, C.-X.; Zhang, W.-Z.; Zhang, N.; Lu, X.-B. 1,3-Dipolar Cycloaddition of Nitrile Imine with Carbon Dioxide: Access to 1,3,4-Oxadiazole-2(3H)-ones. J. Org. Chem. 2017, 82, 7637, DOI: 10.1021/acs.joc.7b00963[ACS Full Text
], [CAS], Google Scholar
3https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVGgt73L&md5=0d77612903d20b3a9087cefd6e7096831,3-Dipolar Cycloaddition of Nitrile Imine with Carbon Dioxide: Access to 1,3,4-Oxadiazole-2(3H)-onesGuo, Chun-Xiao; Zhang, Wen-Zhen; Zhang, Ning; Lu, Xiao-BingJournal of Organic Chemistry (2017), 82 (14), 7637-7642CODEN: JOCEAH; ISSN:0022-3263. (American Chemical Society)Efficient synthesis of 1,3,4-oxadiazole-2(3H)-one was achieved by CsF/18-crown-6 mediated 1,3-dipolar cycloaddn. of nitrile imine and 2.0 MPa of CO2. CsF/18-crown-6 played a key role in enhancing the reactivity of CO2 as a 1,3-dipolarophile. The practical utility of this transition-metal-free approach to 1,3,4-oxadiazole-2(3H)-one is highlighted by the convenient synthesis of a com. herbicide Oxadiazon and a MAO B inhibitor. - 4(a) Zhang, X.; Gross, U.; Seppelt, K. Fluorocarbonate, [FCO2]−: Preparation and Structure. Angew. Chem., Int. Ed. Engl. 1995, 34, 1858– 1860, DOI: 10.1002/anie.199518581[Crossref], [CAS], Google Scholar.4ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXotF2ku78%253D&md5=1cb058cd84cf7b7967dca42bc0b48a22Fluorocarbonate, [FCO2]-: preparation and structureZhang, Xiongzhi; Gross, Udo; Seppelt, KonradAngewandte Chemie, International Edition in English (1995), 34 (17), 1858-60CODEN: ACIEAY; ISSN:0570-0833. (VCH)CO2 reacts with F- salts of Me4N+ or Me3(Me3CCH2)N+ or 1,1,3,3,5,5-hexamethylpiperidinium (pip+) to give [FCO2]- salts. Elemental anal. of pip+[FCO2] and NMR of Me4N[FCO2] confirm the identity of the anion. For the latter, multinuclear NMR, IR and Raman data are given. Ab initio calcns. for enthalpy of reaction of F- with CO, N2O, or H2C:C:CH2 indicate the fluoride affinity of CO2 may be the approx. lower limit for a reaction.(b) Arnold, D. W.; Bradforth, S. E.; Kim, E. H.; Neumark, D. M. Study of halogen-carbon dioxide clusters and the fluoroformyloxyl radical by photodetachment of X–(CO2) (X= I, Cl, Br) and FCO2–. J. Chem. Phys. 1995, 102, 3493, DOI: 10.1063/1.468575[Crossref], [CAS], Google Scholar4bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXjvFGqtbY%253D&md5=2d1eeb7b405437e16aa326c4d58c8ed9Study of halogen-carbon dioxide clusters and the fluoroformyloxyl radical by photodetachment of X-(CO2) (X = I,Cl,Br) and FCO-2Arnold, Don W.; Bradforth, Stephen E.; Kim, Eun H.; Neumark, Daniel M.Journal of Chemical Physics (1995), 102 (9), 3493-509CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Photoelectron spectra have been measured for the anions X-(CO2), with X = I, Br, Cl, and F. The vibrationally resolved spectra show that I-(CO2), Br-(CO2), and Cl-(CO2) are primarily electrostatically bound clusters, although the charge-quadrupole interaction is strong enough to distort the CO2 mol. by as much as 10° [in Cl-(CO2)]. Ab initio calcns. and electrostatic models are used to describe the geometry and bonding of these clusters. The photoelectron spectrum of FCO2- is qual. different and shows transitions to both the ~X 2B2 ground and the ~A2A2 first excited electronic states of the covalently bound FCO2 radical. The previously unobserved ~A 2A2 state is measured to lie 0.579 eV above the ground state. Vibrational frequencies are assigned with the assistance of ab initio calcns. The FCO2 heat of formation is detd. to be ΔfH0298(FCO2)= -85.2±2.8 kcal/mol. While both FCO2- and FCO2 are more strongly bound than the other halide-CO2 clusters, the C-F bonds are very weak relative to C-F bounds found in other halocarbon compds.
- 5(a) Ukai, K.; Aoki, M.; Takaya, J.; Iwasawa, N. Rhodium(I)-Catalyzed Carboxylation of Aryl- and Alkenylboronic Esters with CO2. J. Am. Chem. Soc. 2006, 128, 8706, DOI: 10.1021/ja061232m[ACS Full Text.
], [CAS], Google Scholar
5ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XmtVSjsbo%253D&md5=0c880ab4e3f61a8cf0f5e6ddab46953cRhodium(I)-Catalyzed Carboxylation of Aryl- and Alkenylboronic Esters with CO2Ukai, Kazutoshi; Aoki, Masao; Takaya, Jun; Iwasawa, NobuharuJournal of the American Chemical Society (2006), 128 (27), 8706-8707CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)When the esters of arylboronic acids with 2,2-dimethyl-1,3-propanol were treated with a catalytic amt. of [Rh(OH)(cod)]2 in the presence of 1,3-bis(diphenylphosphino)propane and CsF in dioxane at 60 °C under carbon dioxide atm., the benzoic acid derivs. were obtained in good yields. Reactions of alkenylboronic esters also proceeded under similar conditions to give α,β-unsatd. carboxylic acids. As these boronic esters are now easily available through coupling or direct borylation reactions, this method would be a useful method for the prepn. of various functionalized aryl- and alkenyl carboxylic acids. For example, the rhodium-catalyzed reaction of a boronic acid ester, i.e., 5,5-dimethyl-2-phenyl-1,3,2-dioxaborinane, with carbon dioxide gave benzoic acid.(b) Li, S.; Ma, S. Quadri-Synergetic Effect for Highly Effective Carbon Dioxide Fixation and Its Application to Indoloquinolinone. Adv. Synth. Catal. 2012, 354, 2387– 2394, DOI: 10.1002/adsc.201200469[Crossref], [CAS], Google Scholar.5bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38Xht1Glu7fE&md5=b15476302f7b12a0de48cb35cefd67a2Quadri-Synergetic Effect for Highly Effective Carbon Dioxide Fixation and Its Application to IndoloquinolinoneLi, Suhua; Ma, ShengmingAdvanced Synthesis & Catalysis (2012), 354 (13), 2387-2394, S2387/1-S2387/65CODEN: ASCAF7; ISSN:1615-4150. (Wiley-VCH Verlag GmbH & Co. KGaA)The first copper-catalyzed cyclic anti-nucleometallation-carboxylation of 2-alkynylanilines with carbon dioxide in the presence of dimethylzinc (ZnMe2) and cesium fluoride (CsF) for the effective synthesis of indolyl-3-carboxylic acids and indolodihydropyran-2-one is described. Through a mechanistic study, it is unveiled that the metal ions Cu2+, Zn2+, Cs+ and F- are working together for this CO2-based highly efficient carboxylation.(c) Takaya, J.; Tadami, S.; Ukai, K.; Iwasawa, N. Copper (I)-catalyzed carboxylation of aryl-and alkenylboronic esters. Org. Lett. 2008, 10, 2697– 2700, DOI: 10.1021/ol800829q[ACS Full Text], [CAS], Google Scholar
5chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXmsVKhtbg%253D&md5=80cfe34f459e75c4618d72e8245dac2eCopper(I)-Catalyzed Carboxylation of Aryl- and Alkenylboronic EstersTakaya, Jun; Tadami, Satoshi; Ukai, Kazutoshi; Iwasawa, NobuharuOrganic Letters (2008), 10 (13), 2697-2700CODEN: ORLEF7; ISSN:1523-7060. (American Chemical Society)The copper(I)-catalyzed carboxylation reaction of aryl- and alkenylboronic esters proceeded smoothly under CO2 to give the corresponding carboxylic acids in good yield. This reaction showed wide generality with higher functional group tolerance compared to the corresponding Rh(I)-catalyzed reaction. - 6Denmark, S. E.; Beutner, G. L. Lewis Base Catalysis in Organic Synthesis. Angew. Chem., Int. Ed. 2008, 47, 1560– 1638, DOI: 10.1002/anie.200604943[Crossref], [CAS], Google Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXis1Krtr8%253D&md5=234cf65fa06fdaaabd490440bd2d7ce4Lewis base catalysis in organic synthesisDenmark, Scott E.; Beutner, Gregory L.Angewandte Chemie, International Edition (2008), 47 (9), 1560-1638CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. The legacy of Gilbert Newton Lewis (1875-1946) pervades the lexicon of chem. bonding and reactivity; the power of his concept of donor-acceptor bonding is evident in the eponymous foundations of electron-pair acceptors (Lewis acids) and donors (Lewis bases). His discovery ushered in the use of Lewis acids as reagents and catalysts for org. reactions. However, in recent years, the recognition that Lewis bases can also serve in this capacity has grown enormously. Most importantly, it has become increasingly apparent that the behavior of Lewis bases as agents for promoting chem. reactions is not merely as an electronic complement of the cognate Lewis acids: in fact Lewis bases are capable of enhancing both the electrophilic and nucleophilic character of mols. to which they are bound. This diversity of behavior leads to a remarkable versatility for the catalysis of reactions by Lewis bases.
- 7Murillo, F.; Barroso, J.; de los Santos, M. G.; Ávila, G.; Pan, S.; Fernández-Herrera, M. A.; Merino, G. Revisiting the Formation Mechanism of 1,3,4-Oxadiazole-2(3H)-ones from Hydrazonyl Chloride and Carbon Dioxide. J. Org. Chem. 2018, 83, 13045– 13050, DOI: 10.1021/acs.joc.8b01676[ACS Full Text
], [CAS], Google Scholar
7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhvVKgtb7L&md5=a2af3f4295d724218efd82995fcb7c29Revisiting the Formation Mechanism of 1,3,4-Oxadiazole-2(3H)-ones from Hydrazonyl Chloride and Carbon DioxideMurillo, Fernando; Barroso, Jorge; de los Santos, Maria G.; Avila, Gustavo; Pan, Sudip; Fernandez-Herrera, Maria A.; Merino, GabrielJournal of Organic Chemistry (2018), 83 (21), 13045-13050CODEN: JOCEAH; ISSN:0022-3263. (American Chemical Society)The reaction mechanism for the synthesis of 1,3,4-oxadiazole-2(3H)-ones from hydrazonyl chloride and CO2 in the presence of CsF/18-crown-6 and toluene, is revisited via d. functional theory computations. Although this reaction was earlier classified as a 1,3-dipolar cycloaddn., we found some competing pathways involved therein. The mechanisms including the (F-CO2)- anion and the nitrile imine intermediate are some options. The dimerization of nitrile imine is another competing mechanism in this reaction. Our results show that the most favorable mechanism proceeds via a stepwise pathway without involving any nitrile imine intermediate or the (F-CO2)- anion. The F- anion, resulting from the formation of a complex between 18-crown-6 and Cs+ cation, deprotonates the nitrile imine precursor easily, which acts then as a nucleophilic anion, enhancing the reactivity of CO2 toward it. The mechanism for the reaction with COS, an isoelectronic analog of CO2, is also explored. - 8(a) Wang, X. S.; Lee, Y. J.; Liu, W. R. The nitrilimine-alkene cycloaddition is an ultra rapid click reaction. Chem. Commun. 2014, 50, 3176– 3179, DOI: 10.1039/C3CC48682F[Crossref], [PubMed], [CAS], Google Scholar.8ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXjtFSnsb8%253D&md5=110e7d55cf4cdb059796fad9c551e83dThe nitrilimine-alkene cycloaddition is an ultra rapid click reactionWang, Xiaoshan Shayna; Lee, Yan-Jiun; Liu, Wenshe R.Chemical Communications (Cambridge, United Kingdom) (2014), 50 (24), 3176-3179CODEN: CHCOFS; ISSN:1359-7345. (Royal Society of Chemistry)The transient formation of nitrilimine in aq. conditions is greatly influenced by pH and chloride. In basic conditions (pH 10) with no chloride, a diarylnitrilimine precursor readily ionizes to form diarylnitrilimine that reacts almost instantly with an acrylamide-contg. protein and fluorescently labels it.(b) Molteni, G.; Ponti, A. The Nitrilimine-Alkene Cycloaddition Regioselectivity Rationalized by Density Functional Theory Reactivity Indices. Molecules 2017, 22, 202, DOI: 10.3390/molecules22020202[Crossref], [CAS], Google Scholar.8bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXovVyrsbY%253D&md5=6911064fc37031e100a76f6c48fa66e6The nitrilimine-alkene cycloaddition regioselectivity rationalized by density functional theory reactivity indicesMolteni, Giorgio; Ponti, AlessandroMolecules (2017), 22 (2), 202/1-202/12CODEN: MOLEFW; ISSN:1420-3049. (MDPI AG)Conventional frontier MO theory is not able to satisfactorily explain the regioselectivity outcome of the nitrilimine-alkene cycloaddn. We considered that conceptual d. functional theory (DFT) could be an effective theor. framework to rationalize the regioselectivity of the title reaction. Several nitrilimine-alkene cycloaddns. were analyzed, for which we could find regioselectivity data in the literature. We computed DFT reactivity indexes at the B3LYP/6-311G(2d,p)//B3LYP/6-31G(d,p) and employed the grand potential stabilization criterion to calc. the preferred regioisomer. Exptl. and calcd. regioselectivity agree in the vast majority of cases. It was concluded that predominance of a single regioisomer can be obtained by maximizing (i) the chem. p.d. between nitrilimine and alkene and (ii) the local softness difference between the reactive at. sites within each reactant. Such maximization can be achieved by carefully selecting the substituents on both reactants.(c) Shawali, A. S. Chemoselectivity in 1,3-dipolar cycloaddition reactions of nitrilimines with multifunctionalized dipolarophiles. Curr. Org. Chem. 2014, 18, 598– 614, DOI: 10.2174/1385272819666140201002900[Crossref], [CAS], Google Scholar8chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXmsVWmtrg%253D&md5=ee2c4aa64224c841e816496fa766445eChemoselectivity in 1,3-Dipolar Cycloaddition Reactions of Nitrilimines with Multifunctionalized DipolarophilesShawali, Ahmad S.Current Organic Chemistry (2014), 18 (5), 598-614CODEN: CORCFE; ISSN:1385-2728. (Bentham Science Publishers Ltd.)A review. This review presents a survey of literature reports dealing with both site and peri-selectivities in cycloaddn. reactions of nitrilimines with multifunctionalized dipolarophiles. The literature results covered in this review, during the period from 1970 to mid 2013, demonstrated that chemo-selectivity plays an important role in synthetic design.
- 9(a) Bickelhaupt, F. M.; Houk, K. N. Analyzing Reaction Rates with the Distortion/Interaction-Activation Strain Model. Angew. Chem., Int. Ed. 2017, 56, 10070– 10086, DOI: 10.1002/anie.201701486[Crossref], [CAS], Google Scholar.9ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtFOlt7nL&md5=6e86c7cfaf7b09e58c0a50a589601578Analyzing Reaction Rates with the Distortion/Interaction-Activation Strain ModelBickelhaupt, F. Matthias; Houk, Kendall N.Angewandte Chemie, International Edition (2017), 56 (34), 10070-10086CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)The activation strain or distortion/interaction model is a tool to analyze activation barriers that det. reaction rates. For bimol. 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 mols. The energy required to distort the mols. 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 org. and inorg. chem., including substitutions and eliminations, cycloaddns., and several types of organometallic reactions.Bickelhaupt, F. M.; Houk, K. N. Angew. Chem. 2017, 129, 10204– 10221, DOI: 10.1002/ange.201701486 .(b) Vermeeren, P.; van der Lubbe, S. C. C.; Fonseca Guerra, C.; Bickelhaupt, F. M.; Hamlin, T. A. Understanding Chemical Reactivity Using the Activation Strain Model. Nat. Protoc. 2020, 15, 649– 667, DOI: 10.1038/s41596-019-0265-0[Crossref], [PubMed], [CAS], Google Scholar.9bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXntFOisg%253D%253D&md5=419807df5cfe254a7efa44da869ac4d7Understanding chemical reactivity using the activation strain modelVermeeren, Pascal; van der Lubbe, Stephanie C. C.; Fonseca Guerra, Celia; Bickelhaupt, F. Matthias; Hamlin, Trevor A.Nature Protocols (2020), 15 (2), 649-667CODEN: NPARDW; ISSN:1750-2799. (Nature Research)Understanding chem. 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 chem. 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 mol. 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 methodol. analyzes the change in energy along a reaction coordinate. The use of this methodol. has been proven to be crit. to the understanding of reactions, spanning the realms of the inorg. and org., as well as the supramol. and biochem., fields. This protocol provides step-by-step instructions-starting from the optimization of the stationary points and extending through calcn. of the potential energy surface and anal. of the trend-decisive energy terms-that can serve as a guide for carrying out the anal. of any given reaction of interest within hours to days, depending on the size of the mol. system.(c) Fernández, I.; Bickelhaupt, F. M. The activation strain model and molecular orbital theory: understanding and designing chemical reactions. Chem. Soc. Rev. 2014, 43, 4953– 4967, DOI: 10.1039/C4CS00055B[Crossref], [PubMed], [CAS], Google Scholar9chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhtValsrfN&md5=13f0b1e56b57a8605b030db8c0d5340aThe activation strain model and molecular orbital theory: understanding and designing chemical reactionsFernandez, Israel; Bickelhaupt, F. MatthiasChemical Society Reviews (2014), 43 (14), 4953-4967CODEN: CSRVBR; ISSN:0306-0012. (Royal Society of Chemistry)A review. In this Tutorial Review, we make the point that a true understanding of trends in reactivity (as opposed to measuring or simply computing them) requires a causal reactivity model. To this end, we present and discuss the Activation Strain Model (ASM). The ASM establishes the desired causal relationship between reaction barriers, on one hand, and the properties of reactants and characteristics of reaction mechanisms, on the other hand. In the ASM, the potential energy surface ΔE(ζ) along the reaction coordinate ζ is decompd. into the strain ΔEstrain(ζ) of the reactants that become increasingly deformed as the reaction proceeds, plus the interaction ΔEint(ζ) between these deformed reactants, i.e., ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). The ASM can be used in conjunction with any quantum chem. program. An anal. of the method and its application to problems in org. and organometallic chem. illustrate the power of the ASM as a unifying concept and a tool for rational design of reactants and catalysts.
- 10Bickelhaupt, F. M.; Baerends, E. J. Kohn-Sham Density Functional Theory: Predicting and Understanding Chemistry. In Reviews in Computational Chemistry; Lipkowitz, K. B., Boyd, D. B., Eds.; Wiley-VCH: New York, 2000; Vol. 15, pp 1– 86.
- 11(a) van Meer, R.; Gritsenko, O. V.; Baerends, E. J. Physical Meaning of Virtual Kohn-Sham Orbitals and Orbital Energies: An Ideal Basis for the Description of Molecular Excitations. J. Chem. Theory Comput. 2014, 10, 4432– 4441, DOI: 10.1021/ct500727c[ACS Full Text
], [CAS], Google Scholar
11ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhsFWjsrjL&md5=7142e461c840d366cba04ca734f4278dPhysical Meaning of Virtual Kohn-Sham Orbitals and Orbital Energies: An Ideal Basis for the Description of Molecular Excitationsvan Meer, R.; Gritsenko, O. V.; Baerends, E. J.Journal of Chemical Theory and Computation (2014), 10 (10), 4432-4441CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)In recent years, several benchmark studies on the performance of large sets of functionals in time-dependent d. functional theory (TDDFT) calcns. of excitation energies have been performed. The tested functionals do not approx. exact Kohn-Sham orbitals and orbital energies closely. We highlight the advantages of (close to) exact Kohn-Sham orbitals and orbital energies for a simple description, very often as just a single orbital-to-orbital transition, of mol. excitations. Benchmark calcns. are performed for the statistical av. of orbital potentials (SAOP) functional for the potential [J. Chem. Phys. 2000, 112, 1344; 2001, 114, 652], which approximates the true Kohn-Sham potential much better than LDA, GGA, mGGA, and hybrid potentials do. An accurate Kohn-Sham potential not only performs satisfactorily for calcd. vertical excitation energies of both valence and Rydberg transitions, but also exhibits appealing properties of the KS orbitals including occupied orbital energies close to ionization energies, virtual-occupied orbital energy gaps very close to excitation energies, realistic shapes of virtual orbitals, leading to straightforward interpretation of most excitations as single orbital transitions. We stress that such advantages are completely lost in time-dependent Hartree-Fock and partly in hybrid approaches. Many excitations and excitation energies calcd. with local d., generalized gradient, and hybrid functionals are spurious. There is, with an accurate KS, or even the LDA or GGA potentials, nothing problematic about the "band gap" in mols.: the HOMO-LUMO gap is close to the first excitation energy (the optical gap). - 12(a) Yu, S.; Vermeeren, P.; van Dommelen, K.; Bickelhaupt, F. M.; Hamlin, T. A. Understanding the 1,3-Dipolar Cycloadditions of Allenes. Chem. - Eur. J. 2020, 26, 11529– 11539, DOI: 10.1002/chem.202000857[Crossref], [PubMed], [CAS], Google Scholar.12ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhsFGhs73K&md5=d5177b3c6ce4715505e99e8b9d846dffUnderstanding the 1,3-Dipolar Cycloadditions of AllenesYu, Song; Vermeeren, Pascal; van Dommelen, Kevin; Bickelhaupt, F. Matthias; Hamlin, Trevor A.Chemistry - A European Journal (2020), 26 (50), 11529-11539CODEN: CEUJED; ISSN:0947-6539. (Wiley-VCH Verlag GmbH & Co. KGaA)We have quantum chem. studied the reactivity, site-, and regioselectivity of the 1,3-dipolar cycloaddn. between Me azide and various allenes, including the archetypal allene propadiene, heteroallenes, and cyclic allenes, by using d. functional theory (DFT). The 1,3-dipolar cycloaddn. reactivity of linear (hetero)allenes decreases as the no. of heteroatoms in the allene increases, and formation of the 1,5-adduct is, in all cases, favored over the 1,4-adduct. Both effects find their origin in the strength of the primary orbital interactions. The cycloaddn. reactivity of cyclic allenes was also investigated, and the increased predistortion of allenes, that results upon cyclization, leads to systematically lower activation barriers not due to the expected variations in the strain energy, but instead from the differences in the interaction energy. The geometric predistortion of cyclic allenes enhances the reactivity compared to linear allenes through a unique mechanism that involves a smaller HOMO-LUMO gap, which manifests as more stabilizing orbital interactions.(b) Hamlin, T. A.; Levandowski, B. J.; Narsaria, A. K.; Houk, K. N.; Bickelhaupt, F. M. Bickelhaupt. Structural Distortion of Cycloalkynes Influences Cycloaddition Rates by both Strain and Interaction Energies. Chem. - Eur. J. 2019, 25, 6342– 6348, DOI: 10.1002/chem.201900295[Crossref], [PubMed], [CAS], Google Scholar.12bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXmtFGks78%253D&md5=45443ff96577c491f65e9ac1e983cfd3Structural Distortion of Cycloalkynes Influences Cycloaddition Rates both by Strain and Interaction EnergiesHamlin, Trevor A.; Levandowski, Brian J.; Narsaria, Ayush K.; Houk, Kendall N.; Bickelhaupt, F. MatthiasChemistry - A European Journal (2019), 25 (25), 6342-6348CODEN: CEUJED; ISSN:0947-6539. (Wiley-VCH Verlag GmbH & Co. KGaA)The reactivities of 2-butyne, cycloheptyne, cyclooctyne, and cyclononyne in the 1,3-dipolar cycloaddn. reaction with Me azide were evaluated through DFT calcns. at the M06-2X/6-311++G(d)//M06-2X/6-31+G(d) level of theory. Computed activation free energies for the cycloaddns. 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 anal. 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.(c) Hamlin, T. A.; Svatunek, D.; Yu, S.; Ridder, L.; Infante, I.; Visscher, L.; Bickelhaupt, F. M. Elucidating the Trends in Reactivity of Aza-1,3-Dipolar Cycloadditions. Eur. J. Org. Chem. 2019, 2019, 378– 386, DOI: 10.1002/ejoc.201800572[Crossref], [CAS], Google Scholar12chttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsFKqtrjI&md5=f3243b5ec83ba4fe8cfa4711292358eaElucidating the Trends in Reactivity of Aza-1,3-Dipolar CycloadditionsHamlin, Trevor A.; Svatunek, Dennis; Yu, Song; Ridder, Lars; Infante, Ivan; Visscher, Lucas; Bickelhaupt, F. MatthiasEuropean Journal of Organic Chemistry (2019), 2019 (2-3), 378-386CODEN: EJOCFK; ISSN:1099-0690. (Wiley-VCH Verlag GmbH & Co. KGaA)This report describes a d. functional theory investigation into the reactivities of a series of aza-1,3-dipoles with ethylene at the BP86/TZ2P level. A benchmark study was carried out using QMflows, a newly developed program for automated workflows of quantum chem. calcns. In total, 24 1,3-dipolar cycloaddn. (1,3-DCA) reactions were benchmarked using the highly accurate G3B3 method as a ref. We screened a no. of exchange and correlation functionals, including PBE, OLYP, BP86, BLYP, both with and without explicit dispersion corrections, to assess their accuracies and to det. which of these computationally efficient functionals performed the best for calcg. the energetics for cycloaddn. reactions. The BP86/TZ2P method produced the smallest errors for the activation and reaction enthalpies. Then, to understand the factors controlling the reactivity in these reactions, seven archetypal aza-1,3-dipolar cycloaddns. were investigated using the activation strain model and energy decompn. anal. Our investigations highlight the fact that differences in activation barrier for these 1,3-DCA reactions do not arise from differences in strain energy of the dipole, as previously proposed. Instead, relative reactivities originate from differences in interaction energy. Anal. of the 1,3-dipole-dipolarophile interactions reveals the reactivity trends primarily result from differences in the extent of the primary orbital interactions.
- 13
Differences in the dispersion curves were neglectable and were omitted for clarity.
There is no corresponding record for this reference. - 14Svatunek, D.; Houk, K. N. autoDIAS: a python tool for an automated distortion/interaction activation strain analysis. J. Comput. Chem. 2019, 40, 2509, DOI: 10.1002/jcc.26023[Crossref], [PubMed], [CAS], Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtlejtrfI&md5=52541195d19bc0c9afa1c1da7a34cfc7autoDIAS: a python tool for an automated distortion/interaction activation strain analysisSvatunek, Dennis; Houk, Kendall N.Journal of Computational Chemistry (2019), 40 (28), 2509-2515CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)The distortion/interaction activation strain (DIAS) anal. is a powerful tool for the investigation of energy barriers. However, setup and data anal. of such a calcn. can be cumbersome and requires lengthy intervention of the user. We present autoDIAS, a python tool for the automated setup, performance, and data extn. of the DIAS anal., including automated detection of fragments and relevant geometric parameters. © 2019 Wiley Periodicals, Inc.
- 15Sun, X.; Soini, T. M.; Poater, J.; Hamlin, T. A.; Bickelhaupt, F. M. PyFrag 2019—Automating the exploration and analysis of reaction mechanisms. J. Comput. Chem. 2019, 40, 2227, DOI: 10.1002/jcc.25871[Crossref], [PubMed], [CAS], Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtFWlsrrN&md5=732dcd3dcee95362fd072fd71191236fPyFrag 2019-Automating the exploration and analysis of reaction mechanismsSun, Xiaobo; Soini, Thomas M.; Poater, Jordi; Hamlin, Trevor A.; Bickelhaupt, F. MatthiasJournal of Computational Chemistry (2019), 40 (25), 2227-2233CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)We present a substantial update to the PyFrag 2008 program, which was originally designed to perform a fragment-based activation strain anal. 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 assocd. with the automated computational exploration of reaction mechanisms: it (1) computes the reaction path by carrying out multiple parallel calcns. 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 decompn. 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.
- 16(a) Pracht, P.; Bohle, F.; Grimme, S. Automated exploration of the low-energy chemical space with fast quantum chemical methods. Phys. Chem. Chem. Phys. 2020, 22, 7169– 7192, DOI: 10.1039/C9CP06869D[Crossref], [PubMed], [CAS], Google Scholar.16ahttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXjsVSmtr0%253D&md5=2eb2432e4672d91f5fe1acb28ad2c1a0Automated exploration of the low-energy chemical space with fast quantum chemical methodsPracht, Philipp; Bohle, Fabian; Grimme, StefanPhysical Chemistry Chemical Physics (2020), 22 (14), 7169-7192CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We propose and discuss an efficient scheme for the in silico sampling for parts of the mol. chem. space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodn. ensembles at a quantum chem. level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated min. energy structures normally form the basis of further, mostly DFT computational work, such as the calcn. of spectra or macroscopic properties. By using basic quantum chem. methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chem. methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized mols.) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chem. and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.(b) Grimme, S.; Bannwarth, C.; Dohm, S.; Hansen, A.; Pisarek, J.; Pracht, P.; Neese, F. Fully Automated Quantum-Chemistry-Based Computation of Spin-Spin-Coupled Nuclear Magnetic Resonance Spectra. Angew. Chem., Int. Ed. 2017, 56, 14763– 14769, DOI: 10.1002/anie.201708266[Crossref], [CAS], Google Scholar16bhttps://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhs1emu7jI&md5=d73a66caa9934ea2e7ec995376b39a3eFully Automated Quantum-Chemistry-Based Computation of Spin-Spin-Coupled Nuclear Magnetic Resonance SpectraGrimme, Stefan; Bannwarth, Christoph; Dohm, Sebastian; Hansen, Andreas; Pisarek, Jana; Pracht, Philipp; Seibert, Jakob; Neese, FrankAngewandte Chemie, International Edition (2017), 56 (46), 14763-14769CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)The authors present a composite procedure for the quantum-chem. computation of spin-spin-coupled 1H NMR spectra for general, flexible mols. in soln. that is based on four main steps, namely conformer/rotamer ensemble (CRE) generation by the fast tight-binding method GFN-xTB and a newly developed search algorithm, computation of the relative free energies and NMR parameters, and solving the spin Hamiltonian. In this way the NMR-specific nuclear permutation problem is solved, and the correct spin symmetries were obtained. Energies, shielding consts., and spin-spin couplings are computed at state-of-the-art DFT levels with continuum solvation. A few (in)org. and transition-metal complexes are presented, and very good, unprecedented agreement between the theor. and exptl. spectra was achieved. The approach is routinely applicable to systems with up to 100-150 atoms and may open new avenues for the detailed (conformational) structure elucidation of, for example, natural products or drug mols.
- 17Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M. Gaussian 09, rev. D.01; Gaussian Inc.: Wallingford, CT, 2009.Google ScholarThere is no corresponding record for this reference.
- 18Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215, DOI: 10.1007/s00214-007-0310-x[Crossref], [CAS], Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXltFyltbY%253D&md5=c31d6f319d7c7a45aa9b716220e4a422The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionalsZhao, Yan; Truhlar, Donald G.Theoretical Chemistry Accounts (2008), 120 (1-3), 215-241CODEN: TCACFW; ISSN:1432-881X. (Springer GmbH)We present two new hybrid meta exchange-correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amt. of nonlocal exchange (2X), and it is parametrized only for nonmetals. The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree-Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree-Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochem., four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for mol. excitation energies. We also illustrate the performance of these 17 methods for three databases contg. 40 bond lengths and for databases contg. 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochem., kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chem. and for noncovalent interactions.
- 19Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297– 3305, DOI: 10.1039/b508541a[Crossref], [PubMed], [CAS], Google Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpsFWgu7o%253D&md5=a820fb6055c993b50c405ba0fc62b194Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracyWeigend, Florian; Ahlrichs, ReinhartPhysical Chemistry Chemical Physics (2005), 7 (18), 3297-3305CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 mols. representing (nearly) all elements-except lanthanides-in their common oxidn. states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, d. functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
- 20Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Generalized born solvation model SM12. J. Phys. Chem. B 2009, 113, 6378, DOI: 10.1021/jp810292n[ACS Full Text
], [CAS], Google Scholar
20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXksV2is74%253D&md5=54931a64c70d28445ee53876a8b1a4b9Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface TensionsMarenich, Aleksandr V.; Cramer, Christopher J.; Truhlar, Donald G.Journal of Physical Chemistry B (2009), 113 (18), 6378-6396CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)We present a new continuum solvation model based on the quantum mech. charge d. of a solute mol. interacting with a continuum description of the solvent. The model is called SMD, where the "D" stands for "d." to denote that the full solute electron d. is used without defining partial at. charges. "Continuum" denotes that the solvent is not represented explicitly but rather as a dielec. medium with surface tension at the solute-solvent boundary. SMD is a universal solvation model, where "universal" denotes its applicability to any charged or uncharged solute in any solvent or liq. medium for which a few key descriptors are known (in particular, dielec. const., refractive index, bulk surface tension, and acidity and basicity parameters). The model separates the observable solvation free energy into two main components. The first component is the bulk electrostatic contribution arising from a self-consistent reaction field treatment that involves the soln. of the nonhomogeneous Poisson equation for electrostatics in terms of the integral-equation-formalism polarizable continuum model (IEF-PCM). The cavities for the bulk electrostatic calcn. are defined by superpositions of nuclear-centered spheres. The second component is called the cavity-dispersion-solvent-structure term and is the contribution arising from short-range interactions between the solute and solvent mols. in the first solvation shell. This contribution is a sum of terms that are proportional (with geometry-dependent proportionality consts. called at. surface tensions) to the solvent-accessible surface areas of the individual atoms of the solute. The SMD model has been parametrized with a training set of 2821 solvation data including 112 aq. ionic solvation free energies, 220 solvation free energies for 166 ions in acetonitrile, methanol, and DMSO, 2346 solvation free energies for 318 neutral solutes in 91 solvents (90 nonaq. org. solvents and water), and 143 transfer free energies for 93 neutral solutes between water and 15 org. solvents. The elements present in the solutes are H, C, N, O, F, Si, P, S, Cl, and Br. The SMD model employs a single set of parameters (intrinsic at. Coulomb radii and at. surface tension coeffs.) optimized over six electronic structure methods: M05-2X/MIDI!6D, M05-2X/6-31G*, M05-2X/6-31+G**, M05-2X/cc-pVTZ, B3LYP/6-31G*, and HF/6-31G*. Although the SMD model has been parametrized using the IEF-PCM protocol for bulk electrostatics, it may also be employed with other algorithms for solving the nonhomogeneous Poisson equation for continuum solvation calcns. in which the solute is represented by its electron d. in real space. This includes, for example, the conductor-like screening algorithm. With the 6-31G* basis set, the SMD model achieves mean unsigned errors of 0.6-1.0 kcal/mol in the solvation free energies of tested neutrals and mean unsigned errors of 4 kcal/mol on av. for ions with either Gaussian03 or GAMESS. - 21Goerigk, L.; Hansen, A.; Bauer, C.; Ehrlich, S.; Najibi, A.; Grimme, S. A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions. Phys. Chem. Chem. Phys. 2017, 19, 32184, DOI: 10.1039/C7CP04913G[Crossref], [PubMed], [CAS], Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhslajtLnF&md5=f9393c9e3907336c4da053743797f8dfA look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactionsGoerigk, Lars; Hansen, Andreas; Bauer, Christoph; Ehrlich, Stephan; Najibi, Asim; Grimme, StefanPhysical Chemistry Chemical Physics (2017), 19 (48), 32184-32215CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We present the GMTKN55 benchmark database for general main group thermochem., kinetics and noncovalent interactions. Compared to its popular predecessor GMTKN30, it allows assessment across a larger variety of chem. problems - with 13 new benchmark sets being presented for the first time - and it also provides ref. values of significantly higher quality for most sets. GMTKN55 comprises 1505 relative energies based on 2462 single-point calcns. and it is accessible to the user community via a dedicated website. Herein, we demonstrate the importance of better ref. values, and we re-emphasize the need for London-dispersion corrections in d. functional theory (DFT) treatments of thermochem. problems, including Minnesota methods. We assessed 217 variations of dispersion-cor. and -uncorrected d. functional approxns., and carried out a detailed anal. of 83 of them to identify robust and reliable approaches. Double-hybrid functionals are the most reliable approaches for thermochem. and noncovalent interactions, and they should be used whenever tech. feasible. These are, in particular, DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), and B2GPPLYP-D3(BJ). The best hybrids are ωB97X-V, M052X-D3(0), and ωB97X-D3, but we also recommend PW6B95-D3(BJ) as the best conventional global hybrid. At the meta-generalized-gradient (meta-GGA) level, the SCAN-D3(BJ) method can be recommended. Other meta-GGAs are outperformed by the GGA functionals revPBE-D3(BJ), B97-D3(BJ), and OLYP-D3(BJ). We note that many popular methods, such as B3LYP, are not part of our recommendations. In fact, with our results we hope to inspire a change in the user community's perception of common DFT methods. We also encourage method developers to use GMTKN55 for cross-validation studies of new methodologies.
- 22Ribeiro, R. F.; Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Use of solution-phase vibrational frequencies in continuum models for the free energy of solvation. J. Phys. Chem. B 2011, 115, 14556, DOI: 10.1021/jp205508z[ACS Full Text
], [CAS], Google Scholar
22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhsFSjtr3O&md5=3a164fbab7255d92e1099064e7f72261Use of Solution-Phase Vibrational Frequencies in Continuum Models for the Free Energy of SolvationRibeiro, Raphael F.; Marenich, Aleksandr V.; Cramer, Christopher J.; Truhlar, Donald G.Journal of Physical Chemistry B (2011), 115 (49), 14556-14562CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)We find that vibrational contributions to a solute's free energy are in general insensitive to whether the solute vibrational frequencies are computed in the gas phase or in soln. In most cases, the difference is smaller than the intrinsic error in solvation free energies assocd. with the continuum approxn. to solvation modeling, although care must be taken to avoid spurious results assocd. with limitations in the quantum-mech. harmonic-oscillator approxn. for very low-frequency mol. vibrations. We compute solute vibrational partition functions in aq. and carbon tetrachloride soln. and compare them to gas-phase mol. partition functions computed with the same level of theory and the same quasiharmonic approxn. for the diverse and extensive set of mols. and ions included in the training set of the SMD continuum solvation model, and we find mean unsigned differences in vibrational contributions to the solute free energy of only about 0.2 kcal/mol. On the basis of these results and a review of the theory, we conclude, in contrast to previous work, that using partition functions computed for mols. optimized in soln. is a correct and useful approach for averaging over solute degrees of freedom when computing free energies of solutes in soln., and it is moreover recommended for cases where liq. and gas-phase solute structures differ appreciably or when stationary points present in liq. soln. do not exist in the gas phase, for which we provide some examples. When gas-phase and soln.-phase geometries and frequencies are similar, the use of gas-phase geometries and frequencies is a useful approxn. - 23te Velde, G.; Bickelhaupt, F. M.; Baerends, E. J.; Fonseca Guerra, C.; van Gisbergen, S. J. A.; Snijders, J. G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931, DOI: 10.1002/jcc.1056[Crossref], [CAS], Google Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXjtlGntrw%253D&md5=314e7e942de9b28e664afc5adb2f574fChemistry with ADFTe Velde, G.; Bickelhaupt, F. M.; Baerends, E. J.; Fonseca Guerra, C.; Van Gisbergen, S. J. A.; Snijders, J. G.; Ziegler, T.Journal of Computational Chemistry (2001), 22 (9), 931-967CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)A review with 241 refs. We present the theor. and tech. foundations of the Amsterdam D. Functional (ADF) program with a survey of the characteristics of the code (numerical integration, d. fitting for the Coulomb potential, and STO basis functions). Recent developments enhance the efficiency of ADF (e.g., parallelization, near order-N scaling, QM/MM) and its functionality (e.g., NMR chem. shifts, COSMO solvent effects, ZORA relativistic method, excitation energies, frequency-dependent (hyper)polarizabilities, at. VDD charges). In the Applications section we discuss the phys. model of the electronic structure and the chem. bond, i.e., the Kohn-Sham MO (MO) theory, and illustrate the power of the Kohn-Sham MO model in conjunction with the ADF-typical fragment approach to quant. understand and predict chem. phenomena. We review the "Activation-strain TS interaction" (ATS) model of chem. reactivity as a conceptual framework for understanding how activation barriers of various types of (competing) reaction mechanisms arise and how they may be controlled, for example, in org. chem. or homogeneous catalysis. Finally, we include a brief discussion of exemplary applications in the field of biochem. (structure and bonding of DNA) and of time-dependent d. functional theory (TDDFT) to indicate how this development further reinforces the ADF tools for the anal. of chem. phenomena.
- 24Legault, C. Y. CYLview, 1.0b; Université de Sherbrooke: Canada, Sherbrooke, QC, 2009, www.cylview.org.Google ScholarThere is no corresponding record for this reference.
- 25Fonseca Guerra, C.; Handgraaf, J. W.; Baerends, E. J.; Bickelhaupt, F. M. Voronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysis?. J. Comput. Chem. 2004, 25, 189– 210, DOI: 10.1002/jcc.10351[Crossref], [PubMed], [CAS], Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD3srmvVOktw%253D%253D&md5=eb0605ed105ca7e7f2064a26620fa1cdVoronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysisFonseca Guerra Celia; Handgraaf Jan-Willem; Baerends Evert Jan; Bickelhaupt F MatthiasJournal of computational chemistry (2004), 25 (2), 189-210 ISSN:0192-8651.We present the Voronoi Deformation Density (VDD) method for computing atomic charges. The VDD method does not explicitly use the basis functions but calculates the amount of electronic density that flows to or from a certain atom due to bond formation by spatial integration of the deformation density over the atomic Voronoi cell. We compare our method to the well-known Mulliken, Hirshfeld, Bader, and Weinhold [Natural Population Analysis (NPA)] charges for a variety of biological, organic, and inorganic molecules. The Mulliken charges are (again) shown to be useless due to heavy basis set dependency, and the Bader charges (and often also the NPA charges) are not realistic, yielding too extreme values that suggest much ionic character even in the case of covalent bonds. The Hirshfeld and VDD charges, which prove to be numerically very similar, are to be recommended because they yield chemically meaningful charges. We stress the need to use spatial integration over an atomic domain to get rid of basis set dependency, and the need to integrate the deformation density in order to obtain a realistic picture of the charge rearrangement upon bonding. An asset of the VDD charges is the transparency of the approach owing to the simple geometric partitioning of space. The deformation density based charges prove to conform to chemical experience.
Supporting Information
Supporting Information
ARTICLE SECTIONSThe Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.joc.0c02963.
Additional computational results; and Cartesian coordinates, energies, and number of imaginary frequencies of all stationary points (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.