1
|
Grotjahn R, Purnomo J, Jin D, Lutfi N, Furche F. Chemically Accurate Singlet-Triplet Gaps of Arylcarbenes from Local Hybrid Density Functionals. J Phys Chem A 2024; 128:6046-6060. [PMID: 39012067 DOI: 10.1021/acs.jpca.4c02852] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/17/2024]
Abstract
Singlet-triplet (ST) gaps are key descriptors of carbenes, because their properties and reactivity are strongly spin-dependent. However, the theoretical prediction of ST gaps is challenging and generally thought to require elaborate correlated wave function methods or double-hybrid density functionals. By evaluating two recent test sets of arylcarbenes (AC12 and AC18), we show that local hybrid functionals based on the "common t" local mixing function (LMF) model achieve mean absolute errors below 1 kcal/mol at a computational cost only slightly higher than that of global hybrid functionals. An analysis of correlation contributions to the ST gaps suggests that the accuracy of the common t-LMF model is mainly due to an improved description of nondynamical correlation which, unlike exchange, is not additive in each spin-channel. Although spin-nonadditivity can be achieved using the local spin polarization alone, using the "common", i.e., spin-unresolved, iso-orbital indicator t for constructing the LMF is found to be critical for consistent accuracy in ST gaps of arylcarbenes. The results support the view of LHs as vehicles to improve the description of nondynamical correlation rather than sophisticated exchange mixing approaches.
Collapse
Affiliation(s)
- Robin Grotjahn
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Justin Purnomo
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Dayun Jin
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Nicolas Lutfi
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
2
|
Kaupp M, Wodyński A, Arbuznikov AV, Fürst S, Schattenberg CJ. Toward the Next Generation of Density Functionals: Escaping the Zero-Sum Game by Using the Exact-Exchange Energy Density. Acc Chem Res 2024; 57:1815-1826. [PMID: 38905497 PMCID: PMC11223257 DOI: 10.1021/acs.accounts.4c00209] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2024] [Revised: 05/29/2024] [Accepted: 06/03/2024] [Indexed: 06/23/2024]
Abstract
ConspectusKohn-Sham density functional theory (KS DFT) is arguably the most widely applied electronic-structure method with tens of thousands of publications each year in a wide variety of fields. Its importance and usefulness can thus hardly be overstated. The central quantity that determines the accuracy of KS DFT calculations is the exchange-correlation functional. Its exact form is unknown, or better "unknowable", and therefore the derivation of ever more accurate yet efficiently applicable approximate functionals is the "holy grail" in the field. In this context, the simultaneous minimization of so-called delocalization errors and static correlation errors is the greatest challenge that needs to be overcome as we move toward more accurate yet computationally efficient methods. In many cases, an improvement on one of these two aspects (also often termed fractional-charge and fractional-spin errors, respectively) generates a deterioration in the other one. Here we report on recent notable progress in escaping this so-called "zero-sum-game" by constructing new functionals based on the exact-exchange energy density. In particular, local hybrid and range-separated local hybrid functionals are discussed that incorporate additional terms that deal with static correlation as well as with delocalization errors. Taking hints from other coordinate-space models of nondynamical and strong electron correlations (the B13 and KP16/B13 models), position-dependent functions that cover these aspects in real space have been devised and incorporated into the local-mixing functions determining the position-dependence of exact-exchange admixture of local hybrids as well as into the treatment of range separation in range-separated local hybrids. While initial functionals followed closely the B13 and KP16/B13 frameworks, meanwhile simpler real-space functions based on ratios of semilocal and exact-exchange energy densities have been found, providing a basis for relatively simple and numerically convenient functionals. Notably, the correction terms can either increase or decrease exact-exchange admixture locally in real space (and in interelectronic-distance space), leading even to regions with negative admixture in cases of particularly strong static correlations. Efficient implementations into a fast computer code (Turbomole) using seminumerical integration techniques make such local hybrid and range-separated local hybrid functionals promising new tools for complicated composite systems in many research areas, where simultaneously small delocalization errors and static correlation errors are crucial. First real-world application examples of the new functionals are provided, including stretched bonds, symmetry-breaking and hyperfine coupling in open-shell transition-metal complexes, as well as a reduction of static correlation errors in the computation of nuclear shieldings and magnetizabilities. The newest versions of range-separated local hybrids (e.g., ωLH23tdE) retain the excellent frontier-orbital energies and correct asymptotic exchange-correlation potential of the underlying ωLH22t functional while improving substantially on strong-correlation cases. The form of these functionals can be further linked to the performance of the recent impactful deep-neural-network "black-box" functional DM21, which itself may be viewed as a range-separated local hybrid.
Collapse
Affiliation(s)
- Martin Kaupp
- Institut für Chemie,
Theoretische Chemie/Quantenchemie, Technische
Universität Berlin, Sekr. C7, Strasse des 17. Juni 115, 10623 Berlin, Germany
| | - Artur Wodyński
- Institut für Chemie,
Theoretische Chemie/Quantenchemie, Technische
Universität Berlin, Sekr. C7, Strasse des 17. Juni 115, 10623 Berlin, Germany
| | - Alexei V. Arbuznikov
- Institut für Chemie,
Theoretische Chemie/Quantenchemie, Technische
Universität Berlin, Sekr. C7, Strasse des 17. Juni 115, 10623 Berlin, Germany
| | - Susanne Fürst
- Institut für Chemie,
Theoretische Chemie/Quantenchemie, Technische
Universität Berlin, Sekr. C7, Strasse des 17. Juni 115, 10623 Berlin, Germany
| | - Caspar J. Schattenberg
- Institut für Chemie,
Theoretische Chemie/Quantenchemie, Technische
Universität Berlin, Sekr. C7, Strasse des 17. Juni 115, 10623 Berlin, Germany
| |
Collapse
|
3
|
Grotjahn R. Learning from the 4-(dimethylamino)benzonitrile twist: Two-parameter range-separated local hybrid functional with high accuracy for triplet and charge-transfer excitations. J Chem Phys 2023; 159:174102. [PMID: 37909451 DOI: 10.1063/5.0173701] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 10/06/2023] [Indexed: 11/03/2023] Open
Abstract
The recent ωLH22t range-separated local hybrid (RSLH) is shown to provide outstanding accuracy for the notorious benchmark problem of the two lowest excited-state potential energy curves for the amino group twist in 4-(dimethylamino)benzonitrile (DMABN). However, the design of ωLH22t as a general-purpose functional resulted in less convincing performance for triplet excitations, which is an important advantage of previous LHs. Furthermore, ωLH22t uses 8 empirical parameters to achieve broad accuracy. In this work, the RSLH ωLH23ct-sir is constructed with minimal empiricism by optimizing its local mixing function prefactor and range-separation parameter for only 8 excitation energies. ωLH23ct-sir maintains the excellent performance of ωLH22t for the DMABN twist and charge-transfer benchmarks but significantly improves the errors for triplet excitation energies (0.17 vs 0.24 eV). Additional test calculations for the AE6BH6 thermochemistry test set and large dipole moment and static polarizability test sets confirm that the focus on excitation energies in the optimization of ωLH23ct-sir has not caused any dramatic errors for ground-state properties. Although ωLH23ct-sir cannot replace ωLH22t as a general-purpose functional, it is preferable for problems requiring a universally good description of localized and charge-transfer excitations of both singlet and triplet multiplicity. Current limitations on the application of ωLH23ct-sir and other RSLHs to the study of singlet-triplet gaps of emitters for thermally activated delayed fluorescence are discussed. This work also includes the first systematic analysis of the influence of the local mixing function prefactor and the range-separation parameter in an RSLH on different types of excitations.
Collapse
Affiliation(s)
- Robin Grotjahn
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
4
|
Franzke Y, Holzer C, Andersen JH, Begušić T, Bruder F, Coriani S, Della Sala F, Fabiano E, Fedotov DA, Fürst S, Gillhuber S, Grotjahn R, Kaupp M, Kehry M, Krstić M, Mack F, Majumdar S, Nguyen BD, Parker SM, Pauly F, Pausch A, Perlt E, Phun GS, Rajabi A, Rappoport D, Samal B, Schrader T, Sharma M, Tapavicza E, Treß RS, Voora V, Wodyński A, Yu JM, Zerulla B, Furche F, Hättig C, Sierka M, Tew DP, Weigend F. TURBOMOLE: Today and Tomorrow. J Chem Theory Comput 2023; 19:6859-6890. [PMID: 37382508 PMCID: PMC10601488 DOI: 10.1021/acs.jctc.3c00347] [Citation(s) in RCA: 44] [Impact Index Per Article: 44.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Indexed: 06/30/2023]
Abstract
TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.
Collapse
Affiliation(s)
- Yannick
J. Franzke
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| | - Christof Holzer
- Institute
of Theoretical Solid State Physics, Karlsruhe
Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Josefine H. Andersen
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Tomislav Begušić
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Florian Bruder
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| | - Sonia Coriani
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Fabio Della Sala
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for
Biomolecular Nanotechnologies @UNILE, Istituto
Italiano di Tecnologia, Via Barsanti, 73010 Arnesano, Italy
| | - Eduardo Fabiano
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for
Biomolecular Nanotechnologies @UNILE, Istituto
Italiano di Tecnologia, Via Barsanti, 73010 Arnesano, Italy
| | - Daniil A. Fedotov
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Susanne Fürst
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Sebastian Gillhuber
- Institute
of Inorganic Chemistry, Karlsruhe Institute
of Technology (KIT), Engesserstr. 15, 76131 Karlsruhe, Germany
| | - Robin Grotjahn
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Martin Kaupp
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Max Kehry
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Marjan Krstić
- Institute
of Theoretical Solid State Physics, Karlsruhe
Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Fabian Mack
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Sourav Majumdar
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Brian D. Nguyen
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Shane M. Parker
- Department
of Chemistry, Case Western Reserve University, 10900 Euclid Ave, Cleveland, Ohio 44106 United States
| | - Fabian Pauly
- Institute
of Physics, University of Augsburg, Universitätsstr. 1, 86159 Augsburg, Germany
| | - Ansgar Pausch
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Eva Perlt
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Gabriel S. Phun
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Ahmadreza Rajabi
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Dmitrij Rappoport
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Bibek Samal
- Department
of Chemical Sciences, Tata Institute of
Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Tim Schrader
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Manas Sharma
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Enrico Tapavicza
- Department
of Chemistry and Biochemistry, California
State University, Long Beach, 1250 Bellflower Boulevard, Long
Beach, California 90840-9507, United States
| | - Robert S. Treß
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Vamsee Voora
- Department
of Chemical Sciences, Tata Institute of
Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Artur Wodyński
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Jason M. Yu
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Benedikt Zerulla
- Institute
of Nanotechnology, Karlsruhe Institute of
Technology (KIT), Hermann-von-Helmholtz-Platz
1, 76344 Eggenstein-Leopoldshafen Germany
| | - Filipp Furche
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Christof Hättig
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Marek Sierka
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - David P. Tew
- Physical
and Theoretical Chemistry Laboratory, University
of Oxford, South Parks
Road, Oxford OX1 3QZ, United Kingdom
| | - Florian Weigend
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| |
Collapse
|
5
|
Liu X, McKemmish L, Pérez-Ríos J. The performance of CCSD(T) for the calculation of dipole moments in diatomics. Phys Chem Chem Phys 2023; 25:4093-4104. [PMID: 36651174 DOI: 10.1039/d2cp05060a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This work analyzes the accuracy of the coupled cluster with single, double, and perturbative triple excitation [CCSD(T)] method for predicting dipole moments. In particular, we benchmark CCSD(T) predictions for the equilibrium bond length, vibrational frequency, and dipole moment versus accurate experimental data. As a result, we find that CCSD(T) leads to accurate dipole moments. However, in some cases, it disagrees with the experimental values, and the disagreement can not be satisfactorily explained via relativistic or multi-reference effects. Therefore, our results indicate that benchmark studies for energy and geometry properties do not accurately describe other electron density magnitudes.
Collapse
Affiliation(s)
- Xiangyue Liu
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Laura McKemmish
- School of Chemistry, UNSW Sydney, Sydney, NSW 2052, Australia
| | - Jesús Pérez-Ríos
- Department of Physics and Astronomy, Stony Brook University, Stony Brook 11794, New York, USA. .,Institute for Advanced Computational Science, Stony Brook University, Stony Brook, NY 11794-3800, USA
| |
Collapse
|
6
|
Fürst S, Haasler M, Grotjahn R, Kaupp M. Full Implementation, Optimization, and Evaluation of a Range-Separated Local Hybrid Functional with Wide Accuracy for Ground and Excited States. J Chem Theory Comput 2023; 19:488-502. [PMID: 36625881 DOI: 10.1021/acs.jctc.2c00782] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We report the first full and efficient implementation of range-separated local hybrid functionals (RSLHs) into the TURBOMOLE program package. This enables the computation of ground-state energies and nuclear gradients as well as excitation energies. Regarding the computational effort, RSLHs scale like regular local hybrid functionals (LHs) with system or basis set size and increase timings by a factor of 2-3 in total. An advanced RSLH, ωLH22t, has been optimized for atomization energies and reaction barriers. It is an extension of the recent LH20t local hybrid and is based on short-range PBE and long-range HF exchange-energy densities, a pig2 calibration function to deal with the gauge ambiguity of exchange-energy densities, and reoptimized B95c correlation. ωLH22t has been evaluated for a wide range of ground-state and excited-state quantities. It further improves upon the already successful LH20t functional for the GMTKN55 main-group energetics test suite, and it outperforms any global hybrid while performing close to the top rung-4 functional, ωB97M-V, for these evaluations when augmented by D4 dispersion corrections. ωLH22t performs excellently for transition-metal reactivity and provides good balance between delocalization errors and left-right correlation for mixed-valence systems, with a somewhat larger bias toward localized states compared to LH20t. It approaches the accuracy of the best local hybrids to date for core, valence singlet and triplet, and Rydberg excitation energies while improving strikingly on intra- and intermolecular charge-transfer excitations, comparable to the most successful range-separated hybrids available.
Collapse
Affiliation(s)
- Susanne Fürst
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Matthias Haasler
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
7
|
Wodyński A, Kaupp M. Local Hybrid Functional Applicable to Weakly and Strongly Correlated Systems. J Chem Theory Comput 2022; 18:6111-6123. [PMID: 36170626 DOI: 10.1021/acs.jctc.2c00795] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recent idea (Wodyński, A.; Arbuznikov, A. V.; Kaupp M. J. Chem. Phys. 2021, 155, 144101) to augment local hybrid functionals by a strong-correlation (sc) factor obtained from the adiabatic connection in the spirit of the KP16 model has been extended and applied to generate the accurate sc-corrected local hybrid functional scLH22t. By damping small values of the ratio between nondynamical and dynamical correlation entering the correction factor, it has become possible to avoid double counting of nondynamical correlation for weakly correlated situations and thereby preserve the excellent accuracy of the underlying LH20t local hybrid for such cases almost perfectly. On the other hand, scLH22t improves substantially over LH20t in reducing fractional-spin errors (FSEs), in providing improved spin-restricted bond dissociation curves, and in treating some typical systems with multireference character. The obtained FSEs are similar to those of the KP16/B13 model and slightly larger than for B13, but performance for weakly correlated systems is better than for these two related methods, which are also difficult to use self-consistently. The recent DM21 functional based on the training of a deep neural network still performs somewhat better than scLH22t but allows no physical insights into the origins of reduced FSEs. Examination of local mixing functions (LMFs) for the corrected scLH22t and uncorrected LH20t functionals provides further insights: in weakly correlated situations, the LMF remains essentially unchanged. Strong-correlation effects manifest in a reduction of the LMF values in certain regions of space, even to the extent of producing negative LMF values. It is suggested that this is the mechanism by which also DM21, which may be viewed as a range-separated local hybrid, is able to reduce FSEs.
Collapse
Affiliation(s)
- Artur Wodyński
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
8
|
Grotjahn R, Kaupp M. A Look at Real‐World Transition‐Metal Thermochemistry and Kinetics with Local Hybrid Functionals. Isr J Chem 2022. [DOI: 10.1002/ijch.202200021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Robin Grotjahn
- Technische Universität Berlin Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7 Straße des 17. Juni 135 D-10623 Berlin Germany
| | - Martin Kaupp
- Technische Universität Berlin Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7 Straße des 17. Juni 135 D-10623 Berlin Germany
| |
Collapse
|
9
|
Soyemi A, Szilvási T. Benchmarking Semiempirical QM Methods for Calculating the Dipole Moment of Organic Molecules. J Phys Chem A 2022; 126:1905-1921. [PMID: 35290045 DOI: 10.1021/acs.jpca.1c10144] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The dipole moment is a simple descriptor of the charge distribution and polarity and is important for understanding and predicting various molecular properties. Semiempirical (SE) methods offer a cost-effective way to calculate dipole moment that can be used in high-throughput screening applications although the accuracy of the methods is still in question. Therefore, we have evaluated AM1, GFN0-xTB, GFN1-xTB, GFN2-xTB, PM3, PM6, PM7, B97-3c, HF-3c, and PBEh-3c SE methods, which cover a variety of SE approximations, to directly assess the performance of SE methods in predicting molecular dipole moments and their directions using 7211 organic molecules contained in the QM7b database. We find that B97-3c and PBEh-3c perform best against coupled-cluster reference values yielding dipole moments with a mean absolute error (MAE) of 0.10 and 0.11 D, respectively, which is similar to the MAE of density functional theory (DFT) methods (∼0.1 D) reported earlier. Analysis of the atomic composition shows that all SE methods show good performance for hydrocarbons for which the spread of error was within 1 D of the reference data, while the worst performances are for sulfur-containing compounds for which only B97-3c and PBEh-3c show acceptable performance. We also evaluate the effect of SE optimized geometry, instead of the benchmark DFT geometry, that shows a dramatic drop in the performance of AM1 and PM3 for which the range of error tripled. Based on our overall findings, we highlight that there is an optimal compromise between accuracy and computational cost using GFN2-xTB (MAE: 0.25 D) that is 3 orders of magnitude faster than B97-3c and PBEh-3c. Thus, we recommend using GFN2-xTB for cost-efficient calculation of the dipole moment of organic molecules containing C, H, O, and N atoms, whereas, for sulfur-containing organic molecules, we suggest PBEh-3c.
Collapse
Affiliation(s)
- Ademola Soyemi
- Department of Chemical and Biological Engineering, University of Alabama, Tuscaloosa, Alabama 35487, United States
| | - Tibor Szilvási
- Department of Chemical and Biological Engineering, University of Alabama, Tuscaloosa, Alabama 35487, United States
| |
Collapse
|
10
|
Grotjahn R, Kaupp M. Assessment of hybrid functionals for singlet and triplet excitations: Why do some local hybrid functionals perform so well for triplet excitation energies? J Chem Phys 2021; 155:124108. [PMID: 34598568 DOI: 10.1063/5.0063751] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
The performance of various hybrid density functionals is assessed for 105 singlet and 105 corresponding triplet vertical excitation energies from the QUEST database. The overall lowest mean absolute error is obtained with the local hybrid (LH) functional LH12ct-SsirPW92 with individual errors of 0.11 eV (0.11 eV) for singlet (triplet) n → π* excitations and 0.29 eV (0.17 eV) for π → π* excitations. This is slightly better than with the overall best performing global hybrid M06-2X [n → π*: 0.13 eV (0.17 eV), π → π*: 0.30 eV (0.20 eV)], while most other global and range-separated hybrids and some LHs suffer from the "triplet problem" of time-dependent density functional theory. This is exemplified by correlating the errors for singlet and triplet excitations on a state-by-state basis. The excellent performance of LHs based on a common local mixing function, i.e., an LMF constructed from the spin-summed rather than the spin-resolved semilocal quantities, is systematically investigated by the introduction of a spin-channel interpolation scheme that allows us to continuously modulate the fraction of opposite-spin terms used in the LMF. The correlation of triplet and singlet errors is systematically improved for the n → π* excitations when larger fractions of the opposite-spin-channel are used in the LMF, whereas this effect is limited for the π → π* excitations. This strongly supports a previously made hypothesis that attributes the excellent performance of LHs based on a common LMF to cross-spin-channel nondynamical correlation terms.
Collapse
Affiliation(s)
- Robin Grotjahn
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| |
Collapse
|
11
|
Grotjahn R, Kaupp M. Reliable TDDFT Protocol Based on a Local Hybrid Functional for the Prediction of Vibronic Phosphorescence Spectra Applied to Tris(2,2'-bipyridine)-Metal Complexes. J Phys Chem A 2021; 125:7099-7110. [PMID: 34370482 DOI: 10.1021/acs.jpca.1c05101] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An efficient computational protocol for the prediction of vibrationally resolved phosphorescence spectra is developed and validated for five tris(2,2'-bipyridine)-metal complexes ([M(bpy)3]n+, where M = Zn, Ru, Rh, Os, Ir). The outstanding feature of this protocol is the use of full linear-response time-dependent density functional theory (TDDFT) for the excited-state triplet calculation, i.e., the commonly seen strategies employing the Tamm-Dancoff approximation (TDA) or unrestricted density functional theory (DFT) calculations for the T1 state are not needed. This is achieved by the use of a local hybrid functional (LH12ct-SsirPW92) that features a real-space dependent admixture of exact exchange governed by a local mixing function. The excellent performance of this LH for triplet excitation energies known from previous studies transfers to a remarkable mean absolute error of 0.06 eV for the phosphorescence 0-0 energies investigated herein, while the popular B3PW91 functional gives an error of 0.27 eV in TDDFT and 0.09 eV in unrestricted DFT calculations, respectively. The advantages of the local hybrid are particularly apparent for excited states with a mixed-valence character. The influence of spin-orbit coupling was found to be significant for [Os(bpy)3]2+ red-shifting the 0-0 energy for phosphorescence by 0.17 eV, while the effect is negligible for the other complexes (<0.03 eV). The influence of the basis-set and integration-grid sizes is evaluated, and a computationally lighter protocol is validated that leads to drastic savings in computation time with negligible loss in accuracy.
Collapse
Affiliation(s)
- Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
12
|
Loukhovitski BI, Sharipov AS. Molecular Collision Diameters and Electronic Polarizabilities: Inherent Relationship and Fast Evaluation. J Phys Chem A 2021; 125:5117-5123. [PMID: 34077208 DOI: 10.1021/acs.jpca.1c02201] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The collision diameter σ for a large set of molecular species is related to the static electronic polarizability αel. A remarkable correlation between these quantities conceptually similar to the analogous one previously identified for atoms is revealed. Our recommended model is the function σ(αel) = p1 + p2αel1/3, where p1 = 0.768 Å and p2 = 2.168 are the fitting parameters providing the best overall match to the reference data for collision diameter (181 data points). The obtained correlation allows one to easily find the collision diameter of molecules from the known polarizability and vice versa. These findings can be useful for many applications, where there is a need for inexpensive assessments of the collision diameters or electronic polarizabilities, for example, when developing the transport property databases for modeling of chemically reacting flows.
Collapse
|
13
|
Hait D, Liang YH, Head-Gordon M. Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems. J Chem Phys 2021; 154:074109. [DOI: 10.1063/5.0038694] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Yu Hsuan Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|