1
|
Holzer C, Franzke YJ. Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory. Chemphyschem 2024; 25:e202400120. [PMID: 38456204 DOI: 10.1002/cphc.202400120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Revised: 03/07/2024] [Accepted: 03/08/2024] [Indexed: 03/09/2024]
Abstract
Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.
Collapse
Affiliation(s)
- Christof Holzer
- Karlsruhe Institute of Technology (KIT), Institute of Theoretical Solid State Physics, Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Yannick J Franzke
- Friedrich Schiller University Jena, Otto Schott Institute of Materials Research, Löbdergraben 32, 07743, Jena, Germany
| |
Collapse
|
2
|
Franzke YJ, Bruder F, Gillhuber S, Holzer C, Weigend F. Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods. J Phys Chem A 2024; 128:670-686. [PMID: 38195394 DOI: 10.1021/acs.jpca.3c07093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2024]
Abstract
An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Florian Bruder
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Sebastian Gillhuber
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology (KIT), Engesserstraße 15, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| |
Collapse
|
3
|
Yuan X, Halbert L, Pototschnig JV, Papadopoulos A, Coriani S, Visscher L, Pereira Gomes AS. Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures. J Chem Theory Comput 2024; 20:677-694. [PMID: 38193434 DOI: 10.1021/acs.jctc.3c00812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.
Collapse
Affiliation(s)
- Xiang Yuan
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, F-59000 Lille, France
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Loïc Halbert
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Johann Valentin Pototschnig
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Anastasios Papadopoulos
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Sonia Coriani
- DTU Chemistry─Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
4
|
Bruder F, Franzke YJ, Holzer C, Weigend F. Zero-field splitting parameters within exact two-component theory and modern density functional theory using seminumerical integration. J Chem Phys 2023; 159:194117. [PMID: 37987521 DOI: 10.1063/5.0175758] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Accepted: 10/26/2023] [Indexed: 11/22/2023] Open
Abstract
An efficient implementation of zero-field splitting parameters based on the work of Schmitt et al. [J. Chem. Phys. 134, 194113 (2011)] is presented. Seminumerical integration techniques are used for the two-electron spin-dipole contribution and the response equations of the spin-orbit perturbation. The original formulation is further generalized. First, it is extended to meta-generalized gradient approximations and local hybrid functionals. For these functional classes, the response of the paramagnetic current density is considered in the coupled-perturbed Kohn-Sham equations for the spin-orbit perturbation term. Second, the spin-orbit perturbation is formulated within relativistic exact two-component theory and the screened nuclear spin-orbit (SNSO) approximation. The accuracy of the implementation is demonstrated for transition-metal and diatomic main-group compounds. The efficiency is assessed for Mn and Mo complexes. Here, it is found that coarse integration grids for the seminumerical schemes lead to drastic speedups while introducing clearly negligible errors. In addition, the SNSO approximation substantially reduces the computational demands and leads to very similar results as the spin-orbit mean field Ansatz.
Collapse
Affiliation(s)
- Florian Bruder
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| |
Collapse
|
5
|
Franzke YJ, Holzer C. Exact two-component theory becoming an efficient tool for NMR shieldings and shifts with spin-orbit coupling. J Chem Phys 2023; 159:184102. [PMID: 37937936 DOI: 10.1063/5.0171509] [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/09/2023] [Accepted: 09/04/2023] [Indexed: 11/09/2023] Open
Abstract
We present a gauge-origin invariant exact two-component (X2C) approach within a modern density functional framework, supporting meta-generalized gradient approximations such as TPSS and range-separated hybrid functionals such as CAM-B3LYP. The complete exchange-correlation kernel is applied, including the direct contribution of the field-dependent basis functions and the reorthonormalization contribution from the perturbed overlap matrix. Additionally, the finite nucleus model is available for the electron-nucleus potential and the vector potential throughout. Efficiency is ensured by the diagonal local approximation to the unitary decoupling transformation in X2C as well as the (multipole-accelerated) resolution of the identity approximation for the Coulomb term (MARI-J, RI-J) and the seminumerical exchange approximation. Errors introduced by these approximations are assessed and found to be clearly negligible. The applicability of our implementation to large-scale calculations is demonstrated for a tin pincer-type system as well as low-valent tin and lead complexes. Here, the calculation of the Sn nuclear magnetic resonance shifts for the pincer-type ligand with about 2400 basis functions requires less than 1 h for hybrid density functionals. Further, the impact of spin-orbit coupling on the nucleus-independent chemical shifts and the corresponding ring currents of all-metal aromatic systems is studied.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| |
Collapse
|
6
|
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: 37] [Impact Index Per Article: 37.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
|
7
|
Franzke YJ. Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy. J Chem Theory Comput 2023; 19:2010-2028. [PMID: 36939092 DOI: 10.1021/acs.jctc.2c01248] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2023]
Abstract
The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| |
Collapse
|
8
|
Krätschmer F, Sun X, Gillhuber S, Kucher H, Franzke YJ, Weigend F, Roesky PW. Fully Tin-Coated Coinage Metal Ions: A Pincer-Type Bis-stannylene Ligand for Exclusive Tetrahedral Complexation. Chemistry 2022; 29:e202203583. [PMID: 36533713 DOI: 10.1002/chem.202203583] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 12/16/2022] [Accepted: 12/19/2022] [Indexed: 12/23/2022]
Abstract
The synthesis of a novel bis-stannylene pincer ligand and its complexation with coinage metals (CuI , AgI and AuI ) are described. All coinage metal centres are in tetrahedral coordination environments in the solid state and are exclusively coordinated by four neutral SnII donors. 119 Sn NMR provided information about the behaviour in solution. All of the isolated compounds have photoluminescent properties, and these were investigated at low and elevated temperatures. Compared to the free bis-stannylene ligand, coordination to coinage metals led to an increase in the luminescence intensity. The new compounds were investigated in detail through all-electron relativistic density functional theory (DFT) calculations.
Collapse
Affiliation(s)
- Frederic Krätschmer
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstr. 15, 76131, Karlsruhe, Germany
| | - Xiaofei Sun
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstr. 15, 76131, Karlsruhe, Germany
| | - Sebastian Gillhuber
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstr. 15, 76131, Karlsruhe, Germany
| | - Hannes Kucher
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstr. 15, 76131, Karlsruhe, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032, Marburg, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032, Marburg, Germany
| | - Peter W Roesky
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstr. 15, 76131, Karlsruhe, Germany
| |
Collapse
|
9
|
Bruder F, Franzke YJ, Weigend F. Paramagnetic NMR Shielding Tensors Based on Scalar Exact Two-Component and Spin-Orbit Perturbation Theory. J Phys Chem A 2022; 126:5050-5069. [PMID: 35857421 DOI: 10.1021/acs.jpca.2c03579] [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/30/2022]
Abstract
The temperature-dependent Fermi-contact and pseudocontact terms are important contributions to the paramagnetic NMR shielding tensor. Herein, we augment the scalar-relativistic (local) exact two-component (X2C) framework with spin-orbit perturbation theory including the screened nuclear spin-orbit correction for the EPR hyperfine coupling and g tensor to compute these temperature-dependent terms. The accuracy of this perturbative ansatz is assessed with the self-consistent spin-orbit two-component and four-component treatments serving as reference. This shows that the Fermi-contact and pseudocontact interaction is sufficiently described for paramagnetic NMR shifts; however, larger deviations are found for the EPR spectra and the principle components of the EPR properties of heavy elements. The impact of the perturbative treatment is further compared to that of the density functional approximation and the basis set. Large-scale calculations are routinely possible with the multipole-accelerated resolution of the identity approximation and the seminumerical exchange approximation, as shown for [CeTi6O3(OiPr)9(salicylate)6].
Collapse
Affiliation(s)
- Florian Bruder
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| |
Collapse
|
10
|
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
|
11
|
Holzer C, Franzke YJ. A Local Hybrid Exchange Functional Approximation from First Principles. J Chem Phys 2022; 157:034108. [DOI: 10.1063/5.0100439] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Local hybrid functionals are a more flexible class of density functional approximations allowing for a position-dependent admixture of exact exchange. This additional flexibility, however, comes with a more involved mathematical form and a more complicated design. A common denominator for previously constructed local hybrid funtionals is usage of thermochemical benchmark data to construct these functionals. Herein, we design a local hybrid functional without relying on benchmark data. Instead, we construct it in a more ab initio manner, following the principles of modern meta-generalized gradient approximations and considering theoretical constrains. To achieve this, we make use of the density matrix expansion and a local mixing function based on an approximate correlation length. The accuracy of the developed density functional approximation is assessed for thermochemistry, excitation energies, polarizabilities, magnetizabilities, NMR spin-spincoupling constants, NMR shieldings and shifts, as well as EPR g-tensors and hyperfine coupling constants. Here, the new exchange functional shows a robust performance and is especially well suited for atomization energies, barrier heights, excitation energies, NMR coupling constants, and EPR properties, whereas it looses some ground for the NMR shifts.Therefore, the designed functional is a major step forwards for functionals that have been designed from first principles.
Collapse
Affiliation(s)
- Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruher Institut für Technologie Fakultät für Physik, Germany
| | - Yannick J. Franzke
- Fachbereich Chemie, Philipps-Universität Marburg Fachbereich Chemie, Germany
| |
Collapse
|
12
|
Franzke YJ, Holzer C. Communication: Impact of the current density on paramagnetic NMR properties. J Chem Phys 2022; 157:031102. [DOI: 10.1063/5.0103898] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Meta-generalized gradient approximations (meta-GGAs) and local hybrid functionals generally depend on the kinetic energy density τ. For magnetic properties, this necessitates generalizations to ensure gauge invariance. In most implementations, τ is generalized by incorporating the external magnetic field. However, this introduces artifacts in the response of the density matrix and does not satisfy the iso-orbital constraint. Here, we extend previous approaches based on the current density to paramagnetic NMR shieldings and EPR g-tensors. The impact is assessed for main-group compounds and transition-metal complexes considering 25 density functional approximations. It is shown that the current density leads to substantial improvements-especially for the popular Minnesota and SCAN functional families. Thus, we strongly recommend to use the current density generalized τ in paramagnetic NMR and EPR calculations with meta-GGAs.
Collapse
Affiliation(s)
- Yannick J. Franzke
- Fachbereich Chemie, Philipps-Universität Marburg Fachbereich Chemie, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruher Institut für Technologie Fakultät für Physik, Germany
| |
Collapse
|
13
|
Pausch A, Holzer C, Klopper W. Efficient Calculation of Magnetic Circular Dichroism Spectra Using Spin-Noncollinear Linear-Response Time-Dependent Density Functional Theory in Finite Magnetic Fields. J Chem Theory Comput 2022; 18:3747-3758. [PMID: 35576504 DOI: 10.1021/acs.jctc.2c00232] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Excited-state calculations in finite magnetic fields are presented in the framework of spin-noncollinear linear-response time-dependent density functional theory. To ensure gauge-origin invariance, London atomic orbitals are employed throughout. An efficient implementation into the Turbomole package, which also includes the resolution of the identity approximation, allows for the investigation of excited states of large molecular systems. The implementation is used to investigate the magnetic circular dichroism spectra of sizable organometallic molecules such as a zinc tetraazaporphyrin with two fused naphthalene units, which is a molecule with 57 atoms.
Collapse
Affiliation(s)
- Ansgar Pausch
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Strasse 1, 76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany.,Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| |
Collapse
|
14
|
Franzke YJ, Yu JM. Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance. J Chem Theory Comput 2022; 18:2246-2266. [PMID: 35354319 DOI: 10.1021/acs.jctc.1c01175] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C5Me5)2Y(μ-S)2Mo(μ-S)2Y(C5Me5)2]-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C6Cl5)4]- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Jason M Yu
- Department of Chemistry, University of California─Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
15
|
Franzke YJ, Holzer C, Mack F. NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation. J Chem Theory Comput 2022; 18:1030-1045. [PMID: 34981925 DOI: 10.1021/acs.jctc.1c00999] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
16
|
Schattenberg CJ, Lehmann M, Bühl M, Kaupp M. Systematic Evaluation of Modern Density Functional Methods for the Computation of NMR Shifts of 3d Transition-Metal Nuclei. J Chem Theory Comput 2021; 18:273-292. [PMID: 34968062 DOI: 10.1021/acs.jctc.1c00964] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A wide range of density functionals from all rungs of Jacob's ladder have been evaluated systematically for a set of experimental 3d transition-metal NMR shifts of 70 complexes encompassing 12 × 49Ti, 10 × 51V, 10 × 53Cr, 11 × 55Mn, 9 × 57Fe, 9 × 59Co, and 9 × 61Ni shift values, as well as a diverse range of electronic structure characteristics. The overall 39 functionals evaluated include one LDA, eight GGAs, seven meta-GGAs (including their current-density-functional─CDFT─versions), nine global hybrids, four range-separated hybrids, eight local hybrids, and two double hybrids, and we also include Hartree-Fock and MP2 calculations. While recent evaluations of the same functionals for a very large coupled-cluster-based benchmark of main-group shieldings and shifts achieved in some cases aggregate percentage mean absolute errors clearly below 2%, the best results for the present 3d-nuclei set are in the range between 4 and 5%. Strikingly, the overall best-performing functionals are the recently implemented CDFT versions of two meta-GGAs, namely cM06-L (4.0%) and cVSXC (4.3%), followed by cLH14t-calPBE (4.9%), B3LYP (5.0%), and cLH07t-SVWN (5.1%), i.e., the previously best-performing global hybrid and two local hybrids. A number of further functionals achieve aggregate deviations in the range 5-6%. Range-separated hybrids offer no particular advantage over global hybrids. Due to the overall poor performance of Hartree-Fock theory for all systems except the titanium complexes, MP2 and double-hybrid functionals are unsuitable for these 3d-nucleus shifts and provide large errors. Global hybrid functionals with larger EXX admixtures, such as BHLYP or M06-2X, also perform poorly, and some other highly parametrized global hybrids also are unsuitable. For many functionals depending on local kinetic energy τ, their CDFT variants perform much better than their "non-CDFT" versions. This holds notably also for the above-mentioned M06-L and VSXC, while the effect is small for τ-dependent local hybrids and can even be somewhat detrimental to the agreement with experiment for a few other cases. The separation between well-performing and more poorly performing functionals is mainly determined by their results for the most critical nuclei 55Mn, 57Fe, and 59Co. Here either moderate exact-exchange admixtures or CDFT versions of meta-GGAs are beneficial for the accuracy. The overall deviations of the better-performing global or local hybrids are then typically dominated by the 53Cr shifts, where triplet instabilities appear to disfavor exact-exchange admixture. Further detailed analyses help to pinpoint specific nuclei and specific types of complexes that are challenges for a given functional.
Collapse
Affiliation(s)
- Caspar Jonas Schattenberg
- Institut für Chemie, Theoretische Chemie/Quantenchemie Sekretariat C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Morten Lehmann
- Institut für Chemie, Theoretische Chemie/Quantenchemie Sekretariat C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Michael Bühl
- School of Chemistry, University of St. Andrews, Purdie Building, North Haugh, St Andrews KY16 9ST, Fife, U.K
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie Sekretariat C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623, Berlin, Germany
| |
Collapse
|
17
|
McGrady JE, Weigend F, Dehnen S. Electronic structure and bonding in endohedral Zintl clusters. Chem Soc Rev 2021; 51:628-649. [PMID: 34931207 DOI: 10.1039/d1cs00775k] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Endohedral Zintl clusters-multi-metallic anionic molecules in which a d-block or f-block metal atom is enclosed by p-block (semi)metal atoms-are very topical in contemporary inorganic chemistry. Not only do they provide insight into the embryonic states of intermetallic compounds and show promise in catalytic applications, they also shed light on the nature of chemical bonding between metal atoms. Over the past two decades, a plethora of endohedral Zintl clusters have been synthesized, revealing a fascinating diversity of molecular architectures. Many different perspectives on the bonding in them have emerged in the literature, sometimes complementary and sometimes conflicting, and there has been no concerted effort to classify the entire family based on a small number of unifying principles. A closer look, however, reveals distinct patterns in structure and bonding that reflect the extent to which valence electrons are shared between the endohedral atom and the cluster shell. We show that there is a much more uniform relationship between the total valence electron count and the structure and bonding patterns of these clusters than previously anticipated. All of the p-block (semi)metal shells can be placed on a ladder of total valence electron count that ranges between 4n+2 (closo deltahedra), 5n (closed, three-bonded polyhedra) and 6n (crown-like structures). Although some structural isomerism can occur for a given electron count, the presence of a central metal cation imposes a preference for rather regular and approximately spherical structures which maximise electrostatic interactions between the metal and the shell. In cases where the endohedral metal has relatively accessible valence electrons (from the d or f shells), it can also contribute its valence electrons to the total electron count of the cluster shell, raising the effective electron count and often altering the structural preferences. The electronic situation in any given cluster is considered from different perspectives, some more physical and some more chemical, in a way that highlights the important point that, in the end, they explain the same situation. This article provides a unifying perspective of bonding that captures the structural diversity across this diverse family of multimetallic clusters.
Collapse
Affiliation(s)
- John E McGrady
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, OX1 3QZ, UK.
| | - Florian Weigend
- Fachbereich Chemie and Wissenschaftliches Zentrum für Materialwissenschaften, Philipps University Marburg, Hans-Meerwein-Straße 4, 35043 Marburg, Germany.
| | - Stefanie Dehnen
- Fachbereich Chemie and Wissenschaftliches Zentrum für Materialwissenschaften, Philipps University Marburg, Hans-Meerwein-Straße 4, 35043 Marburg, Germany.
| |
Collapse
|
18
|
Franzke YJ, Yu JM. Hyperfine Coupling Constants in Local Exact Two-Component Theory. J Chem Theory Comput 2021; 18:323-343. [PMID: 34928142 DOI: 10.1021/acs.jctc.1c01027] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within the one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, that is, the derivatives of the unitary decoupling transformation are considered. This requires the solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened nuclear spin-orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12-17 transition-metal compounds. The error introduced by DLU is negligible, and the DLU-X2C Hamiltonian accurately reproduces its four-component "fully" relativistic parent results. Functionals with a large amount of Hartree-Fock exchange such as CAM-QTP-02 and ωB97X-D are generally favorable. The pure density functional r2SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-ζ bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C6Cl5)4]- with more than 4700 primitive basis functions and four rare-earth single-molecule magnets: [La(OAr*)3]-, [Lu(NR2)3]-, [Lu(OAr*)3]-, and [TbPc2]-. Here, the results with the spin-orbit DLU-X2C Hamiltonian are in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
19
|
Gillhuber S, Franzke YJ, Weigend F. Paramagnetic NMR Shielding Tensors and Ring Currents: Efficient Implementation and Application to Heavy Element Compounds. J Phys Chem A 2021; 125:9707-9723. [PMID: 34723533 DOI: 10.1021/acs.jpca.1c07793] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present an efficient implementation of paramagnetic NMR shielding tensors and shifts in a nonrelativistic and scalar-relativistic density functional theory framework. For the latter, we make use of the scalar exact two-component Hamiltonian in its local approximation, and generally we apply the well established (multipole-accelerated) resolution of the identity approximation and the seminumerical exchange approximation. The perturbed density matrix of a paramagnetic NMR shielding calculation is further used to study the magnetically induced current density and ring currents of open-shell systems as illustrated for [U@Bi12]3-. [U@Bi12]3- features delocalized highest occupied molecular orbitals and sustains a net diatropic ring current of ca. 18 nA/T through the Bi12 torus similar to the all-metal aromatic heavy-element cluster [Th@Bi12]4-.
Collapse
Affiliation(s)
- Sebastian Gillhuber
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| |
Collapse
|
20
|
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
|
21
|
Peters B, Stuhrmann G, Mack F, Weigend F, Dehnen S. Highly Soluble Supertetrahedra upon Selective Partial Butylation of Chalcogenido Metalate Clusters in Ionic Liquids. Angew Chem Int Ed Engl 2021; 60:17622-17628. [PMID: 33974339 PMCID: PMC8362043 DOI: 10.1002/anie.202104867] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 05/10/2021] [Indexed: 11/07/2022]
Abstract
Supertetrahedral clusters have been reported in two generally different types so far: one type possessing an organic ligand shell, no or low charges, and high solubility, while the other cluster type is ligand-free with usually high charges and low or no solubility in common solvents. The latter is a tremendous disadvantage regarding further use of the clusters in solution. However, as organic substituents usually broaden the HOMO-LUMO gaps, which cannot be overcompensated by the (limited) cluster sizes, a full organic shielding comes along with drawbacks regarding opto-electronic properties. We therefore sought to find a way of generating soluble clusters with a minimum number of organic substituents. Here, we present the synthesis and full characterization of two salts of [Sn10 O4 S16 (SBu)4 ]4- that are high soluble in CH2 Cl2 or CH3 CN, which includes first NMR and mass spectra obtained from solutions of such salts with mostly inorganic supertetrahedral clusters. The optical absorption properties of this new class of compounds indicates nearly unaffected band gaps. The synthetic approach and the spectroscopic findings were rationalized and explained by means of high-level quantum chemical studies.
Collapse
Affiliation(s)
- Bertram Peters
- Fachbereich Chemie and Wissenschaftliches Zentrum für MaterialwissenschaftenPhilipps-Universität MarburgHans-Meerwein-Straße 435043MarburgGermany
| | - Gina Stuhrmann
- Fachbereich Chemie and Wissenschaftliches Zentrum für MaterialwissenschaftenPhilipps-Universität MarburgHans-Meerwein-Straße 435043MarburgGermany
| | - Fabian Mack
- Institut für Physikalische ChemieKarlsruher Institut für Technologie (KIT)Kaiserstr. 1276131KarlsruheGermany
| | - Florian Weigend
- Fachbereich Chemie and Wissenschaftliches Zentrum für MaterialwissenschaftenPhilipps-Universität MarburgHans-Meerwein-Straße 435043MarburgGermany
| | - Stefanie Dehnen
- Fachbereich Chemie and Wissenschaftliches Zentrum für MaterialwissenschaftenPhilipps-Universität MarburgHans-Meerwein-Straße 435043MarburgGermany
| |
Collapse
|
22
|
Holzer C, Franzke YJ, Kehry M. Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties. J Chem Theory Comput 2021; 17:2928-2947. [PMID: 33914504 DOI: 10.1021/acs.jctc.1c00203] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.
Collapse
Affiliation(s)
- Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany.,Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Max Kehry
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|