1
|
Sheldon C, Paier J, Usvyat D, Sauer J. Hybrid RPA:DFT Approach for Adsorption on Transition Metal Surfaces: Methane and Ethane on Platinum (111). J Chem Theory Comput 2024; 20:2219-2227. [PMID: 38330551 PMCID: PMC10938501 DOI: 10.1021/acs.jctc.3c01308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 01/17/2024] [Accepted: 01/18/2024] [Indexed: 02/10/2024]
Abstract
The hybrid QM:QM approach is extended to adsorption on transition metal surfaces. The random phase approximation (RPA) as the high-level method is applied to cluster models and, using the subtractive scheme, embedded in periodic models which are treated with density functional theory (DFT) that is the low-level method. The PBE functional, both without dispersion and augmented with the many-body dispersion (MBD), is employed. Adsorption of methane and ethane on the Pt(111) surface is studied. For methane in a 2 × 2 surface cell, the hybrid RPA:PBE and RPA:PBE+MBD results, -14.3 and -16.0 kJ mol-1, respectively, are in close agreement with the periodic RPA value of -13.8 kJ mol-1 at significantly reduced computational cost (factor of ∼50). For methane and ethane, the RPA:PBE results (-14.3 and -17.8 kJ mol-1, respectively) indicate underbinding relative to energies derived from experimental desorption barriers for relevant loadings (-15.6 ± 1.6 and -27.2 ± 2.9 kJ mol-1, respectively), whereas the hybrid RPA:PBE+MBD results (-16.0 and -24.9 kJ mol-1, respectively) agree with the experiment well within experimental uncertainty limits (deviation of -0.4 ± 1.5 and +2.3 ± 2.9 kJ mol-1, respectively). Finding a cluster that adequately and robustly represents the adsorbate at the bulk surface is important for the success of the RPA-based QM:QM scheme for metals.
Collapse
Affiliation(s)
- Christopher Sheldon
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4, Berlin 14195, Germany
| | - Joachim Paier
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
- Lehrstuhl
für Theoretische Chemie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Egerlandstrasse 3, Erlangen 91058, Germany
| | - Denis Usvyat
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
| | - Joachim Sauer
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
| |
Collapse
|
2
|
Hellgren M, Baguet L. Strengths and limitations of the adiabatic exact-exchange kernel for total energy calculations. J Chem Phys 2023; 158:2889488. [PMID: 37158324 DOI: 10.1063/5.0146423] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 04/24/2023] [Indexed: 05/10/2023] Open
Abstract
We investigate the adiabatic approximation to the exact-exchange kernel for calculating correlation energies within the adiabatic-connection fluctuation-dissipation framework of time-dependent density functional theory. A numerical study is performed on a set of systems having bonds of different character (H2 and N2 molecules, H-chain, H2-dimer, solid-Ar, and the H2O-dimer). We find that the adiabatic kernel can be sufficient in strongly bound covalent systems, yielding similar bond lengths and binding energies. However, for non-covalent systems, the adiabatic kernel introduces significant errors around equilibrium geometry, systematically overestimating the interaction energy. The origin of this behavior is investigated by studying a model dimer composed of one-dimensional, closed-shell atoms, interacting via soft-Coulomb potentials. The kernel is shown to exhibit a strong frequency dependence at small to intermediate atomic separation that affects both the low-energy spectrum and the exchange-correlation hole obtained from the corresponding diagonal of the two-particle density matrix.
Collapse
Affiliation(s)
- Maria Hellgren
- Sorbonne Université, MNHN, UMR CNRS 7590, IMPMC, 4 place Jussieu, 75005 Paris, France
| | - Lucas Baguet
- CEA, DAM, DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière en Conditions Extrêmes, 91680 Bruyères-le-Châtel, France
| |
Collapse
|
3
|
Neiss C, Fauser S, Görling A. Geometries and vibrational frequencies with Kohn-Sham methods using σ-functionals for the correlation energy. J Chem Phys 2023; 158:044107. [PMID: 36725500 DOI: 10.1063/5.0129524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022] Open
Abstract
Recently, Kohn-Sham (KS) methods with new correlation functionals, called σ-functionals, have been introduced. Technically, σ-functionals are closely related to the well-known random phase approximation (RPA); formally, σ-functionals are rooted in perturbation theory along the adiabatic connection. If employed in a post-self-consistent field manner in a Gaussian basis set framework, then, σ-functional methods are computationally very efficient. Moreover, for main group chemistry, σ-functionals are highly accurate and can compete with high-level wave-function methods. For reaction and transition state energies, e.g., chemical accuracy of 1 kcal/mol is reached. Here, we show how to calculate first derivatives of the total energy with respect to nuclear coordinates for methods using σ-functionals and then carry out geometry optimizations for test sets of main group molecules, transition metal compounds, and non-covalently bonded systems. For main group molecules, we additionally calculate vibrational frequencies. σ-Functional methods are found to yield very accurate geometries and vibrational frequencies for main group molecules superior not only to those from conventional KS methods but also to those from RPA methods. For geometries of transition metal compounds, not surprisingly, best geometries are found for RPA methods, while σ-functional methods yield somewhat less good results. This is attributed to the fact that in the optimization of σ-functionals, transition metal compounds could not be represented well due to the lack of reliable reference data. For non-covalently bonded systems, σ-functionals yield geometries of the same quality as the RPA or as conventional KS schemes combined with dispersion corrections.
Collapse
Affiliation(s)
- Christian Neiss
- Lehrstuhl für Theoretische Chemie, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Steffen Fauser
- Lehrstuhl für Theoretische Chemie, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
4
|
Drontschenko V, Graf D, Laqua H, Ochsenfeld C. Efficient Method for the Computation of Frozen-Core Nuclear Gradients within the Random Phase Approximation. J Chem Theory Comput 2022; 18:7359-7372. [PMID: 36331398 DOI: 10.1021/acs.jctc.2c00774] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
A method for the evaluation of analytical frozen-core gradients within the random phase approximation is presented. We outline an efficient way to evaluate the response of the density of active electrons arising only when introducing the frozen-core approximation and constituting the main difficulty, together with the response of the standard Kohn-Sham density. The general framework allows to extend the outlined procedure to related electron correlation methods in the atomic orbital basis that require the evaluation of density responses, such as second-order Møller-Plesset perturbation theory or coupled cluster variants. By using Cholesky decomposed densities─which reintroduce the occupied index in the time-determining steps─we are able to achieve speedups of 20-30% (depending on the size of the basis set) by using the frozen-core approximation, which is of similar magnitude as for molecular orbital formulations. We further show that the errors introduced by the frozen-core approximation are practically insignificant for molecular geometries.
Collapse
Affiliation(s)
- Viktoria Drontschenko
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany
| | - Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany
| | - Henryk Laqua
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany.,Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| |
Collapse
|
5
|
Förster A. Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies. J Chem Theory Comput 2022; 18:5948-5965. [PMID: 36150190 PMCID: PMC9558381 DOI: 10.1021/acs.jctc.2c00366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
With increasing interelectronic distance, the screening
of the
electron–electron interaction by the presence of other electrons
becomes the dominant source of electron correlation. This effect is
described by the random phase approximation (RPA) which is therefore
a promising method for the calculation of weak interactions. The success
of the RPA relies on the cancellation of errors, which can be traced
back to the violation of the crossing symmetry of the 4-point vertex,
leading to strongly overestimated total correlation energies. By the
addition of second-order screened exchange (SOSEX) to the correlation
energy, this issue is substantially reduced. In the adiabatic connection
(AC) SOSEX formalism, one of the two electron–electron interaction
lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A
related SOSEX expression in which both electron–electron interaction
lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the
evaluation of this correlation energy expression does not require
an expensive numerical frequency integration and is therefore advantageous
from a computational perspective. We compare the accuracy of the statically
screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical
reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for
charged excitations and noncovalent interactions where they lead to
major improvements over RPA.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands
| |
Collapse
|
6
|
Erhard J, Fauser S, Trushin E, Görling A. Scaled σ-functionals for the Kohn-Sham correlation energy with scaling functions from the homogeneous electron gas. J Chem Phys 2022; 157:114105. [PMID: 36137780 DOI: 10.1063/5.0101641] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. σ-Functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct random phase approximation (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigenvalues σ of the Kohn-Sham response matrix that enters the coupling constant and frequency integration in the dRPA is replaced by another function optimized with the help of reference sets of atomization, reaction, transition state, and non-covalent interaction energies. σ-Functionals are highly accurate and yield chemical accuracy of 1 kcal/mol in reaction or transition state energies, in main group chemistry. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron number, such as ionizations or electron attachments. This problem is attributed to unphysical self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are determined with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addition yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts.
Collapse
Affiliation(s)
- Jannis Erhard
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Steffen Fauser
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
7
|
Śmiga S, Della Sala F, Gori-Giorgi P, Fabiano E. Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit. J Chem Theory Comput 2022; 18:5936-5947. [PMID: 36094908 PMCID: PMC9558377 DOI: 10.1021/acs.jctc.2c00352] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Adiabatic connection
models (ACMs), which interpolate between the
limits of weak and strong interaction, are powerful tools to build
accurate exchange–correlation functionals. If the exact weak-interaction
expansion from the second-order perturbation theory is included, a
self-consistent implementation of these functionals is challenging
and still absent in the literature. In this work, we fill this gap
by presenting a fully self-consistent-field (SCF) implementation of
some popular ACM functionals. While using second-order perturbation
theory at weak interactions, we have also introduced new generalized
gradient approximations (GGAs), beyond the usual point-charge-plus-continuum
model, for the first two leading terms at strong interactions, which
are crucial to ensure robustness and reliability. We then assess the
SCF–ACM functionals for molecular systems and for prototypical
strong-correlation problems. We find that they perform well for both
the total energy and the electronic density and that the impact of
SCF orbitals is directly connected to the accuracy of the ACM functional
form. For the H2 dissociation, the SCF–ACM functionals
yield significant improvements with respect to standard functionals
also thanks to the use of the new GGAs for the strong-coupling functionals.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziądzka 5, 87-100 Toruń, Poland
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| | - Paola Gori-Giorgi
- Department of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| |
Collapse
|
8
|
Everhart LM, Derteano JA, Bates JE. Tension between predicting accurate ground state correlation energies and excitation energies from adiabatic approximations in TDDFT. J Chem Phys 2022; 156:084116. [DOI: 10.1063/5.0080382] [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
The connection between the adiabatic excitation energy of time-dependent density functional theory and the ground state correlation energy from the adiabatic connection fluctuation–dissipation theorem (ACFDT) is explored in the limiting case of one excited state. An exact expression is derived for any adiabatic Hartree-exchange–correlation kernel that connects the excitation energy and the potential contribution to correlation. The resulting formula is applied to the asymmetric Hubbard dimer, a system where this limit is exact. Results from a hierarchy of approximations to the kernel, including the random phase approximation (RPA) with and without exchange and the adiabatically exact (AE) approximation, are compared to the exact ones. At full coupling, the numerical results indicate a tension between predicting an accurate excitation energy and an accurate potential contribution to correlation. The AE approximation is capable of making accurate predictions of both quantities, but only in parts of the parameter space that classify as weakly correlated, while RPA tends to be unable to accurately predict these properties simultaneously anywhere. For a strongly correlated dimer, the AE approximation greatly overestimates the excitation energy yet continues to yield an accurate ground state correlation energy due to its accurate prediction of the adiabatic connection integrand. If similar trends hold for real systems, the development of correlation kernels will be important for applications of the ACFDT in systems with large potential contributions to correlation.
Collapse
Affiliation(s)
- Lucas M. Everhart
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| | - Julio A. Derteano
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| | - Jefferson E. Bates
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| |
Collapse
|
9
|
Siecińska S, Śmiga S, Grabowski I, Della Sala F, Fabiano E. Boosting the OEP2-sc method with spin-component scaling. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2037771] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Sylwia Siecińska
- Faculty of Physics, Astronomy and Informatics, Institute of Physics, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Szymon Śmiga
- Faculty of Physics, Astronomy and Informatics, Institute of Physics, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Ireneusz Grabowski
- Faculty of Physics, Astronomy and Informatics, Institute of Physics, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Italy
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Arnesano, Italy
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Italy
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Arnesano, Italy
| |
Collapse
|
10
|
Vieira D. Readdressing molecular dissociation within the Kohn–Sham formalism of density-functional theory: simple models and a different point of view. Mol Phys 2021. [DOI: 10.1080/00268976.2021.2008037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Daniel Vieira
- Departamento de Física, Programa de Pós-Graduação em Física, Universidade do Estado de Santa Catarina, Joinville, SC, Brazil
| |
Collapse
|
11
|
Di Sabatino S, Loos PF, Romaniello P. Scrutinizing GW-Based Methods Using the Hubbard Dimer. Front Chem 2021; 9:751054. [PMID: 34778206 PMCID: PMC8586429 DOI: 10.3389/fchem.2021.751054] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Accepted: 09/28/2021] [Indexed: 12/03/2022] Open
Abstract
Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.
Collapse
Affiliation(s)
- S. Di Sabatino
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS and ETSF, Toulouse, France
| | - P.-F. Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| | - P. Romaniello
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS and ETSF, Toulouse, France
| |
Collapse
|
12
|
Fauser S, Trushin E, Neiss C, Görling A. Chemical accuracy with σ-functionals for the Kohn-Sham correlation energy optimized for different input orbitals and eigenvalues. J Chem Phys 2021; 155:134111. [PMID: 34624971 DOI: 10.1063/5.0059641] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recently, a new type of orbital-dependent functional for the Kohn-Sham (KS) correlation energy, σ-functionals, was introduced. Technically, σ-functionals are closely related to the well-known direct random phase approximation (dRPA). Within the dRPA, a function of the eigenvalues σ of the frequency-dependent KS response function is integrated over purely imaginary frequencies. In σ-functionals, this function is replaced by one that is optimized with respect to reference sets of atomization, reaction, transition state, and non-covalent interaction energies. The previously introduced σ-functional uses input orbitals and eigenvalues from KS calculations with the generalized gradient approximation (GGA) exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). Here, σ-functionals using input orbitals and eigenvalues from the meta-GGA TPSS and the hybrid-functionals PBE0 and B3LYP are presented and tested. The number of reference sets taken into account in the optimization of the σ-functionals is larger than in the first PBE based σ-functional and includes sets with 3d-transition metal compounds. Therefore, also a reparameterized PBE based σ-functional is introduced. The σ-functionals based on PBE0 and B3LYP orbitals and eigenvalues reach chemical accuracy for main group chemistry. For the 10 966 reactions from the highly accurate W4-11RE reference set, the B3LYP based σ-functional exhibits a mean average deviation of 1.03 kcal/mol compared to 1.08 kcal/mol for the coupled cluster singles doubles perturbative triples method if the same valence quadruple zeta basis set is used. For 3d-transition metal chemistry, accuracies of about 2 kcal/mol are reached. The computational effort for the post-self-consistent evaluation of the σ-functional is lower than that of a preceding PBE0 or B3LYP calculation for typical systems.
Collapse
Affiliation(s)
- Steffen Fauser
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91 058 Erlangen, Germany
| | - Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91 058 Erlangen, Germany
| | - Christian Neiss
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91 058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91 058 Erlangen, Germany
| |
Collapse
|
13
|
Burton HGA, Marut C, Daas TJ, Gori-Giorgi P, Loos PF. Variations of the Hartree-Fock fractional-spin error for one electron. J Chem Phys 2021; 155:054107. [PMID: 34364354 DOI: 10.1063/5.0056968] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Fractional-spin errors are inherent in all current approximate density functionals, including Hartree-Fock theory, and their origin has been related to strong static correlation effects. The conventional way to encode fractional-spin calculations is to construct an ensemble density that scales between the high-spin and low-spin densities. In this article, we explore the variation of the Hartree-Fock fractional-spin (or ghost-interaction) error in one-electron systems using restricted and unrestricted ensemble densities and the exact generalized Hartree-Fock representation. By considering the hydrogen atom and H+ 2 cation, we analyze how the unrestricted and generalized Hartree-Fock schemes minimize this error by localizing the electrons or rotating the spin coordinates. We also reveal a clear similarity between the Coulomb hole of He-like ions and the density depletion near the nucleus induced by the fractional-spin error in the unpolarized hydrogen atom. Finally, we analyze the effect of the fractional-spin error on the Møller-Plesset adiabatic connection, excited states, and functional- and density-driven errors.
Collapse
Affiliation(s)
- Hugh G A Burton
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Clotilde Marut
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Timothy J Daas
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
14
|
Ruan S, Ren X, Gould T, Ruzsinszky A. Self-Interaction-Corrected Random Phase Approximation. J Chem Theory Comput 2021; 17:2107-2115. [PMID: 33689324 DOI: 10.1021/acs.jctc.0c01079] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The short-range correlation energy of the random phase approximation (RPA) is too negative and is often corrected by local or nonlocal methods. These beyond-RPA corrections usually lead to a mixed performance for thermodynamics and dissociation properties. RPA+ is an additive correction based on density functional approximations that often gives realistic total energies for atoms or solids. RPA+ adds a moderate correction to the ionization energies/electron affinities of RPA but does not yield an improvement beyond RPA for atomization energies of molecules. This incompleteness results in severely underestimated atomization energies just like in RPA. Exchange-correlation kernels within the Dyson equation could simultaneously improve atomization, ionization energies, and electron affinities, but their implementation is computationally less feasible in localized basis set codes. In preceding work ( Phys. Rev. A 100, 2019022515), two of the authors proposed a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities. gRPA+ was found to yield a large improvement of ionization energies and electron affinities of light atoms over RPA, and a smaller improvement over RPA+. Within this work, we investigate to what extent this improvement transfers to atomization energies, ionization energies, and electron affinities of molecules, using a modified gRPA+ (mgRPA+) method that can be applied in codes with localized basis functions. We thereby aim to understand the applicability of beyond-RPA corrections based on density functional approximations.
Collapse
Affiliation(s)
- Shiqi Ruan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Xinguo Ren
- Institute of Physics, Chinese Academy of Sciences, 3rd South Street 8, Beijing 100190, China
| | - Tim Gould
- Qld Micro- and Nanotechnology Center, Griffith University, Nathan, Qld 4111, Australia
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
15
|
Trushin E, Thierbach A, Görling A. Toward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energy. J Chem Phys 2021; 154:014104. [PMID: 33412877 DOI: 10.1063/5.0026849] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce new functionals for the Kohn-Sham correlation energy that are based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem and are named σ-functionals. Like in the well-established direct random phase approximation (dRPA), σ-functionals require as input exclusively eigenvalues σ of the frequency-dependent KS response function. In the new functionals, functions of σ replace the σ-dependent dRPA expression in the coupling-constant and frequency integrations contained in the ACFD theorem. We optimize σ-functionals with the help of reference sets for atomization, reaction, transition state, and non-covalent interaction energies. The optimized functionals are to be used in a post-self-consistent way using orbitals and eigenvalues from conventional Kohn-Sham calculations employing the exchange-correlation functional of Perdew, Burke, and Ernzerhof. The accuracy of the presented approach is much higher than that of dRPA methods and is comparable to that of high-level wave function methods. Reaction and transition state energies from σ-functionals exhibit accuracies close to 1 kcal/mol and thus approach chemical accuracy. For the 10 966 reactions of the W4-11RE reference set, the mean absolute deviation is 1.25 kcal/mol compared to 3.21 kcal/mol in the dRPA case. Non-covalent binding energies are accurate to a few tenths of a kcal/mol. The presented approach is highly efficient, and the post-self-consistent calculation of the total energy requires less computational time than a density-functional calculation with a hybrid functional and thus can be easily carried out routinely. σ-Functionals can be implemented in any existing dRPA code with negligible programming effort.
Collapse
Affiliation(s)
- Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Adrian Thierbach
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
16
|
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
|
17
|
Erhard J, Fauser S, Kalaß S, Moerman E, Trushin E, Görling A. Lieb-Oxford bound and pair correlation functions for density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem. Faraday Discuss 2020; 224:79-97. [PMID: 32935700 DOI: 10.1039/d0fd00047g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Compliance with the Lieb-Oxford bound for the indirect Coulomb energy and for the exchange-correlation energy is investigated for a number of density-functional methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem to treat correlation. Furthermore, the correlation contribution to the pair density resulting from these methods is compared with highly accurate reference values for the helium atom and for the hydrogen molecule at several bond distances. For molecules, the Lieb-Oxford bound is obeyed by all considered methods. For the homogeneous electron gas, it is violated by all methods for low electron densities. The simplest considered ACFD method, the direct random phase approximation (dRPA), violates the Lieb-Oxford bound much earlier than more advanced ACFD methods that, in addition to the simple Hartree kernel, take into account the exchange kernel and an approximate correlation kernel in the calculation of the correlation energy. While the dRPA yields quite poor correlation contributions to the pair density, those from more advanced ACFD methods are physically reasonable but still leave room for improvements, particularly in the case of the stretched hydrogen molecule.
Collapse
Affiliation(s)
- Jannis Erhard
- Lehrstuhl für Theoretische Chemie, Egerlandstraße 3, 91058 Erlangen, Germany.
| | | | | | | | | | | |
Collapse
|
18
|
Thierbach A, Görling A. Analytic energy gradients for the self-consistent direct random phase approximation. J Chem Phys 2020; 153:134113. [DOI: 10.1063/5.0021809] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Adrian Thierbach
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
19
|
Thierbach A, Görling A. Analytic energy gradients for the exact exchange Kohn–Sham method. J Chem Phys 2020; 152:114113. [DOI: 10.1063/1.5142711] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Adrian Thierbach
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
20
|
Mezei PD, Kállay M. Construction of a Range-Separated Dual-Hybrid Direct Random Phase Approximation. J Chem Theory Comput 2019; 15:6678-6687. [PMID: 31693355 DOI: 10.1021/acs.jctc.9b00891] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Blending the good performance of the global hybrid PBE0 functional at short-range and the dual-hybrid dRPA75 functional at long range, we propose a new range-separated direct random phase approximation (dRPA75rs), which considerably improves on the accuracy of the calculated reaction energies and barrier heights compared to the parent approaches and provides a good description of noncovalent interactions without any dispersion correction. We also combine the new scheme with spin-component scaling (SCS-dRPA75rs), which enables the accurate calculation of energy differences for processes involving electron pair breaking, such as atomization. The new method scaling as the fourth power of the system size shows a balanced performance on a broad test set involving radicals, transition metal atoms, and heavy atoms, which makes it competitive with the best double-hybrid functionals based on the second-order perturbation theory. According to the results for the homogeneous electron gas, our dRPA75rs method expectedly gives errors for metallic systems similar to the dRPA approach with an additional error cancellation in the case of partial spin polarization.
Collapse
Affiliation(s)
- Pál D Mezei
- Department of Chemistry , University of Basel , Basel 4056 , Switzerland
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science , Budapest University of Technology and Economics , P.O. Box 91, H-1521 Budapest , Hungary
| |
Collapse
|
21
|
Mezei PD, Ruzsinszky A, Kállay M. Reducing the Many-Electron Self-Interaction Error in the Second-Order Screened Exchange Method. J Chem Theory Comput 2019; 15:6607-6616. [DOI: 10.1021/acs.jctc.9b00672] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pál D. Mezei
- Department of Chemistry, University of Basel, Basel 4056, Switzerland
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
22
|
Chehaibou B, Badawi M, Bučko T, Bazhirov T, Rocca D. Computing RPA Adsorption Enthalpies by Machine Learning Thermodynamic Perturbation Theory. J Chem Theory Comput 2019; 15:6333-6342. [DOI: 10.1021/acs.jctc.9b00782] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Bilal Chehaibou
- Université de Lorraine, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
| | - Michael Badawi
- Université de Lorraine, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
| | - Tomáš Bučko
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Mlynská Dolina, Ilkovičova 6, SK-84215 Bratislava, Slovakia
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská cesta 9, SK-84236 Bratislava, Slovakia
| | - Timur Bazhirov
- Exabyte Inc., San Francisco, California 94103, United States
| | - Dario Rocca
- Université de Lorraine, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
| |
Collapse
|
23
|
Thierbach A, Schmidtel D, Görling A. Robust and accurate hybrid random-phase-approximation methods. J Chem Phys 2019; 151:144117. [DOI: 10.1063/1.5120587] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrian Thierbach
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Daniel Schmidtel
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
24
|
Mordovina U, Reinhard TE, Theophilou I, Appel H, Rubio A. Self-Consistent Density-Functional Embedding: A Novel Approach for Density-Functional Approximations. J Chem Theory Comput 2019; 15:5209-5220. [PMID: 31490684 PMCID: PMC6785802 DOI: 10.1021/acs.jctc.9b00063] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2019] [Indexed: 11/29/2022]
Abstract
In the present work, we introduce a self-consistent density-functional embedding technique, which leaves the realm of standard energy-functional approaches in density functional theory and targets directly the density-to-potential mapping that lies at its heart. Inspired by the density matrix embedding theory, we project the full system onto a set of small interacting fragments that can be solved accurately. Based on the rigorous relation of density and potential in density functional theory, we then invert the fragment densities to local potentials. Combining these results in a continuous manner provides an update for the Kohn-Sham potential of the full system, which is then used to update the projection. We benchmark our approach for molecular bond stretching in one and two dimensions and show that, in these cases, the scheme converges to accurate approximations for densities and Kohn-Sham potentials. We demonstrate that the known steps and peaks of the exact exchange-correlation potential are reproduced by our method with remarkable accuracy.
Collapse
Affiliation(s)
- Uliana Mordovina
- Max
Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Teresa E. Reinhard
- Max
Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Iris Theophilou
- Max
Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Heiko Appel
- Max
Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Angel Rubio
- Max
Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
- Center
for Computational Quantum Physics (CCQ), Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United
States
| |
Collapse
|
25
|
Zhang IY, Xu X. Simultaneous Attenuation of Both Self-Interaction Error and Nondynamic Correlation Error in Density Functional Theory: A Spin-Pair Distinctive Adiabatic-Connection Approximation. J Phys Chem Lett 2019; 10:2617-2623. [PMID: 31046289 DOI: 10.1021/acs.jpclett.9b00946] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We present a spin-pair distinctive algorithm in the context of adiabatic-connection fluctuation-dissipation (ACFD) theorem, which enables to quantify the self-interaction error (SIE) and the nondynamic/strong correlation error (NCE) in the direct random-phase approximation (dRPA). Using this knowledge, we propose a spin-component scaled dRPA (scsRPA) correlation model with simultaneous attenuation of both the SIE and the NCE. Along with the exact exchange, scsRPA is shown to present a comprehensive improvement over dRPA, as well as the well-established PBE and PBE0 functionals, for bonding energies of pronounced multireference characters and transition-metal complexes of strongly correlated systems, while consistently providing an accurate description for reaction energies, reaction barriers, and noncovalent bond interactions of weakly correlated systems.
Collapse
Affiliation(s)
- Igor Ying Zhang
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Collaborative Innovation Center of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry , Fudan University , Shanghai 200433 , China
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Collaborative Innovation Center of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry , Fudan University , Shanghai 200433 , China
| |
Collapse
|
26
|
Assessment of a range-separated orbital-optimised random-phase approximation electron correlation method. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2363-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
27
|
Bates JE, Sengupta N, Sensenig J, Ruzsinszky A. Adiabatic Connection without Coupling Constant Integration. J Chem Theory Comput 2018; 14:2979-2990. [DOI: 10.1021/acs.jctc.8b00067] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jefferson E. Bates
- Department of Chemistry, Appalachian State University, Boone, North Carolina 28607, United States
| | - Niladri Sengupta
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Jonathon Sensenig
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
28
|
Intermolecular dispersion energies from coupled exact-exchange Kohn-Sham excitation energies and vectors. COMPUT THEOR CHEM 2018. [DOI: 10.1016/j.comptc.2018.02.019] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
29
|
Huang C, Chi YC. Directly patching high-level exchange-correlation potential based on fully determined optimized effective potentials. J Chem Phys 2017; 147:244111. [PMID: 29289130 DOI: 10.1063/1.5003663] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The key element in Kohn-Sham (KS) density functional theory is the exchange-correlation (XC) potential. We recently proposed the exchange-correlation potential patching (XCPP) method with the aim of directly constructing high-level XC potential in a large system by patching the locally computed, high-level XC potentials throughout the system. In this work, we investigate the patching of the exact exchange (EXX) and the random phase approximation (RPA) correlation potentials. A major challenge of XCPP is that a cluster's XC potential, obtained by solving the optimized effective potential equation, is only determined up to an unknown constant. Without fully determining the clusters' XC potentials, the patched system's XC potential is "uneven" in the real space and may cause non-physical results. Here, we developed a simple method to determine this unknown constant. The performance of XCPP-RPA is investigated on three one-dimensional systems: H20, H10Li8, and the stretching of the H19-H bond. We investigated two definitions of EXX: (i) the definition based on the adiabatic connection and fluctuation dissipation theorem (ACFDT) and (ii) the Hartree-Fock (HF) definition. With ACFDT-type EXX, effective error cancellations were observed between the patched EXX and the patched RPA correlation potentials. Such error cancellations were absent for the HF-type EXX, which was attributed to the fact that for systems with fractional occupation numbers, the integral of the HF-type EXX hole is not -1. The KS spectra and band gaps from XCPP agree reasonably well with the benchmarks as we make the clusters large.
Collapse
Affiliation(s)
- Chen Huang
- Department of Scientific Computing, Florida State University, Tallahassee, Florida 32306-4120, USA
| | - Yu-Chieh Chi
- Department of Scientific Computing, Florida State University, Tallahassee, Florida 32306-4120, USA
| |
Collapse
|
30
|
Dixit A, Claudot J, Lebègue S, Rocca D. Improving the Efficiency of Beyond-RPA Methods within the Dielectric Matrix Formulation: Algorithms and Applications to the A24 and S22 Test Sets. J Chem Theory Comput 2017; 13:5432-5442. [DOI: 10.1021/acs.jctc.7b00837] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Anant Dixit
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandœuvre-lès-Nancy, France
| | - Julien Claudot
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandœuvre-lès-Nancy, France
| | - Sébastien Lebègue
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandœuvre-lès-Nancy, France
- CNRS, CRM2, UMR 7036, 54506 Vandœuvre-lès-Nancy, France
| | - Dario Rocca
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandœuvre-lès-Nancy, France
- CNRS, CRM2, UMR 7036, 54506 Vandœuvre-lès-Nancy, France
| |
Collapse
|
31
|
Heßelmann A. Low scaling random-phase approximation electron correlation method including exchange interactions using localised orbitals. J Chem Phys 2017; 146:174110. [DOI: 10.1063/1.4981817] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andreas Heßelmann
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
32
|
Bates JE, Mezei PD, Csonka GI, Sun J, Ruzsinszky A. Reference Determinant Dependence of the Random Phase Approximation in 3d Transition Metal Chemistry. J Chem Theory Comput 2016; 13:100-109. [DOI: 10.1021/acs.jctc.6b00900] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- J. E. Bates
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - P. D. Mezei
- Department
of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, H-1521 Budapest, Hungary
| | - G. I. Csonka
- Department
of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, H-1521 Budapest, Hungary
| | - J. Sun
- Department
of Physics, University of Texas El Paso, El Paso, Texas 79968, United States
| | - A. Ruzsinszky
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
33
|
Erhard J, Bleiziffer P, Görling A. Power Series Approximation for the Correlation Kernel Leading to Kohn-Sham Methods Combining Accuracy, Computational Efficiency, and General Applicability. PHYSICAL REVIEW LETTERS 2016; 117:143002. [PMID: 27740821 DOI: 10.1103/physrevlett.117.143002] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/24/2016] [Indexed: 05/27/2023]
Abstract
A power series approximation for the correlation kernel of time-dependent density-functional theory is presented. Using this approximation in the adiabatic-connection fluctuation-dissipation (ACFD) theorem leads to a new family of Kohn-Sham methods. The new methods yield reaction energies and barriers of unprecedented accuracy and enable a treatment of static (strong) correlation with an accuracy of high-level multireference configuration interaction methods but are single-reference methods allowing for a black-box-like handling of static correlation. The new methods exhibit a better scaling of the computational effort with the system size than rivaling wave-function-based electronic structure methods. Moreover, the new methods do not suffer from the problem of singularities in response functions plaguing previous ACFD methods and therefore are applicable to any type of electronic system.
Collapse
Affiliation(s)
- Jannis Erhard
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstrasse 3, D-91058 Erlangen, Germany
| | - Patrick Bleiziffer
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstrasse 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstrasse 3, D-91058 Erlangen, Germany
| |
Collapse
|
34
|
Zhang IY, Rinke P, Perdew JP, Scheffler M. Towards Efficient Orbital-Dependent Density Functionals for Weak and Strong Correlation. PHYSICAL REVIEW LETTERS 2016; 117:133002. [PMID: 27715089 DOI: 10.1103/physrevlett.117.133002] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/15/2016] [Indexed: 06/06/2023]
Abstract
We present a new paradigm for the design of exchange-correlation functionals in density-functional theory. Electron pairs are correlated explicitly by means of the recently developed second order Bethe-Goldstone equation (BGE2) approach. Here we propose a screened BGE2 (sBGE2) variant that efficiently regulates the coupling of a given electron pair. sBGE2 correctly dissociates H_{2} and H_{2}^{+}, a problem that has been regarded as a great challenge in density-functional theory for a long time. The sBGE2 functional is then taken as a building block for an orbital-dependent functional, termed ZRPS, which is a natural extension of the PBE0 hybrid functional. While worsening the good performance of sBGE2 in H_{2} and H_{2}^{+}, ZRPS yields a remarkable and consistent improvement over other density functionals across various chemical environments from weak to strong correlation.
Collapse
Affiliation(s)
- Igor Ying Zhang
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Patrick Rinke
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
- Department of Applied Physics, Aalto University, P.O. Box 11100, Aalto FI-00076, Finland
| | - John P Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Matthias Scheffler
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
- Department of Chemistry and Biochemistry and Materials Department, University of California-Santa Barbara, Santa Barbara, California 93106-5050, USA
| |
Collapse
|
35
|
Dixit A, Ángyán JG, Rocca D. Improving the accuracy of ground-state correlation energies within a plane-wave basis set: The electron-hole exchange kernel. J Chem Phys 2016; 145:104105. [DOI: 10.1063/1.4962352] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Anant Dixit
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandoeuvre-lès-Nancy, France
| | - János G. Ángyán
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, CRM2, UMR 7036, 54506 Vandoeuvre-lès-Nancy, France
| | - Dario Rocca
- Université de Lorraine, CRM2, UMR 7036, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, CRM2, UMR 7036, 54506 Vandoeuvre-lès-Nancy, France
| |
Collapse
|
36
|
Pykal M, Jurečka P, Karlický F, Otyepka M. Modelling of graphene functionalization. Phys Chem Chem Phys 2016; 18:6351-72. [DOI: 10.1039/c5cp03599f] [Citation(s) in RCA: 168] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
This perspective describes the available theoretical methods and models for simulating graphene functionalization based on quantum and classical mechanics.
Collapse
Affiliation(s)
- Martin Pykal
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| | - František Karlický
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| |
Collapse
|
37
|
Pavlíková Přecechtělová J, Bahmann H, Kaupp M, Ernzerhof M. Design of exchange-correlation functionals through the correlation factor approach. J Chem Phys 2015; 143:144102. [DOI: 10.1063/1.4932074] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|
38
|
Patrick CE, Thygesen KS. Adiabatic-connection fluctuation-dissipation DFT for the structural properties of solids—The renormalized ALDA and electron gas kernels. J Chem Phys 2015; 143:102802. [DOI: 10.1063/1.4919236] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Christopher E. Patrick
- Center for Atomic-Scale Materials Design (CAMD), Department of Physics, Technical University of Denmark, DK—2800 Kongens Lyngby, Denmark
| | - Kristian S. Thygesen
- Center for Atomic-Scale Materials Design (CAMD), Department of Physics, Technical University of Denmark, DK—2800 Kongens Lyngby, Denmark
| |
Collapse
|
39
|
Bleiziffer P, Krug M, Görling A. Self-consistent Kohn-Sham method based on the adiabatic-connection fluctuation-dissipation theorem and the exact-exchange kernel. J Chem Phys 2015; 142:244108. [DOI: 10.1063/1.4922517] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Patrick Bleiziffer
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Marcel Krug
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
40
|
Kállay M. Linear-scaling implementation of the direct random-phase approximation. J Chem Phys 2015; 142:204105. [DOI: 10.1063/1.4921542] [Citation(s) in RCA: 78] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
van Aggelen H, Yang Y, Yang W. Exchange-correlation energy from pairing matrix fluctuation and the particle-particle random phase approximation. J Chem Phys 2015; 140:18A511. [PMID: 24832319 DOI: 10.1063/1.4865816] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Despite their unmatched success for many applications, commonly used local, semi-local, and hybrid density functionals still face challenges when it comes to describing long-range interactions, static correlation, and electron delocalization. Density functionals of both the occupied and virtual orbitals are able to address these problems. The particle-hole (ph-) Random Phase Approximation (RPA), a functional of occupied and virtual orbitals, has recently known a revival within the density functional theory community. Following up on an idea introduced in our recent communication [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013)], we formulate more general adiabatic connections for the correlation energy in terms of pairing matrix fluctuations described by the particle-particle (pp-) propagator. With numerical examples of the pp-RPA, the lowest-order approximation to the pp-propagator, we illustrate the potential of density functional approximations based on pairing matrix fluctuations. The pp-RPA is size-extensive, self-interaction free, fully anti-symmetric, describes the strong static correlation limit in H2, and eliminates delocalization errors in H2(+) and other single-bond systems. It gives surprisingly good non-bonded interaction energies--competitive with the ph-RPA--with the correct R(-6) asymptotic decay as a function of the separation R, which we argue is mainly attributable to its correct second-order energy term. While the pp-RPA tends to underestimate absolute correlation energies, it gives good relative energies: much better atomization energies than the ph-RPA, as it has no tendency to underbind, and reaction energies of similar quality. The adiabatic connection in terms of pairing matrix fluctuation paves the way for promising new density functional approximations.
Collapse
Affiliation(s)
- Helen van Aggelen
- Department of Inorganic and Physical Chemistry, Ghent University, Ghent, Belgium
| | - Yang Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry and Department of Physics, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
42
|
Heßelmann A. Molecular Excitation Energies from Time-Dependent Density Functional Theory Employing Random-Phase Approximation Hessians with Exact Exchange. J Chem Theory Comput 2015; 11:1607-20. [DOI: 10.1021/acs.jctc.5b00024] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andreas Heßelmann
- Lehrstuhl für Theoretische
Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
43
|
Hollett JW. Non-pairwise additivity of the leading-order dispersion energy. J Chem Phys 2015; 142:084105. [PMID: 25725710 DOI: 10.1063/1.4908134] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The leading-order (i.e., dipole-dipole) dispersion energy is calculated for one-dimensional (1D) and two-dimensional (2D) infinite lattices, and an infinite 1D array of infinitely long lines, of doubly occupied locally harmonic wells. The dispersion energy is decomposed into pairwise and non-pairwise additive components. By varying the force constant and separation of the wells, the non-pairwise additive contribution to the dispersion energy is shown to depend on the overlap of density between neighboring wells. As well separation is increased, the non-pairwise additivity of the dispersion energy decays. The different rates of decay for 1D and 2D lattices of wells is explained in terms of a Jacobian effect that influences the number of nearest neighbors. For an array of infinitely long lines of wells spaced 5 bohrs apart, and an inter-well spacing of 3 bohrs within a line, the non-pairwise additive component of the leading-order dispersion energy is -0.11 kJ mol(-1) well(-1), which is 7% of the total. The polarizability of the wells and the density overlap between them are small in comparison to that of the atomic densities that arise from the molecular density partitioning used in post-density-functional theory (DFT) damped dispersion corrections, or DFT-D methods. Therefore, the nonadditivity of the leading-order dispersion observed here is a conservative estimate of that in molecular clusters.
Collapse
Affiliation(s)
- Joshua W Hollett
- Department of Chemistry, University of Winnipeg, Winnipeg, Manitoba R3B 2G3, Canada and Department of Chemistry, University of Manitoba, Winnipeg, Manitoba R3T 2N2, Canada
| |
Collapse
|
44
|
The Ring and Exchange-Ring Approximations Based on Kohn–Sham Reference States. Top Curr Chem (Cham) 2014. [DOI: 10.1007/128_2014_557] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
|
45
|
Bleiziffer P, Schmidtel D, Görling A. Stability conditions for exact-exchange Kohn-Sham methods and their relation to correlation energies from the adiabatic-connection fluctuation-dissipation theorem. J Chem Phys 2014; 141:204107. [DOI: 10.1063/1.4901924] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
46
|
Přecechtělová J, Bahmann H, Kaupp M, Ernzerhof M. Communication: A non-empirical correlation factor model for the exchange-correlation energy. J Chem Phys 2014; 141:111102. [DOI: 10.1063/1.4896057] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jana Přecechtělová
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
- Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Hilke Bahmann
- Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Matthias Ernzerhof
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| |
Collapse
|
47
|
Grabowski I, Fabiano E, Teale AM, Śmiga S, Buksztel A, Sala FD. Orbital-dependent second-order scaled-opposite-spin correlation functionals in the optimized effective potential method. J Chem Phys 2014; 141:024113. [DOI: 10.1063/1.4887097] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Ireneusz Grabowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland
| | - Eduardo Fabiano
- National Nanotechnology Laboratory, Istituto Nanoscienze–CNR, Via per Arnesano, I-73100 Lecce, Italy
| | - Andrew M. Teale
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland
| | - Adam Buksztel
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland
| | - Fabio Della Sala
- National Nanotechnology Laboratory, Istituto Nanoscienze–CNR, Via per Arnesano, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano (LE), Italy
| |
Collapse
|
48
|
Lu D. Evaluation of model exchange-correlation kernels in the adiabatic connection fluctuation-dissipation theorem for inhomogeneous systems. J Chem Phys 2014; 140:18A520. [DOI: 10.1063/1.4867538] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Olsen T, Thygesen KS. Static correlation beyond the random phase approximation: Dissociating H2 with the Bethe-Salpeter equation and time-dependent GW. J Chem Phys 2014; 140:164116. [DOI: 10.1063/1.4871875] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Bates JE, Furche F. Communication: Random phase approximation renormalized many-body perturbation theory. J Chem Phys 2013; 139:171103. [DOI: 10.1063/1.4827254] [Citation(s) in RCA: 81] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
|