1
|
Lee J, Bertels LW, Small DW, Head-Gordon M. Kohn-Sham Density Functional Theory with Complex, Spin-Restricted Orbitals: Accessing a New Class of Densities without the Symmetry Dilemma. PHYSICAL REVIEW LETTERS 2019; 123:113001. [PMID: 31573235 DOI: 10.1103/physrevlett.123.113001] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2019] [Indexed: 06/10/2023]
Abstract
We show that using complex, spin-restricted orbitals in Kohn-Sham (KS) density functional theory allows one to access a new class of densities that is not accessible by either spin-restricted (RKS) or spin-unrestricted (UKS) orbitals. We further show that the real part of a complex RKS (CRKS) density matrix can be nonidempotent when the imaginary part of the density matrix is not zero. Using CRKS orbitals shows significant improvements in the triplet-singlet gaps of a benchmark set, called TS12, for well-established, widely used density functionals. Moreover, it was shown that RKS and UKS yield qualitatively wrong charge densities and spin densities, respectively, leading to worse energetics. We demonstrate that representative modern density functionals show surprisingly no improvement even with a qualitatively more accurate density from CRKS orbitals. To this end, our work not only provides a way to escape the symmetry dilemma whenever there exists a CRKS solution, but also suggests a new route to design better approximate density functionals.
Collapse
Affiliation(s)
- Joonho Lee
- College of Chemistry, University of California, Berkeley, California 94720, USA
| | - Luke W Bertels
- College of Chemistry, University of California, Berkeley, California 94720, USA
| | - David W Small
- College of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- College of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
2
|
Li Manni G, Carlson RK, Luo S, Ma D, Olsen J, Truhlar DG, Gagliardi L. Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2015; 10:3669-80. [PMID: 26588512 DOI: 10.1021/ct500483t] [Citation(s) in RCA: 299] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a new theoretical framework, called Multiconfiguration Pair-Density Functional Theory (MC-PDFT), which combines multiconfigurational wave functions with a generalization of density functional theory (DFT). A multiconfigurational self-consistent-field (MCSCF) wave function with correct spin and space symmetry is used to compute the total electronic density, its gradient, the on-top pair density, and the kinetic and Coulomb contributions to the total electronic energy. We then use a functional of the total density, its gradient, and the on-top pair density to calculate the remaining part of the energy, which we call the on-top-density-functional energy in contrast to the exchange-correlation energy of Kohn-Sham DFT. Because the on-top pair density is an element of the two-particle density matrix, this goes beyond the Hohenberg-Kohn theorem that refers only to the one-particle density. To illustrate the theory, we obtain first approximations to the required new type of density functionals by translating conventional density functionals of the spin densities using a simple prescription, and we perform post-SCF density functional calculations using the total density, density gradient, and on-top pair density from the MCSCF calculations. Double counting of dynamic correlation or exchange does not occur because the MCSCF energy is not used. The theory is illustrated by applications to the bond energies and potential energy curves of H2, N2, F2, CaO, Cr2, and NiCl and the electronic excitation energies of Be, C, N, N(+), O, O(+), Sc(+), Mn, Co, Mo, Ru, N2, HCHO, C4H6, c-C5H6, and pyrazine. The method presented has a computational cost and scaling similar to MCSCF, but a quantitative accuracy, even with the present first approximations to the new types of density functionals, that is comparable to much more expensive multireference perturbation theory methods.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Rebecca K Carlson
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Sijie Luo
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Dongxia Ma
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University , Langelandsgade 140, DK-8000, Aarhus C, Denmark
| | - Donald G Truhlar
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota , Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota , Minneapolis, Minnesota 55455, United States
| |
Collapse
|
5
|
Seth M, Ziegler T. Calculation of excitation energies of open-shell molecules with spatially degenerate ground states. II. Transformed reference via intermediate configuration Kohn-Sham time dependent density functional theory oscillator strengths and magnetic circular dichroism C terms. J Chem Phys 2006; 124:144105. [PMID: 16626178 DOI: 10.1063/1.2187004] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An extension of the transformed reference via an intermediate configuration Kohn-Sham time dependent density functional theory (TRICKS-TDDFT) method for calculating the transition energies of molecules with spatially degenerate ground states is proposed that enables oscillator strengths to also be evaluated. The oscillator strengths are calculated starting from a description of the degenerate ground state and the excited states of interest in terms of linear combinations of Slater determinants based upon the F-vectors obtained in the TRICKS-TDDFT calculation. This approach for calculating oscillator strengths can also be applied to several other properties that involve excited states. An example of such a property, the C term of magnetic circular dichroism (MCD), is presented. The new method is illustrated through example calculations of the absorbance and MCD spectra of several octahedral and C(4v) d(5) transition metal complexes. The calculated transition energies and oscillator strengths are somewhat too low and too high, respectively, but these errors are within the range that would normally be expected for TDDFT calculations of this type. The calculations help to resolve some previously unexplained problems with the MCD spectra of the C(4v) complexes.
Collapse
Affiliation(s)
- Michael Seth
- Department of Chemistry, University of Calgary, University Drive 2500, Calgary, Alberta T2N-1N4, Canada
| | | |
Collapse
|
6
|
Abstract
This essay provides a perspective on several issues in valence bond theory: the physical significance of semilocal bonding orbitals, the capability of valence bond concepts to explain systems with multireferences character, the use of valence bond theory to provide analytical representations of potential energy surfaces for chemical dynamics by the method of semiempirical valence bond potential energy surfaces (an early example of specific reaction parameters), by multiconfiguration molecular mechanics, by the combined valence bond-molecular mechanics method, and by the use of valence bond states as coupled diabatic states for describing electronically nonadiabatic processes (photochemistry). The essay includes both ab initio and semiempirical approaches.
Collapse
Affiliation(s)
- Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, USA.
| |
Collapse
|
8
|
Seth M, Ziegler T. Calculation of excitation energies of open-shell molecules with spatially degenerate ground states. I. Transformed reference via an intermediate configuration Kohn-Sham density-functional theory and applications to d1 and d2 systems with octahedral and tetrahedral symmetries. J Chem Phys 2005; 123:144105. [PMID: 16238372 DOI: 10.1063/1.2047553] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method for calculating the UV-vis spectra of molecules with spatially degenerate ground states using time-dependent density-functional theory (TDDFT) is proposed. The new transformed reference via an intermediate configuration Kohn-Sham TDDFT (TRICKS-TDDFT) method avoids the difficulties caused by the multireference nature of spatially degenerate states by rather than utilizing the ground state instead taking a nondegenerate excited state with desirable properties as the reference for the TDDFT calculation. The scope and practical application of the method are discussed. Like all open-shell TDDFT calculations this method at times suffers from the inability to produce transitions to states that are eigenfunctions of the total spin operator. A technique for alleviating this difficulty to some extent is proposed. The applicability and accuracy of the TRICKS-TDDFT method is demonstrated through example calculations of several d(1) and d(2) transition metal complexes with tetrahedral and octahedral symmetries. For the most part, the results of these calculations are similar in quality to to those obtained from standard TDDFT calculations.
Collapse
Affiliation(s)
- Michael Seth
- Department of Chemistry, University of Calgary, Calgary, AB T2N-1N4, Canada
| | | |
Collapse
|
10
|
Srinivas S, Jellinek J. Structural and electronic properties of small beryllium clusters: A theoretical study. J Chem Phys 2004; 121:7243-52. [PMID: 15473792 DOI: 10.1063/1.1791071] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
Geometric structures and electronic properties of small beryllium clusters (Be(n), 2< or = n< or =9) are investigated within the gradient-corrected density functional theory. The computations are performed with the Becke exchange and Perdew-Wang correlation functionals. Both low and high multiplicity states are considered. A predominance of higher multiplicity states among the low-energy isomers of the larger clusters is found. An analysis of the variations in the structural and electronic properties with cluster size is presented, and the results are compared with those of earlier studies.
Collapse
Affiliation(s)
- Sudha Srinivas
- Department of Physics, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | | |
Collapse
|