Benchmark Density Functional Theory Approach for the Calculation of Bond Dissociation Energies of the M–O2 Bond: A Key Step in Water Splitting ReactionsClick to copy article linkArticle link copied!
- Naveen Kosar*Naveen Kosar*Email: [email protected]. Phone: +92-343-7628474.Department of Chemistry, University of Management and Technology (UMT), C11, Johar Town, Lahore 54770, PakistanMore by Naveen Kosar
- Khurshid AyubKhurshid AyubDepartment of Chemistry, COMSATS University Islamabad, Abbottabad Campus, Abbottabad 22060, PakistanMore by Khurshid Ayub
- Mazhar Amjad GilaniMazhar Amjad GilaniDepartment of Chemistry, COMSATS University Islamabad, Lahore Campus, Lahore 54600, PakistanMore by Mazhar Amjad Gilani
- Shabbir MuhammadShabbir MuhammadDepartment of Chemistry, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi ArabiaMore by Shabbir Muhammad
- Tariq Mahmood*Tariq Mahmood*Email: [email protected]Department of Chemistry, COMSATS University Islamabad, Abbottabad Campus, Abbottabad 22060, PakistanMore by Tariq Mahmood
Abstract
A very fascinating aspect in quantum chemical research is to determine the accurate and cost-effective methods for the calculation of electronic and structural properties through a benchmark study. The current study focuses on the performance evaluation of density functional theory methods for the accurate measurement of bond dissociation energies (BDEs) of chemically important M–O2 bonds in water splitting reactions. The BDE measurement has got noteworthy attention due to its importance in all areas of chemistry. For BDE measurements of M–O2 bonds in five metal complexes with oxygen molecules, 14 density functionals (DFs) are chosen from seven classes of DFs with two series of mixed basis sets. A combination of pseudopotential and Pople basis sets [LANL2DZ & 6-31G(d) and SDD & 6-31+G(d)] are used as a series of mixed basis sets. The B3LYP-GD3BJ functional with LANL2DZ & 6-31G(d) gives outstanding results due to low deviations, error, and the best Pearson’s correlation (R) between the experimental and theoretical data. Our study suggested an efficient, low-cost, precise, and accurate B3LYP-GD3BJ/LANL2DZ & 6-31G(d) level of theory for BDE of the M–O2 bond, which may be useful for chemists working in the field of energy generation and utilization.
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.
1. Introduction
Figure 1
Figure 1. Catalytic water oxidation and oxygen molecule evolution mechanism by Ru complexes.
2. Results and Discussion
2.1. Evaluation of DFs with LANL2DZ & 6-31G(d) Basis Sets
Figure 2
Figure 2. RMSD of different DFs with two series of basis sets for M–O2 BDEs.
Figure 3
Figure 3. SD of different DFs with two series of basis sets for M–O2 BDEs.
Figure 4
Figure 4. MAE of different DFs with two series of basis sets for M–O2 BDEs.
classes of DFT | DFs | rmsd | SD | R | MAE |
---|---|---|---|---|---|
LDA | LSDA | 9.99 | 8.17 | 0.57 | –7.09 |
GGA | BP86 | 7.12 | 6.47 | 0.52 | –0.37 |
meta-GGA | TPSSTPSS | 7.21 | 7.41 | 0.52 | –1.26 |
H-GGA | B3LYP | 4.36 | 1.55 | 0.81 | –1.29 |
B3PW91 | 4.87 | 3.11 | 0.34 | –1.08 | |
PBE0 | 6.67 | 6.22 | 0.45 | –0.70 | |
MPWPW91 | 7.28 | 6.43 | 0.53 | –0.10 | |
B97 | 12.20 | 9.15 | 0.54 | 3.41 | |
GH meta-GGA | M05-2X | 17.54 | 16.33 | 0.54 | 5.51 |
M05 | 17.15 | 12.35 | 0.77 | 7.50 | |
M06-2X | 17.50 | 13.82 | 0.56 | 6.37 | |
M06 | 5.74 | 5.48 | 0.38 | –1.56 | |
RS H-GGA | CAM-B3LYP | 14.79 | 12.68 | 0.53 | 4.37 |
H-GGA-D | B3LYP-GD3BJ | 4.12 | 1.23 | 0.88 | –3.16 |
All values are given in kcal/mol, except R which is presented as fraction of 1.0.
Figure 5
Figure 5. Pearson’s correlation (R) of B3LYP-GD3Bj with SDD & 6-31+G(d) basis set for BDE calculation of the M–O2 bond.
2.2. Evaluation of DFs with SDD & 6-31+G(d) Basis Sets
classes of DFT | DFs | rmsd | SD | R | MAE |
---|---|---|---|---|---|
LDA | LSDA | 16.66 | 19.86 | 0.02 | –12.45 |
GGA | BP86 | 14.36 | 27.08 | 0.48 | 2.37 |
meta-GGA | TPSSTPSS | 7.85 | 18.97 | 0.03 | –6.59 |
H-GGA | B3LYP | 17.22 | 27.95 | 0.44 | 4.29 |
B3PW91 | 14.92 | 24.03 | 0.42 | 3.75 | |
PBE0 | 6.97 | 18.78 | 0.02 | –6.08 | |
MPWPW91 | 6.25 | 18.93 | –0.02 | –5.49 | |
B97 | 16.33 | 29.31 | 0.40 | 3.73 | |
GH meta-GGA | M05-2X | 5.06 | 28.22 | 0.04 | –4.44 |
M05 | 6.20 | 9.57 | 0.39 | 0.05 | |
M06-2X | 18.65 | 16.65 | 0.82 | 6.49 | |
M06 | 3.98 | 23.99 | 0.13 | –2.20 | |
RS H-GGA | CAM-B3LYP | 4.70 | 22.51 | 0.04 | –4.08 |
H-GGA-D | B3LYP-GD3BJ | 6.62 | 2.47 | 0.35 | –5.81 |
All values are given in kcal/mol, except R which is presented as fraction of 1.0.
Figure 6
Figure 6. Structures of transition-metal complexes having oxygen molecule with known experimental BDEs of M–O2 bond. In all the complexes, O2 show side on the binding model with selected transition-metal complexes.
3. Conclusions


4. Computational Methodologies Adopted in the Current Work
classes of DFT | DFs |
---|---|
local density approximation (LDA) | LSDA |
generalized gradient approximation (GGA) | BP86 |
meta generalized gradient approximation (meta-GGA) | TPSSTPSS |
hybrid generalized gradient approximation (H-GGA) | B3LYP |
B3PW91 | |
PBE0 | |
MPWPW91 | |
B97 | |
global hybrid generalized gradient approximation (GH meta-GGA) | M05-2X |
M05 | |
M06-2X | |
M06 | |
range separated hybrid generalized gradient approximation (RS H-GGA) | CAM-B3LYP |
dispersion corrected hybrid GGA (H-GGA-D) | B3LYP-GD3BJ |
Figure 7
Figure 7. Modal reaction for M–O2 bond dissociation of the rhenium-based complex with oxygen molecule. The rest of the complexes follow the same pattern.
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsomega.2c01331.
Rmsd; SD; Pearson’s correlation (R); MAE; and experimental BDEs of all selected transition-metal complexes with oxygen molecules (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
The author from the King Khalid University extends his appreciation to Deanship of Scientific Research at the King Khalid University for funding this work through Large Groups RGP.2/194/43. Additionally, Higher Education Commission of Pakistan and COMSATS University Islamabad, Abbottabad Campus are also acknowledged for financial and technical support.
References
This article references 55 other publications.
- 1Shahbaz, M.; Mutascu, M.; Azim, P. Environmental Kuznets Curve in Romania and the Role of Energy Consumption. Renewable Sustainable Energy Rev. 2013, 18, 165– 173, DOI: 10.1016/j.rser.2012.10.012Google ScholarThere is no corresponding record for this reference.
- 2Bard, A. J.; Fox, M. A. Artificial Photosynthesis: Solar Splitting of Water to Hydrogen and Oxygen. Acc. Chem. Res. 1995, 28, 141– 145, DOI: 10.1021/ar00051a007Google Scholar2https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXjvFOksrc%253D&md5=62af0ea101f65e4fb2932b80cb29e6ebArtificial Photosynthesis: Solar Splitting of Water to Hydrogen and OxygenBard, Allen J.; Fox, Marye AnneAccounts of Chemical Research (1995), 28 (3), 141-5CODEN: ACHRE4; ISSN:0001-4842. (American Chemical Society)A review with 33 refs. of history and progress of water splitting, semiconductor solid state photovoltaic based systems, semiconductor electrode liq. junction systems, semiconductor particle systems, sensitized semiconductor systems, homogeneous and microheterogeneous systems, and future prospects.
- 3Chow, J.; Kopp, R. J.; Portney, P. R. Energy Resources and Global Development. Science 2003, 302, 1528– 1531, DOI: 10.1126/science.1091939Google Scholar3https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD3srntF2gsw%253D%253D&md5=38a0f15c28ab1b0135f4c53c76049f17Energy resources and global developmentChow Jeffrey; Kopp Raymond J; Portney Paul RScience (New York, N.Y.) (2003), 302 (5650), 1528-31 ISSN:.In order to address the economic and environmental consequences of our global energy system, we consider the availability and consumption of energy resources. Problems arise from our dependence on combustible fuels, the environmental risks associated with their extraction, and the environmental damage caused by their emissions. Yet no primary energy source, be it renewable or nonrenewable, is free of environmental or economic limitations. As developed and developing economies continue to grow, conversion to and adoption of environmentally benign energy technology will depend on political and economic realities.
- 4Navarro Yerga, R. M.; Álvarez Galván, M. C.; del Valle, F.; Villoria de la Mano, J. A.; Fierro, J. L. G. Water Splitting on Semiconductor Catalysts under Visible-Light Irradiation. ChemSusChem 2009, 2, 471– 485, DOI: 10.1002/cssc.200900018Google Scholar4https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXnvV2nt7k%253D&md5=0faa86f3527feeb18239f291fd6e74ecWater Splitting on Semiconductor Catalysts under Visible-Light IrradiationNavarro Yerga, Rufino M.; Alvarez Galvan, M. Consuelo; del Valle, F.; Villoria de la Mano, Jose A.; Fierro, Jose L. G.ChemSusChem (2009), 2 (6), 471-485CODEN: CHEMIZ; ISSN:1864-5631. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. The conversion of solar energy into hydrogen via the water-splitting process, assisted by photo-semiconductor catalysts, is one of the most promising technologies for the future because large quantities of hydrogen can potentially be generated in a clean and sustainable manner. This paper provides an overview of the principles, approaches, and research progress on solar hydrogen prodn. via the water-splitting reaction on photo-semiconductor catalysts. It presents a survey of the advances made over the last decades in the development of catalysts for photochem. water splitting under visible-light irradn. The paper also analyzes the energy requirements and main factors that det. the activity of photocatalysts in the conversion of water into hydrogen and oxygen using sunlight. Remarkable progress has been made since the pioneering work by Fujishima and Honda in 1972, but his development of photocatalysts with improved efficiencies for hydrogen prodn. from water using solar energy still faces major challenges. Research strategies and approaches adopted in the search for active and efficient photocatalysts, for example through new materials and synthesis methods, are presented and analyzed.
- 5Barton, E. E.; Rampulla, D. M.; Bocarsly, A. B. Selective Solar-Driven Reduction of CO 2 to Methanol Using a Catalyzed p-GaP Based Photoelectrochemical Cell. J. Am. Chem. Soc. 2008, 130, 6342– 6344, DOI: 10.1021/ja0776327Google Scholar5https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXltFens7w%253D&md5=ba5b46e11004ee2431e64b139112d411Selective Solar-Driven Reduction of CO2 to Methanol Using a Catalyzed p-GaP Based Photoelectrochemical CellBarton, Emily E.; Rampulla, David M.; Bocarsly, Andrew B.Journal of the American Chemical Society (2008), 130 (20), 6342-6344CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)With rising atm. CO2 levels there is interest in artificial photosynthetic schemes for converting this greenhouse gas into fuel and small orgs. Photoelectrochem. schemes for activating the inert CO2 mol. operate at high overpotentials and thus do not convert light energy to chem. energy. A selective conversion of CO2 to MeOH at a p-GaP semiconductor electrode with a homogeneous pyridinium ion catalyst, driving the reaction with light energy to yield faradaic efficiencies near 100% at potentials well below the std. potential, is described.
- 6Rashatasakhon, P.; Ozdemir, A. D.; Willis, J.; Padwa, A. Six- versus Five-Membered Ring Formation in Radical Cyclizations of 7-Bromo-Substituted Hexahydroindolinones. Org. Lett. 2004, 6, 917– 920, DOI: 10.1021/ol036347zGoogle Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXht1WitL8%253D&md5=6ec21c4ff38a7d5d1aa702455c0f8aa4Six- versus Five-Membered Ring Formation in Radical Cyclizations of 7-Bromo-Substituted HexahydroindolinonesRashatasakhon, Paitoon; Ozdemir, Ayse Daut; Willis, Jerremey; Padwa, AlbertOrganic Letters (2004), 6 (6), 917-920CODEN: ORLEF7; ISSN:1523-7060. (American Chemical Society)Radical cyclization of N-allyl-7-bromo-3a-methyl-hexahydroindol-2-one affords a six-membered ring product that prevails over the isomeric five-membered compd. The former product is generated through two reaction pathways: (a) 6-endo-trig ring closure and (b) rearrangement of an intermediate methylenecyclopentyl radical obtained by 5-exo-trig cyclization.
- 7Service, R. F. Turning Over a New Leaf. Science 2011, 334, 925– 927, DOI: 10.1126/science.334.6058.925Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhs1eit7bE&md5=547ccba7aeb313c77a9f8262fc5a1ff0Turning over a new leafService, Robert F.Science (Washington, DC, United States) (2011), 334 (6058), 925-927CODEN: SCIEAS; ISSN:0036-8075. (American Association for the Advancement of Science)There is no expanded citation for this reference.
- 8Kurz, P.; Berggren, G.; Anderlund, M. F.; Styring, S. Oxygen Evolving Reactions Catalysed by Synthetic Manganese Complexes: A Systematic Screening. Dalton Trans. 2007, 38, 4258, DOI: 10.1039/b710761gGoogle ScholarThere is no corresponding record for this reference.
- 9Zong, R.; Thummel, R. P. A New Family of Ru Complexes for Water Oxidation. J. Am. Chem. Soc. 2005, 127, 12802– 12803, DOI: 10.1021/ja054791mGoogle Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpt1Cnsb0%253D&md5=2ace4686bbdedbe6af1c7ffb439fce24A New Family of Ru Complexes for Water OxidationZong, Ruifa; Thummel, Randolph P.Journal of the American Chemical Society (2005), 127 (37), 12802-12803CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Reaction of the bridging ligand 3,6-bis-[6'-(1'',8''-naphthyrid-2''-yl)-pyrid-2'-yl]pyridazine (1) with [Ru(DMSO)4Cl2] in aq. ethanol followed by excess 4-substituted pyridine (4-R-py) in the presence of triethylamine provides a series of three well-organized dinuclear complexes characterized by 1H NMR, MS, and X-ray. Mononuclear analogs are prepd. from 4-tert-butyl-2,6-di(1',8'-naphthyrid-2'-yl)pyridine (5) and characterized in a similar fashion. All six complexes show electronic absorption and redox properties consistent with the electron donor/acceptor ability of the axial 4-R-py ligand. When an acetonitrile soln. of the catalyst is added to an aq. Ce(IV)-CF3SO3H soln. (pH = 1.0) at 24 °C, oxygen evolution is obsd. for both mono and dinuclear systems. Turnover nos. range from 50 to 3200 with the best results being found when the axial ligand is 4-methylpyridine (mononuclear TN = 580 and dinuclear TN = 3200).
- 10Polyansky, D. E.; Muckerman, J. T.; Rochford, J.; Zong, R.; Thummel, R. P.; Fujita, E. Water Oxidation by a Mononuclear Ruthenium Catalyst: Characterization of the Intermediates. J. Am. Chem. Soc. 2011, 133, 14649– 14665, DOI: 10.1021/ja203249eGoogle Scholar10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhtV2isbjP&md5=6a2129651f23f68e39edcdd2c5e58ec0Water Oxidation by a Mononuclear Ruthenium Catalyst: Characterization of the IntermediatesPolyansky, Dmitry E.; Muckerman, James T.; Rochford, Jonathan; Zong, Ruifa; Thummel, Randolph P.; Fujita, EtsukoJournal of the American Chemical Society (2011), 133 (37), 14649-14665CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)A detailed characterization of intermediates in water oxidn. catalyzed by a mononuclear Ru polypyridyl complex [RuII-OH2]2+ (Ru = Ru complex with one 4-t-butyl-2,6-di-(1',8'-naphthyrid-2'-yl)-pyridine ligand and two 4-picoline ligands) has been carried out using electrochem., UV-vis and resonance Raman spectroscopy, pulse radiolysis, stopped flow, and electrospray ionization mass spectrometry (ESI-MS) with H218O labeling expts. and theor. calcns. The results reveal a no. of intriguing properties of intermediates such as [RuIV=O]2+ and [RuIV-OO]2+. At pH > 2.9, two consecutive proton-coupled one-electron steps take place at the potential of the [RuIII-OH]2+/[RuII-OH2]2+ couple, which is equal to or higher than the potential of the [RuIV=O]2+/[RuIII-OH]2+ couple (i.e., the observation of a two-electron oxidn. in cyclic voltammetry). At pH 1, the rate const. of the first one-electron oxidn. by Ce(IV) is k1 = 2 × 104 M-1 s-1. While pH-independent oxidn. of [RuIV=O]2+ takes place at 1420 mV vs NHE, bulk electrolysis of [RuII-OH2]2+ at 1260 mV vs NHE at pH 1 (0.1 M triflic acid) and 1150 mV at pH 6 (10 mM sodium phosphate) yielded a red colored soln. with a Coulomb count corresponding to a net four-electron oxidn. ESI-MS with labeling expts. clearly indicates that this species has an O-O bond. This species required an addnl. oxidn. to liberate an oxygen mol., and without any addnl. oxidant it completely decompd. slowly to form [RuII-OOH]+ over 2 wk. While there remains some conflicting evidence, we have assigned this species as 1[RuIV-η2-OO]2+ based on our electrochem., spectroscopic, and theor. observations alongside a previously reported anal.
- 11Baran, J. D.; Grönbeck, H.; Hellman, A. Analysis of Porphyrines as Catalysts for Electrochemical Reduction of O 2 and Oxidation of H 2 O. J. Am. Chem. Soc. 2014, 136, 1320– 1326, DOI: 10.1021/ja4060299Google Scholar11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXksVag&md5=e249eaac0eeb7649475c68e94850177fAnalysis of Porphyrins as Catalysts for Electrochemical Reduction of O2 and Oxidation of H2OBaran, Jakub D.; Groenbeck, Henrik; Hellman, AndersJournal of the American Chemical Society (2014), 136 (4), 1320-1326CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Bioinspired structures are promising as improved catalysts for various redox reactions. One example is metal hangman-porphyrins (MHP), which recently were suggested for O redn./evolution reaction (ORR/OER). The unique properties of the MHPs are attributed to both the hangman scaffold and the C6F5 side groups. Herein, the OER/ORR over various transition metal MHPs was studied by d. functional theory calcns. within an electrochem. framework. A comparison of the reaction landscape for MHP, metal porphyrins (MP) and metal tetrafluorophenyloporphyrins (MTFPP), allow for a disentanglement of the different roles of the hangman motif and the side groups. In agreement with exptl. studies, Fe and Co are the best MHP metal centers to catalyze these reactions. The addn. of the three-dimensional moiety as hangman scaffold does not break the apparently universal energy relation between *OH and *OOH intermediates. However, the hangman motif is found to stabilize the O intermediate, whereas addn. of C6F5 groups reduces the binding energy of all reaction intermediates. The authors' results indicate that the combination of these two effects allow new design possibilities for macromol. systems with enhanced catalytic OER/ORR activity.
- 12Mahmood, T.; Kosar, N.; Ayub, K. DFT Study of Acceleration of Electrocyclization in Photochromes under Radical Cationic Conditions: Comparison with Recent Experimental Data. Tetrahedron 2017, 73, 3521– 3528, DOI: 10.1016/j.tet.2017.05.031Google Scholar12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXnvFGhsbc%253D&md5=539b2d64d314ab526c84d46409eea38aDFT study of acceleration of electrocyclization in photochromes under radical cationic conditions: Comparison with recent experimental dataMahmood, Tariq; Kosar, Naveen; Ayub, KhurshidTetrahedron (2017), 73 (25), 3521-3528CODEN: TETRAB; ISSN:0040-4020. (Elsevier Ltd.)Radical cation formation is proposed for the rapid cyclization of 1, 2-bis[5-phenyl-2-methylthien-3-yl]cyclopentene and oligothiophene functionalized dimethyldihydropyrenes (DMDHP). D. functional theory calcns. have been performed to rationalize the effect of a radical cation on the activation barrier of different classes of electrocyclic photochromes (DHP, dithienylethene, dihydroazulene and fulgide). For exact comparative anal., the activation barrier of neutral (singlet) analogs at the same level of theory are also calcd. In addn., the concerted nature and aromaticity of transition states were investigated with the help of synchronicity (Sy.) and nuclear independent chem. shift values NICS(0) calcns., resp., for both the radical cation and neutral systems. In case of the radical cation, thermal return of CPD to DHP, the activation barrier is very low (ΔH = 3.13 kcal mol-1, ΔG = 4.01 kcal mol-1) as compared to the neutral analog (ΔH = 20.6 kcal mol-1, ΔG = 20.98 kcal mol-1), which is consistent with exptl. observations. Similarly for dithenylethenes, radical cation formation has a large impact on the activation barrier (ΔH = 19.44 kcal mol-1, ΔG = 22.29 kcal mol-1). However, radical cation formation has almost negligible impact on the activation barrier of VHF-DHA and fulgide isomerization. The significant difference has been obsd. for synchronicity and NICS(0) values of all types of photochromes under radical cation conditions as compared to the neutral system.
- 13Chiatti, F.; Delle Piane, M.; Ugliengo, P.; Corno, M. Water at Hydroxyapatite Surfaces: The Effect of Coverage and Surface Termination as Investigated by All-Electron B3LYP-D* Simulations. Theor. Chem. Acc. 2016, 135, 54, DOI: 10.1007/s00214-016-1818-8Google ScholarThere is no corresponding record for this reference.
- 14Civalleri, B.; Maschio, L.; Ugliengo, P.; Zicovich-Wilson, C. M. Role of Dispersive Interactions in the CO Adsorption on MgO(001): Periodic B3LYP Calculations Augmented with an Empirical Dispersion Term. Phys. Chem. Chem. Phys. 2010, 12, 6382, DOI: 10.1039/c001192dGoogle Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXntF2htrY%253D&md5=4e2b5f13bf5d1d4742a54ddda90d47b9Role of dispersive interactions in the CO adsorption on MgO(001): periodic B3LYP calculations augmented with an empirical dispersion termCivalleri, Bartolomeo; Maschio, Lorenzo; Ugliengo, Piero; Zicovich-Wilson, Claudio M.Physical Chemistry Chemical Physics (2010), 12 (24), 6382-6386CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Empirically dispersion cor. B3LYP method (i.e. B3LYP-D) is demonstrated to give excellent results for structure, adsorption energy and vibrational frequency shift for the CO mol. adsorbed on the MgO(001) surface, a system considered a challenge for c.d. functional methods. A periodic approach was adopted to model the interaction using a three-layer slab model. For the B3LYP-D* method an interaction energy of -13.1 kJ mol-1 is computed at low-coverage in very good agreement with exptl. evidence (-12.6 kJ mol-1) as well as a pos. CO vibrational shift of 10 cm-1 to be compared with the exptl. value of 14 cm-1.
- 15Lousada, C. M.; Johansson, A. J.; Brinck, T.; Jonsson, M. Reactivity of Metal Oxide Clusters with Hydrogen Peroxide and Water – a DFT Study Evaluating the Performance of Different Exchange–Correlation Functionals. Phys. Chem. Chem. Phys. 2013, 15, 5539, DOI: 10.1039/c3cp44559cGoogle Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXkt1Gnsrw%253D&md5=aae4ff6dd76e2c755a726dd9bbf6c2f0Reactivity of metal oxide clusters with hydrogen peroxide and water - a DFT study evaluating the performance of different exchange-correlation functionalsLousada, Claudio M.; Johansson, Adam Johannes; Brinck, Tore; Jonsson, MatsPhysical Chemistry Chemical Physics (2013), 15 (15), 5539-5552CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We have performed a d. functional theory (DFT) investigation of the interactions of H2O2, H2O and HO radicals with clusters of ZrO2, TiO2 and Y2O3. Different modes of H2O adsorption onto the clusters were studied. In almost all the cases the dissociative adsorption is more exothermic than mol. adsorption. At the surfaces where H2O has undergone dissociative adsorption, the adsorption of H2O2 and the transition state for its decompn. are mediated by hydrogen bonding with the surface HO groups. Using the functionals B3LYP, B3LYP-D and M06 with clusters of 26 and 8 units of ZrO2, the M06 functional performed better than B3LYP in describing the reaction of decompn. of H2O2 and the adsorption of H2O. Addnl., we investigated clusters of the type (ZrO2)2, (TiO2)2 and (Y2O3) and the performance of the functionals B3LYP, B3LYP-D, B3LYP*, M06, M06-L, PBE0, PBE and PWPW91 in describing H2O2, H2O and HO√ adsorption and the energy barrier for decompn. of H2O2. The trends obtained for HO√ adsorption onto the clusters are discussed in terms of the ionization energy of the metal cation present in the oxide. In order to correctly account for the existence of an energy barrier for the decompn. of H2O2, the functional used must include Hartree-Fock exchange. Using minimal cluster models, the best performance in describing the energy barrier for H2O2 decompn. was obtained with the M06 and PBE0 functionals - the av. abs. deviations from expts. are 6 kJ mol-1 and 5 kJ mol-1 resp. With the M06 functional and a larger monoclinic (ZrO2)8 cluster model, the performance is in excellent agreement with exptl. data. For the different oxides, PBE0 was found to be the most effective functional in terms of performance and computational time cost.
- 16Minkin, V. I.; Starikov, A. G.; Starikova, A. A. Computational Insight into Magnetic Behavior and Properties of the Transition Metal Complexes with Redox-Active Ligands: A DFT Approach. Pure Appl. Chem. 2018, 90, 811– 824, DOI: 10.1515/pac-2017-0803Google Scholar16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXptV2mu7g%253D&md5=b4217a1239e16b5e8a56df8ffe7129ceComputational insight into magnetic behavior and properties of the transition metal complexes with redox-active ligands: a DFT approachMinkin, Vladimir I.; Starikov, Andrey G.; Starikova, Alyona A.Pure and Applied Chemistry (2018), 90 (5), 811-824CODEN: PACHAS; ISSN:0033-4545. (Walter de Gruyter, Inc.)Various aspects related to the use of DFT method for the study of magnetic, geometry and energetic properties of transition metal complexes with redox-active ligands are considered. Particular attention is given to the correct choice of model compds. and methodol. of the calcns.
- 17Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104, DOI: 10.1063/1.3382344Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXkvVyks7o%253D&md5=2bca89d904579d5565537a0820dc2ae8A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-PuGrimme, Stefan; Antony, Jens; Ehrlich, Stephan; Krieg, HelgeJournal of Chemical Physics (2010), 132 (15), 154104/1-154104/19CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The method of dispersion correction as an add-on to std. Kohn-Sham d. functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coeffs. and cutoff radii that are both computed from first principles. The coeffs. for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination nos. (CN). They are used to interpolate between dispersion coeffs. of atoms in different chem. environments. The method only requires adjustment of two global parameters for each d. functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of at. forces. Three-body nonadditivity terms are considered. The method has been assessed on std. benchmark sets for inter- and intramol. noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean abs. deviations for the S22 benchmark set of noncovalent interactions for 11 std. d. functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C6 coeffs. also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in mols. and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems. (c) 2010 American Institute of Physics.
- 18Starikov, A. G.; Ivanov, D. G.; Starikova, A. A.; Minkin, V. I. Dispersion Interactions in Oligomerization of Metal Diketonates: A DFT Evaluation. Chem. Pap. 2018, 72, 829– 839, DOI: 10.1007/s11696-017-0225-5Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVShsrjP&md5=dc816f1d25691017e5e3f96a18662106Dispersion interactions in oligomerization of metal diketonates: a DFT evaluationStarikov, Andrey G.; Ivanov, Dmitrii G.; Starikova, Alyona A.; Minkin, Vladimir I.Chemical Papers (2018), 72 (4), 829-839CODEN: CHPAEG; ISSN:1336-9075. (Springer International Publishing AG)The structures of the oligomers of Co(II), Ni(II), Cu(II) and Zn(II) β-diketonates and their stability with respect to dissocn. into the monomeric components were studied with the use of d. functional theory [DFT B3LYP/6-311++G(d,p)] method taking addnl. into account corrections for dispersion interactions in CAM and D3BJ approxns. The calcd. geometries of the oligomers well match those exptl. detd., but stabilization energies obtained at the B3LYP-CAM and B3LYP-D3BJ approxns. are significantly overestimated.
- 19Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157– 167, DOI: 10.1021/ar700111aGoogle Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXksV2iug%253D%253D&md5=95aa78f44294d3eed013ff9719cafe76Density Functionals with Broad Applicability in ChemistryZhao, Yan; Truhlar, Donald G.Accounts of Chemical Research (2008), 41 (2), 157-167CODEN: ACHRE4; ISSN:0001-4842. (American Chemical Society)A review. Although d. functional theory is widely used in the computational chem. community, the most popular d. functional, B3LYP, has some serious shortcomings: (i) it is better for main-group chem. than for transition metals; (ii) it systematically underestimates reaction barrier heights; (iii) it is inaccurate for interactions dominated by medium-range correlation energy, such as van der Waals attraction, arom.-arom. stacking, and alkane isomerization energies. We have developed a variety of databases for testing and designing new d. functionals. We used these data to design new d. functionals, called M06-class (and, earlier, M05-class) functionals, for which we enforced some fundamental exact constraints such as the uniform-electron-gas limit and the absence of self-correlation energy. Our M06-class functionals depend on spin-up and spin-down electron densities (i.e., spin densities), spin d. gradients, spin kinetic energy densities, and, for nonlocal (also called hybrid) functionals, Hartree-Fock exchange. We have developed four new functionals that overcome the above-mentioned difficulties: (a) M06, a hybrid meta functional, is a functional with good accuracy "across-the-board" for transition metals, main group thermochem., medium-range correlation energy, and barrier heights; (b) M06-2X, another hybrid meta functional, is not good for transition metals but has excellent performance for main group chem., predicts accurate valence and Rydberg electronic excitation energies, and is an excellent functional for arom.-arom. stacking interactions; (c) M06-L is not as accurate as M06 for barrier heights but is the most accurate functional for transition metals and is the only local functional (no Hartree-Fock exchange) with better across-the-board av. performance than B3LYP; this is very important because only local functionals are affordable for many demanding applications on very large systems; (d) M06-HF has good performance for valence, Rydberg, and charge transfer excited states with minimal sacrifice of ground-state accuracy. In this Account, we compared the performance of the M06-class functionals and one M05-class functional (M05-2X) to that of some popular functionals for diverse databases and their performance on several difficult cases. The tests include barrier heights, conformational energy, and the trend in bond dissocn. energies of Grubbs' ruthenium catalysts for olefin metathesis. Based on these tests, we recommend (1) the M06-2X, BMK, and M05-2X functionals for main-group thermochem. and kinetics, (2) M06-2X and M06 for systems where main-group thermochem., kinetics, and noncovalent interactions are all important, (3) M06-L and M06 for transition metal thermochem., (4) M06 for problems involving multireference rearrangements or reactions where both org. and transition-metal bonds are formed or broken, (5) M06-2X, M05-2X, M06-HF, M06, and M06-L for the study of noncovalent interactions, (6) M06-HF when the use of full Hartree-Fock exchange is important, for example, to avoid the error of self-interaction at long-range, (7) M06-L when a local functional is required, because a local functional has much lower cost for large systems.
- 20Lonsdale, R.; Harvey, J. N.; Mulholland, A. J. Inclusion of Dispersion Effects Significantly Improves Accuracy of Calculated Reaction Barriers for Cytochrome P450 Catalyzed Reactions. J. Phys. Chem. Lett. 2010, 1, 3232– 3237, DOI: 10.1021/jz101279nGoogle Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhtleitbvK&md5=9bc8d74cf3976d4da46535ec7d926182Inclusion of Dispersion Effects Significantly Improves Accuracy of Calculated Reaction Barriers for Cytochrome P450 Catalyzed ReactionsLonsdale, Richard; Harvey, Jeremy N.; Mulholland, Adrian J.Journal of Physical Chemistry Letters (2010), 1 (21), 3232-3237CODEN: JPCLCD; ISSN:1948-7185. (American Chemical Society)Prediction of cytochrome P 450 reactivity is of great importance to the development of new medicinal compds. D. functional theory (DFT) has proven itself as a useful tool in the characterization of the elusive reactive species, compd. I, and of the mechanisms of substrate oxidn. B3LYP is the most widely used d. functional in the study of P 450s; however, a major drawback of B3LYP is its inaccurate treatment of dispersion, leading to discrepancies between expt. and theory in some systems. Recent work has shown that an added empirical dispersion correction to B3LYP (B3LYP-D) yields more promising results for similar systems. In the present work, two previously studied systems, camphor hydroxylation and alkene oxidn., have been recalcd. using B3LYP-D. Our work shows that inclusion of dispersion has a significant effect on the energies and geometries of transition states and encounter complexes; furthermore, an improved agreement with exptl. data is obsd.
- 21Stein, T.; Kronik, L.; Baer, R. Reliable Prediction of Charge Transfer Excitations in Molecular Complexes Using Time-Dependent Density Functional Theory. J. Am. Chem. Soc. 2009, 131, 2818– 2820, DOI: 10.1021/ja8087482Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhvVSns7s%253D&md5=b01ac7f2580fdb9495f7e51be1910d57Reliable prediction of charge transfer excitations in molecular complexes using time-dependent density functional theoryStein, Tamar; Kronik, Leeor; Baer, RoiJournal of the American Chemical Society (2009), 131 (8), 2818-2820CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)We show how charge transfer excitations at mol. complexes can be calcd. quant. using time-dependent d. functional theory. Predictive power is obtained from range-sepd. hybrid functionals using nonempirical tuning of the range-splitting parameter. Excellent performance of this approach is obtained for a series of complexes composed of various arom. donors and the tetracyanoethylene acceptor, paving the way to systematic nonempirical quant. studies of charge-transfer excitations in real systems.
- 22Hirao, H. Which DFT Functional Performs Well in the Calculation of Methylcobalamin? Comparison of the B3LYP and BP86 Functionals and Evaluation of the Impact of Empirical Dispersion Correction. J. Phys. Chem. A 2011, 115, 9308– 9313, DOI: 10.1021/jp2052807Google Scholar22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXpsFOgtr4%253D&md5=37163d022b89bd1f37e6a1d8a97029a9Which DFT Functional Performs Well in the Calculation of Methylcobalamin? Comparison of the B3LYP and BP86 Functionals and Evaluation of the Impact of Empirical Dispersion CorrectionHirao, HajimeJournal of Physical Chemistry A (2011), 115 (33), 9308-9313CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Controversy remains regarding the suitable d. functionals for the calcn. of vitamin B12 systems that contain cobalt. To identify the optimum functionals, geometry optimization calcns. were performed on a full-size model of methylcobalamin (MeCbl) using the B3LYP, B3LYP-D, BP86, and BP86-D methods in conjunction with the 6-31G* basis set. Single-point energy evaluations were also performed with the 6-311+G(2d,p) basis set. Consistent with previous studies, the BP86-optimized geometry showed fairly good agreement with the exptl. geometry. Various factors that may influence the homolytic bond dissocn. energy (BDE) of the Co-C bond of MeCbl were systematically evaluated with these methods. Our anal. demonstrated that dispersion was the largest correction term that influenced the magnitude of BDE. Previous studies have shown that B3LYP significantly underestimates BDE, whereas BP86 gives BDE values that are fairly close to the exptl. values (36-37 kcal/mol). The same trend in the relative magnitudes of the BDEs was obsd. in the present calcns. However, BP86 underestimated the BDE for a full model of MeCbl. When the amt. of Hartree-Fock exchange in the B3LYP functional was reduced to 15% and the dispersion correction was made (i.e., B3LYP*-D), the calcd. BDE was in good accord with exptl. values. B3P86-D also performed well. A detailed anal. was undertaken to det. which atoms in cobalamin have large dispersion interactions with a Me fragment of MeCbl.
- 23Minenkov, Y.; Occhipinti, G.; Jensen, V. R. Metal–Phosphine Bond Strengths of the Transition Metals: A Challenge for DFT. J. Phys. Chem. A 2009, 113, 11833– 11844, DOI: 10.1021/jp902940cGoogle Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhtV2mtb7L&md5=c0df6248a1cdde1cf246c445f7dbf4a3Metal-phosphine bond strengths of the transition metals: A challenge for DFTMinenkov, Yury; Occhipinti, Giovanni; Jensen, Vidar R.Journal of Physical Chemistry A (2009), 113 (43), 11833-11844CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Previous promising tests of the new M06 family of functionals in predicting ruthenium-metal phosphine bond dissocn. energies have been extended to a series of phosphine complexes of chromium, molybdenum, nickel, and ruthenium for which relevant exptl. data are available. In addn. to the M06 family of functionals, bond dissocn. enthalpies have been calcd. using a selection of d. functionals and hybrid functionals based on the generalized gradient approxn. (GGA), and with or without an empirical term (i.e., DFT-D) accounting for long-range dispersion. For the ruthenium complexes, second-order Moller-Plesset perturbation theory (MP2) has also been applied. Electrostatic and nonelectrostatic solvent effects have been estd. using the polarizable continuum model (PCM), allowing for comparison with exptl. data obtained for dissocn. reactions in org. solvents. Whereas the GGA and hybrid-GGA functionals grossly underestimate the abs. metal-phosphine bond enthalpies, with mean unsigned errors (MUEs) for a set of 10 phosphine dissocn. reactions in the range 13-27 kcal/mol, the recently developed DFT-based methods for inclusion of attractive noncovalent interactions and dispersion (the DFT-D and M06 functionals) dramatically improve upon the situation. The best agreement with expt. is obsd. for BLYP-D (MUE = 2.2 kcal/mol), and with the exception for M06-2X, all these methods provide MUEs well below 5 kcal/mol, which should be sufficient for a broad range of applications. The improvements in predicted relative bond enthalpies are less convincing, however. In several cases the GGA and hybrid-GGA functionals are better at reproducing substitution effects than the DFT-D and M06 methods.
- 24Grimme, S. Accurate Description of van Der Waals Complexes by Density Functional Theory Including Empirical Corrections. J. Comput. Chem. 2004, 25, 1463– 1473, DOI: 10.1002/jcc.20078Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmtFKgt78%253D&md5=16e36ed7a1c098813d7d947ac72b9078Accurate description of van der Waals complexes by density functional theory including empirical correctionsGrimme, StefanJournal of Computational Chemistry (2004), 25 (12), 1463-1473CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)An empirical method to account for van der Waals interactions in practical calcns. in the framework of the d. functional theory (termed DFT-D) was tested for a wide variety of mol. complexes. As in previous schemes, the dispersive energy was described by damped interat. potentials of the form C6R-6. The use of pure, gradient-cor. d. functionals (BLYP and PBE), together with the resoln.-of-the-identity (RI) approxn. for the Coulomb operator, allows very efficient computations for large systems. In contrast to the previous work, extended AO basis sets of polarized TZV or QZV quality were employed, which reduced the basis set superposition error to a negligible extend. By using a global scaling factor for the at. C6 coeffs., the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes was found to be insignificant if steep damping functions were employed. The method was applied to a total of 29 complexes of atoms and small mols. (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene·H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theor. or exptl. results for binding energies and intermol. distances is obtained. For stacked arom. systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to std. MP2 treatments that systematically over-bind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chem. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calcns. of the interaction energies.
- 25Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787– 1799, DOI: 10.1002/jcc.20495Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XhtFenu7bO&md5=0b4aa16bebc3a0a2ec175d4b161ab0e4Semiempirical GGA-type density functional constructed with a long-range dispersion correctionGrimme, StefanJournal of Computational Chemistry (2006), 27 (15), 1787-1799CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)A new d. functional (DF) of the generalized gradient approxn. (GGA) type for general chem. applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C6·R-6. A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common d. functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on std. thermochem. benchmark sets, for 40 noncovalently bound complexes, including large stacked arom. mols. and group II element clusters, and for the computation of mol. geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for std. functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean abs. deviation of only 3.8 kcal mol-1. The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the av. CCSD(T) accuracy. The basic strategy in the development to restrict the d. functional description to shorter electron correlation lengths scales and to describe situations with medium to large interat. distances by damped C6·R-6 terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chem. method for large systems where dispersion forces are of general importance.
- 26Johnson, E. R.; Becke, A. D. A Post-Hartree-Fock Model of Intermolecular Interactions: Inclusion of Higher-Order Corrections. J. Chem. Phys. 2006, 124, 174104, DOI: 10.1063/1.2190220Google Scholar26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XksFChurY%253D&md5=99cb745d5a698e2c43e8d1a2bc732d3eA post-Hartree-Fock model of intermolecular interactions: Inclusion of higher-order correctionsJohnson, Erin R.; Becke, Axel D.Journal of Chemical Physics (2006), 124 (17), 174104/1-174104/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We have previously demonstrated that the dipole moment of the exchange hole can be used to derive intermol. C6 dispersion coeffs. [J. Chem. Phys. 122, 154104 (2005)]. This was subsequently the basis for a novel post-Hartree-Fock model of intermol. interactions [J. Chem. Phys. 123, 024101 (2005)]. In the present work, the model is extended to include higher-order dispersion coeffs. C8 and C10. The extended model performs very well for prediction of intermonomer sepns. and binding energies of 45 van der Waals complexes. In particular, it performs twice as well as basis-set extrapolated MP2 theory for dispersion-bound complexes, with minimal computational cost.
- 27Fey, N.; Ridgway, B. M.; Jover, J.; McMullin, C. L.; Harvey, J. N. Organometallic Reactivity: The Role of Metal–Ligand Bond Energies from a Computational Perspective. Dalton Trans. 2011, 40, 11184, DOI: 10.1039/c1dt10909jGoogle Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhtlejt7vP&md5=8d85df3f5cd06663f9f5d04b747a962bOrganometallic reactivity: the role of metal-ligand bond energies from a computational perspectiveFey, Natalie; Ridgway, Benjamin M.; Jover, Jesus; McMullin, Claire L.; Harvey, Jeremy N.Dalton Transactions (2011), 40 (42), 11184-11191CODEN: DTARAF; ISSN:1477-9226. (Royal Society of Chemistry)A review; the assocn. and dissocn. of ligands plays a vital role in detg. the reactivity of organometallic catalysts. Computational studies with d. functional theory often fail to reproduce exptl. metal-ligand bond energies, but recently functionals which better capture dispersion effects have been developed. Here we explore their application and discuss future challenges for computational studies of organometallic catalysis.
- 28Frisch, 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.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N., Staroverov, V. N., Normand, R. K. J., Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Rega, M. C. N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Gomperts, J. J. R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J. D. J. F. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford CT, 2009.Google ScholarThere is no corresponding record for this reference.
- 29Mo, Y.; Tian, G.; Car, R.; Staroverov, V. N.; Scuseria, G. E.; Tao, J. Performance of a Nonempirical Density Functional on Molecules and Hydrogen-Bonded Complexes. J. Chem. Phys. 2016, 145, 234306, DOI: 10.1063/1.4971853Google Scholar29https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XitFamur%252FN&md5=394568c0db0ff1d048c1ecfd073f7401Performance of a nonempirical density functional on molecules and hydrogen-bonded complexesMo, Yuxiang; Tian, Guocai; Car, Roberto; Staroverov, Viktor N.; Scuseria, Gustavo E.; Tao, JianminJournal of Chemical Physics (2016), 145 (23), 234306/1-234306/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recently, Tao and Mo derived a meta-generalized gradient approxn. functional based on a model exchange-correlation hole. In this work, the performance of this functional is assessed on std. test sets, using the 6-311++G(3df,3pd) basis set. These test sets include 223 G3/99 enthalpies of formation, 99 atomization energies, 76 barrier heights, 58 electron affinities, 8 proton affinities, 96 bond lengths, 82 harmonic vibrational frequencies, 10 hydrogen-bonded mol. complexes, and 22 at. excitation energies. Our calcns. show that the Tao-Mo functional can achieve high accuracy for most properties considered, relative to the local spin-d. approxn., Perdew-Burke-Ernzerhof, and Tao-Perdew-Staroverov-Scuseria functionals. In particular, it yields the best accuracy for proton affinities, harmonic vibrational frequencies, hydrogen-bond dissocn. energies and bond lengths, and at. excitation energies. (c) 2016 American Institute of Physics.
- 30Zhao, Y.; Truhlar, D. G. A New Local Density Functional for Main-Group Thermochemistry, Transition Metal Bonding, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2006, 125, 194101, DOI: 10.1063/1.2370993Google Scholar30https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xht1CnsrzL&md5=9106de7f4660e3d76691b03dfb530956A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactionsZhao, Yan; Truhlar, Donald G.Journal of Chemical Physics (2006), 125 (19), 194101/1-194101/18CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a new local d. functional, called M06-L, for main-group and transition element thermochem., thermochem. kinetics, and noncovalent interactions. The functional is designed to capture the main dependence of the exchange-correlation energy on local spin d., spin d. gradient, and spin kinetic energy d., and it is parametrized to satisfy the uniform-electron-gas limit and to have good performance for both main-group chem. and transition metal chem. The M06-L functional and 14 other functionals have been comparatively assessed against 22 energetic databases. Among the tested functionals, which include the popular B3LYP, BLYP, and BP86 functionals as well as our previous M05 functional, the M06-L functional gives the best overall performance for a combination of main-group thermochem., thermochem. kinetics, and organometallic, inorganometallic, biol., and noncovalent interactions. It also does very well for predicting geometries and vibrational frequencies. Because of the computational advantages of local functionals, the present functional should be very useful for many applications in chem., esp. for simulations on moderate-sized and large systems and when long time scales must be addressed.
- 31Tao, J.; Perdew, J. P.; Staroverov, V. N.; Scuseria, G. E. Climbing the Density Functional Ladder: Nonempirical Meta–Generalized Gradient Approximation Designed for Molecules and Solids. Phys. Rev. Lett. 2003, 91, 146401, DOI: 10.1103/PhysRevLett.91.146401Google Scholar31https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXnvVKmur4%253D&md5=ee0673dbddc2c14c6bdc2eddf4945d04Climbing the Density Functional Ladder: Nonempirical Meta-Generalized Gradient Approximation Designed for Molecules and SolidsTao, Jianmin; Perdew, John P.; Staroverov, Viktor N.; Scuseria, Gustavo E.Physical Review Letters (2003), 91 (14), 146401/1-146401/4CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)The electron d., its gradient, and the Kohn-Sham orbital kinetic energy d. are the local ingredients of a meta-generalized gradient approxn. (meta-GGA). We construct a meta-GGA d. functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The exchange and correlation terms respect two paradigms: one- or two-electron densities and slowly varying densities, and so describe both mols. and solids with high accuracy, as shown by extensive numerical tests. This functional completes the third rung of "Jacob's ladder" of approxns., above the local spin d. and GGA rungs.
- 32Zhao, Y.; Lynch, B. J.; Truhlar, D. G. Development and Assessment of a New Hybrid Density Functional Model for Thermochemical Kinetics. J. Phys. Chem. A 2004, 108, 2715– 2719, DOI: 10.1021/jp049908sGoogle Scholar32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXhvF2nu7o%253D&md5=204a63ab21ed0a0ca8c89bee27ed90cfDevelopment and Assessment of a New Hybrid Density Functional Model for Thermochemical KineticsZhao, Yan; Lynch, Benjamin J.; Truhlar, Donald G.Journal of Physical Chemistry A (2004), 108 (14), 2715-2719CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)A new hybrid Hartree-Fock-d. functional model called the Becke88-Becke95 1-parameter model for kinetics (BB1K) was optimized against a database of three forward barrier heights, three reverse barrier heights, and three energies of reaction for the reactions in the BH6 representative barrier height database. We then assessed the newly developed BB1K method against a saddle point geometries database, a database of 42 barrier heights, the AE6 representative atomization energy database, a mol. geometries data set, and a set of 13 zero point energies. The results show that BB1K can give excellent saddle point geometries and barrier heights, and its performance for calcg. atomization energies is 40% better than MPW1K. Using a mean mean unsigned error criterion that equally wts. the errors in barrier heights and in bond energies, the new BB1K method outperforms all other DFT and hybrid DFT methods by a large margin, and we therefore conclude that it is the best d. functional-type method for thermochem. kinetics.
- 33Becke, A. D. New Mixing of Hartree–Fock and Local Density-functional Theories. J. Chem. Phys. 1993, 98, 1372– 1377, DOI: 10.1063/1.464304Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3sXhtlagt7o%253D&md5=02a9289532197e7445af684bb29f4e11A new mixing of Hartree-Fock and local-density-functional theoriesBecke, Axel D.Journal of Chemical Physics (1993), 98 (2), 1372-7CODEN: JCPSA6; ISSN:0021-9606.Previous attempts to combine the Hartree-Fock theory with the local d.-functional theory were unsuccessful in applications to mol. bonding. A new coupling is presented of these two theories that maintains their simplicity and computational efficiency, and yet greatly improves their predictive power. Very encouraging results of tests on atomization energies, ionization potentials, and proton affinities are reported, and the potential for future development is discussed.
- 34Schmider, H. L.; Becke, A. D. Optimized Density Functionals from the Extended G2 Test Set. J. Chem. Phys. 1998, 108, 9624– 9631, DOI: 10.1063/1.476438Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXjsFCrtLs%253D&md5=e57a83f02bc6ee46e6b1275e308ed7e1Optimized density functionals from the extended G2 test setSchmider, Hartmut L.; Becke, Axel D.Journal of Chemical Physics (1998), 108 (23), 9624-9631CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A recently suggested procedure for the systematic optimization of gradient-cor. exchange-correlation functionals [A. D. Becke, J. Chem. Phys. 107, 8554 (1997)] has been applied to the extended G2 test set [L. A. Curtiss et al., J. Chem. Phys. 106, 1063 (1997)], which consists of the std. heats of formation of 148 mols. The limit of reprodn. of the exptl. data in this test set is found to be 1.78 kcal/mol mean abs. error, with a max. of 8.89 kcal/mol error for the ozone mol. This compares rather well with previous results for G2 theory itself (1.58 and 8.2 kcal/mol, resp.). We show that fair stability can be obtained by our optimization procedure.
- 35Wilson, P. J.; Bradley, T. J.; Tozer, D. J. Hybrid Exchange-Correlation Functional Determined from Thermochemical Data and Ab Initio Potentials. J. Chem. Phys. 2001, 115, 9233– 9242, DOI: 10.1063/1.1412605Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXotlyisL0%253D&md5=14c531ba16b99b613ac4213e29e0f8a8Hybrid exchange-correlation functional determined from thermochemical data and ab initio potentialsWilson, Philip J.; Bradley, Thomas J.; Tozer, David J.Journal of Chemical Physics (2001), 115 (20), 9233-9242CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Multiplicative potentials, appropriate for adding to the non-multiplicative fractional orbital exchange term in the Kohn-Sham equations, are detd. from correlated ab initio electron densities. The potentials are examd. graphically and are used in conjunction with conventional thermochem. data to det. a new hybrid exchange-correlation functional, denoted B97-2. Calcns. using B97-2 are compared with those from (a) the B97-1 functional [J. Chem. Phys. 109, 6264 (1998)], which has the same functional form and fraction of orbital exchange, but was fitted to just thermochem. data; and (b) the widely used B3LYP functional [J. Chem. Phys. 98, 5648 (1993)]. B97-2 atomization energies are close to those from B97-1; total electronic energies and ionization potentials are less accurate, but remain an improvement over B3LYP. Mol. structures from all three functionals are comparable. Static isotropic polarizabilities improve from B3LYP to B97-1 to B97-2; the B97-2 functional underestimates exptl. values, which is consistent with the neglect of zero-point vibrational corrections. NMR shielding consts. - detd. as the conventional second deriv. of the electronic energy - improve from B3LYP to B97-1 to B97-2. Shieldings detd. directly from these DFT electron densities using the recently proposed MKS approach [Chem. Phys. Lett. 337, 341 (2001)] are two to three times more accurate than the conventional shieldings, and exhibit an analogous improvement across the three functionals. Classical reaction barriers for sixteen chem. reactions improve significantly from B3LYP to B97-1 to B97-2. The introduction of multiplicative potentials into semi-empirical hybrid functional development therefore appears beneficial.
- 36Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158– 6170, DOI: 10.1063/1.478522Google Scholar36https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXitVCmt7Y%253D&md5=cad4185c69f9232753497f5203d6dc9fToward reliable density functional methods without adjustable parameters: the PBE0 modelAdamo, Carlo; Barone, VincenzoJournal of Chemical Physics (1999), 110 (13), 6158-6170CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present an anal. of the performances of a parameter free d. functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amt. of exact exchange. The results obtained for structural, thermodn., kinetic and spectroscopic (magnetic, IR and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chem. and condensed matter physics.
- 37Zhao, Y.; Schultz, N. E.; Truhlar, D. G. Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput. 2006, 2, 364– 382, DOI: 10.1021/ct0502763Google Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XhtVGqurg%253D&md5=19477e4012ae45e2f0b980eda666c82aDesign of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent InteractionsZhao, Yan; Schultz, Nathan E.; Truhlar, Donald G.Journal of Chemical Theory and Computation (2006), 2 (2), 364-382CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present a new hybrid meta exchange-correlation functional, called M05-2X, for thermochem., thermochem. kinetics, and noncovalent interactions. We also provide a full discussion of the new M05 functional, previously presented in a short communication. The M05 functional was parametrized including both metals and nonmetals, whereas M05-2X is a high-nonlocality functional with double the amt. of nonlocal exchange (2X) that is parametrized only for nonmetals. In particular, M05 was parametrized against 35 data values, and M05-2X is parametrized against 34 data values. Both functionals, along with 28 other functionals, have been comparatively assessed against 234 data values: the MGAE109/3 main-group atomization energy database, the IP13/3 ionization potential database, the EA13/3 electron affinity database, the HTBH38/4 database of barrier height for hydrogen-transfer reactions, five noncovalent databases, two databases involving metal-metal and metal-ligand bond energies, a dipole moment database, a database of four alkyl bond dissocn. energies of alkanes and ethers, and three total energies of one-electron systems. We also tested the new functionals and 12 others for eight hydrogen-bonding and stacking interaction energies in nucleobase pairs, and we tested M05 and M05-2X and 19 other functionals for the geometry, dipole moment, and binding energy of HCN-BF3, which has recently been shown to be a very difficult case for d. functional theory. We tested eight functionals for four more alkyl bond dissocn. energies, and we tested 12 functionals for several addnl. bond energies with varying amts. of multireference character. On the basis of all the results for 256 data values in 18 databases in the present study, we recommend M05-2X, M05, PW6B95, PWB6K, and MPWB1K for general-purpose applications in thermochem., kinetics, and noncovalent interactions involving nonmetals and we recommend M05 for studies involving both metallic and nonmetallic elements. The M05 functional, essentially uniquely among the functionals with broad applicability to chem., also performs well not only for main-group thermochem. and radical reaction barrier heights but also for transition-metal-transition-metal interactions. The M05-2X functional has the best performance for thermochem. kinetics, noncovalent interactions (esp. weak interaction, hydrogen bonding, π···π stacking, and interactions energies of nucleobases), and alkyl bond dissocn. energies and the best composite results for energetics, excluding metals.
- 38Zhao, 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 Function. Theor. Chem. Acc. 2008, 120, 215– 241, DOI: 10.1007/s00214-007-0310-xGoogle Scholar38https://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.
- 39Zhao, Y.; Truhlar, D. G. Comparative DFT Study of van Der Waals Complexes: Rare-Gas Dimers, Alkaline-Earth Dimers, Zinc Dimer, and Zinc-Rare-Gas Dimers. J. Phys. Chem. A 2006, 110, 5121– 5129, DOI: 10.1021/jp060231dGoogle Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XisFWju7g%253D&md5=a57949711cc5b6992f7a2d6f1a6aa37fComparative DFT Study of van der Waals Complexes: Rare-Gas Dimers, Alkaline-Earth Dimers, Zinc Dimer, and Zinc-Rare-Gas DimersZhao, Yan; Truhlar, Donald G.Journal of Physical Chemistry A (2006), 110 (15), 5121-5129CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Recent interest in the application of d. functional theory prompted us to test various functionals for the van der Waals interactions in the rare-gas dimers, the alk.-earth metal dimers, zinc dimer, and zinc-rare-gas dimers. In the present study, we report such tests for 18 DFT functionals, including both some very recent functionals and some well-established older ones. We draw the following conclusions based on the mean errors in binding energies and complex geometries: (1) B97-1 gives the best performance for predicting the geometry of rare-gas dimers, whereas M05-2X and B97-1 give the best energetics for rare-gas dimers; (2) PWB6K gives the best performance for the prediction of the geometry of the alk.-earth metal dimers, zinc dimers, and zinc-rare-gas dimers. M05-2X gives the best energetics for the metal dimers, whereas B97-1 gives the best energetics for the zinc-rare-gas dimers; (3) The M05 functional is unique in providing good accuracy for both covalent transition-metal dimers and van der Waals metal dimers; (4) The combined mean percentage unsigned error in geometries and energetics shows that M05-2X and MPWB1K are the overall best methods for the prediction of van der Waals interactions in metal and rare-gas van der Waals dimers.
- 40Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456– 1465, DOI: 10.1002/jcc.21759Google Scholar40https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXjsF2isL0%253D&md5=370c4fe3164f548718b4bfcf22d1c753Effect of the damping function in dispersion corrected density functional theoryGrimme, Stefan; Ehrlich, Stephan; Goerigk, LarsJournal of Computational Chemistry (2011), 32 (7), 1456-1465CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)It is shown by an extensive benchmark on mol. energy data that the math. form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a std. "zero-damping" formula and rational damping to finite values for small interat. distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coeffs. is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interat. forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramol. dispersion in four representative mol. structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermol. distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of cor. GGAs for non-covalent interactions. According to the thermodn. benchmarks BJ-damping is more accurate esp. for medium-range electron correlation problems and only small and practically insignificant double-counting effects are obsd. It seems to provide a phys. correct short-range behavior of correlation/dispersion even with unmodified std. functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying d. functional. © 2011 Wiley Periodicals, Inc.; J. Comput. Chem., 2011.
- 41Chai, J.-D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom–Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615, DOI: 10.1039/b810189bGoogle Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhtlCksbfO&md5=c7848f8bf050e11972d438aaebd68fdfLong-range corrected hybrid density functionals with damped atom-atom dispersion correctionsChai, Jeng-Da; Head-Gordon, MartinPhysical Chemistry Chemical Physics (2008), 10 (44), 6615-6620CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We report re-optimization of a recently proposed long-range cor. (LC) hybrid d. functional [J.-D. Chai and M. Head-Gordon, J. Chem. Phys., 2008, 128, 084106] to include empirical atom-atom dispersion corrections. The resulting functional, ωB97X-D yields satisfactory accuracy for thermochem., kinetics, and non-covalent interactions. Tests show that for non-covalent systems, ωB97X-D shows slight improvement over other empirical dispersion-cor. d. functionals, while for covalent systems and kinetics it performs noticeably better. Relative to our previous functionals, such as ωB97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions.
- 42Ditchfield, R.; Hehre, W. J.; Pople, J. A. Self-Consistent Molecular-Orbital Methods. IX. An Extended Gaussian-Type Basis for Molecular-Orbital Studies of Organic Molecules. J. Chem. Phys. 1971, 54, 724– 728, DOI: 10.1063/1.1674902Google Scholar42https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE3MXksFOiuw%253D%253D&md5=abce439752b61fad3aa748515ec78c71Self-consistent molecular-orbital methods. IX. Extended Gaussian-type basis for molecular-orbital studies of organic moleculesDitchfield, R.; Hehre, Warren J.; Pople, John A.Journal of Chemical Physics (1971), 54 (2), 724-8CODEN: JCPSA6; ISSN:0021-9606.An extended basis set of at. functions expressed as fixed linear combinations of Gaussian functions is presented for H and the first-row atoms C to F. In this set. described as 4-31 G, each inner shell is represented by a single basis function taken as a sum of 4 Gaussians, and each valence orbital is split into inner and outer parts described by 3 and 1 Gaussian function, resp. The expansion coeffs. and Gaussian exponents are detd. by minimizing the total calcd. energy of the at. ground state. This basis set is then used in single-determinant MO studies of a group of small polyat. mols. Optimization of valence-shell scaling factors shows that considerable rescaling of at. functions occurs in mols., the largest effects being obsd. for H and C. However, the range of optimum scale factors for each atom is small enough to allow the selection of a std. mol. set. The use of this std. basis gives theoretical equil. geometries in reasonable agreement with expt.
- 43Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations. Potentials for the Transition Metal Atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270– 283, DOI: 10.1063/1.448799Google Scholar43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL2MXhtlyju70%253D&md5=29271d2a54b5c81acd19762c570e64d7Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms scandium to mercuryHay, P. Jeffrey; Wadt, Willard R.Journal of Chemical Physics (1985), 82 (1), 270-83CODEN: JCPSA6; ISSN:0021-9606.Ab initio effective core potentials (ECP's) were generated to replace the Coulomb, exchange, and core-orthogonality effects of the chem. inert core electron in the transition metal atoms Sc to Hg. For the second and third transition series relative ECP's were generated which also incorporate the mass-velocity and Darwin relativistic effects into the potential. The ab initio ECP's should facilitate valence electron calcns. on mols. contg. transition-metal atoms with accuracies approaching all-electron calcns. at a fraction of the computational cost. Analytic fits to the potentials are presented for use in multicenter integral evaluation. Gaussian orbital valence basis sets are developed for the (3d,4s,4p), (4d,5s,5p), and (5d,6s,6p) orbitals of the first, second, and third transition series atoms, resp. All-electron and valence-electron at. excitation energies are also compared for the low-lying states of Sc-Hg, and the valence-electron calcns. reproduce the all-electron excitation energies (typically within a few tenths of an eV).
- 44Andrae, D.; Häuβermann, U.; Dolg, M.; Stoll, H.; Preuβ, H. Energy-Adjustedab Initio Pseudopotentials for the Second and Third Row Transition Elements. Theor. Chim. Acta 1990, 77, 123– 141, DOI: 10.1007/BF01114537Google Scholar44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3cXkt12ntLo%253D&md5=8203c6bc6149924cbd4b23b3063715e1Energy-adjusted ab initio pseudopotentials for the second and third row transition elementsAndrae, D.; Haeussermann, U.; Dolg, M.; Stoll, H.; Preuss, H.Theoretica Chimica Acta (1990), 77 (2), 123-41CODEN: TCHAAM; ISSN:0040-5744.Nonrelativistic and quasirelativistic ab initio pseudopotentials substituting the M(z-28)+-core orbitals of the second row transition elements and M(z-60)+-core orbitals of the third row transition elements, resp., and optimized (8s7p6d)/[6s5p3d]-GTO valence basis sets for use in mol. calcns. have been generated. Addnl., corresponding spin-orbit operators have also been derived. At. excitation and ionization energies from numerical HF as well as from SCF pseudopotential calcns. using the derived basis sets differ in most cases by less than 0.1 eV from corresponding numerical all-electron results. Spin-orbit splittings for low-lying states are in reasonable agreement with corresponding all-electron Dirac-Fock (DF) results.
- 45Schwerdtfeger, P. The Pseudopotential Approximation in Electronic Structure Theory. ChemPhysChem 2011, 12, 3143– 3155, DOI: 10.1002/cphc.201100387Google Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhsFKrtLjJ&md5=111797c515719980dddada3d8d319185The Pseudopotential Approximation in Electronic Structure TheorySchwerdtfeger, PeterChemPhysChem (2011), 12 (17), 3143-3155CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. A short review is presented on one of the most successful theories for electronic structure calcns., the pseudopotential approxn., originally introduced by Hans G. A. Hellmann in 1934. Recent developments in relativistic quantum theory allow for the accurate adjustment of pseudopotential parameters to valence spectra, producing results for properties of atoms, mols., and the solid-state in excellent agreement with more accurate all-electron results if a small-core definition is used. Thus the relativistic pseudopotential approxn. is now the most widely applied method for systems contg. heavy elements.
- 46Lynch, B. J.; Zhao, Y.; Truhlar, D. G. Effectiveness of Diffuse Basis Functions for Calculating Relative Energies by Density Functional Theory. J. Phys. Chem. A 2003, 107, 1384– 1388, DOI: 10.1021/jp021590lGoogle Scholar46https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXpvVajtg%253D%253D&md5=ab05489f5339a0733ee43e0fbed3cb24Effectiveness of Diffuse Basis Functions for Calculating Relative Energies by Density Functional TheoryLynch, Benjamin J.; Zhao, Yan; Truhlar, Donald G.Journal of Physical Chemistry A (2003), 107 (9), 1384-1388CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)The addn. of diffuse functions to a double-ζ basis set is more important than increasing to a triple-ζ basis when calcg. reaction energies, reaction barrier heights, and conformational energies with d. functional theory, in particular with the modified Perdew-Wang d. functional. Diffuse basis functions are vital to describe the relative energies between reactants, products, and transition states in isogyric reactions, and they provide enormous improvement in accuracy for conformational equil., using 1, 2-ethanediol and butadiene as examples. As a byproduct of the present study, the authors present a 1-parameter hybrid d. functional method optimized for sugars and sugar-like mols.; this is called MPW1S.
- 47Lanci, M. P.; Roth, J. P. Oxygen Isotope Effects upon Reversible O 2 -Binding Reactions: Characterizing Mononuclear Superoxide and Peroxide Structures. J. Am. Chem. Soc. 2006, 128, 16006– 16007, DOI: 10.1021/ja0669326Google Scholar47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xht1GjsLfI&md5=24588ed1a1a25ecb204314f907b41f2cOxygen Isotope Effects upon Reversible O2-Binding Reactions: Characterizing Mononuclear Superoxide and Peroxide StructuresLanci, Michael P.; Roth, Justine P.Journal of the American Chemical Society (2006), 128 (50), 16006-16007CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Identifying intermediates in catalytic oxidn. reactions requires the development of new probes of structure and mechanism. Reported here are proof-of-concept studies of oxygen (18O) isotope effects upon reversible O2-binding reactions of classic inorg. compds. It is shown that the 18O equil. isotope effects may be used to differentiate structures where O2 is bound as a side-on peroxide ligand vs. an end-on superoxide ligand. The application of 18O equil. isotope effects to the interpretation of 18O kinetic isotope effects and the study of O2 activation mechanisms is also discussed.
- 48Vaska, L.; Chen, L. S.; Senoff, C. V. Oxygen-Carrying Iridium Complexes: Kinetics, Mechanism, and Thermodynamics. Science 1971, 174, 587– 589, DOI: 10.1126/science.174.4009.587Google Scholar48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE3MXlsFGhtLg%253D&md5=ede51c21e64aedcdc03bc26f4d1f1888Oxygen-carrying iridium complexes. Kinetics, mechanism, and thermodynamicsVaska, L.; Chen, Loomis S.; Senoff, C. V.Science (Washington, DC, United States) (1971), 174 (4009), 587-9CODEN: SCIEAS; ISSN:0036-8075.The rates of oxygenation and deoxygenation of a series of Ir complexes increase and decrease, resp., with increasing electron-releasing tendency of the anionic ligands (A) attached to the metal atom in the O-carrying compds. [IrA(CO)(Ph3P)2]. Calcd. heats of oxygenation (-ΔH02), related to Ir-O2 bond energies, are proportional to the previously reported O-O bond lengths in the O adducts, [O2IrA(CO)(Ph3P)2].
- 49Saracini, C.; Liakos, D. G.; Zapata Rivera, J. E.; Neese, F.; Meyer, G. J.; Karlin, K. D. Excitation Wavelength Dependent O 2 Release from Copper(II)–Superoxide Compounds: Laser Flash-Photolysis Experiments and Theoretical Studies. J. Am. Chem. Soc. 2014, 136, 1260– 1263, DOI: 10.1021/ja4115314Google Scholar49https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXotVCisw%253D%253D&md5=55e083b2fdc6b05a0a2d3dd1c42ebba2Excitation Wavelength Dependent O2 Release from Copper(II)-Superoxide Compounds: Laser Flash-Photolysis Experiments and Theoretical StudiesSaracini, Claudio; Liakos, Dimitrios G.; Zapata Rivera, Jhon E.; Neese, Frank; Meyer, Gerald J.; Karlin, Kenneth D.Journal of the American Chemical Society (2014), 136 (4), 1260-1263CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Irradn. of the Cu(II)-superoxide synthetic complexes [(TMG3tren)-CuII(O2)]+ (1) and [(PV-TMPA)-CuII(O2)]+ (2) with visible light resulted in direct photogeneration of O2 gas at low temp. (from -40 to -70° C for 1 and from -125 to -135° C for 2) in 2-methyltetrahydrofuran (MeTHF) solvent (TMG3tren = 1,1,1-tris-{2-[N2-1,1,3,3-tetramethylguanidino]-ethyl}-amine; PV-TMPA = bis(pyrid-2-ylmethyl){[6-(pivalamido)pyrid-2-yl]methyl}amine) . The yield of O2 release was wavelength dependent: λexc = 436 nm, .vphi. = 0.29 (for 1), .vphi. = 0.11 (for 2), and λexc = 683 nm, .vphi. = 0.035 (for 1), .vphi. = 0.078 (for 2), which was followed by fast O2-recombination with [(TMG3tren)-CuI]+ (3) and [(PV-TMPA)-CuI]+ (4). Enthalpic barriers for O2 rebinding to the copper-(I) center (∼10 kJ mol-1) and for O2 dissocn. from the superoxide compd. 1 (45 kJ mol-1) were detd. TD-DFT studies, carried out for 1, support the exptl. results confirming the dissociative character of the excited states formed upon blue- or red-light laser excitation.
- 50Mardirossian, N.; Head-Gordon, M. Thirty Years of Density Functional Theory in Computational Chemistry: An Overview and Extensive Assessment of 200 Density Functionals. Mol. Phys. 2017, 115, 2315– 2372, DOI: 10.1080/00268976.2017.1333644Google Scholar50https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVCltb3O&md5=ba27d707ee3f5fcdd949644d3d2cbd5eThirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionalsMardirossian, Narbe; Head-Gordon, MartinMolecular Physics (2017), 115 (19), 2315-2372CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)In the past 30 years, Kohn-Sham d. functional theory has emerged as the most popular electronic structure method in computational chem. To assess the ever-increasing no. of approx. exchange-correlation functionals, this review benchmarks a total of 200 d. functionals on a mol. database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochem., and barrier heights. In addn., the evolution of non-empirical and semi-empirical d. functional design is reviewed, and guidelines are provided for the proper and effective use of d. functionals. The most promising functional considered is ωB97M-V, a range-sepd. hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chem. applications, and the principal remaining limitations are assocd. with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.
- 51Kosar, N.; Ayub, K.; Mahmood, T. Accurate Theoretical Method for Homolytic Cleavage of C Sn Bond: A Benchmark Approach. Comput. Theor. Chem. 2018, 1140, 134– 144, DOI: 10.1016/j.comptc.2018.08.003Google Scholar51https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsFCku7vM&md5=519271b95bdd51e97f01e8889944f3f9Accurate theoretical method for homolytic cleavage of C-Sn bond: A benchmark approachKosar, Naveen; Ayub, Khurshid; Mahmood, TariqComputational & Theoretical Chemistry (2018), 1140 (), 134-144CODEN: CTCOA5; ISSN:2210-271X. (Elsevier B.V.)Stille coupling is a well-known cross coupling reaction, where the rate detn. step is the dissocn. of carbon stannous (C-Sn) bond. The organotin compds. are also used as precursors in the manufg. of tin oxide films, solar cells, gas sensors, flat panel display technol. and low emission glass materials, etc. The reactivity of organotin compds. has direct relationship with the homolytic cleavage of C-Sn bond. Therefore, accurate detn. of C-Sn bond has direct relevance in understanding many phenomena. The current benchmark study is aimed at finding out the accurate theor. method for the homolytic cleavage (bond dissocn. energy) of C-Sn bond. In this regard, nineteen DFs from eight different classes of DFT with two effective core potential basis sets (LANL2DZ and SDD) and two Karlsruhe basis sets (def2-SVP and def2-TZVP) are selected for the BDE calcn. of C-Sn bond. Ten structurally diverse organotin compds. with exptl. known BDE of C-Sn bond are selected from the literature. The statistical [root mean square deviation (RMSD), std. deviation (SD), Pearson's correlation (R) and mean abs. error (MAE)] results are obtained by the comparison of theor. data with the exptl. BDE values of C-Sn bond of selected organotin compds. Among all DFT classes, GGA-D class is a batter class and BLYP-D3 functional of this class is selected as the best functional for homolytic bond dissocn. energy (BDE) calcn. of C-Sn bond. This functional with SDD basis set shows remarkable performance in reproducing the BDE of C-Sn bond with more accuracy. The SD, RMSD, R and MAE are 4.11 kcal mol-1 is 3.9 kcal mol-1 is 0.963 and -0.01 kcal mol-1, resp.
- 52Kosar, N.; Mahmood, T.; Ayub, K. Role of Dispersion Corrected Hybrid GGA Class in Accurately Calculating the Bond Dissociation Energy of Carbon Halogen Bond: A Benchmark Study. J. Mol. Struct. 2017, 1150, 447– 458, DOI: 10.1016/j.molstruc.2017.08.104Google Scholar52https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhsVyrtbbE&md5=1a6a3793f854353bd3dae05928e18607Role of dispersion corrected hybrid GGA class in accurately calculating the bond dissociation energy of carbon halogen bond: A benchmark studyKosar, Naveen; Mahmood, Tariq; Ayub, KhurshidJournal of Molecular Structure (2017), 1150 (), 447-458CODEN: JMOSB4; ISSN:0022-2860. (Elsevier B.V.)Benchmark study has been carried out to find a cost effective and accurate method for bond dissocn. energy (BDE) of carbon halogen (C-X) bond. BDE of C-X bond plays a vital role in chem. reactions, particularly for kinetic barrier and thermochem. etc. The compds. with C-X bond used for the benchmark study are important reactants in org., inorg. and bioorg. chem. Exptl. data of C-X bond dissocn. energy is compared with theor. results. The statistical anal. tools such as root mean square deviation (RMSD), std. deviation (SD), Pearson's correlation (R) and mean abs. error (MAE) are used for comparison. Overall, thirty-one d. functionals from eight different classes of d. functional theory (DFT) along with Pople and Dunning basis sets are evaluated. Among different classes of DFT, the dispersion cor. range sepd. hybrid GGA class along with 6-31G(d), 6-311G(d), aug-cc-pVDZ and aug-cc-pVTZ basis sets performed best for bond dissocn. energy calcn. of C-X bond. The ωB97XD showed the best performance with less deviations (RMSD, SD), mean abs. error (MAE) and a significant Pearson's correlation (R) when compared to exptl. data. The ωB97XD along with Pople basis set 6-311g(d) has RMSD, SD, R and MAE of 3.14 kcal mol-1, 3.05 kcal mol-1, 0.97 and -1.07 kcal mol-1, resp.
- 53Van Voorhis, T.; Head-Gordon, M. Benchmark Variational Coupled Cluster Doubles Results. J. Chem. Phys. 2000, 113, 8873– 8879, DOI: 10.1063/1.1319643Google Scholar53https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXnvFaitrY%253D&md5=6d6dadea9a39b990c3d2f5109b066b89Benchmark variational coupled cluster doubles resultsVan Voorhis, Troy; Head-Gordon, MartinJournal of Chemical Physics (2000), 113 (20), 8873-8879CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present the first application of the Rayleigh-Ritz variational procedure to the coupled cluster doubles trial function. The variational approach is applied to the potential surface of H4, the double dissocn. of water and the dissocn. of N2, and the results are compared to std. coupled cluster doubles calcns. It is found that the variational approach gives a greatly improved description of strongly correlated systems, where the std. approach is known to fail. Some examn. of the basis set dependence of the results is presented.
- 54Mardirossian, N.; Parkhill, J. A.; Head-Gordon, M. Benchmark Results for Empirical Post-GGA Functionals: Difficult Exchange Problems and Independent Tests. Phys. Chem. Chem. Phys. 2011, 13, 19325, DOI: 10.1039/c1cp21635jGoogle Scholar54https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhtlKju73J&md5=5b6c30e1d90d38ddd4d3af2e0a8263e2Benchmark results for empirical post-GGA functionals: Difficult exchange problems and independent testsMardirossian, Narbe; Parkhill, John A.; Head-Gordon, MartinPhysical Chemistry Chemical Physics (2011), 13 (43), 19325-19337CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)A review. Many of the most promising new d. functionals have improved the treatment of non-local exchange effects with the help of semi-empirical information and more sophisticated recipes for combining Hartree-Fock and local exchange approxns. In order to quantify recent advancements and identify directions for improvement, we have examd. a broad spectrum of test problems. We evaluate the performance of several new hybrid d. functionals (ωB97, ωB97X, ωB97X-D, LRC-ωPBEh, M06, M06-2X, and M06-HF) on a variety of chem. problems, some sensitive to the treatment of exact exchange (which we have hoped to systematically improve) and some which require a balanced treatment of correlation. Since all of the functionals under consideration are parameterized with ground-state thermochem. data, the benchmark aims to det. the applicability of the new d. functionals to cases that have not been considered in the optimization of the semi-empirical parameters. The first class of benchmarks includes the excitation energies of 21 mols. (83 states) primarily from a recent benchmark conducted by Tozer and co-workers, with some addnl. refs. from data made available from the groups of Thiel and Truhlar. We briefly examine the conformational preferences of a small peptide and complete our study with two recently published sets of data that have shown large, systematic errors in simple alkane thermochem. While our results indicate that the more general hybrids currently under development perform well for problems outside of their parameterization and improve over the std. hybrid d. functionals in an essentially systematic way, there is still a significant self-interaction error in the more difficult cases. Functionals based on a range-sepn. of exchange and functionals depending on the kinetic-energy d. both perform comparably, and there is evidence for complementary strengths.
- 55Gomes, J.; Zimmerman, P. M.; Head-Gordon, M.; Bell, A. T. Accurate Prediction of Hydrocarbon Interactions with Zeolites Utilizing Improved Exchange-Correlation Functionals and QM/MM Methods: Benchmark Calculations of Adsorption Enthalpies and Application to Ethene Methylation by Methanol. J. Phys. Chem. C 2012, 116, 15406– 15414, DOI: 10.1021/jp303321sGoogle Scholar55https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XoslGhs7c%253D&md5=edf14af043baa6493a5f71e1e9d9f11dAccurate Prediction of Hydrocarbon Interactions with Zeolites Utilizing Improved Exchange-Correlation Functionals and QM/MM Methods: Benchmark Calculations of Adsorption Enthalpies and Application to Ethene Methylation by MethanolGomes, Joseph; Zimmerman, Paul M.; Head-Gordon, Martin; Bell, Alexis T.Journal of Physical Chemistry C (2012), 116 (29), 15406-15414CODEN: JPCCCK; ISSN:1932-7447. (American Chemical Society)The adsorption enthalpy of light hydrocarbon mols. in both acidic and neutral zeolite MFI has been investigated with a range of computational methods. The role of cluster model size and d. functional theory methodol. is examd. by comparison with high quality ab initio wave function theory MP2 results and exptl. detd. heats of adsorption. The commonly applied B3LYP functional performs poorly in benchmark studies due to an inadequate description of intermol. interactions. The functionals ωB97X-D and M06-2X predict adsorption enthalpies consistent with exptl. values over three classes of adsorbates. A hybrid quantum mechanics/mol. mechanics (QM/MM) method is required to converge calcd. thermochem. properties with respect to cluster model size in a manner that is computationally efficient. The accuracy of both QM and QM/MM methods is highly sensitive to choice of level of theory and cluster size, requiring the use of large basis sets, large cluster models, and a d. functional capable of capturing intermol. interactions for achieving the desired chem. accuracy of 2 kcal/mol with respect to exptl. detd. adsorption enthalpies and activation barriers. The computational effort for performing QM/MM simulations is considerably lower than that of similar quality QM results, and allows for the chem. accurate simulation of chem. reactions occurring in zeolites in a manner that is computationally cost efficient without sacrificing accuracy. The resulting QM/MM procedure is applied to study the reaction of ethene methylation by methanol.
Cited By
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by ACS Publications if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
This article is cited by 13 publications.
- Taghrid S. AlOmar, Abdur Rauf, Umer Rashid, Sehrish Sarfaraz, Khurshid Ayub, Fahad Hussain, Najla Almasoud, Abdulaziz S. AlOmar, Gauhar Rehman, Zubair Ahmad, Naveed Muhammad, Zafar Ali Shah. Molecular docking, DFT studies, and anti-inflammatory evaluation of peshawaraquinone isolated from
Fernandoa adenophylla. Journal of Biomolecular Structure and Dynamics 2024, 42
(20)
, 10604-10616. https://doi.org/10.1080/07391102.2023.2258402
- Aurélien Moncomble, Jean‐Paul Cornard. Orthogonal effect of exchange and correlation parameters in
density functional theory
to compute geometric and spectroscopic quantities. Journal of Computational Chemistry 2024, 45
(24)
, 2071-2079. https://doi.org/10.1002/jcc.27382
- Yu-Lun Wu, Jing Du, Hai-Ying Zhang, Ming-Jun Hou, Qiao-Yue Li, Wei-Chao Chen, Kui-Zhan Shao, Bo Zhu, Chao Qin, Xin-Long Wang, Zhong-Min Su. Dimensional regulation in gigantic molybdenum blue wheels featuring {(W)Mo5} motifs for enhanced proton conductivity. Nano Research 2024, 17
(9)
, 8261-8268. https://doi.org/10.1007/s12274-024-6868-y
- M. Guin, M. Afzal, B. Jana, S. Halder, S. Chatterjee, S. Konar. Co(II) Coordination Compound: Structural and Computational Insights via Crystal Structure, DFT, MEP, NBO and Hirshfeld Surface Analyses. Journal of Structural Chemistry 2024, 65
(9)
, 1875-1885. https://doi.org/10.1134/S0022476624090154
- Tiago da Silva Arouche, Julio Cesar Mendes Lobato, Rosivaldo dos Santos Borges, Mozaniel Santana de Oliveira, Antonio Maia de Jesus Chaves Neto. Molecular interactions of the Omicron, Kappa, and Delta SARS-CoV-2 spike proteins with quantum dots of graphene oxide. Journal of Molecular Modeling 2024, 30
(7)
https://doi.org/10.1007/s00894-024-05996-z
- Batoul Alipour, Mohamad Zaman Kassaee. Comparison of Cu3, Cu5, and Cu7 clusters as potential antioxidants: A theoretical quest. Journal of Molecular Modeling 2024, 30
(5)
https://doi.org/10.1007/s00894-024-05933-0
- Inyang Oyo-Ita, Victory C. Nsofor, Ibtehaj F. Alshdoukhi, Hewa Y. Abdullah, N. Sfina, Bassey B. Asuquo, Terkumbur E. Gber, Adedapo S. Adeyinka, Muyiwa M. Orosun, Hitler Louis. Adsorption, excitation analysis, and sensor properties of heteroatoms (S, P, Si) encapsulated gallium nitride nanotube for hexanol application: A computational approach. Materials Today Communications 2024, 38 , 107679. https://doi.org/10.1016/j.mtcomm.2023.107679
- Tiago da Silva Arouche, Julio Cesar Mendes Lobato, Rosivaldo Borges, Mozaniel Santana Oliveira, Antonio Maia de Jesus Chaves Neto. Molecular Interactions of the Omicron, Kappa, and Delta SARS-CoV-2 Spike Proteins with Quantum Dots of Graphene Oxide. 2024https://doi.org/10.21203/rs.3.rs-3846852/v1
- Ismail Badran, Kotaybah Hashlamoun, Nashaat N. Nassar. Bond dissociation energies of the fifth‐row elements (
InI
): A quantum theoretical benchmark study. International Journal of Quantum Chemistry 2023, 123
(23)
https://doi.org/10.1002/qua.27222
- Ganesh Masilamani, Gamidi Rama Krishna, Sashi Debnath, Anjan Bedi. Origin of Optoelectronic Contradictions in 3,4-Cycloalkyl[c]-chalcogenophenes: A Computational Study. Polymers 2023, 15
(21)
, 4240. https://doi.org/10.3390/polym15214240
- Kiran Shehzadi, Ayesha Tariq, Muhammad Zubair, Tariq Mahmood, Naveen Kosar, Idris Karakaya, Nasir Rasool, Komal Rizwan. Synthesis of pyridine and furan based arylated ketones through palladium catalyst with DFT study of their static and frequency dependent NLO response. Inorganic Chemistry Communications 2023, 151 , 110566. https://doi.org/10.1016/j.inoche.2023.110566
- Muhammad Mohsin, Tehmeena Ishaq, Ijaz Ahmad Bhatti, Maryam Maryam, Asim Jilani, Ammar A. Melaibari, Nidal H. Abu-Hamdeh. Semiconductor Nanomaterial Photocatalysts for Water-Splitting Hydrogen Production: The Holy Grail of Converting Solar Energy to Fuel. Nanomaterials 2023, 13
(3)
, 546. https://doi.org/10.3390/nano13030546
- Muhammad Hanif, Naveen Kosar, Tariq Mahmood, Mian Muhammad, Farhat Ullah, Muhammad Nawaz Tahir, Alany Ingrid Ribeiro, Ezzat Khan. Schiff Bases Derived from 2‐Amino‐6‐methylbenzothiazole, 2‐Amino‐5‐chloropyridine and 4‐Chlorobenzaldehyde: Structure, Computational Studies and Evaluation of Biological Activity. ChemistrySelect 2022, 7
(38)
https://doi.org/10.1002/slct.202203386
- Nelson H. Morgon, Srijit Biswas, Surajit Duari, Aguinaldo R. de Souza. Solvent Effects in the Regioselective N-Functionalization of Tautomerizable Heterocycles Catalyzed by Methyl Trifluoromethanesulfonate: A Density Functional Theory Study with Implicit Solvent Model. Computation 2022, 10
(10)
, 172. https://doi.org/10.3390/computation10100172
Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.
Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.
The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.
Recommended Articles
Abstract
Figure 1
Figure 1. Catalytic water oxidation and oxygen molecule evolution mechanism by Ru complexes.
Figure 2
Figure 2. RMSD of different DFs with two series of basis sets for M–O2 BDEs.
Figure 3
Figure 3. SD of different DFs with two series of basis sets for M–O2 BDEs.
Figure 4
Figure 4. MAE of different DFs with two series of basis sets for M–O2 BDEs.
Figure 5
Figure 5. Pearson’s correlation (R) of B3LYP-GD3Bj with SDD & 6-31+G(d) basis set for BDE calculation of the M–O2 bond.
Figure 6
Figure 6. Structures of transition-metal complexes having oxygen molecule with known experimental BDEs of M–O2 bond. In all the complexes, O2 show side on the binding model with selected transition-metal complexes.
Figure 7
Figure 7. Modal reaction for M–O2 bond dissociation of the rhenium-based complex with oxygen molecule. The rest of the complexes follow the same pattern.
References
This article references 55 other publications.
- 1Shahbaz, M.; Mutascu, M.; Azim, P. Environmental Kuznets Curve in Romania and the Role of Energy Consumption. Renewable Sustainable Energy Rev. 2013, 18, 165– 173, DOI: 10.1016/j.rser.2012.10.012There is no corresponding record for this reference.
- 2Bard, A. J.; Fox, M. A. Artificial Photosynthesis: Solar Splitting of Water to Hydrogen and Oxygen. Acc. Chem. Res. 1995, 28, 141– 145, DOI: 10.1021/ar00051a0072https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXjvFOksrc%253D&md5=62af0ea101f65e4fb2932b80cb29e6ebArtificial Photosynthesis: Solar Splitting of Water to Hydrogen and OxygenBard, Allen J.; Fox, Marye AnneAccounts of Chemical Research (1995), 28 (3), 141-5CODEN: ACHRE4; ISSN:0001-4842. (American Chemical Society)A review with 33 refs. of history and progress of water splitting, semiconductor solid state photovoltaic based systems, semiconductor electrode liq. junction systems, semiconductor particle systems, sensitized semiconductor systems, homogeneous and microheterogeneous systems, and future prospects.
- 3Chow, J.; Kopp, R. J.; Portney, P. R. Energy Resources and Global Development. Science 2003, 302, 1528– 1531, DOI: 10.1126/science.10919393https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD3srntF2gsw%253D%253D&md5=38a0f15c28ab1b0135f4c53c76049f17Energy resources and global developmentChow Jeffrey; Kopp Raymond J; Portney Paul RScience (New York, N.Y.) (2003), 302 (5650), 1528-31 ISSN:.In order to address the economic and environmental consequences of our global energy system, we consider the availability and consumption of energy resources. Problems arise from our dependence on combustible fuels, the environmental risks associated with their extraction, and the environmental damage caused by their emissions. Yet no primary energy source, be it renewable or nonrenewable, is free of environmental or economic limitations. As developed and developing economies continue to grow, conversion to and adoption of environmentally benign energy technology will depend on political and economic realities.
- 4Navarro Yerga, R. M.; Álvarez Galván, M. C.; del Valle, F.; Villoria de la Mano, J. A.; Fierro, J. L. G. Water Splitting on Semiconductor Catalysts under Visible-Light Irradiation. ChemSusChem 2009, 2, 471– 485, DOI: 10.1002/cssc.2009000184https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXnvV2nt7k%253D&md5=0faa86f3527feeb18239f291fd6e74ecWater Splitting on Semiconductor Catalysts under Visible-Light IrradiationNavarro Yerga, Rufino M.; Alvarez Galvan, M. Consuelo; del Valle, F.; Villoria de la Mano, Jose A.; Fierro, Jose L. G.ChemSusChem (2009), 2 (6), 471-485CODEN: CHEMIZ; ISSN:1864-5631. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. The conversion of solar energy into hydrogen via the water-splitting process, assisted by photo-semiconductor catalysts, is one of the most promising technologies for the future because large quantities of hydrogen can potentially be generated in a clean and sustainable manner. This paper provides an overview of the principles, approaches, and research progress on solar hydrogen prodn. via the water-splitting reaction on photo-semiconductor catalysts. It presents a survey of the advances made over the last decades in the development of catalysts for photochem. water splitting under visible-light irradn. The paper also analyzes the energy requirements and main factors that det. the activity of photocatalysts in the conversion of water into hydrogen and oxygen using sunlight. Remarkable progress has been made since the pioneering work by Fujishima and Honda in 1972, but his development of photocatalysts with improved efficiencies for hydrogen prodn. from water using solar energy still faces major challenges. Research strategies and approaches adopted in the search for active and efficient photocatalysts, for example through new materials and synthesis methods, are presented and analyzed.
- 5Barton, E. E.; Rampulla, D. M.; Bocarsly, A. B. Selective Solar-Driven Reduction of CO 2 to Methanol Using a Catalyzed p-GaP Based Photoelectrochemical Cell. J. Am. Chem. Soc. 2008, 130, 6342– 6344, DOI: 10.1021/ja07763275https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXltFens7w%253D&md5=ba5b46e11004ee2431e64b139112d411Selective Solar-Driven Reduction of CO2 to Methanol Using a Catalyzed p-GaP Based Photoelectrochemical CellBarton, Emily E.; Rampulla, David M.; Bocarsly, Andrew B.Journal of the American Chemical Society (2008), 130 (20), 6342-6344CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)With rising atm. CO2 levels there is interest in artificial photosynthetic schemes for converting this greenhouse gas into fuel and small orgs. Photoelectrochem. schemes for activating the inert CO2 mol. operate at high overpotentials and thus do not convert light energy to chem. energy. A selective conversion of CO2 to MeOH at a p-GaP semiconductor electrode with a homogeneous pyridinium ion catalyst, driving the reaction with light energy to yield faradaic efficiencies near 100% at potentials well below the std. potential, is described.
- 6Rashatasakhon, P.; Ozdemir, A. D.; Willis, J.; Padwa, A. Six- versus Five-Membered Ring Formation in Radical Cyclizations of 7-Bromo-Substituted Hexahydroindolinones. Org. Lett. 2004, 6, 917– 920, DOI: 10.1021/ol036347z6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXht1WitL8%253D&md5=6ec21c4ff38a7d5d1aa702455c0f8aa4Six- versus Five-Membered Ring Formation in Radical Cyclizations of 7-Bromo-Substituted HexahydroindolinonesRashatasakhon, Paitoon; Ozdemir, Ayse Daut; Willis, Jerremey; Padwa, AlbertOrganic Letters (2004), 6 (6), 917-920CODEN: ORLEF7; ISSN:1523-7060. (American Chemical Society)Radical cyclization of N-allyl-7-bromo-3a-methyl-hexahydroindol-2-one affords a six-membered ring product that prevails over the isomeric five-membered compd. The former product is generated through two reaction pathways: (a) 6-endo-trig ring closure and (b) rearrangement of an intermediate methylenecyclopentyl radical obtained by 5-exo-trig cyclization.
- 7Service, R. F. Turning Over a New Leaf. Science 2011, 334, 925– 927, DOI: 10.1126/science.334.6058.9257https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhs1eit7bE&md5=547ccba7aeb313c77a9f8262fc5a1ff0Turning over a new leafService, Robert F.Science (Washington, DC, United States) (2011), 334 (6058), 925-927CODEN: SCIEAS; ISSN:0036-8075. (American Association for the Advancement of Science)There is no expanded citation for this reference.
- 8Kurz, P.; Berggren, G.; Anderlund, M. F.; Styring, S. Oxygen Evolving Reactions Catalysed by Synthetic Manganese Complexes: A Systematic Screening. Dalton Trans. 2007, 38, 4258, DOI: 10.1039/b710761gThere is no corresponding record for this reference.
- 9Zong, R.; Thummel, R. P. A New Family of Ru Complexes for Water Oxidation. J. Am. Chem. Soc. 2005, 127, 12802– 12803, DOI: 10.1021/ja054791m9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpt1Cnsb0%253D&md5=2ace4686bbdedbe6af1c7ffb439fce24A New Family of Ru Complexes for Water OxidationZong, Ruifa; Thummel, Randolph P.Journal of the American Chemical Society (2005), 127 (37), 12802-12803CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Reaction of the bridging ligand 3,6-bis-[6'-(1'',8''-naphthyrid-2''-yl)-pyrid-2'-yl]pyridazine (1) with [Ru(DMSO)4Cl2] in aq. ethanol followed by excess 4-substituted pyridine (4-R-py) in the presence of triethylamine provides a series of three well-organized dinuclear complexes characterized by 1H NMR, MS, and X-ray. Mononuclear analogs are prepd. from 4-tert-butyl-2,6-di(1',8'-naphthyrid-2'-yl)pyridine (5) and characterized in a similar fashion. All six complexes show electronic absorption and redox properties consistent with the electron donor/acceptor ability of the axial 4-R-py ligand. When an acetonitrile soln. of the catalyst is added to an aq. Ce(IV)-CF3SO3H soln. (pH = 1.0) at 24 °C, oxygen evolution is obsd. for both mono and dinuclear systems. Turnover nos. range from 50 to 3200 with the best results being found when the axial ligand is 4-methylpyridine (mononuclear TN = 580 and dinuclear TN = 3200).
- 10Polyansky, D. E.; Muckerman, J. T.; Rochford, J.; Zong, R.; Thummel, R. P.; Fujita, E. Water Oxidation by a Mononuclear Ruthenium Catalyst: Characterization of the Intermediates. J. Am. Chem. Soc. 2011, 133, 14649– 14665, DOI: 10.1021/ja203249e10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhtV2isbjP&md5=6a2129651f23f68e39edcdd2c5e58ec0Water Oxidation by a Mononuclear Ruthenium Catalyst: Characterization of the IntermediatesPolyansky, Dmitry E.; Muckerman, James T.; Rochford, Jonathan; Zong, Ruifa; Thummel, Randolph P.; Fujita, EtsukoJournal of the American Chemical Society (2011), 133 (37), 14649-14665CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)A detailed characterization of intermediates in water oxidn. catalyzed by a mononuclear Ru polypyridyl complex [RuII-OH2]2+ (Ru = Ru complex with one 4-t-butyl-2,6-di-(1',8'-naphthyrid-2'-yl)-pyridine ligand and two 4-picoline ligands) has been carried out using electrochem., UV-vis and resonance Raman spectroscopy, pulse radiolysis, stopped flow, and electrospray ionization mass spectrometry (ESI-MS) with H218O labeling expts. and theor. calcns. The results reveal a no. of intriguing properties of intermediates such as [RuIV=O]2+ and [RuIV-OO]2+. At pH > 2.9, two consecutive proton-coupled one-electron steps take place at the potential of the [RuIII-OH]2+/[RuII-OH2]2+ couple, which is equal to or higher than the potential of the [RuIV=O]2+/[RuIII-OH]2+ couple (i.e., the observation of a two-electron oxidn. in cyclic voltammetry). At pH 1, the rate const. of the first one-electron oxidn. by Ce(IV) is k1 = 2 × 104 M-1 s-1. While pH-independent oxidn. of [RuIV=O]2+ takes place at 1420 mV vs NHE, bulk electrolysis of [RuII-OH2]2+ at 1260 mV vs NHE at pH 1 (0.1 M triflic acid) and 1150 mV at pH 6 (10 mM sodium phosphate) yielded a red colored soln. with a Coulomb count corresponding to a net four-electron oxidn. ESI-MS with labeling expts. clearly indicates that this species has an O-O bond. This species required an addnl. oxidn. to liberate an oxygen mol., and without any addnl. oxidant it completely decompd. slowly to form [RuII-OOH]+ over 2 wk. While there remains some conflicting evidence, we have assigned this species as 1[RuIV-η2-OO]2+ based on our electrochem., spectroscopic, and theor. observations alongside a previously reported anal.
- 11Baran, J. D.; Grönbeck, H.; Hellman, A. Analysis of Porphyrines as Catalysts for Electrochemical Reduction of O 2 and Oxidation of H 2 O. J. Am. Chem. Soc. 2014, 136, 1320– 1326, DOI: 10.1021/ja406029911https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXksVag&md5=e249eaac0eeb7649475c68e94850177fAnalysis of Porphyrins as Catalysts for Electrochemical Reduction of O2 and Oxidation of H2OBaran, Jakub D.; Groenbeck, Henrik; Hellman, AndersJournal of the American Chemical Society (2014), 136 (4), 1320-1326CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Bioinspired structures are promising as improved catalysts for various redox reactions. One example is metal hangman-porphyrins (MHP), which recently were suggested for O redn./evolution reaction (ORR/OER). The unique properties of the MHPs are attributed to both the hangman scaffold and the C6F5 side groups. Herein, the OER/ORR over various transition metal MHPs was studied by d. functional theory calcns. within an electrochem. framework. A comparison of the reaction landscape for MHP, metal porphyrins (MP) and metal tetrafluorophenyloporphyrins (MTFPP), allow for a disentanglement of the different roles of the hangman motif and the side groups. In agreement with exptl. studies, Fe and Co are the best MHP metal centers to catalyze these reactions. The addn. of the three-dimensional moiety as hangman scaffold does not break the apparently universal energy relation between *OH and *OOH intermediates. However, the hangman motif is found to stabilize the O intermediate, whereas addn. of C6F5 groups reduces the binding energy of all reaction intermediates. The authors' results indicate that the combination of these two effects allow new design possibilities for macromol. systems with enhanced catalytic OER/ORR activity.
- 12Mahmood, T.; Kosar, N.; Ayub, K. DFT Study of Acceleration of Electrocyclization in Photochromes under Radical Cationic Conditions: Comparison with Recent Experimental Data. Tetrahedron 2017, 73, 3521– 3528, DOI: 10.1016/j.tet.2017.05.03112https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXnvFGhsbc%253D&md5=539b2d64d314ab526c84d46409eea38aDFT study of acceleration of electrocyclization in photochromes under radical cationic conditions: Comparison with recent experimental dataMahmood, Tariq; Kosar, Naveen; Ayub, KhurshidTetrahedron (2017), 73 (25), 3521-3528CODEN: TETRAB; ISSN:0040-4020. (Elsevier Ltd.)Radical cation formation is proposed for the rapid cyclization of 1, 2-bis[5-phenyl-2-methylthien-3-yl]cyclopentene and oligothiophene functionalized dimethyldihydropyrenes (DMDHP). D. functional theory calcns. have been performed to rationalize the effect of a radical cation on the activation barrier of different classes of electrocyclic photochromes (DHP, dithienylethene, dihydroazulene and fulgide). For exact comparative anal., the activation barrier of neutral (singlet) analogs at the same level of theory are also calcd. In addn., the concerted nature and aromaticity of transition states were investigated with the help of synchronicity (Sy.) and nuclear independent chem. shift values NICS(0) calcns., resp., for both the radical cation and neutral systems. In case of the radical cation, thermal return of CPD to DHP, the activation barrier is very low (ΔH = 3.13 kcal mol-1, ΔG = 4.01 kcal mol-1) as compared to the neutral analog (ΔH = 20.6 kcal mol-1, ΔG = 20.98 kcal mol-1), which is consistent with exptl. observations. Similarly for dithenylethenes, radical cation formation has a large impact on the activation barrier (ΔH = 19.44 kcal mol-1, ΔG = 22.29 kcal mol-1). However, radical cation formation has almost negligible impact on the activation barrier of VHF-DHA and fulgide isomerization. The significant difference has been obsd. for synchronicity and NICS(0) values of all types of photochromes under radical cation conditions as compared to the neutral system.
- 13Chiatti, F.; Delle Piane, M.; Ugliengo, P.; Corno, M. Water at Hydroxyapatite Surfaces: The Effect of Coverage and Surface Termination as Investigated by All-Electron B3LYP-D* Simulations. Theor. Chem. Acc. 2016, 135, 54, DOI: 10.1007/s00214-016-1818-8There is no corresponding record for this reference.
- 14Civalleri, B.; Maschio, L.; Ugliengo, P.; Zicovich-Wilson, C. M. Role of Dispersive Interactions in the CO Adsorption on MgO(001): Periodic B3LYP Calculations Augmented with an Empirical Dispersion Term. Phys. Chem. Chem. Phys. 2010, 12, 6382, DOI: 10.1039/c001192d14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXntF2htrY%253D&md5=4e2b5f13bf5d1d4742a54ddda90d47b9Role of dispersive interactions in the CO adsorption on MgO(001): periodic B3LYP calculations augmented with an empirical dispersion termCivalleri, Bartolomeo; Maschio, Lorenzo; Ugliengo, Piero; Zicovich-Wilson, Claudio M.Physical Chemistry Chemical Physics (2010), 12 (24), 6382-6386CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Empirically dispersion cor. B3LYP method (i.e. B3LYP-D) is demonstrated to give excellent results for structure, adsorption energy and vibrational frequency shift for the CO mol. adsorbed on the MgO(001) surface, a system considered a challenge for c.d. functional methods. A periodic approach was adopted to model the interaction using a three-layer slab model. For the B3LYP-D* method an interaction energy of -13.1 kJ mol-1 is computed at low-coverage in very good agreement with exptl. evidence (-12.6 kJ mol-1) as well as a pos. CO vibrational shift of 10 cm-1 to be compared with the exptl. value of 14 cm-1.
- 15Lousada, C. M.; Johansson, A. J.; Brinck, T.; Jonsson, M. Reactivity of Metal Oxide Clusters with Hydrogen Peroxide and Water – a DFT Study Evaluating the Performance of Different Exchange–Correlation Functionals. Phys. Chem. Chem. Phys. 2013, 15, 5539, DOI: 10.1039/c3cp44559c15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXkt1Gnsrw%253D&md5=aae4ff6dd76e2c755a726dd9bbf6c2f0Reactivity of metal oxide clusters with hydrogen peroxide and water - a DFT study evaluating the performance of different exchange-correlation functionalsLousada, Claudio M.; Johansson, Adam Johannes; Brinck, Tore; Jonsson, MatsPhysical Chemistry Chemical Physics (2013), 15 (15), 5539-5552CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We have performed a d. functional theory (DFT) investigation of the interactions of H2O2, H2O and HO radicals with clusters of ZrO2, TiO2 and Y2O3. Different modes of H2O adsorption onto the clusters were studied. In almost all the cases the dissociative adsorption is more exothermic than mol. adsorption. At the surfaces where H2O has undergone dissociative adsorption, the adsorption of H2O2 and the transition state for its decompn. are mediated by hydrogen bonding with the surface HO groups. Using the functionals B3LYP, B3LYP-D and M06 with clusters of 26 and 8 units of ZrO2, the M06 functional performed better than B3LYP in describing the reaction of decompn. of H2O2 and the adsorption of H2O. Addnl., we investigated clusters of the type (ZrO2)2, (TiO2)2 and (Y2O3) and the performance of the functionals B3LYP, B3LYP-D, B3LYP*, M06, M06-L, PBE0, PBE and PWPW91 in describing H2O2, H2O and HO√ adsorption and the energy barrier for decompn. of H2O2. The trends obtained for HO√ adsorption onto the clusters are discussed in terms of the ionization energy of the metal cation present in the oxide. In order to correctly account for the existence of an energy barrier for the decompn. of H2O2, the functional used must include Hartree-Fock exchange. Using minimal cluster models, the best performance in describing the energy barrier for H2O2 decompn. was obtained with the M06 and PBE0 functionals - the av. abs. deviations from expts. are 6 kJ mol-1 and 5 kJ mol-1 resp. With the M06 functional and a larger monoclinic (ZrO2)8 cluster model, the performance is in excellent agreement with exptl. data. For the different oxides, PBE0 was found to be the most effective functional in terms of performance and computational time cost.
- 16Minkin, V. I.; Starikov, A. G.; Starikova, A. A. Computational Insight into Magnetic Behavior and Properties of the Transition Metal Complexes with Redox-Active Ligands: A DFT Approach. Pure Appl. Chem. 2018, 90, 811– 824, DOI: 10.1515/pac-2017-080316https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXptV2mu7g%253D&md5=b4217a1239e16b5e8a56df8ffe7129ceComputational insight into magnetic behavior and properties of the transition metal complexes with redox-active ligands: a DFT approachMinkin, Vladimir I.; Starikov, Andrey G.; Starikova, Alyona A.Pure and Applied Chemistry (2018), 90 (5), 811-824CODEN: PACHAS; ISSN:0033-4545. (Walter de Gruyter, Inc.)Various aspects related to the use of DFT method for the study of magnetic, geometry and energetic properties of transition metal complexes with redox-active ligands are considered. Particular attention is given to the correct choice of model compds. and methodol. of the calcns.
- 17Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104, DOI: 10.1063/1.338234417https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXkvVyks7o%253D&md5=2bca89d904579d5565537a0820dc2ae8A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-PuGrimme, Stefan; Antony, Jens; Ehrlich, Stephan; Krieg, HelgeJournal of Chemical Physics (2010), 132 (15), 154104/1-154104/19CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The method of dispersion correction as an add-on to std. Kohn-Sham d. functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coeffs. and cutoff radii that are both computed from first principles. The coeffs. for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination nos. (CN). They are used to interpolate between dispersion coeffs. of atoms in different chem. environments. The method only requires adjustment of two global parameters for each d. functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of at. forces. Three-body nonadditivity terms are considered. The method has been assessed on std. benchmark sets for inter- and intramol. noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean abs. deviations for the S22 benchmark set of noncovalent interactions for 11 std. d. functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C6 coeffs. also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in mols. and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems. (c) 2010 American Institute of Physics.
- 18Starikov, A. G.; Ivanov, D. G.; Starikova, A. A.; Minkin, V. I. Dispersion Interactions in Oligomerization of Metal Diketonates: A DFT Evaluation. Chem. Pap. 2018, 72, 829– 839, DOI: 10.1007/s11696-017-0225-518https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVShsrjP&md5=dc816f1d25691017e5e3f96a18662106Dispersion interactions in oligomerization of metal diketonates: a DFT evaluationStarikov, Andrey G.; Ivanov, Dmitrii G.; Starikova, Alyona A.; Minkin, Vladimir I.Chemical Papers (2018), 72 (4), 829-839CODEN: CHPAEG; ISSN:1336-9075. (Springer International Publishing AG)The structures of the oligomers of Co(II), Ni(II), Cu(II) and Zn(II) β-diketonates and their stability with respect to dissocn. into the monomeric components were studied with the use of d. functional theory [DFT B3LYP/6-311++G(d,p)] method taking addnl. into account corrections for dispersion interactions in CAM and D3BJ approxns. The calcd. geometries of the oligomers well match those exptl. detd., but stabilization energies obtained at the B3LYP-CAM and B3LYP-D3BJ approxns. are significantly overestimated.
- 19Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157– 167, DOI: 10.1021/ar700111a19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXksV2iug%253D%253D&md5=95aa78f44294d3eed013ff9719cafe76Density Functionals with Broad Applicability in ChemistryZhao, Yan; Truhlar, Donald G.Accounts of Chemical Research (2008), 41 (2), 157-167CODEN: ACHRE4; ISSN:0001-4842. (American Chemical Society)A review. Although d. functional theory is widely used in the computational chem. community, the most popular d. functional, B3LYP, has some serious shortcomings: (i) it is better for main-group chem. than for transition metals; (ii) it systematically underestimates reaction barrier heights; (iii) it is inaccurate for interactions dominated by medium-range correlation energy, such as van der Waals attraction, arom.-arom. stacking, and alkane isomerization energies. We have developed a variety of databases for testing and designing new d. functionals. We used these data to design new d. functionals, called M06-class (and, earlier, M05-class) functionals, for which we enforced some fundamental exact constraints such as the uniform-electron-gas limit and the absence of self-correlation energy. Our M06-class functionals depend on spin-up and spin-down electron densities (i.e., spin densities), spin d. gradients, spin kinetic energy densities, and, for nonlocal (also called hybrid) functionals, Hartree-Fock exchange. We have developed four new functionals that overcome the above-mentioned difficulties: (a) M06, a hybrid meta functional, is a functional with good accuracy "across-the-board" for transition metals, main group thermochem., medium-range correlation energy, and barrier heights; (b) M06-2X, another hybrid meta functional, is not good for transition metals but has excellent performance for main group chem., predicts accurate valence and Rydberg electronic excitation energies, and is an excellent functional for arom.-arom. stacking interactions; (c) M06-L is not as accurate as M06 for barrier heights but is the most accurate functional for transition metals and is the only local functional (no Hartree-Fock exchange) with better across-the-board av. performance than B3LYP; this is very important because only local functionals are affordable for many demanding applications on very large systems; (d) M06-HF has good performance for valence, Rydberg, and charge transfer excited states with minimal sacrifice of ground-state accuracy. In this Account, we compared the performance of the M06-class functionals and one M05-class functional (M05-2X) to that of some popular functionals for diverse databases and their performance on several difficult cases. The tests include barrier heights, conformational energy, and the trend in bond dissocn. energies of Grubbs' ruthenium catalysts for olefin metathesis. Based on these tests, we recommend (1) the M06-2X, BMK, and M05-2X functionals for main-group thermochem. and kinetics, (2) M06-2X and M06 for systems where main-group thermochem., kinetics, and noncovalent interactions are all important, (3) M06-L and M06 for transition metal thermochem., (4) M06 for problems involving multireference rearrangements or reactions where both org. and transition-metal bonds are formed or broken, (5) M06-2X, M05-2X, M06-HF, M06, and M06-L for the study of noncovalent interactions, (6) M06-HF when the use of full Hartree-Fock exchange is important, for example, to avoid the error of self-interaction at long-range, (7) M06-L when a local functional is required, because a local functional has much lower cost for large systems.
- 20Lonsdale, R.; Harvey, J. N.; Mulholland, A. J. Inclusion of Dispersion Effects Significantly Improves Accuracy of Calculated Reaction Barriers for Cytochrome P450 Catalyzed Reactions. J. Phys. Chem. Lett. 2010, 1, 3232– 3237, DOI: 10.1021/jz101279n20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhtleitbvK&md5=9bc8d74cf3976d4da46535ec7d926182Inclusion of Dispersion Effects Significantly Improves Accuracy of Calculated Reaction Barriers for Cytochrome P450 Catalyzed ReactionsLonsdale, Richard; Harvey, Jeremy N.; Mulholland, Adrian J.Journal of Physical Chemistry Letters (2010), 1 (21), 3232-3237CODEN: JPCLCD; ISSN:1948-7185. (American Chemical Society)Prediction of cytochrome P 450 reactivity is of great importance to the development of new medicinal compds. D. functional theory (DFT) has proven itself as a useful tool in the characterization of the elusive reactive species, compd. I, and of the mechanisms of substrate oxidn. B3LYP is the most widely used d. functional in the study of P 450s; however, a major drawback of B3LYP is its inaccurate treatment of dispersion, leading to discrepancies between expt. and theory in some systems. Recent work has shown that an added empirical dispersion correction to B3LYP (B3LYP-D) yields more promising results for similar systems. In the present work, two previously studied systems, camphor hydroxylation and alkene oxidn., have been recalcd. using B3LYP-D. Our work shows that inclusion of dispersion has a significant effect on the energies and geometries of transition states and encounter complexes; furthermore, an improved agreement with exptl. data is obsd.
- 21Stein, T.; Kronik, L.; Baer, R. Reliable Prediction of Charge Transfer Excitations in Molecular Complexes Using Time-Dependent Density Functional Theory. J. Am. Chem. Soc. 2009, 131, 2818– 2820, DOI: 10.1021/ja808748221https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhvVSns7s%253D&md5=b01ac7f2580fdb9495f7e51be1910d57Reliable prediction of charge transfer excitations in molecular complexes using time-dependent density functional theoryStein, Tamar; Kronik, Leeor; Baer, RoiJournal of the American Chemical Society (2009), 131 (8), 2818-2820CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)We show how charge transfer excitations at mol. complexes can be calcd. quant. using time-dependent d. functional theory. Predictive power is obtained from range-sepd. hybrid functionals using nonempirical tuning of the range-splitting parameter. Excellent performance of this approach is obtained for a series of complexes composed of various arom. donors and the tetracyanoethylene acceptor, paving the way to systematic nonempirical quant. studies of charge-transfer excitations in real systems.
- 22Hirao, H. Which DFT Functional Performs Well in the Calculation of Methylcobalamin? Comparison of the B3LYP and BP86 Functionals and Evaluation of the Impact of Empirical Dispersion Correction. J. Phys. Chem. A 2011, 115, 9308– 9313, DOI: 10.1021/jp205280722https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXpsFOgtr4%253D&md5=37163d022b89bd1f37e6a1d8a97029a9Which DFT Functional Performs Well in the Calculation of Methylcobalamin? Comparison of the B3LYP and BP86 Functionals and Evaluation of the Impact of Empirical Dispersion CorrectionHirao, HajimeJournal of Physical Chemistry A (2011), 115 (33), 9308-9313CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Controversy remains regarding the suitable d. functionals for the calcn. of vitamin B12 systems that contain cobalt. To identify the optimum functionals, geometry optimization calcns. were performed on a full-size model of methylcobalamin (MeCbl) using the B3LYP, B3LYP-D, BP86, and BP86-D methods in conjunction with the 6-31G* basis set. Single-point energy evaluations were also performed with the 6-311+G(2d,p) basis set. Consistent with previous studies, the BP86-optimized geometry showed fairly good agreement with the exptl. geometry. Various factors that may influence the homolytic bond dissocn. energy (BDE) of the Co-C bond of MeCbl were systematically evaluated with these methods. Our anal. demonstrated that dispersion was the largest correction term that influenced the magnitude of BDE. Previous studies have shown that B3LYP significantly underestimates BDE, whereas BP86 gives BDE values that are fairly close to the exptl. values (36-37 kcal/mol). The same trend in the relative magnitudes of the BDEs was obsd. in the present calcns. However, BP86 underestimated the BDE for a full model of MeCbl. When the amt. of Hartree-Fock exchange in the B3LYP functional was reduced to 15% and the dispersion correction was made (i.e., B3LYP*-D), the calcd. BDE was in good accord with exptl. values. B3P86-D also performed well. A detailed anal. was undertaken to det. which atoms in cobalamin have large dispersion interactions with a Me fragment of MeCbl.
- 23Minenkov, Y.; Occhipinti, G.; Jensen, V. R. Metal–Phosphine Bond Strengths of the Transition Metals: A Challenge for DFT. J. Phys. Chem. A 2009, 113, 11833– 11844, DOI: 10.1021/jp902940c23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhtV2mtb7L&md5=c0df6248a1cdde1cf246c445f7dbf4a3Metal-phosphine bond strengths of the transition metals: A challenge for DFTMinenkov, Yury; Occhipinti, Giovanni; Jensen, Vidar R.Journal of Physical Chemistry A (2009), 113 (43), 11833-11844CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Previous promising tests of the new M06 family of functionals in predicting ruthenium-metal phosphine bond dissocn. energies have been extended to a series of phosphine complexes of chromium, molybdenum, nickel, and ruthenium for which relevant exptl. data are available. In addn. to the M06 family of functionals, bond dissocn. enthalpies have been calcd. using a selection of d. functionals and hybrid functionals based on the generalized gradient approxn. (GGA), and with or without an empirical term (i.e., DFT-D) accounting for long-range dispersion. For the ruthenium complexes, second-order Moller-Plesset perturbation theory (MP2) has also been applied. Electrostatic and nonelectrostatic solvent effects have been estd. using the polarizable continuum model (PCM), allowing for comparison with exptl. data obtained for dissocn. reactions in org. solvents. Whereas the GGA and hybrid-GGA functionals grossly underestimate the abs. metal-phosphine bond enthalpies, with mean unsigned errors (MUEs) for a set of 10 phosphine dissocn. reactions in the range 13-27 kcal/mol, the recently developed DFT-based methods for inclusion of attractive noncovalent interactions and dispersion (the DFT-D and M06 functionals) dramatically improve upon the situation. The best agreement with expt. is obsd. for BLYP-D (MUE = 2.2 kcal/mol), and with the exception for M06-2X, all these methods provide MUEs well below 5 kcal/mol, which should be sufficient for a broad range of applications. The improvements in predicted relative bond enthalpies are less convincing, however. In several cases the GGA and hybrid-GGA functionals are better at reproducing substitution effects than the DFT-D and M06 methods.
- 24Grimme, S. Accurate Description of van Der Waals Complexes by Density Functional Theory Including Empirical Corrections. J. Comput. Chem. 2004, 25, 1463– 1473, DOI: 10.1002/jcc.2007824https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmtFKgt78%253D&md5=16e36ed7a1c098813d7d947ac72b9078Accurate description of van der Waals complexes by density functional theory including empirical correctionsGrimme, StefanJournal of Computational Chemistry (2004), 25 (12), 1463-1473CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)An empirical method to account for van der Waals interactions in practical calcns. in the framework of the d. functional theory (termed DFT-D) was tested for a wide variety of mol. complexes. As in previous schemes, the dispersive energy was described by damped interat. potentials of the form C6R-6. The use of pure, gradient-cor. d. functionals (BLYP and PBE), together with the resoln.-of-the-identity (RI) approxn. for the Coulomb operator, allows very efficient computations for large systems. In contrast to the previous work, extended AO basis sets of polarized TZV or QZV quality were employed, which reduced the basis set superposition error to a negligible extend. By using a global scaling factor for the at. C6 coeffs., the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes was found to be insignificant if steep damping functions were employed. The method was applied to a total of 29 complexes of atoms and small mols. (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene·H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theor. or exptl. results for binding energies and intermol. distances is obtained. For stacked arom. systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to std. MP2 treatments that systematically over-bind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chem. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calcns. of the interaction energies.
- 25Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787– 1799, DOI: 10.1002/jcc.2049525https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XhtFenu7bO&md5=0b4aa16bebc3a0a2ec175d4b161ab0e4Semiempirical GGA-type density functional constructed with a long-range dispersion correctionGrimme, StefanJournal of Computational Chemistry (2006), 27 (15), 1787-1799CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)A new d. functional (DF) of the generalized gradient approxn. (GGA) type for general chem. applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C6·R-6. A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common d. functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on std. thermochem. benchmark sets, for 40 noncovalently bound complexes, including large stacked arom. mols. and group II element clusters, and for the computation of mol. geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for std. functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean abs. deviation of only 3.8 kcal mol-1. The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the av. CCSD(T) accuracy. The basic strategy in the development to restrict the d. functional description to shorter electron correlation lengths scales and to describe situations with medium to large interat. distances by damped C6·R-6 terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chem. method for large systems where dispersion forces are of general importance.
- 26Johnson, E. R.; Becke, A. D. A Post-Hartree-Fock Model of Intermolecular Interactions: Inclusion of Higher-Order Corrections. J. Chem. Phys. 2006, 124, 174104, DOI: 10.1063/1.219022026https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XksFChurY%253D&md5=99cb745d5a698e2c43e8d1a2bc732d3eA post-Hartree-Fock model of intermolecular interactions: Inclusion of higher-order correctionsJohnson, Erin R.; Becke, Axel D.Journal of Chemical Physics (2006), 124 (17), 174104/1-174104/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We have previously demonstrated that the dipole moment of the exchange hole can be used to derive intermol. C6 dispersion coeffs. [J. Chem. Phys. 122, 154104 (2005)]. This was subsequently the basis for a novel post-Hartree-Fock model of intermol. interactions [J. Chem. Phys. 123, 024101 (2005)]. In the present work, the model is extended to include higher-order dispersion coeffs. C8 and C10. The extended model performs very well for prediction of intermonomer sepns. and binding energies of 45 van der Waals complexes. In particular, it performs twice as well as basis-set extrapolated MP2 theory for dispersion-bound complexes, with minimal computational cost.
- 27Fey, N.; Ridgway, B. M.; Jover, J.; McMullin, C. L.; Harvey, J. N. Organometallic Reactivity: The Role of Metal–Ligand Bond Energies from a Computational Perspective. Dalton Trans. 2011, 40, 11184, DOI: 10.1039/c1dt10909j27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhtlejt7vP&md5=8d85df3f5cd06663f9f5d04b747a962bOrganometallic reactivity: the role of metal-ligand bond energies from a computational perspectiveFey, Natalie; Ridgway, Benjamin M.; Jover, Jesus; McMullin, Claire L.; Harvey, Jeremy N.Dalton Transactions (2011), 40 (42), 11184-11191CODEN: DTARAF; ISSN:1477-9226. (Royal Society of Chemistry)A review; the assocn. and dissocn. of ligands plays a vital role in detg. the reactivity of organometallic catalysts. Computational studies with d. functional theory often fail to reproduce exptl. metal-ligand bond energies, but recently functionals which better capture dispersion effects have been developed. Here we explore their application and discuss future challenges for computational studies of organometallic catalysis.
- 28Frisch, 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.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N., Staroverov, V. N., Normand, R. K. J., Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Rega, M. C. N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Gomperts, J. J. R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J. D. J. F. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford CT, 2009.There is no corresponding record for this reference.
- 29Mo, Y.; Tian, G.; Car, R.; Staroverov, V. N.; Scuseria, G. E.; Tao, J. Performance of a Nonempirical Density Functional on Molecules and Hydrogen-Bonded Complexes. J. Chem. Phys. 2016, 145, 234306, DOI: 10.1063/1.497185329https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XitFamur%252FN&md5=394568c0db0ff1d048c1ecfd073f7401Performance of a nonempirical density functional on molecules and hydrogen-bonded complexesMo, Yuxiang; Tian, Guocai; Car, Roberto; Staroverov, Viktor N.; Scuseria, Gustavo E.; Tao, JianminJournal of Chemical Physics (2016), 145 (23), 234306/1-234306/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recently, Tao and Mo derived a meta-generalized gradient approxn. functional based on a model exchange-correlation hole. In this work, the performance of this functional is assessed on std. test sets, using the 6-311++G(3df,3pd) basis set. These test sets include 223 G3/99 enthalpies of formation, 99 atomization energies, 76 barrier heights, 58 electron affinities, 8 proton affinities, 96 bond lengths, 82 harmonic vibrational frequencies, 10 hydrogen-bonded mol. complexes, and 22 at. excitation energies. Our calcns. show that the Tao-Mo functional can achieve high accuracy for most properties considered, relative to the local spin-d. approxn., Perdew-Burke-Ernzerhof, and Tao-Perdew-Staroverov-Scuseria functionals. In particular, it yields the best accuracy for proton affinities, harmonic vibrational frequencies, hydrogen-bond dissocn. energies and bond lengths, and at. excitation energies. (c) 2016 American Institute of Physics.
- 30Zhao, Y.; Truhlar, D. G. A New Local Density Functional for Main-Group Thermochemistry, Transition Metal Bonding, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2006, 125, 194101, DOI: 10.1063/1.237099330https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xht1CnsrzL&md5=9106de7f4660e3d76691b03dfb530956A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactionsZhao, Yan; Truhlar, Donald G.Journal of Chemical Physics (2006), 125 (19), 194101/1-194101/18CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a new local d. functional, called M06-L, for main-group and transition element thermochem., thermochem. kinetics, and noncovalent interactions. The functional is designed to capture the main dependence of the exchange-correlation energy on local spin d., spin d. gradient, and spin kinetic energy d., and it is parametrized to satisfy the uniform-electron-gas limit and to have good performance for both main-group chem. and transition metal chem. The M06-L functional and 14 other functionals have been comparatively assessed against 22 energetic databases. Among the tested functionals, which include the popular B3LYP, BLYP, and BP86 functionals as well as our previous M05 functional, the M06-L functional gives the best overall performance for a combination of main-group thermochem., thermochem. kinetics, and organometallic, inorganometallic, biol., and noncovalent interactions. It also does very well for predicting geometries and vibrational frequencies. Because of the computational advantages of local functionals, the present functional should be very useful for many applications in chem., esp. for simulations on moderate-sized and large systems and when long time scales must be addressed.
- 31Tao, J.; Perdew, J. P.; Staroverov, V. N.; Scuseria, G. E. Climbing the Density Functional Ladder: Nonempirical Meta–Generalized Gradient Approximation Designed for Molecules and Solids. Phys. Rev. Lett. 2003, 91, 146401, DOI: 10.1103/PhysRevLett.91.14640131https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXnvVKmur4%253D&md5=ee0673dbddc2c14c6bdc2eddf4945d04Climbing the Density Functional Ladder: Nonempirical Meta-Generalized Gradient Approximation Designed for Molecules and SolidsTao, Jianmin; Perdew, John P.; Staroverov, Viktor N.; Scuseria, Gustavo E.Physical Review Letters (2003), 91 (14), 146401/1-146401/4CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)The electron d., its gradient, and the Kohn-Sham orbital kinetic energy d. are the local ingredients of a meta-generalized gradient approxn. (meta-GGA). We construct a meta-GGA d. functional for the exchange-correlation energy that satisfies exact constraints without empirical parameters. The exchange and correlation terms respect two paradigms: one- or two-electron densities and slowly varying densities, and so describe both mols. and solids with high accuracy, as shown by extensive numerical tests. This functional completes the third rung of "Jacob's ladder" of approxns., above the local spin d. and GGA rungs.
- 32Zhao, Y.; Lynch, B. J.; Truhlar, D. G. Development and Assessment of a New Hybrid Density Functional Model for Thermochemical Kinetics. J. Phys. Chem. A 2004, 108, 2715– 2719, DOI: 10.1021/jp049908s32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXhvF2nu7o%253D&md5=204a63ab21ed0a0ca8c89bee27ed90cfDevelopment and Assessment of a New Hybrid Density Functional Model for Thermochemical KineticsZhao, Yan; Lynch, Benjamin J.; Truhlar, Donald G.Journal of Physical Chemistry A (2004), 108 (14), 2715-2719CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)A new hybrid Hartree-Fock-d. functional model called the Becke88-Becke95 1-parameter model for kinetics (BB1K) was optimized against a database of three forward barrier heights, three reverse barrier heights, and three energies of reaction for the reactions in the BH6 representative barrier height database. We then assessed the newly developed BB1K method against a saddle point geometries database, a database of 42 barrier heights, the AE6 representative atomization energy database, a mol. geometries data set, and a set of 13 zero point energies. The results show that BB1K can give excellent saddle point geometries and barrier heights, and its performance for calcg. atomization energies is 40% better than MPW1K. Using a mean mean unsigned error criterion that equally wts. the errors in barrier heights and in bond energies, the new BB1K method outperforms all other DFT and hybrid DFT methods by a large margin, and we therefore conclude that it is the best d. functional-type method for thermochem. kinetics.
- 33Becke, A. D. New Mixing of Hartree–Fock and Local Density-functional Theories. J. Chem. Phys. 1993, 98, 1372– 1377, DOI: 10.1063/1.46430433https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3sXhtlagt7o%253D&md5=02a9289532197e7445af684bb29f4e11A new mixing of Hartree-Fock and local-density-functional theoriesBecke, Axel D.Journal of Chemical Physics (1993), 98 (2), 1372-7CODEN: JCPSA6; ISSN:0021-9606.Previous attempts to combine the Hartree-Fock theory with the local d.-functional theory were unsuccessful in applications to mol. bonding. A new coupling is presented of these two theories that maintains their simplicity and computational efficiency, and yet greatly improves their predictive power. Very encouraging results of tests on atomization energies, ionization potentials, and proton affinities are reported, and the potential for future development is discussed.
- 34Schmider, H. L.; Becke, A. D. Optimized Density Functionals from the Extended G2 Test Set. J. Chem. Phys. 1998, 108, 9624– 9631, DOI: 10.1063/1.47643834https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXjsFCrtLs%253D&md5=e57a83f02bc6ee46e6b1275e308ed7e1Optimized density functionals from the extended G2 test setSchmider, Hartmut L.; Becke, Axel D.Journal of Chemical Physics (1998), 108 (23), 9624-9631CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A recently suggested procedure for the systematic optimization of gradient-cor. exchange-correlation functionals [A. D. Becke, J. Chem. Phys. 107, 8554 (1997)] has been applied to the extended G2 test set [L. A. Curtiss et al., J. Chem. Phys. 106, 1063 (1997)], which consists of the std. heats of formation of 148 mols. The limit of reprodn. of the exptl. data in this test set is found to be 1.78 kcal/mol mean abs. error, with a max. of 8.89 kcal/mol error for the ozone mol. This compares rather well with previous results for G2 theory itself (1.58 and 8.2 kcal/mol, resp.). We show that fair stability can be obtained by our optimization procedure.
- 35Wilson, P. J.; Bradley, T. J.; Tozer, D. J. Hybrid Exchange-Correlation Functional Determined from Thermochemical Data and Ab Initio Potentials. J. Chem. Phys. 2001, 115, 9233– 9242, DOI: 10.1063/1.141260535https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXotlyisL0%253D&md5=14c531ba16b99b613ac4213e29e0f8a8Hybrid exchange-correlation functional determined from thermochemical data and ab initio potentialsWilson, Philip J.; Bradley, Thomas J.; Tozer, David J.Journal of Chemical Physics (2001), 115 (20), 9233-9242CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Multiplicative potentials, appropriate for adding to the non-multiplicative fractional orbital exchange term in the Kohn-Sham equations, are detd. from correlated ab initio electron densities. The potentials are examd. graphically and are used in conjunction with conventional thermochem. data to det. a new hybrid exchange-correlation functional, denoted B97-2. Calcns. using B97-2 are compared with those from (a) the B97-1 functional [J. Chem. Phys. 109, 6264 (1998)], which has the same functional form and fraction of orbital exchange, but was fitted to just thermochem. data; and (b) the widely used B3LYP functional [J. Chem. Phys. 98, 5648 (1993)]. B97-2 atomization energies are close to those from B97-1; total electronic energies and ionization potentials are less accurate, but remain an improvement over B3LYP. Mol. structures from all three functionals are comparable. Static isotropic polarizabilities improve from B3LYP to B97-1 to B97-2; the B97-2 functional underestimates exptl. values, which is consistent with the neglect of zero-point vibrational corrections. NMR shielding consts. - detd. as the conventional second deriv. of the electronic energy - improve from B3LYP to B97-1 to B97-2. Shieldings detd. directly from these DFT electron densities using the recently proposed MKS approach [Chem. Phys. Lett. 337, 341 (2001)] are two to three times more accurate than the conventional shieldings, and exhibit an analogous improvement across the three functionals. Classical reaction barriers for sixteen chem. reactions improve significantly from B3LYP to B97-1 to B97-2. The introduction of multiplicative potentials into semi-empirical hybrid functional development therefore appears beneficial.
- 36Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158– 6170, DOI: 10.1063/1.47852236https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXitVCmt7Y%253D&md5=cad4185c69f9232753497f5203d6dc9fToward reliable density functional methods without adjustable parameters: the PBE0 modelAdamo, Carlo; Barone, VincenzoJournal of Chemical Physics (1999), 110 (13), 6158-6170CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present an anal. of the performances of a parameter free d. functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amt. of exact exchange. The results obtained for structural, thermodn., kinetic and spectroscopic (magnetic, IR and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chem. and condensed matter physics.
- 37Zhao, Y.; Schultz, N. E.; Truhlar, D. G. Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput. 2006, 2, 364– 382, DOI: 10.1021/ct050276337https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XhtVGqurg%253D&md5=19477e4012ae45e2f0b980eda666c82aDesign of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent InteractionsZhao, Yan; Schultz, Nathan E.; Truhlar, Donald G.Journal of Chemical Theory and Computation (2006), 2 (2), 364-382CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present a new hybrid meta exchange-correlation functional, called M05-2X, for thermochem., thermochem. kinetics, and noncovalent interactions. We also provide a full discussion of the new M05 functional, previously presented in a short communication. The M05 functional was parametrized including both metals and nonmetals, whereas M05-2X is a high-nonlocality functional with double the amt. of nonlocal exchange (2X) that is parametrized only for nonmetals. In particular, M05 was parametrized against 35 data values, and M05-2X is parametrized against 34 data values. Both functionals, along with 28 other functionals, have been comparatively assessed against 234 data values: the MGAE109/3 main-group atomization energy database, the IP13/3 ionization potential database, the EA13/3 electron affinity database, the HTBH38/4 database of barrier height for hydrogen-transfer reactions, five noncovalent databases, two databases involving metal-metal and metal-ligand bond energies, a dipole moment database, a database of four alkyl bond dissocn. energies of alkanes and ethers, and three total energies of one-electron systems. We also tested the new functionals and 12 others for eight hydrogen-bonding and stacking interaction energies in nucleobase pairs, and we tested M05 and M05-2X and 19 other functionals for the geometry, dipole moment, and binding energy of HCN-BF3, which has recently been shown to be a very difficult case for d. functional theory. We tested eight functionals for four more alkyl bond dissocn. energies, and we tested 12 functionals for several addnl. bond energies with varying amts. of multireference character. On the basis of all the results for 256 data values in 18 databases in the present study, we recommend M05-2X, M05, PW6B95, PWB6K, and MPWB1K for general-purpose applications in thermochem., kinetics, and noncovalent interactions involving nonmetals and we recommend M05 for studies involving both metallic and nonmetallic elements. The M05 functional, essentially uniquely among the functionals with broad applicability to chem., also performs well not only for main-group thermochem. and radical reaction barrier heights but also for transition-metal-transition-metal interactions. The M05-2X functional has the best performance for thermochem. kinetics, noncovalent interactions (esp. weak interaction, hydrogen bonding, π···π stacking, and interactions energies of nucleobases), and alkyl bond dissocn. energies and the best composite results for energetics, excluding metals.
- 38Zhao, 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 Function. Theor. Chem. Acc. 2008, 120, 215– 241, DOI: 10.1007/s00214-007-0310-x38https://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.
- 39Zhao, Y.; Truhlar, D. G. Comparative DFT Study of van Der Waals Complexes: Rare-Gas Dimers, Alkaline-Earth Dimers, Zinc Dimer, and Zinc-Rare-Gas Dimers. J. Phys. Chem. A 2006, 110, 5121– 5129, DOI: 10.1021/jp060231d39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XisFWju7g%253D&md5=a57949711cc5b6992f7a2d6f1a6aa37fComparative DFT Study of van der Waals Complexes: Rare-Gas Dimers, Alkaline-Earth Dimers, Zinc Dimer, and Zinc-Rare-Gas DimersZhao, Yan; Truhlar, Donald G.Journal of Physical Chemistry A (2006), 110 (15), 5121-5129CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Recent interest in the application of d. functional theory prompted us to test various functionals for the van der Waals interactions in the rare-gas dimers, the alk.-earth metal dimers, zinc dimer, and zinc-rare-gas dimers. In the present study, we report such tests for 18 DFT functionals, including both some very recent functionals and some well-established older ones. We draw the following conclusions based on the mean errors in binding energies and complex geometries: (1) B97-1 gives the best performance for predicting the geometry of rare-gas dimers, whereas M05-2X and B97-1 give the best energetics for rare-gas dimers; (2) PWB6K gives the best performance for the prediction of the geometry of the alk.-earth metal dimers, zinc dimers, and zinc-rare-gas dimers. M05-2X gives the best energetics for the metal dimers, whereas B97-1 gives the best energetics for the zinc-rare-gas dimers; (3) The M05 functional is unique in providing good accuracy for both covalent transition-metal dimers and van der Waals metal dimers; (4) The combined mean percentage unsigned error in geometries and energetics shows that M05-2X and MPWB1K are the overall best methods for the prediction of van der Waals interactions in metal and rare-gas van der Waals dimers.
- 40Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456– 1465, DOI: 10.1002/jcc.2175940https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXjsF2isL0%253D&md5=370c4fe3164f548718b4bfcf22d1c753Effect of the damping function in dispersion corrected density functional theoryGrimme, Stefan; Ehrlich, Stephan; Goerigk, LarsJournal of Computational Chemistry (2011), 32 (7), 1456-1465CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)It is shown by an extensive benchmark on mol. energy data that the math. form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a std. "zero-damping" formula and rational damping to finite values for small interat. distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coeffs. is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interat. forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramol. dispersion in four representative mol. structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermol. distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of cor. GGAs for non-covalent interactions. According to the thermodn. benchmarks BJ-damping is more accurate esp. for medium-range electron correlation problems and only small and practically insignificant double-counting effects are obsd. It seems to provide a phys. correct short-range behavior of correlation/dispersion even with unmodified std. functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying d. functional. © 2011 Wiley Periodicals, Inc.; J. Comput. Chem., 2011.
- 41Chai, J.-D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom–Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615, DOI: 10.1039/b810189b41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhtlCksbfO&md5=c7848f8bf050e11972d438aaebd68fdfLong-range corrected hybrid density functionals with damped atom-atom dispersion correctionsChai, Jeng-Da; Head-Gordon, MartinPhysical Chemistry Chemical Physics (2008), 10 (44), 6615-6620CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We report re-optimization of a recently proposed long-range cor. (LC) hybrid d. functional [J.-D. Chai and M. Head-Gordon, J. Chem. Phys., 2008, 128, 084106] to include empirical atom-atom dispersion corrections. The resulting functional, ωB97X-D yields satisfactory accuracy for thermochem., kinetics, and non-covalent interactions. Tests show that for non-covalent systems, ωB97X-D shows slight improvement over other empirical dispersion-cor. d. functionals, while for covalent systems and kinetics it performs noticeably better. Relative to our previous functionals, such as ωB97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions.
- 42Ditchfield, R.; Hehre, W. J.; Pople, J. A. Self-Consistent Molecular-Orbital Methods. IX. An Extended Gaussian-Type Basis for Molecular-Orbital Studies of Organic Molecules. J. Chem. Phys. 1971, 54, 724– 728, DOI: 10.1063/1.167490242https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE3MXksFOiuw%253D%253D&md5=abce439752b61fad3aa748515ec78c71Self-consistent molecular-orbital methods. IX. Extended Gaussian-type basis for molecular-orbital studies of organic moleculesDitchfield, R.; Hehre, Warren J.; Pople, John A.Journal of Chemical Physics (1971), 54 (2), 724-8CODEN: JCPSA6; ISSN:0021-9606.An extended basis set of at. functions expressed as fixed linear combinations of Gaussian functions is presented for H and the first-row atoms C to F. In this set. described as 4-31 G, each inner shell is represented by a single basis function taken as a sum of 4 Gaussians, and each valence orbital is split into inner and outer parts described by 3 and 1 Gaussian function, resp. The expansion coeffs. and Gaussian exponents are detd. by minimizing the total calcd. energy of the at. ground state. This basis set is then used in single-determinant MO studies of a group of small polyat. mols. Optimization of valence-shell scaling factors shows that considerable rescaling of at. functions occurs in mols., the largest effects being obsd. for H and C. However, the range of optimum scale factors for each atom is small enough to allow the selection of a std. mol. set. The use of this std. basis gives theoretical equil. geometries in reasonable agreement with expt.
- 43Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations. Potentials for the Transition Metal Atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270– 283, DOI: 10.1063/1.44879943https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL2MXhtlyju70%253D&md5=29271d2a54b5c81acd19762c570e64d7Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms scandium to mercuryHay, P. Jeffrey; Wadt, Willard R.Journal of Chemical Physics (1985), 82 (1), 270-83CODEN: JCPSA6; ISSN:0021-9606.Ab initio effective core potentials (ECP's) were generated to replace the Coulomb, exchange, and core-orthogonality effects of the chem. inert core electron in the transition metal atoms Sc to Hg. For the second and third transition series relative ECP's were generated which also incorporate the mass-velocity and Darwin relativistic effects into the potential. The ab initio ECP's should facilitate valence electron calcns. on mols. contg. transition-metal atoms with accuracies approaching all-electron calcns. at a fraction of the computational cost. Analytic fits to the potentials are presented for use in multicenter integral evaluation. Gaussian orbital valence basis sets are developed for the (3d,4s,4p), (4d,5s,5p), and (5d,6s,6p) orbitals of the first, second, and third transition series atoms, resp. All-electron and valence-electron at. excitation energies are also compared for the low-lying states of Sc-Hg, and the valence-electron calcns. reproduce the all-electron excitation energies (typically within a few tenths of an eV).
- 44Andrae, D.; Häuβermann, U.; Dolg, M.; Stoll, H.; Preuβ, H. Energy-Adjustedab Initio Pseudopotentials for the Second and Third Row Transition Elements. Theor. Chim. Acta 1990, 77, 123– 141, DOI: 10.1007/BF0111453744https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3cXkt12ntLo%253D&md5=8203c6bc6149924cbd4b23b3063715e1Energy-adjusted ab initio pseudopotentials for the second and third row transition elementsAndrae, D.; Haeussermann, U.; Dolg, M.; Stoll, H.; Preuss, H.Theoretica Chimica Acta (1990), 77 (2), 123-41CODEN: TCHAAM; ISSN:0040-5744.Nonrelativistic and quasirelativistic ab initio pseudopotentials substituting the M(z-28)+-core orbitals of the second row transition elements and M(z-60)+-core orbitals of the third row transition elements, resp., and optimized (8s7p6d)/[6s5p3d]-GTO valence basis sets for use in mol. calcns. have been generated. Addnl., corresponding spin-orbit operators have also been derived. At. excitation and ionization energies from numerical HF as well as from SCF pseudopotential calcns. using the derived basis sets differ in most cases by less than 0.1 eV from corresponding numerical all-electron results. Spin-orbit splittings for low-lying states are in reasonable agreement with corresponding all-electron Dirac-Fock (DF) results.
- 45Schwerdtfeger, P. The Pseudopotential Approximation in Electronic Structure Theory. ChemPhysChem 2011, 12, 3143– 3155, DOI: 10.1002/cphc.20110038745https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhsFKrtLjJ&md5=111797c515719980dddada3d8d319185The Pseudopotential Approximation in Electronic Structure TheorySchwerdtfeger, PeterChemPhysChem (2011), 12 (17), 3143-3155CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. A short review is presented on one of the most successful theories for electronic structure calcns., the pseudopotential approxn., originally introduced by Hans G. A. Hellmann in 1934. Recent developments in relativistic quantum theory allow for the accurate adjustment of pseudopotential parameters to valence spectra, producing results for properties of atoms, mols., and the solid-state in excellent agreement with more accurate all-electron results if a small-core definition is used. Thus the relativistic pseudopotential approxn. is now the most widely applied method for systems contg. heavy elements.
- 46Lynch, B. J.; Zhao, Y.; Truhlar, D. G. Effectiveness of Diffuse Basis Functions for Calculating Relative Energies by Density Functional Theory. J. Phys. Chem. A 2003, 107, 1384– 1388, DOI: 10.1021/jp021590l46https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXpvVajtg%253D%253D&md5=ab05489f5339a0733ee43e0fbed3cb24Effectiveness of Diffuse Basis Functions for Calculating Relative Energies by Density Functional TheoryLynch, Benjamin J.; Zhao, Yan; Truhlar, Donald G.Journal of Physical Chemistry A (2003), 107 (9), 1384-1388CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)The addn. of diffuse functions to a double-ζ basis set is more important than increasing to a triple-ζ basis when calcg. reaction energies, reaction barrier heights, and conformational energies with d. functional theory, in particular with the modified Perdew-Wang d. functional. Diffuse basis functions are vital to describe the relative energies between reactants, products, and transition states in isogyric reactions, and they provide enormous improvement in accuracy for conformational equil., using 1, 2-ethanediol and butadiene as examples. As a byproduct of the present study, the authors present a 1-parameter hybrid d. functional method optimized for sugars and sugar-like mols.; this is called MPW1S.
- 47Lanci, M. P.; Roth, J. P. Oxygen Isotope Effects upon Reversible O 2 -Binding Reactions: Characterizing Mononuclear Superoxide and Peroxide Structures. J. Am. Chem. Soc. 2006, 128, 16006– 16007, DOI: 10.1021/ja066932647https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xht1GjsLfI&md5=24588ed1a1a25ecb204314f907b41f2cOxygen Isotope Effects upon Reversible O2-Binding Reactions: Characterizing Mononuclear Superoxide and Peroxide StructuresLanci, Michael P.; Roth, Justine P.Journal of the American Chemical Society (2006), 128 (50), 16006-16007CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Identifying intermediates in catalytic oxidn. reactions requires the development of new probes of structure and mechanism. Reported here are proof-of-concept studies of oxygen (18O) isotope effects upon reversible O2-binding reactions of classic inorg. compds. It is shown that the 18O equil. isotope effects may be used to differentiate structures where O2 is bound as a side-on peroxide ligand vs. an end-on superoxide ligand. The application of 18O equil. isotope effects to the interpretation of 18O kinetic isotope effects and the study of O2 activation mechanisms is also discussed.
- 48Vaska, L.; Chen, L. S.; Senoff, C. V. Oxygen-Carrying Iridium Complexes: Kinetics, Mechanism, and Thermodynamics. Science 1971, 174, 587– 589, DOI: 10.1126/science.174.4009.58748https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE3MXlsFGhtLg%253D&md5=ede51c21e64aedcdc03bc26f4d1f1888Oxygen-carrying iridium complexes. Kinetics, mechanism, and thermodynamicsVaska, L.; Chen, Loomis S.; Senoff, C. V.Science (Washington, DC, United States) (1971), 174 (4009), 587-9CODEN: SCIEAS; ISSN:0036-8075.The rates of oxygenation and deoxygenation of a series of Ir complexes increase and decrease, resp., with increasing electron-releasing tendency of the anionic ligands (A) attached to the metal atom in the O-carrying compds. [IrA(CO)(Ph3P)2]. Calcd. heats of oxygenation (-ΔH02), related to Ir-O2 bond energies, are proportional to the previously reported O-O bond lengths in the O adducts, [O2IrA(CO)(Ph3P)2].
- 49Saracini, C.; Liakos, D. G.; Zapata Rivera, J. E.; Neese, F.; Meyer, G. J.; Karlin, K. D. Excitation Wavelength Dependent O 2 Release from Copper(II)–Superoxide Compounds: Laser Flash-Photolysis Experiments and Theoretical Studies. J. Am. Chem. Soc. 2014, 136, 1260– 1263, DOI: 10.1021/ja411531449https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXotVCisw%253D%253D&md5=55e083b2fdc6b05a0a2d3dd1c42ebba2Excitation Wavelength Dependent O2 Release from Copper(II)-Superoxide Compounds: Laser Flash-Photolysis Experiments and Theoretical StudiesSaracini, Claudio; Liakos, Dimitrios G.; Zapata Rivera, Jhon E.; Neese, Frank; Meyer, Gerald J.; Karlin, Kenneth D.Journal of the American Chemical Society (2014), 136 (4), 1260-1263CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Irradn. of the Cu(II)-superoxide synthetic complexes [(TMG3tren)-CuII(O2)]+ (1) and [(PV-TMPA)-CuII(O2)]+ (2) with visible light resulted in direct photogeneration of O2 gas at low temp. (from -40 to -70° C for 1 and from -125 to -135° C for 2) in 2-methyltetrahydrofuran (MeTHF) solvent (TMG3tren = 1,1,1-tris-{2-[N2-1,1,3,3-tetramethylguanidino]-ethyl}-amine; PV-TMPA = bis(pyrid-2-ylmethyl){[6-(pivalamido)pyrid-2-yl]methyl}amine) . The yield of O2 release was wavelength dependent: λexc = 436 nm, .vphi. = 0.29 (for 1), .vphi. = 0.11 (for 2), and λexc = 683 nm, .vphi. = 0.035 (for 1), .vphi. = 0.078 (for 2), which was followed by fast O2-recombination with [(TMG3tren)-CuI]+ (3) and [(PV-TMPA)-CuI]+ (4). Enthalpic barriers for O2 rebinding to the copper-(I) center (∼10 kJ mol-1) and for O2 dissocn. from the superoxide compd. 1 (45 kJ mol-1) were detd. TD-DFT studies, carried out for 1, support the exptl. results confirming the dissociative character of the excited states formed upon blue- or red-light laser excitation.
- 50Mardirossian, N.; Head-Gordon, M. Thirty Years of Density Functional Theory in Computational Chemistry: An Overview and Extensive Assessment of 200 Density Functionals. Mol. Phys. 2017, 115, 2315– 2372, DOI: 10.1080/00268976.2017.133364450https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVCltb3O&md5=ba27d707ee3f5fcdd949644d3d2cbd5eThirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionalsMardirossian, Narbe; Head-Gordon, MartinMolecular Physics (2017), 115 (19), 2315-2372CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)In the past 30 years, Kohn-Sham d. functional theory has emerged as the most popular electronic structure method in computational chem. To assess the ever-increasing no. of approx. exchange-correlation functionals, this review benchmarks a total of 200 d. functionals on a mol. database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochem., and barrier heights. In addn., the evolution of non-empirical and semi-empirical d. functional design is reviewed, and guidelines are provided for the proper and effective use of d. functionals. The most promising functional considered is ωB97M-V, a range-sepd. hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chem. applications, and the principal remaining limitations are assocd. with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.
- 51Kosar, N.; Ayub, K.; Mahmood, T. Accurate Theoretical Method for Homolytic Cleavage of C Sn Bond: A Benchmark Approach. Comput. Theor. Chem. 2018, 1140, 134– 144, DOI: 10.1016/j.comptc.2018.08.00351https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsFCku7vM&md5=519271b95bdd51e97f01e8889944f3f9Accurate theoretical method for homolytic cleavage of C-Sn bond: A benchmark approachKosar, Naveen; Ayub, Khurshid; Mahmood, TariqComputational & Theoretical Chemistry (2018), 1140 (), 134-144CODEN: CTCOA5; ISSN:2210-271X. (Elsevier B.V.)Stille coupling is a well-known cross coupling reaction, where the rate detn. step is the dissocn. of carbon stannous (C-Sn) bond. The organotin compds. are also used as precursors in the manufg. of tin oxide films, solar cells, gas sensors, flat panel display technol. and low emission glass materials, etc. The reactivity of organotin compds. has direct relationship with the homolytic cleavage of C-Sn bond. Therefore, accurate detn. of C-Sn bond has direct relevance in understanding many phenomena. The current benchmark study is aimed at finding out the accurate theor. method for the homolytic cleavage (bond dissocn. energy) of C-Sn bond. In this regard, nineteen DFs from eight different classes of DFT with two effective core potential basis sets (LANL2DZ and SDD) and two Karlsruhe basis sets (def2-SVP and def2-TZVP) are selected for the BDE calcn. of C-Sn bond. Ten structurally diverse organotin compds. with exptl. known BDE of C-Sn bond are selected from the literature. The statistical [root mean square deviation (RMSD), std. deviation (SD), Pearson's correlation (R) and mean abs. error (MAE)] results are obtained by the comparison of theor. data with the exptl. BDE values of C-Sn bond of selected organotin compds. Among all DFT classes, GGA-D class is a batter class and BLYP-D3 functional of this class is selected as the best functional for homolytic bond dissocn. energy (BDE) calcn. of C-Sn bond. This functional with SDD basis set shows remarkable performance in reproducing the BDE of C-Sn bond with more accuracy. The SD, RMSD, R and MAE are 4.11 kcal mol-1 is 3.9 kcal mol-1 is 0.963 and -0.01 kcal mol-1, resp.
- 52Kosar, N.; Mahmood, T.; Ayub, K. Role of Dispersion Corrected Hybrid GGA Class in Accurately Calculating the Bond Dissociation Energy of Carbon Halogen Bond: A Benchmark Study. J. Mol. Struct. 2017, 1150, 447– 458, DOI: 10.1016/j.molstruc.2017.08.10452https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhsVyrtbbE&md5=1a6a3793f854353bd3dae05928e18607Role of dispersion corrected hybrid GGA class in accurately calculating the bond dissociation energy of carbon halogen bond: A benchmark studyKosar, Naveen; Mahmood, Tariq; Ayub, KhurshidJournal of Molecular Structure (2017), 1150 (), 447-458CODEN: JMOSB4; ISSN:0022-2860. (Elsevier B.V.)Benchmark study has been carried out to find a cost effective and accurate method for bond dissocn. energy (BDE) of carbon halogen (C-X) bond. BDE of C-X bond plays a vital role in chem. reactions, particularly for kinetic barrier and thermochem. etc. The compds. with C-X bond used for the benchmark study are important reactants in org., inorg. and bioorg. chem. Exptl. data of C-X bond dissocn. energy is compared with theor. results. The statistical anal. tools such as root mean square deviation (RMSD), std. deviation (SD), Pearson's correlation (R) and mean abs. error (MAE) are used for comparison. Overall, thirty-one d. functionals from eight different classes of d. functional theory (DFT) along with Pople and Dunning basis sets are evaluated. Among different classes of DFT, the dispersion cor. range sepd. hybrid GGA class along with 6-31G(d), 6-311G(d), aug-cc-pVDZ and aug-cc-pVTZ basis sets performed best for bond dissocn. energy calcn. of C-X bond. The ωB97XD showed the best performance with less deviations (RMSD, SD), mean abs. error (MAE) and a significant Pearson's correlation (R) when compared to exptl. data. The ωB97XD along with Pople basis set 6-311g(d) has RMSD, SD, R and MAE of 3.14 kcal mol-1, 3.05 kcal mol-1, 0.97 and -1.07 kcal mol-1, resp.
- 53Van Voorhis, T.; Head-Gordon, M. Benchmark Variational Coupled Cluster Doubles Results. J. Chem. Phys. 2000, 113, 8873– 8879, DOI: 10.1063/1.131964353https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXnvFaitrY%253D&md5=6d6dadea9a39b990c3d2f5109b066b89Benchmark variational coupled cluster doubles resultsVan Voorhis, Troy; Head-Gordon, MartinJournal of Chemical Physics (2000), 113 (20), 8873-8879CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present the first application of the Rayleigh-Ritz variational procedure to the coupled cluster doubles trial function. The variational approach is applied to the potential surface of H4, the double dissocn. of water and the dissocn. of N2, and the results are compared to std. coupled cluster doubles calcns. It is found that the variational approach gives a greatly improved description of strongly correlated systems, where the std. approach is known to fail. Some examn. of the basis set dependence of the results is presented.
- 54Mardirossian, N.; Parkhill, J. A.; Head-Gordon, M. Benchmark Results for Empirical Post-GGA Functionals: Difficult Exchange Problems and Independent Tests. Phys. Chem. Chem. Phys. 2011, 13, 19325, DOI: 10.1039/c1cp21635j54https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhtlKju73J&md5=5b6c30e1d90d38ddd4d3af2e0a8263e2Benchmark results for empirical post-GGA functionals: Difficult exchange problems and independent testsMardirossian, Narbe; Parkhill, John A.; Head-Gordon, MartinPhysical Chemistry Chemical Physics (2011), 13 (43), 19325-19337CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)A review. Many of the most promising new d. functionals have improved the treatment of non-local exchange effects with the help of semi-empirical information and more sophisticated recipes for combining Hartree-Fock and local exchange approxns. In order to quantify recent advancements and identify directions for improvement, we have examd. a broad spectrum of test problems. We evaluate the performance of several new hybrid d. functionals (ωB97, ωB97X, ωB97X-D, LRC-ωPBEh, M06, M06-2X, and M06-HF) on a variety of chem. problems, some sensitive to the treatment of exact exchange (which we have hoped to systematically improve) and some which require a balanced treatment of correlation. Since all of the functionals under consideration are parameterized with ground-state thermochem. data, the benchmark aims to det. the applicability of the new d. functionals to cases that have not been considered in the optimization of the semi-empirical parameters. The first class of benchmarks includes the excitation energies of 21 mols. (83 states) primarily from a recent benchmark conducted by Tozer and co-workers, with some addnl. refs. from data made available from the groups of Thiel and Truhlar. We briefly examine the conformational preferences of a small peptide and complete our study with two recently published sets of data that have shown large, systematic errors in simple alkane thermochem. While our results indicate that the more general hybrids currently under development perform well for problems outside of their parameterization and improve over the std. hybrid d. functionals in an essentially systematic way, there is still a significant self-interaction error in the more difficult cases. Functionals based on a range-sepn. of exchange and functionals depending on the kinetic-energy d. both perform comparably, and there is evidence for complementary strengths.
- 55Gomes, J.; Zimmerman, P. M.; Head-Gordon, M.; Bell, A. T. Accurate Prediction of Hydrocarbon Interactions with Zeolites Utilizing Improved Exchange-Correlation Functionals and QM/MM Methods: Benchmark Calculations of Adsorption Enthalpies and Application to Ethene Methylation by Methanol. J. Phys. Chem. C 2012, 116, 15406– 15414, DOI: 10.1021/jp303321s55https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XoslGhs7c%253D&md5=edf14af043baa6493a5f71e1e9d9f11dAccurate Prediction of Hydrocarbon Interactions with Zeolites Utilizing Improved Exchange-Correlation Functionals and QM/MM Methods: Benchmark Calculations of Adsorption Enthalpies and Application to Ethene Methylation by MethanolGomes, Joseph; Zimmerman, Paul M.; Head-Gordon, Martin; Bell, Alexis T.Journal of Physical Chemistry C (2012), 116 (29), 15406-15414CODEN: JPCCCK; ISSN:1932-7447. (American Chemical Society)The adsorption enthalpy of light hydrocarbon mols. in both acidic and neutral zeolite MFI has been investigated with a range of computational methods. The role of cluster model size and d. functional theory methodol. is examd. by comparison with high quality ab initio wave function theory MP2 results and exptl. detd. heats of adsorption. The commonly applied B3LYP functional performs poorly in benchmark studies due to an inadequate description of intermol. interactions. The functionals ωB97X-D and M06-2X predict adsorption enthalpies consistent with exptl. values over three classes of adsorbates. A hybrid quantum mechanics/mol. mechanics (QM/MM) method is required to converge calcd. thermochem. properties with respect to cluster model size in a manner that is computationally efficient. The accuracy of both QM and QM/MM methods is highly sensitive to choice of level of theory and cluster size, requiring the use of large basis sets, large cluster models, and a d. functional capable of capturing intermol. interactions for achieving the desired chem. accuracy of 2 kcal/mol with respect to exptl. detd. adsorption enthalpies and activation barriers. The computational effort for performing QM/MM simulations is considerably lower than that of similar quality QM results, and allows for the chem. accurate simulation of chem. reactions occurring in zeolites in a manner that is computationally cost efficient without sacrificing accuracy. The resulting QM/MM procedure is applied to study the reaction of ethene methylation by methanol.
Supporting Information
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsomega.2c01331.
Rmsd; SD; Pearson’s correlation (R); MAE; and experimental BDEs of all selected transition-metal complexes with oxygen molecules (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.