1
|
Mehta N, Martin JML. On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment. J Comput Chem 2024; 45:1017-1032. [PMID: 38216516 DOI: 10.1002/jcc.27294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 12/05/2023] [Accepted: 12/10/2023] [Indexed: 01/14/2024]
Abstract
Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| |
Collapse
|
2
|
Wappett DA, Goerigk L. Exploring CPS-Extrapolated DLPNO-CCSD(T 1) Reference Values for Benchmarking DFT Methods on Enzymatically Catalyzed Reactions. J Phys Chem A 2024; 128:62-72. [PMID: 38124376 DOI: 10.1021/acs.jpca.3c05086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
Domain-based local pair natural orbital coupled-cluster singles doubles with perturbative triples [DLPNO-CCSD(T)] is regularly used to calculate reliable benchmark reference values at a computational cost significantly lower than that of canonical CCSD(T). Recent work has shown that even greater accuracy can be obtained at only a small additional cost through extrapolation to the complete PNO space (CPS) limit. Herein, we test two levels of CPS extrapolation, CPS(5,6), which approximates the accuracy of standard TightPNO, and CPS(6,7), which surpasses it, as benchmark values to test density functional approximations (DFAs) on a small set of organic and transition-metal-dependent enzyme active site models. Between the different reference levels of theory, there are changes in the magnitudes of the absolute deviations for all functionals, but these are small and there is minimal impact on the relative rankings of the tested DFAs. The differences are more significant for the metalloenzymes than the organic enzymes, so we repeat the tests on our entire ENZYMES22 set of organic enzyme active site models [Wappett, D. A.; Goerigk, L. J. Phys. Chem. A 2019, 123, 7057-7074] to confirm that using the CPS extrapolations for the reference values has negligible impact on the benchmarking outcomes. This means that we can particularly recommend CPS(5,6) as an alternative to standard TightPNO settings for calculating reference values, increasing the applicability of DLPNO-CCSD(T) in benchmarking reaction energies and barrier heights of larger models of organic enzymes. DLPNO-CCSD(T1)/CPS(6,7) energies for ENZYMES22 are finally presented as updated reference values for the set, reflecting the recent improvements in the method.
Collapse
Affiliation(s)
- Dominique A Wappett
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
3
|
Wappett D, Goerigk L. Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set. J Chem Theory Comput 2023; 19:8365-8383. [PMID: 37943578 PMCID: PMC10688432 DOI: 10.1021/acs.jctc.3c00558] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Revised: 10/13/2023] [Accepted: 10/16/2023] [Indexed: 11/10/2023]
Abstract
We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.
Collapse
Affiliation(s)
- Dominique
A. Wappett
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| |
Collapse
|
4
|
de Oliveira MT, Alves JMA, Vrech NL, Braga AAC, Barboza CA. A comprehensive benchmark investigation of quantum chemical methods for carbocations. Phys Chem Chem Phys 2023; 25:1903-1922. [PMID: 36541431 DOI: 10.1039/d2cp04603b] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.
Collapse
Affiliation(s)
- Marcelo T de Oliveira
- Department of Chemistry and Physics, La Trobe Institute of Molecular Sciences, La Trobe University, Melbourne, Victoria 3086, Australia. .,Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Júlia M A Alves
- Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Natália L Vrech
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Ataualpa A C Braga
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Cristina A Barboza
- Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland.,Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, Warsaw, Poland
| |
Collapse
|
5
|
Santra G, Calinsky R, Martin JML. Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights. J Phys Chem A 2022; 126:5492-5505. [PMID: 35930677 PMCID: PMC9393870 DOI: 10.1021/acs.jpca.2c03922] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
To better understand the thermochemical kinetics and
mechanism
of a specific chemical reaction, an accurate estimation of barrier
heights (forward and reverse) and reaction energies is vital. Because
of the large size of reactants and transition state structures involved
in real-life mechanistic studies (e.g., enzymatically catalyzed reactions),
density functional theory remains the workhorse for such calculations.
In this paper, we have assessed the performance of 91 density functionals
for modeling the reaction energies and barrier heights on a large
and chemically diverse data set (BH9) composed of 449 organic chemistry
reactions. We have shown that range-separated hybrid functionals perform
better than the global hybrids for BH9 barrier heights and reaction
energies. Except for the PBE-based range-separated nonempirical double
hybrids, range separation of the exchange term helps improve the performance
for barrier heights and reaction energies. The 16-parameter Berkeley
double hybrid, ωB97M(2), performs remarkably well for both properties.
However, our minimally empirical range-separated double hybrid functionals
offer marginally better accuracy than ωB97M(2) for BH9 barrier
heights and reaction energies.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Rivka Calinsky
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
6
|
Santra G, Martin JML. Do Double-Hybrid Functionals Benefit from Regularization in the PT2 Term? Observations from an Extensive Benchmark. J Phys Chem Lett 2022; 13:3499-3506. [PMID: 35417181 PMCID: PMC9036584 DOI: 10.1021/acs.jpclett.2c00718] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We put to the test a recent suggestion [Shee, J., et al. J. Phys. Chem. Lett. 2021, 12 (50), 12084-12097] that MP2 regularization might improve the performance of double-hybrid density functionals. Using the very large and chemically diverse GMTKN55 benchmark, we find that κ-regularization is indeed beneficial at lower percentages of Hartree-Fock exchange, especially if spin-component scaling is not applied [such as in B2GP-PLYP or ωB97M(2)]. This benefit dwindles for DSD and DOD functionals and vanishes entirely in the ∼70% HF exchange region optimal for them.
Collapse
|
7
|
Alipour M, Izadkhast T. Do any types of double-hybrid models render the correct order of excited state energies in inverted singlet–triplet emitters? J Chem Phys 2022; 156:064302. [DOI: 10.1063/5.0077722] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| | - Tahereh Izadkhast
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| |
Collapse
|
8
|
Alipour M, Karimi N. Spin-Opposite-Scaled Range-Separated Exchange Double-Hybrid Models (SOS-RSX-DHs): Marriage Between DH and RSX/SOS-RSX Is Not Always a Happy Match. J Chem Theory Comput 2021; 17:4077-4091. [PMID: 34085815 DOI: 10.1021/acs.jctc.1c00271] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The range-separated version of double-hybrid density functional theory (DH-DFT) with a remarkable efficiency for both ground-state and excited-state characteristics has recently come into spotlight. In this work, based on theoretical arguments, several variants of spin-opposite-scaled range-separated exchange double-hybrid models (SOS-RSX-DHs) have been proposed and validated. More specifically, we first extend the RSX-DHs to design some other related models. Next, the SOS version of the resulting approximations is constructed and thoroughly evaluated using standard benchmark compilations of various properties. It is shown that although there are properties for which the RSX-DH and SOS-RSX-DH frameworks are rival, there are still some problems particularly prone to the self-interaction error issues where our proposed models seem to be beneficial. Furthermore, some of the presented models devoid of any additional corrections can also surpass the recently proposed approximations from different rungs of "Jacob's Ladder". Nonetheless, perusing the results of different methods and detailed comparisons with the predecessors discloses that all things may not necessarily be well with the RSX and SOS-RSX schemes, where the parent DHs as well as their SOS counterparts can still come into play.
Collapse
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| | - Niloofar Karimi
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| |
Collapse
|
9
|
Wappett DA, Goerigk L. A guide to benchmarking enzymatically catalysed reactions: the importance of accurate reference energies and the chemical environment. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02770-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
10
|
Mehta N, Fellowes T, White JM, Goerigk L. CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions? J Chem Theory Comput 2021; 17:2783-2806. [PMID: 33881869 DOI: 10.1021/acs.jctc.1c00006] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| | - Thomas Fellowes
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia.,Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Victoria 3010, Australia
| | - Jonathan M White
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia.,Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
11
|
Magalhães RP, Fernandes HS, Sousa SF. Modelling Enzymatic Mechanisms with QM/MM Approaches: Current Status and Future Challenges. Isr J Chem 2020. [DOI: 10.1002/ijch.202000014] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Rita P. Magalhães
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| | - Henriques S. Fernandes
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| | - Sérgio F. Sousa
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| |
Collapse
|
12
|
Woller T, Banerjee A, Sylvetsky N, Santra G, Deraet X, De Proft F, Martin JML, Alonso M. Performance of Electronic Structure Methods for the Description of Hückel-Möbius Interconversions in Extended π-Systems. J Phys Chem A 2020; 124:2380-2397. [PMID: 32093467 PMCID: PMC7307915 DOI: 10.1021/acs.jpca.9b10880] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Revised: 02/21/2020] [Indexed: 12/25/2022]
Abstract
Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. DFT remains the workhorse for modeling such extended macrocycles, when taking into account their size and huge conformational flexibility. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, π···π stacking, steric effects, ring strain, and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wave function methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a stronger degree of static correlation than the Hückel and twisted-Hückel structures, and as a result the relative energies of singly twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between the extrapolated MP2/CBS energies estimated from cc-pV{T,Q}Z, diffuse-augmented heavy-aug-cc-pV{T,Q}Z and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n-particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol-1 with the CCSD(T) relative energies. Regarding DFT methods, the range-separated double hybrids, such as ωB97M(2) and B2GP-PLYP, outperform other functionals with RMSDs of 0.6 and 0.8 kcal mol-1, respectively. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended π-systems, the errors drop down to 1.9 kcal mol-1 for the revised revDOD-PBEP86-NL, which eliminates the same-spin correlation energy. Minnesota meta-GGA functionals with high fractions of exact exchange (M06-2X and M08-HX) also perform reasonably well, outperforming more robust and significantly less empirically parametrized functionals like SCAN0-D3.
Collapse
Affiliation(s)
- Tatiana Woller
- Department
of General Chemistry (ALGC), Faculty of Science and Bio-engineering
Sciences, Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
- Laboratoire
de Chimie Théorique (LCT), Sorbonne
Université, CNRS, F-75005 Paris, France
| | - Ambar Banerjee
- Department
of Organic Chemistry, Weizmann Institute
of Science, 76100 Reḥovot, Israel
| | - Nitai Sylvetsky
- Department
of Organic Chemistry, Weizmann Institute
of Science, 76100 Reḥovot, Israel
| | - Golokesh Santra
- Department
of Organic Chemistry, Weizmann Institute
of Science, 76100 Reḥovot, Israel
| | - Xavier Deraet
- Department
of General Chemistry (ALGC), Faculty of Science and Bio-engineering
Sciences, Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
| | - Frank De Proft
- Department
of General Chemistry (ALGC), Faculty of Science and Bio-engineering
Sciences, Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
| | - Jan M. L. Martin
- Department
of Organic Chemistry, Weizmann Institute
of Science, 76100 Reḥovot, Israel
| | - Mercedes Alonso
- Department
of General Chemistry (ALGC), Faculty of Science and Bio-engineering
Sciences, Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
| |
Collapse
|
13
|
How does SCAN compare to PBE in the framework of parameter-free spin-opposite-scaled double-hybrids? Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2019.136898] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
14
|
Alipour M, Izadkhast T. Appraising spin-state energetics in transition metal complexes using double-hybrid models: accountability of SOS0-PBESCAN0-2(a) as a promising paradigm. Phys Chem Chem Phys 2020; 22:9388-9404. [DOI: 10.1039/d0cp00807a] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Through a comprehensive survey, reliable double-hybrid models have been validated and proposed for spin-state energetics in transition metal complexes.
Collapse
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry
- College of Sciences
- Shiraz University
- Shiraz 71946-84795
- Iran
| | - Tahereh Izadkhast
- Department of Chemistry
- College of Sciences
- Shiraz University
- Shiraz 71946-84795
- Iran
| |
Collapse
|
15
|
Martin JML, Santra G. Empirical Double‐Hybrid Density Functional Theory: A ‘Third Way’ in Between WFT and DFT. Isr J Chem 2019. [DOI: 10.1002/ijch.201900114] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Jan M. L. Martin
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| | - Golokesh Santra
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| |
Collapse
|
16
|
Wappett DA, Goerigk L. Toward a Quantum-Chemical Benchmark Set for Enzymatically Catalyzed Reactions: Important Steps and Insights. J Phys Chem A 2019; 123:7057-7074. [DOI: 10.1021/acs.jpca.9b05088] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
17
|
Mehta N, Casanova-Páez M, Goerigk L. Semi-empirical or non-empirical double-hybrid density functionals: which are more robust? Phys Chem Chem Phys 2018; 20:23175-23194. [PMID: 30062343 DOI: 10.1039/c8cp03852j] [Citation(s) in RCA: 89] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Parkville, Australia.
| | | | | |
Collapse
|
18
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
19
|
Najibi A, Goerigk L. A Comprehensive Assessment of the Effectiveness of Orbital Optimization in Double-Hybrid Density Functionals in the Treatment of Thermochemistry, Kinetics, and Noncovalent Interactions. J Phys Chem A 2018; 122:5610-5624. [DOI: 10.1021/acs.jpca.8b04058] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Asim Najibi
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
20
|
Roch LM, Baldridge KK. General optimization procedure towards the design of a new family of minimal parameter spin-component-scaled double-hybrid density functional theory. Phys Chem Chem Phys 2017; 19:26191-26200. [PMID: 28930316 DOI: 10.1039/c7cp04125j] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A general optimization procedure towards the development and implementation of a new family of minimal parameter spin-component-scaled double-hybrid (mSD) density functional theory (DFT) is presented. The nature of the proposed exchange-correlation functional establishes a methodology with minimal empiricism. This new family of double-hybrid (DH) density functionals is demonstrated using the PBEPBE functional, illustrating the optimization procedure to the mSD-PBEPBE method, and the performance characteristics shown for a set of non-covalent complexes covering a broad regime of weak interactions. With only two parameters, mSD-PBEPBE and its cost-effective counterpart, RI-mSD-PBEPBE, show a mean absolute error of ca. 0.4 kcal mol-1 averaged over 66 weak interacting systems. Following a successive 2D-grid refinement for a CBS extrapolation of the coefficients, the optimization procedure can be recommended for the design and implementation of a variety of additional DH methods using any of the plethora of currently available functionals.
Collapse
Affiliation(s)
- Loïc M Roch
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland.
| | | |
Collapse
|
21
|
On the opposite-spin to same-spin ratio of absolute and interaction MP2 correlation energy in parameter-free spin-opposite-scaled double hybrids. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.07.023] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
22
|
Schwabe T, Goerigk L. Time-Dependent Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling. J Chem Theory Comput 2017; 13:4307-4323. [PMID: 28763220 DOI: 10.1021/acs.jctc.7b00386] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
For the first time, we combine time-dependent double-hybrid density functional approximations (TD-DHDFAs) for the calculation of electronic excitation energies with the concepts of spin-component and spin-opposite scaling (SCS/SOS) of electron-pair contributions to their nonlocal correlation components. Different flavors of this idea, ranging from standard SCS parameters to fully fitted parameter sets, are presented and tested on six different parent DHDFAs. For cross-validation, we assess those methods on three benchmark sets that cover small- to medium-sized chromophores (up to 78 atoms) and different excitation types. For this purpose, we also introduce new CC3 reference values for the popular Gordon benchmark set that we recommend using in future studies. Our results confirm that already the (unscaled) parent TD-DHDFAs are accurate and outperform some wave function methods. Further introduction of SCS/SOS eliminates extreme outliers, reduces deviation spans from reference values by up to 0.5 eV, aligns the performance of the Tamm-Dancoff approximation (TDA) to that of full TD calculations, and also enables a more balanced description of different excitation types. The best-performing TD-based methods in our cross validation have mean absolute deviations as low as 0.14 eV compared to the time- and resource-intensive CC3 approach. A very important finding is that we also obtained SOS variants with excellent performance, contrary to wave function based methods. This opens a future pathway to highly efficient methods for the optimization of excited-state geometries, particularly when paired with computing strategies such as the Laplace transform. We recommend our SCS- and SOS-based variants for further testing and subsequent applications.
Collapse
Affiliation(s)
- Tobias Schwabe
- Center for Bioinformatics and Physical Chemistry Institute, University of Hamburg , Bundesstraße 43, D-20146 Hamburg, Germany
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne , Parkville, Victoria 3010, Australia
| |
Collapse
|
23
|
Alipour M, Karimi N. Dissecting the accountability of parameterized and parameter-free single-hybrid and double-hybrid functionals for photophysical properties of TADF-based OLEDs. J Chem Phys 2017. [PMID: 28641443 DOI: 10.1063/1.4986777] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Organic light emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) emitters are an attractive category of materials that have witnessed a booming development in recent years. In the present contribution, we scrutinize the accountability of parameterized and parameter-free single-hybrid (SH) and double-hybrid (DH) functionals through the two formalisms, full time-dependent density functional theory (TD-DFT) and Tamm-Dancoff approximation (TDA), for the estimation of photophysical properties like absorption energy, emission energy, zero-zero transition energy, and singlet-triplet energy splitting of TADF molecules. According to our detailed analyses on the performance of SHs based on TD-DFT and TDA, the TDA-based parameter-free SH functionals, PBE0 and TPSS0, with one-third of exact-like exchange turned out to be the best performers in comparison to other functionals from various rungs to reproduce the experimental data of the benchmarked set. Such affordable SH approximations can thus be employed to predict and design the TADF molecules with low singlet-triplet energy gaps for OLED applications. From another perspective, considering this point that both the nonlocal exchange and correlation are essential for a more reliable description of large charge-transfer excited states, applicability of the functionals incorporating these terms, namely, parameterized and parameter-free DHs, has also been evaluated. Perusing the role of exact-like exchange, perturbative-like correlation, solvent effects, and other related factors, we find that the parameterized functionals B2π-PLYP and B2GP-PLYP and the parameter-free models PBE-CIDH and PBE-QIDH have respectable performance with respect to others. Lastly, besides the recommendation of reliable computational protocols for the purpose, hopefully this study can pave the way toward further developments of other SHs and DHs for theoretical explorations in the field of OLEDs technology.
Collapse
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, College of Sciences, Shiraz University, Shiraz, Iran
| | - Niloofar Karimi
- Department of Chemistry, College of Sciences, Shiraz University, Shiraz, Iran
| |
Collapse
|