1
|
Vacondio S, Varsano D, Ruini A, Ferretti A. Going Beyond the GW Approximation Using the Time-Dependent Hartree-Fock Vertex. J Chem Theory Comput 2024; 20:4718-4737. [PMID: 38772396 DOI: 10.1021/acs.jctc.4c00100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2024]
Abstract
The time-dependent Hartree-Fock (TDHF) vertex of many-body perturbation theory (MBPT) makes it possible to extend TDHF theory to charged excitations. Here we assess its performance by applying it to spherical atoms in their neutral electronic configuration. On a theoretical level, we recast the TDHF vertex as a reducible vertex, highlighting the emergence of a self-energy expansion purely in orders of the bare Coulomb interaction; then, on a numerical level, we present results for polarizabilities, ionization energies (IEs), and photoemission satellites. We confirm the superiority of THDF over simpler methods such as the random phase approximation for the prediction of atomic polarizabilities. We then find that the TDHF vertex reliably provides better IEs than GW and low-order self-energies do in the light-atom, few-electron regime; its performance degrades in heavier, many-electron atoms instead, where an expansion in orders of an unscreened Coulomb interaction becomes less justified. New relevant features are introduced in the satellite spectrum by the TDHF vertex, but the experimental spectra are not fully reproduced due to a missing account of nonlinear effects connected to hole relaxation. We also explore various truncations of the self-energy given by the TDHF vertex, but do not find them to be more convenient than low-order approximations such as GW and second Born (2B), suggesting that vertex corrections should be carried out consistently both in the self-energy and in the polarizability.
Collapse
Affiliation(s)
- Simone Vacondio
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università degli Studi di Modena e Reggio Emilia, Via G. Campi 213/a, 41125 Modena, Italy
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| | - Daniele Varsano
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| | - Alice Ruini
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università degli Studi di Modena e Reggio Emilia, Via G. Campi 213/a, 41125 Modena, Italy
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| | - Andrea Ferretti
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| |
Collapse
|
2
|
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
|
3
|
Vacondio S, Varsano D, Ruini A, Ferretti A. Numerically Precise Benchmark of Many-Body Self-Energies on Spherical Atoms. J Chem Theory Comput 2022; 18:3703-3717. [PMID: 35561415 PMCID: PMC9202310 DOI: 10.1021/acs.jctc.2c00048] [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
![]()
We investigate the
performance of beyond-GW approaches in many-body
perturbation theory by addressing atoms described within the spherical
approximation via a dedicated numerical treatment based on B-splines
and spherical harmonics. We consider the GW, second Born (2B), and
GW + second order screened exchange (GW+SOSEX) self-energies and use
them to obtain ionization potentials from the quasi-particle equation
(QPE) solved perturbatively on top of independent-particle calculations.
We also solve the linearized Sham–Schlüter equation
(LSSE) and compare the resulting xc potentials against exact data.
We find that the LSSE provides consistent starting points for the
QPE but does not present any practical advantage in the present context.
Still, the features of the xc potentials obtained with it shed light
on possible strategies for the inclusion of beyond-GW diagrams in
the many-body self-energy. Our findings show that solving the QPE
with the GW+SOSEX self-energy on top of a PBE or PBE0 solution is
a viable scheme to go beyond GW in finite systems, even in the atomic
limit. However, GW shows a comparable performance if one agrees to
use a hybrid starting point. We also obtain promising results with
the 2B self-energy on top of Hartree–Fock, suggesting that
the full time-dependent Hartree–Fock vertex may be another
viable beyond-GW scheme for finite systems.
Collapse
Affiliation(s)
- S Vacondio
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università di Modena e Reggio Emilia, Via G. Campi 213/a, Modena 41121, Italy.,Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| | - D Varsano
- Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| | - A Ruini
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università di Modena e Reggio Emilia, Via G. Campi 213/a, Modena 41121, Italy.,Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| | - A Ferretti
- Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| |
Collapse
|
4
|
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
|
5
|
Yu JM, Nguyen BD, Tsai J, Hernandez DJ, Furche F. Selfconsistent random phase approximation methods. J Chem Phys 2021; 155:040902. [PMID: 34340391 DOI: 10.1063/5.0056565] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This Perspective reviews recent efforts toward selfconsistent calculations of ground-state energies within the random phase approximation (RPA) in the (generalized) Kohn-Sham (KS) density functional theory context. Since the RPA correlation energy explicitly depends on the non-interacting KS potential, an additional condition to determine the energy as a functional of the density is necessary. This observation leads to the concept of functional selfconsistency (FSC), which requires that the KS density equals the interacting density defined as the functional derivative of the ground-state energy with respect to the external potential. While all existing selfconsistent RPA schemes violate FSC, the recent generalized KS semicanonical projected RPA (GKS-spRPA) method takes a step toward satisfying it. This leads to systematic improvements in densities, binding energy curves, reference state stability, and molecular properties compared to non-selfconsistent RPA as well as optimized effective potential RPA. GKS-spRPA orbital energies accurately approximate valence and core ionization potentials, and even electron affinities of non-valence bound anions. The computational cost and performance of GKS-spRPA are compared to those of related selfconsistent schemes, including GW and orbital optimization methods, and limitations are discussed. Large differences between KS and interacting densities observed in the absence of FSC and the well-rounded performance of GKS-spRPA suggest that the KS potential as a density functional should be defined via the FSC condition for explicitly potential-dependent density functionals.
Collapse
Affiliation(s)
- Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Brian D Nguyen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Jeffrey Tsai
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Devin J Hernandez
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
6
|
Riemelmoser S, Kaltak M, Kresse G. Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation. J Chem Phys 2021; 154:154103. [PMID: 33887939 DOI: 10.1063/5.0045400] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have, therefore, been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham bandgaps and dielectric constants, and comment on the issue of self-consistency.
Collapse
Affiliation(s)
- Stefan Riemelmoser
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
| | - Merzuk Kaltak
- VASP Software GmbH, Sensengasse 8/17, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
Jana S, Samal P. Screened hybrid meta-GGA exchange–correlation functionals for extended systems. Phys Chem Chem Phys 2019; 21:3002-3015. [DOI: 10.1039/c8cp06715e] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Screened Hartree–Fock exchange integrated with semilocal exchange–correlation functionals often proficiently predict several solid-state properties.
Collapse
Affiliation(s)
- Subrata Jana
- School of Physical Sciences
- National Institute of Science Education and Research
- HBNI
- Bhubaneswar 752050
- India
| | - Prasanjit Samal
- School of Physical Sciences
- National Institute of Science Education and Research
- HBNI
- Bhubaneswar 752050
- India
| |
Collapse
|
13
|
Tran F, Blaha P. Importance of the Kinetic Energy Density for Band Gap Calculations in Solids with Density Functional Theory. J Phys Chem A 2017; 121:3318-3325. [PMID: 28402113 PMCID: PMC5423078 DOI: 10.1021/acs.jpca.7b02882] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2017] [Indexed: 11/29/2022]
Abstract
Recently, exchange-correlation potentials in density functional theory were developed with the goal of providing improved band gaps in solids. Among them, the semilocal potentials are particularly interesting for large systems since they lead to calculations that are much faster than with hybrid functionals or methods like GW. We present an exhaustive comparison of semilocal exchange-correlation potentials for band gap calculations on a large test set of solids, and particular attention is paid to the potential HLE16 proposed by Verma and Truhlar. It is shown that the most accurate potential is the modified Becke-Johnson potential, which, most noticeably, is much more accurate than all other semilocal potentials for strongly correlated systems. This can be attributed to its additional dependence on the kinetic energy density. It is also shown that the modified Becke-Johnson potential is at least as accurate as the hybrid functionals and more reliable for solids with large band gaps.
Collapse
Affiliation(s)
- Fabien Tran
- Institute of Materials Chemistry, Vienna
University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Peter Blaha
- Institute of Materials Chemistry, Vienna
University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| |
Collapse
|
14
|
Klimeš J. Lattice energies of molecular solids from the random phase approximation with singles corrections. J Chem Phys 2016; 145:094506. [DOI: 10.1063/1.4962188] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jiří Klimeš
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-18223 Prague 8, Czech Republic and Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
15
|
Śmiga S, Della Sala F, Buksztel A, Grabowski I, Fabiano E. Accurate Kohn-Sham ionization potentials from scaled-opposite-spin second-order optimized effective potential methods. J Comput Chem 2016; 37:2081-90. [DOI: 10.1002/jcc.24436] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2016] [Revised: 05/26/2016] [Accepted: 06/04/2016] [Indexed: 01/25/2023]
Affiliation(s)
- Szymon Śmiga
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT); via Arnesano Lecce 73100
- Institute of Physics; Faculty of Physics, Astronomy and Informatics; Nicolaus Copernicus University; Grudziadzka 5 Torun 87-100 Poland
- Center for Biomolecular Nanotechnologies @UNILE; Istituto Italiano Di Tecnologia; via Barsanti Arnesano I-73010 Italy
| | - Fabio Della Sala
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT); via Arnesano Lecce 73100
- Center for Biomolecular Nanotechnologies @UNILE; Istituto Italiano Di Tecnologia; via Barsanti Arnesano I-73010 Italy
| | - Adam Buksztel
- Institute of Physics; Faculty of Physics, Astronomy and Informatics; Nicolaus Copernicus University; Grudziadzka 5 Torun 87-100 Poland
| | - Ireneusz Grabowski
- Institute of Physics; Faculty of Physics, Astronomy and Informatics; Nicolaus Copernicus University; Grudziadzka 5 Torun 87-100 Poland
| | - Eduardo Fabiano
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT); via Arnesano Lecce 73100
- Center for Biomolecular Nanotechnologies @UNILE; Istituto Italiano Di Tecnologia; via Barsanti Arnesano I-73010 Italy
| |
Collapse
|
16
|
Pela RR, Marques M, Teles LK. Comparing LDA-1/2, HSE03, HSE06 and G₀W₀ approaches for band gap calculations of alloys. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:505502. [PMID: 26609566 DOI: 10.1088/0953-8984/27/50/505502] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
It has long been known that the local density approximation and the generalized gradient approximation do not furnish reliable band gaps, and one needs to go beyond these approximations to reliably describe these properties. Among alternatives are the use of hybrid functionals (HSE03 and HSE06 being popular), the GW approximation or the recently proposed LDA-1/2 method. In this work, we compare rigorously the performance of these four methods in describing the band gaps of alloys, employing the generalized quasi-chemical approach to treat the disorder of the alloy and to obtain judiciously the band gap for the entire compositional range. Zincblende InGaAs and InGaN were chosen as prototypes due to their importance in optoelectronic applications. The comparison between these four approaches was guided both by the agreement between the predicted band gap and the experimental one, and by the demanded computational effort (time and memory). We observed that the HSE06 method provided the most accurate results (in comparison with experiments), whereas, surprisingly, the LDA-1/2 method gave the best compromise between accuracy and computational resources. Due to its low computational cost and good accuracy, we decided to double the supercell used to describe the alloys, and employing LDA-1/2 we observed that the bowing parameter changed remarkably, only agreeing with the measured one for the larger supercell, where LDA-1/2 plays an important role.
Collapse
Affiliation(s)
- R R Pela
- Grupo de Materiais Semicondutores e Nanotecnologia, Departamento de Física, Instituto Tecnológico de Aeronáutica, 12228-900 São José dos Campos, SP, Brazil
| | | | | |
Collapse
|
17
|
Klimeš J, Kaltak M, Maggio E, Kresse G. Singles correlation energy contributions in solids. J Chem Phys 2015; 143:102816. [PMID: 26374009 DOI: 10.1063/1.4929346] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The random phase approximation to the correlation energy often yields highly accurate results for condensed matter systems. However, ways how to improve its accuracy are being sought and here we explore the relevance of singles contributions for prototypical solid state systems. We set out with a derivation of the random phase approximation using the adiabatic connection and fluctuation dissipation theorem, but contrary to the most commonly used derivation, the density is allowed to vary along the coupling constant integral. This yields results closely paralleling standard perturbation theory. We re-derive the standard singles of Görling-Levy perturbation theory [A. Görling and M. Levy, Phys. Rev. A 50, 196 (1994)], highlight the analogy of our expression to the renormalized singles introduced by Ren and coworkers [Phys. Rev. Lett. 106, 153003 (2011)], and introduce a new approximation for the singles using the density matrix in the random phase approximation. We discuss the physical relevance and importance of singles alongside illustrative examples of simple weakly bonded systems, including rare gas solids (Ne, Ar, Xe), ice, adsorption of water on NaCl, and solid benzene. The effect of singles on covalently and metallically bonded systems is also discussed.
Collapse
Affiliation(s)
- Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Merzuk Kaltak
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Emanuele Maggio
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria
| |
Collapse
|
18
|
Tran F, Blaha P, Schwarz K. How close are the Slater and Becke-Roussel potentials in solids? J Chem Theory Comput 2015; 11:4717-26. [PMID: 26500460 PMCID: PMC4606396 DOI: 10.1021/acs.jctc.5b00675] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2015] [Indexed: 11/30/2022]
Abstract
The Becke–Roussel (BR) potential [Phys. Rev. A 1989, 39, 3761] was proposed as an approximation to the Slater potential, which is the Coulomb potential generated by the exact exchange hole. In the present work, a detailed comparison between the Slater and BR potentials in solids is presented. It is shown that the two potentials usually lead to very similar results for the electronic structure; however, in a few cases, e.g., Si, Ge, or strongly correlated systems like NiO, the fundamental band gap or magnetic properties can differ markedly. Such differences should not be neglected when the computationally expensive Slater potential is replaced by the cheap semilocal BR potential in approximations to the exact-exchange Kohn–Sham potential, such as the one proposed by Becke and Johnson [J. Chem. Phys. 2006, 124, 221101].
Collapse
Affiliation(s)
- Fabien Tran
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Peter Blaha
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| | - Karlheinz Schwarz
- Institute of Materials Chemistry, Vienna University of Technology, Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
| |
Collapse
|
19
|
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
|
20
|
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
|
21
|
Fukazawa T, Akai H. Optimized effective potential method and application to static RPA correlation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:115502. [PMID: 25739010 DOI: 10.1088/0953-8984/27/11/115502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The optimized effective potential (OEP) method is a promising technique for calculating the ground state properties of a system within the density functional theory. However, it is not widely used as its computational cost is rather high and, also, some ambiguity remains in the theoretical framework. In order to overcome these problems, we first introduced a method that accelerates the OEP scheme in a static RPA-level correlation functional. Second, the Krieger-Li-Iafrate (KLI) approximation is exploited to solve the OEP equation. Although seemingly too crude, this approximation did not reduce the accuracy of the description of the magnetic transition metals (Fe, Co, and Ni) examined here, the magnetic properties of which are rather sensitive to correlation effects. Finally, we reformulated the OEP method to render it applicable to the direct RPA correlation functional and other, more precise, functionals. Emphasis is placed on the following three points of the discussion: (i) level-crossing at the Fermi surface is taken into account; (ii) eigenvalue variations in a Kohn-Sham functional are correctly treated; and (iii) the resultant OEP equation is different from those reported to date.
Collapse
Affiliation(s)
- Taro Fukazawa
- 1-3 Machikaneyama-cho, Toyonaka, Osaka, 560-8531, Japan
| | | |
Collapse
|
22
|
Natan A. Fock-exchange for periodic structures in the real-space formalism and the KLI approximation. Phys Chem Chem Phys 2015; 17:31510-5. [DOI: 10.1039/c5cp01093d] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The calculation of Fock-exchange interaction is an important task in the computation of molecule and solid properties.
Collapse
Affiliation(s)
- Amir Natan
- Department of Physical Electronics
- Tel-Aviv University
- Tel-Aviv
- Israel
| |
Collapse
|