1
|
Li H, Mansoori Kermani M, Ottochian A, Crescenzi O, Janesko BG, Truhlar DG, Scalmani G, Frisch MJ, Ciofini I, Adamo C. Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry? J Am Chem Soc 2024; 146:6721-6732. [PMID: 38413362 DOI: 10.1021/jacs.3c12713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree-Fock exchange. These results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.
Collapse
Affiliation(s)
- Hanwei Li
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
| | - Maryam Mansoori Kermani
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Alistar Ottochian
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
| | - Orlando Crescenzi
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di Monte Sant'Angelo, Via Cinthia, Napoli 80126, Italy
| | - Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas 76129, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | | | | | - Ilaria Ciofini
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, Paris F-75005, France
| |
Collapse
|
2
|
Zope RR, Yamamoto Y, Baruah T. How well do one-electron self-interaction-correction methods perform for systems with fractional electrons? J Chem Phys 2024; 160:084102. [PMID: 38385511 DOI: 10.1063/5.0182773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/28/2024] [Indexed: 02/23/2024] Open
Abstract
Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (zσ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew-Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of H2+ and He2+ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li+ to F- at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.
Collapse
Affiliation(s)
- Rajendra R Zope
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Yoh Yamamoto
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
3
|
Shi Y, Shi Y, Wasserman A. Stretching Bonds without Breaking Symmetries in Density Functional Theory. J Phys Chem Lett 2024; 15:826-833. [PMID: 38232318 DOI: 10.1021/acs.jpclett.3c03073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
Kohn-Sham density functional theory (KS-DFT) stands out among electronic structure methods due to its balance of accuracy and computational efficiency. However, to achieve chemically accurate energies, standard density functional approximations in KS-DFT often need to break underlying symmetries, a long-standing "symmetry dilemma". By employing fragment spin densities as the main variables in calculations (rather than total molecular densities, as in KS-DFT), we present an embedding framework in which this symmetry dilemma is understood and partially resolved. The spatial overlap between fragment densities is used as the main ingredient to construct a simple, physically motivated approximation to a universal functional of the fragment densities. This "overlap approximation" is shown to significantly improve semilocal KS-DFT binding energies of molecules without artificially breaking either charge or spin symmetries. The approach is shown to be applicable to covalently bonded molecules and to systems of the "strongly correlated" type.
Collapse
Affiliation(s)
- Yuming Shi
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, United States
| | - Yi Shi
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Adam Wasserman
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, United States
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
4
|
Briccolani-Bandini L, Brémond E, Pagliai M, Cardini G, Ciofini I, Adamo C. Concerted versus stepwise proton transfer reactions in the [2, 2'-bipyridyl]-3-3'-diol molecule: A static and dynamic ab-initio investigation. J Comput Chem 2023; 44:2308-2318. [PMID: 37584183 DOI: 10.1002/jcc.27198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Revised: 07/09/2023] [Accepted: 07/11/2023] [Indexed: 08/17/2023]
Abstract
The double proton transfer (PT) reaction has been investigated in the [2,2'-bipyridyl]-3-3'-diol, a complex molecule where the proton movements is coupled to significant rearrangement of the electronic structure. Moreover, the reaction could be concerted, that is the two protons are exchanged simultaneously, or stepwise, where the two protons are transferred sequentially. To this end, a static exploration of the potential energy surface (PES) was carried together with the analysis of the free-energy surface (FES), both surfaces being evaluated at density functional theory level and different exchange-correlation functionals. While the concerted mechanism has been clearly discharged, the characteristics of the stepwise PT significantly depends on the chosen functionals, some suggesting a clear stepwise mechanism characterized by a stable reaction intermediates and two transitions states, whereas other approaches propend for a asynchronous PT, with a single TS. These features appear on both PES and FES, albeit some differences appears due to their different nature.
Collapse
Affiliation(s)
| | - Eric Brémond
- ITODYS, CNRS, Université Paris Cité, Paris, France
| | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| | - Ilaria Ciofini
- CNRS, i-CLeHS, UMR 8060, Chimie ParisTech, PSL University, Paris, France
| | - Carlo Adamo
- CNRS, i-CLeHS, UMR 8060, Chimie ParisTech, PSL University, Paris, France
- Institut Universitaire de France, Paris, France
| |
Collapse
|
5
|
Huan Lew-Yee JF, Piris M, Del Campo JM. Outstanding improvement in removing the delocalization error by global natural orbital functional. J Chem Phys 2023; 158:084110. [PMID: 36859086 DOI: 10.1063/5.0137378] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
This work assesses the performance of the recently proposed global natural orbital functional (GNOF) against the charge delocalization error. GNOF provides a good balance between static and dynamic electronic correlations leading to accurate total energies while preserving spin, even for systems with a highly multi-configurational character. Several analyses were applied to the functional, namely, (i) how the charge is distributed in super-systems of two fragments, (ii) the stability of ionization potentials while increasing the system size, and (iii) potential energy curves of a neutral and charged diatomic system. GNOF was found to practically eliminate the charge delocalization error in many of the studied systems or greatly improve the results obtained previously with PNOF7.
Collapse
Affiliation(s)
- Juan Felipe Huan Lew-Yee
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City C.P. 04510, Mexico
| | - Mario Piris
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain; Euskal Herriko Unibertsitatea (UPV/EHU), PK 1072, 20080 Donostia, Euskadi, Spain; and Basque Foundation for Science (IKERBASQUE), 48009 Bilbao, Euskadi, Spain
| | - Jorge M Del Campo
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City C.P. 04510, Mexico
| |
Collapse
|
6
|
Lew-Yee JFH, M. del Campo J. Charge delocalization error in Piris Natural Orbital Functionals. J Chem Phys 2022; 157:104113. [DOI: 10.1063/5.0102310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Piris Natural Orbital Functionals (PNOF) have been recognized as a low-scaling alternative to study strong correlated systems. In this work, we address the performance of the fifth functional (PNOF5) and the seventh functional (PNOF7) to deal with another common problem, the charge delocalization error. The effects of this problem can be observed in charged systems of repeated well-separated fragments, where the energy should be the sum of the charged and neutral fragments, regardless of how the charge is distributed. In practice, an energetic overstabilization of fractional charged fragments leads to a preference for having the charge delocalized throughout the system. To establish the performance of PNOF functionals regarding charge delocalization error, charged chains of helium atoms and the W4-17-MR set molecules were used as base fragments and their energy, charge distribution and correlation regime were studied. It was found that PNOF5 prefers localized charge distributions, while PNOF7 improves the treatment of interpair static correlation and tends to the correct energetic limit for several cases, although a preference for delocalized charge distributions may arise in highly strong correlation regimes. Overall, it is concluded that PNOF functionals can simultaneously deal with static correlation and charge delocalization errors, resulting in a promising choice to study charge-related problems.
Collapse
Affiliation(s)
- Juan Felipe Huan Lew-Yee
- Departamento de Física y Química Teórica, Universidad Nacional Autónoma de México Facultad de Química, Mexico
| | - Jorge M. del Campo
- Departamento de Física y Química Teórica, Universidad Nacional Autónoma de México, Mexico
| |
Collapse
|
7
|
Mutual Information in Conjugate Spaces for Neutral Atoms and Ions. ENTROPY 2022; 24:e24020233. [PMID: 35205526 PMCID: PMC8871140 DOI: 10.3390/e24020233] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Revised: 01/18/2022] [Accepted: 01/28/2022] [Indexed: 12/10/2022]
Abstract
The discrepancy among one-electron and two-electron densities for diverse N-electron atomss, enclosing neutral systems (with nuclear charge Z=N) and charge-one ions (|N−Z|=1), is quantified by means of mutual information, I, and Quantum Similarity Index, QSI, in the conjugate spaces position/momentum. These differences can be interpreted as a measure of the electron correlation of the system. The analysis is carried out by considering systems with a nuclear charge up to Z=103 and singly charged ions (cations and anions) as far as N=54. The interelectronic correlation, for any given system, is quantified through the comparison of its double-variable electron pair density and the product of the respective one-particle densities. An in-depth study along the Periodic Table reveals the importance, far beyond the weight of the systems considered, of their shell structure.
Collapse
|
8
|
Zhou C, Hermes MR, Wu D, Bao JJ, Pandharkar R, King DS, Zhang D, Scott TR, Lykhin AO, Gagliardi L, Truhlar DG. Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory. Chem Sci 2022; 13:7685-7706. [PMID: 35865899 PMCID: PMC9261488 DOI: 10.1039/d2sc01022d] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Accepted: 06/01/2022] [Indexed: 12/25/2022] Open
Abstract
Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals. This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.![]()
Collapse
Affiliation(s)
- Chen Zhou
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Matthew R. Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Dihua Wu
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Jie J. Bao
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Riddhish Pandharkar
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Daniel S. King
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Dayou Zhang
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Thais R. Scott
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Aleksandr O. Lykhin
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| |
Collapse
|
9
|
Mosquera MA. Density Functional Calculations Based on the Exponential Ansatz. J Phys Chem A 2021; 125:8751-8763. [PMID: 34582684 DOI: 10.1021/acs.jpca.1c07113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
This work explores the application of the singles-based exponential ansatz to density functional calculations. In contrast to the standard approach where Kohn-Sham (KS) orbitals are determined prior to computing molecular quantities of interest, we consider the single-reference Hartree-Fock wave function as a starting point. Applying the exponential ansatz to this single reference gives an auxiliary wave function that is employed to calculate the electronic properties of the system. This wave function is determined self-consistently through the standard KS Hamiltonian but evaluated over the Hartree-Fock molecular orbital basis. By using spin-symmetry breaking, we recover size-consistent results free of unphysical fractional charges in the dissociation limit. Our method shows consistency with standard KS density functional calculations when the system geometry is similar to the equilibrium one or in repulsive configurations. For moderately long distances between atoms, not at dissociation, because of self-interaction the exponential ansatz may give instabilities in the form of large cluster amplitudes. To avoid these, this work introduces a relatively simple regularization method that preserves size-consistency and penalizes high amplitudes of the cluster operator, whereas the results remain physically meaningful. We also present the time-dependent extension of our theory and show that it can feature quantum states where multiple electrons are excited.
Collapse
Affiliation(s)
- Martín A Mosquera
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, Montana 59718, United States
| |
Collapse
|
10
|
Abstract
Kohn-Sham density functional theory with the available exchange-correlation functionals is less accurate for strongly correlated systems, which require a multiconfigurational description as a zero-order function, than for weakly correlated systems, and available functionals of the spin densities do not accurately predict energies for many strongly correlated systems when one uses multiconfigurational wave functions with spin symmetry. Furthermore, adding a correlation functional to a multiconfigurational reference energy can lead to double counting of electron correlation. Multiconfiguration pair-density functional theory (MC-PDFT) overcomes both obstacles, the second by calculating the quantum mechanical part of the electronic energy entirely by a functional, and the first by using a functional of the total density and the on-top pair density rather than the spin densities. This allows one to calculate the energy of strongly correlated systems efficiently with a pair-density functional and a suitable multiconfigurational reference function. This article reviews MC-PDFT and related background information.
Collapse
Affiliation(s)
- Prachi Sharma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA;
| | - Jie J Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA;
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA;
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, USA;
| |
Collapse
|
11
|
Graf D, Ochsenfeld C. A range-separated generalized Kohn-Sham method including a long-range nonlocal random phase approximation correlation potential. J Chem Phys 2020; 153:244118. [PMID: 33380112 DOI: 10.1063/5.0031310] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Based on our recently published range-separated random phase approximation (RPA) functional [Kreppel et al., "Range-separated density-functional theory in combination with the random phase approximation: An accuracy benchmark," J. Chem. Theory Comput. 16, 2985-2994 (2020)], we introduce self-consistent minimization with respect to the one-particle density matrix. In contrast to the range-separated RPA methods presented so far, the new method includes a long-range nonlocal RPA correlation potential in the orbital optimization process, making it a full-featured variational generalized Kohn-Sham (GKS) method. The new method not only improves upon all other tested RPA schemes including the standard post-GKS range-separated RPA for the investigated test cases covering general main group thermochemistry, kinetics, and noncovalent interactions but also significantly outperforms the popular G0W0 method in estimating the ionization potentials and fundamental gaps considered in this work using the eigenvalue spectra obtained from the GKS Hamiltonian.
Collapse
Affiliation(s)
- Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| |
Collapse
|
12
|
Bao JJ, Zhou C, Truhlar DG. Compressed-State Multistate Pair-Density Functional Theory. J Chem Theory Comput 2020; 16:7444-7452. [PMID: 33141587 DOI: 10.1021/acs.jctc.0c00908] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Multiconfiguration pair-density functional theory (MC-PDFT) is a multireference method that can be used to calculate excited states. However, MC-PDFT potential energy surfaces have the wrong topology at conical intersections because the last step of MC-PDFT is not a diagonalization of a model-space Hamiltonian matrix, as done in, for example, multistate second-order perturbation theory (MS-CASPT2). We have previously proposed methods that solve this problem by diagonalizing a model-space effective Hamiltonian matrix, where the diagonal elements are MC-PDFT energies for intermediate states, and the off-diagonal elements are evaluated by wave function theory. One previous method is called variational multistate PDFT (VMS-PDFT), whose intermediate states maximize the trace of the effective Hamiltonian, namely, the sum of the MC-PDFT energies of the model-space states; the VMS-PDFT is very robust but is more computationally expensive than another method, extended multistate PDFT (XMS-PDFT), in which the transformation to intermediate states is accomplished without needing any density functional evaluations. However, although VMS-PDFT was accurate in all cases tested, XMS-PDFT was accurate in only some of them. In the present paper, we propose a new method, called compressed-state multistate PDFT (CMS-PDFT), that is as efficient as XMS-PDFT and as accurate as VMS-PDFT. The new method maximizes the trace of the classical Coulomb energy of the intermediate states such that the electron densities of the intermediate states are compressed. We show that CMS-PDFT performs robustly even where XMS-PDFT fails.
Collapse
Affiliation(s)
- Jie J Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Chen Zhou
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
13
|
Brémond É, Pérez-Jiménez ÁJ, Sancho-García JC, Adamo C. Range-separated hybrid and double-hybrid density functionals: A quest for the determination of the range-separation parameter. J Chem Phys 2020; 152:244124. [PMID: 32610956 DOI: 10.1063/5.0010976] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We recently derived a new and simple route to the determination of the range-separation parameter in range-separated exchange hybrid and double-hybrid density functionals by imposing an additional constraint to the exchange-correlation energy to recover the total energy of the hydrogen atom [Brémond et al., J. Chem. Phys. 15, 201102 (2019)]. Here, we thoroughly assess this choice by statistically comparing the derived values of the range-separation parameters to the ones obtained using the optimal tuning (OT) approach. We show that both approaches closely agree, thus, confirming the reliability of ours. We demonstrate that it provides very close performances in the computation of properties particularly prone to the one- and many-electron self-interaction errors (i.e., ionization potentials). Our approach arises as an alternative to the OT procedure, conserving the accuracy and efficiency of a standard Kohn-Sham approach to density-functional theory computation.
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), FRE 2027, F-75005 Paris, France
| |
Collapse
|
14
|
Mostafanejad M, Liebenthal MD, DePrince AE. Global Hybrid Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2020; 16:2274-2283. [PMID: 32101416 DOI: 10.1021/acs.jctc.9b01178] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A global hybrid extension of multiconfiguration pair-density functional theory (MC-PDFT) is developed. Using a linear decomposition of the electron-electron repulsion term, a fraction λ of the nonlocal exchange interaction, obtained from variational two-electron reduced-density matrix (v2RDM)-driven complete active-space self-consistent field (CASSCF) theory, is combined with its local counterpart, obtained from an on-top pair-density functional. The resulting scheme (called λ-MC-PDFT) inherits the benefits of MC-PDFT (e.g., its simplicity and the resolution of the symmetry dilemma) and, when combined with the v2RDM approach to CASSCF, requires only polynomially scaling computational effort. As a result, λ-MC-PDFT can efficiently describe static and dynamical correlation effects in strongly correlated systems. The efficacy of the approach is assessed for several challenging multiconfigurational problems, including the dissociation of molecular nitrogen, the double dissociation of a water molecule, and the 1,3-dipolar cycloadditions of ozone to ethylene and ozone to acetylene in the O3ADD6 benchmark set.
Collapse
Affiliation(s)
- Mohammad Mostafanejad
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Marcus Dante Liebenthal
- Department of Chemistry and Biochemistry, Ithaca College, Ithaca, New York 14850, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
15
|
Yu LJ, Dale SG, Chan B, Karton A. Benchmark study of DFT and composite methods for bond dissociation energies in argon compounds. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2019.110676] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
16
|
Bao JL, Carter EA. Surface-Plasmon-Induced Ammonia Decomposition on Copper: Excited-State Reaction Pathways Revealed by Embedded Correlated Wavefunction Theory. ACS NANO 2019; 13:9944-9957. [PMID: 31393708 DOI: 10.1021/acsnano.9b05030] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Ammonia is a promising hydrogen storage medium; however, its decomposition via conventional thermal catalysis requires a significant amount of thermal energy input in order to overcome the reaction barriers. Here, we use embedded correlated wavefunction (ECW) theory to quantify reaction pathways and energetics for ammonia decomposition (N-H bond dissociation and N2 and H2 associative desorption) on copper (Cu) nanoparticles using a Cu (111) surface model. We predict that surface plasmon excitations will be able to facilitate ammonia decomposition by substantially reducing the effective barriers along excited-state pathways. We estimate the reductions in reaction barriers for breaking the first N-H bond and for recombinative desorption of surface-bound nitrogen and hydrogen atoms to be approximately 1.7, 0.8, and 0.5 eV, respectively. Further, by using the experimental N2 desorption barrier as a reference, we compare the accuracy of various theoretical methods, including plane-wave Kohn-Sham density functional theory calculations with commonly used exchange-correlation functionals, embedded complete active space second-order perturbation theory, and embedded multiconfiguration pair-density functional theory. This work offers further confirmation that the ECW theoretical framework is the most robust for treating highly correlated local electronic structures of solids.
Collapse
|
17
|
Gaggioli CA, Stoneburner SJ, Cramer CJ, Gagliardi L. Beyond Density Functional Theory: The Multiconfigurational Approach To Model Heterogeneous Catalysis. ACS Catal 2019. [DOI: 10.1021/acscatal.9b01775] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Carlo Alberto Gaggioli
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Samuel J. Stoneburner
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
18
|
Brémond É, Pérez-Jiménez ÁJ, Sancho-García JC, Adamo C. Range-separated hybrid density functionals made simple. J Chem Phys 2019; 150:201102. [PMID: 31153220 DOI: 10.1063/1.5097164] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this communication, we present a new and simple route to derive range-separated exchange (RSX) hybrid and double hybrid density functionals in a nonempirical fashion. In line with our previous developments [Brémond et al., J. Chem. Theory Comput. 14, 4052 (2018)], we show that by imposing an additional physical constraint to the exchange-correlation energy, i.e., by enforcing to reproduce the total energy of the hydrogen atom, we are able to generalize the nonempirical determination of the range-separation parameter to a family of RSX hybrid density functionals. The success of the resulting models is illustrated by an accurate modeling of several molecular systems and properties, like ionization potentials, particularly prone to the one- and many-electron self-interaction errors.
Collapse
Affiliation(s)
- Éric Brémond
- Université Paris Diderot, Sorbonne Paris Cité, ITODYS, UMR CNRS 7086, 15 rue J.-A. de Baïf, F-75013 Paris, France
| | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), FRE 2027, F-75005 Paris, France
| |
Collapse
|
19
|
Bao JL, Verma P, Truhlar DG. How well can density functional theory and pair-density functional theory predict the correct atomic charges for dissociation and accurate dissociation energetics of ionic bonds? Phys Chem Chem Phys 2018; 20:23072-23078. [PMID: 30167614 DOI: 10.1039/c8cp04280b] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The accuracy of density functional theory (DFT) is often judged by predicted dissociation energies, but one should also consider charge densities as illustrated here for dissociation of heteronuclear diatomic molecules, including ionic bonds for which local density functionals yield erroneous results. Some hybrid density functionals with 100% exact exchange in Kohn-Sham DFT and the local functionals in multiconfiguration pair-density functional theory give relatively acurate dissociation energies for NaCl, and they correctly yield uncharged dissociated atoms.
Collapse
Affiliation(s)
- Junwei Lucas Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, USA.
| | | | | |
Collapse
|
20
|
von Szentpály L. Eliminating symmetry problems in electronegativity equalization and correcting self-interaction errors in conceptual DFT. J Comput Chem 2018; 39:1949-1969. [DOI: 10.1002/jcc.25356] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Revised: 04/16/2018] [Accepted: 04/16/2018] [Indexed: 11/09/2022]
Affiliation(s)
- László von Szentpály
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55; Stuttgart D-70569 Germany
| |
Collapse
|
21
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
22
|
Brémond É, Savarese M, Pérez-Jiménez ÁJ, Sancho-García JC, Adamo C. Range-Separated Double-Hybrid Functional from Nonempirical Constraints. J Chem Theory Comput 2018; 14:4052-4062. [PMID: 29923721 DOI: 10.1021/acs.jctc.8b00261] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
On the basis of our previous developments in the field of nonempirical double hybrids, we present here a new exchange-correlation functional based on a range-separated model for the exchange part and integrating a nonlocal perturbative correction to the electron correlation contribution. Named RSX-QIDH, the functional is free from any kind of empirical parametrization. Its range-separation parameter is set to recover the total energy of the hydrogen atom, thus eliminating the self-interaction error for this one-electron system. Subsequent tests on some relevant benchmark data sets confirm that the self-interaction error is particularly low for RSX-QIDH. This new functional provides also correct dissociation profiles for charged rare-gas dimers and very accurate ionization potentials directly from Kohn-Sham orbital energies. Above all, these good results are not obtained at the expense of other properties. Indeed, further tests on standard benchmarks show that RSX-QIDH is competitive with the more empirical ωB97X-2 double hybrid and outperforms the parent LC-PBE long-range corrected hybrid, thus underlining the important role of the nonlocal perturbative correlation.
Collapse
Affiliation(s)
- Éric Brémond
- Université Paris Diderot, Sorbonne Paris Cité, ITODYS, UMR CNRS 7086, 15 rue J.-A. de Baïf , F-75013 Paris , France
| | - Marika Savarese
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , I-16163 Genoa , Italy
| | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, 11, rue Pierre et Marie Curie , F-75005 Paris , France.,Institut Universitaire de France, 103 Boulevard Saint Michel , F-75005 Paris , France
| |
Collapse
|
23
|
Withanage KPK, Trepte K, Peralta JE, Baruah T, Zope R, Jackson KA. On the Question of the Total Energy in the Fermi–Löwdin Orbital Self-Interaction Correction Method. J Chem Theory Comput 2018; 14:4122-4128. [DOI: 10.1021/acs.jctc.8b00344] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kushantha P. K. Withanage
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| | - Kai Trepte
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| | - Juan E. Peralta
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| | - Tunna Baruah
- Physics Department, University of Texas-El Paso, El Paso, Texas 79968, United States
| | - Rajendra Zope
- Physics Department, University of Texas-El Paso, El Paso, Texas 79968, United States
| | - Koblar A. Jackson
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| |
Collapse
|
24
|
Bao JL, Gagliardi L, Truhlar DG. Self-Interaction Error in Density Functional Theory: An Appraisal. J Phys Chem Lett 2018; 9:2353-2358. [PMID: 29624392 DOI: 10.1021/acs.jpclett.8b00242] [Citation(s) in RCA: 90] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Self-interaction error (SIE) is considered to be one of the major sources of error in most approximate exchange-correlation functionals for Kohn-Sham density-functional theory (KS-DFT), and it is large with all local exchange-correlation functionals and with some hybrid functionals. In this work, we consider systems conventionally considered to be dominated by SIE. For these systems, we demonstrate that by using multiconfiguration pair-density functional theory (MC-PDFT), the error of a translated local density-functional approximation is significantly reduced (by a factor of 3) when using an MCSCF density and on-top density, as compared to using KS-DFT with the parent functional; the error in MC-PDFT with local on-top functionals is even lower than the error in some popular KS-DFT hybrid functionals. Density-functional theory, either in MC-PDFT form with local on-top functionals or in KS-DFT form with some functionals having 50% or more nonlocal exchange, has smaller errors for SIE-prone systems than does CASSCF, which has no SIE.
Collapse
Affiliation(s)
- Junwei Lucas Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
25
|
Wang Y, Liu J, Li J, He X. Fragment-based quantum mechanical calculation of protein-protein binding affinities. J Comput Chem 2018; 39:1617-1628. [PMID: 29707784 DOI: 10.1002/jcc.25236] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2017] [Revised: 03/02/2018] [Accepted: 04/01/2018] [Indexed: 12/13/2022]
Abstract
The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method has been successfully utilized for efficient linear-scaling quantum mechanical (QM) calculation of protein energies. In this work, we applied the EE-GMFCC method for calculation of binding affinity of Endonuclease colicin-immunity protein complex. The binding free energy changes between the wild-type and mutants of the complex calculated by EE-GMFCC are in good agreement with experimental results. The correlation coefficient (R) between the predicted binding energy changes and experimental values is 0.906 at the B3LYP/6-31G*-D level, based on the snapshot whose binding affinity is closest to the average result from the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) calculation. The inclusion of the QM effects is important for accurate prediction of protein-protein binding affinities. Moreover, the self-consistent calculation of PB solvation energy is required for accurate calculations of protein-protein binding free energies. This study demonstrates that the EE-GMFCC method is capable of providing reliable prediction of relative binding affinities for protein-protein complexes. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Yaqian Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Jinfeng Liu
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.,Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Jinjin Li
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.,National Engineering Research Centre for Nanotechnology, Shanghai, 200241, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| |
Collapse
|