Transcript
Page 1: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

Efficient methods for computing exchange-correlation potentials for

orbital-dependent functionals

Viktor N. Staroverov

Department of Chemistry, The University of Western Ontario, London, Ontario, Canada

IWCSE 2013, Taiwan National University, Taipei, October 14β€’17, 2013

Page 2: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

2

Orbital-dependent functionals

𝐸XC [𝜌 ]=∫ 𝑓 ( {πœ™π‘– })  π‘‘ 𝐫

β€’ More flexible than LDA and GGAs (can satisfy more exact constraints)

β€’ Needed for accurate description of molecular properties

Kohn-Sham orbitals

Page 3: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

3

Examples

β€’ Exact exchange

β€’ Hybrids (B3LYP, PBE0, etc.)

β€’ Meta-GGAs (TPSS, M06, etc.)

𝐸Xexact [𝜌 ]=βˆ’

14βˆ‘π‘– , 𝑗=1

𝑁

βˆ«π‘‘ π«βˆ«π‘‘ 𝐫 β€² πœ‘π‘– (𝐫 )πœ‘ π‘—βˆ— (𝐫 )πœ‘π‘–

βˆ— (𝐫 β€²)πœ‘ 𝑗 (𝐫 β€² )  |π«βˆ’π«β€²| 

 

same expression as in the Hartreeβ€’Fock theory

Page 4: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

4

The challenge

𝑣 XC (𝐫 )=𝛿 𝐸X C

❑ [{πœ™π‘– }]  π›ΏπœŒ (𝐫)  

=?

Kohnβ€’Sham potentials corresponding to orbital-dependent functionals

cannot be evaluated in closed form

Page 5: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

5

Optimized effective potential (OEP)method

𝛿𝐸 total❑  

𝛿𝑣 XC (𝐫)  =0

Find as the solution to the minimization problem

OEP = functional derivative of the functional

Page 6: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

6

Computing the OEP

Expand the Kohnβ€’Sham orbitals:

Expand the OEP:

𝑣 X C (𝐫 )=βˆ‘π‘˜=1

π‘š

π‘π‘˜ 𝑓 π‘˜(𝐫)

πœ™π‘– (𝐫 )=βˆ‘π‘˜=1

𝑛

π‘π‘˜π‘– πœ’π‘˜(𝐫 )

Minimize the total energy with respect to {} and {}

orbital basis functions

auxiliary basis functions

Page 7: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

7

Attempts to obtain OEP-X in finite basis sets

size

Page 8: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

8

I. First approximation to the OEP:An orbital-averaged potential (OAP)

οΏ½Μ‚οΏ½XCπœ™ 𝑖 (𝐫 )=𝛿 𝐸X C

❑ [{πœ™π‘– }]  

π›Ώπœ™π‘–βˆ—(𝐫 )  

Define operator such that

The OAP is a weighted average:

𝑣 XC (𝐫 )=βˆ‘π‘–=1

𝑁

πœ™π‘–βˆ— (𝐫 ) οΏ½Μ‚οΏ½XCπœ™ 𝑖 (𝐫 )

βˆ‘π‘–=1

𝑁

πœ™π‘–βˆ—(𝐫 )πœ™π‘– (𝐫 )

Page 9: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

9

Example: Slater potential

Fock exchange operator:

Slater potential:

οΏ½Μ‚οΏ½ πœ™ 𝑖 (𝐫 ) ≑𝛿𝐸 X

exact

π›Ώπœ™π‘–βˆ—(𝐫 )

𝑣 S (𝐫 )= 1𝜌 (𝐫 ) βˆ‘π‘–=1

𝑁

πœ™π‘–βˆ—(𝐫) οΏ½Μ‚οΏ½ πœ™ 𝑖(𝐫)

Page 10: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

10

Calculation of orbital-averaged potentials

β€’ by definition (hard, functional specific)

β€’ by inverting the Kohnβ€’Sham equations (easy, general)

Page 11: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

11

Kohnβ€’Sham inversion

𝜏𝐿

𝜌+𝑣+𝑣H +𝑣XC= 1

πœŒβˆ‘π‘–=1

𝑁

πœ– 𝑖|πœ™π‘–|2

[βˆ’ 12βˆ‡2+𝑣+𝑣H +𝑣XC ]πœ™π‘–=πœ– π‘–πœ™ 𝑖

Kohnβ€’Sham equations:

multiply by ,sum over i,divide by

Page 12: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

12

LDA-X potential via Kohn-Sham inversion

Page 13: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

13

PBE-XC potential via Kohnβ€’Sham inversion

Page 14: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

14

A. P. Gaiduk,I. G. Ryabinkin, VNS,JCTC 9, 3959 (2013)

Removal of oscillations

Page 15: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

15

Kohnβ€’Sham inversion for orbital-specific potentials

𝜏𝐿

𝜌+𝑣+𝑣H +𝑣XC= 1

πœŒβˆ‘π‘–=1

𝑁

πœ– 𝑖|πœ™π‘–|2

[βˆ’ 12βˆ‡2+𝑣+𝑣H +οΏ½Μ‚οΏ½XC ]πœ™π‘–=πœ– π‘–πœ™ 𝑖

Generalized Kohnβ€’Sham equations:

same manipulations

Page 16: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

16

Example: Slater potential through Kohnβ€’Sham inversion

𝑣 S (𝐫 )=

14βˆ‡

2

𝜌 (𝐫 ) βˆ’πœ (𝐫)+βˆ‘π‘–=1

𝑁

πœ– π‘–βˆ¨πœ™ 𝑖❑(𝐫 )|2

𝜌 (𝐫 )βˆ’π‘£ (𝐫 ) βˆ’π‘£H (𝐫)

𝜏=12βˆ‘π‘–=1

𝑁

ΒΏβˆ‡πœ™ π‘–βˆ¨ΒΏ2=𝜏𝐿+14βˆ‡2𝜌 ΒΏ

where

Page 17: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

17

Slater potential via Kohnβ€’Sham inversion

Page 18: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

18

OAPs constructed by Kohnβ€’Sham inversion

Page 19: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

19

Correlation potentials via Kohnβ€’Sham inversion

Page 20: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

20

Kohnβ€’Sham inversion for a fixed set of Hartreeβ€’Fock orbitals

𝑣 XOEP β‰ˆπ‘£ X

model=βˆ’πœ 𝐿

HF+βˆ‘π‘–=1

𝑁

πœ– π‘–ΒΏπœ™ 𝑖HF |2

𝜌HF βˆ’π‘£βˆ’π‘£HHF

Slater potential:

𝑣 SHF=

βˆ’πœπΏHF +βˆ‘

𝑖=1

𝑁

πœ– 𝑖HF ΒΏπœ™ 𝑖

HF |2

𝜌HF βˆ’π‘£βˆ’π‘£HHF

But if , then

Page 21: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

21

Dependence of KS inversion on orbital energies

Page 22: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

22

II. Assumption that the OEP and HF orbitals are the same

The assumption

leads to the eigenvalue-consistent orbital-averaged potential (ECOAP)

πœ™π‘–=πœ™π‘–HF

𝑣 XECOAP=𝑣S

HF + 1𝜌HF βˆ‘

𝑖=1

𝑁

(πœ– π‘–βˆ’πœ–ΒΏΒΏ 𝑖HF)|πœ™ 𝑖HF|2 ΒΏ

Page 23: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

23

ECOAP KLI LHF

Page 24: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

24

Calculated exact-exchange (EXX) energies

, mEh

KLI ELP=LHF=CEDA ECOAP

m.a.v. 2.88 2.84 2.47

Sample: 12 atoms from He to BaBasis set: UGBS

A. A. Kananenka, S. V. Kohut, A. P. Gaiduk, I. G. Ryabinkin, VNS, JCP 139, 074112 (2013)

Page 25: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

25

III. Hartreeβ€’Fock exchange-correlation (HFXC) potential

An HFXC potential is the which reproduces a HF density within the Kohnβ€’Sham scheme:

𝜌 (𝐫 )=βˆ‘π‘–=1

𝑁

|πœ™π‘– (𝐫 )|2=ΒΏβˆ‘

𝑖=1

𝑁

|πœ™ 𝑖HF (𝐫 )|2

=𝜌HF (𝐫 )¿

[βˆ’ 12βˆ‡2+𝑣 (𝐫 )+𝑣H (𝐫 )+𝑣XC (𝐫 )]πœ™ 𝑖(𝐫 )=πœ–π‘–πœ™π‘– (𝐫 )

That is, is such that

Page 26: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

26

Inverting the Kohn–Sham equations

𝜏𝐿

𝜌+𝑣+𝑣H +𝑣XC= 1

πœŒβˆ‘π‘–=1

𝑁

πœ– 𝑖|πœ™π‘–|2

[βˆ’ 12βˆ‡2+𝑣+𝑣H +𝑣XC ]πœ™π‘–=πœ– π‘–πœ™ 𝑖

Kohnβ€’Sham equations:

local ionizationpotential

multiply by ,sum over i,divide by

Page 27: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

27

Inverting the Hartree–Fock equations

𝜏 𝐿HF

𝜌HF +𝑣+𝑣H +𝑣SHF= 1

𝜌HF βˆ‘π‘–=1

𝑁

πœ–π‘–HF|πœ™ 𝑖

HF|2

Hartreeβ€’Fock equations:

Slater potential builtwith HF orbitals

[βˆ’ 12βˆ‡2+𝑣+𝑣H +𝐾 ]πœ™π‘–

HF=πœ– 𝑖HF πœ™π‘–

HF

same manipulations

Page 28: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

28

Closed-form expression for the HFXC potential

𝑣 XCHF=𝑣S

HF + 1πœŒβˆ‘π‘–=1

𝑁

πœ– π‘–βˆ¨πœ™ 𝑖 |2 βˆ’

1𝜌HF βˆ‘

𝑖=1

𝑁

πœ– 𝑖HF|πœ™π‘–

HF|2+ 𝜏

HF

𝜌HF βˆ’πœπœŒ

, but , , and

We treat this expression as a model potential within the Kohnβ€’Sham SCF scheme.

Here

Computational cost: same as KLI and Beckeβ€’Johnson (BJ)

Page 29: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

29

HFXC potentials are practically exact OEPs!

Numerical OEP: Engel et al.

Page 30: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

30

Page 31: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

31

Page 32: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

32

HFXC potentials can be easily computed for molecules

Numerical OEP: Makmal et al.

Page 33: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

33

Energies from exchange potentials

, mEh

KLI LHF BJ Basis-set OEP HFXC

m.a.v. 1.74 1.66 5.30 0.12 0.05

Sample: 12 atoms from Li to CdBasis set for LHF, BJ, OEP (aux=orb), HFXC: UGBS KLI and true OEP values are from Engel et al.

I. G. Ryabinkin, A. A. Kananenka, VNS, PRL 139, 013001 (2013)

Page 34: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

34

Virial energy discrepancies

, mEh

KLI LHF BJ Basis-set OEP HFXC

m.a.v. 438.0 629.2 1234.1 1.76 2.76

where

𝐸 vir= ∫ 𝑣X (𝐫 ) [3 𝜌 (𝐫 )+𝐫 β‹…βˆ‡ 𝜌 (𝐫) ]𝑑𝐫

For exact OEPs,

𝐸 vir βˆ’πΈEXX=0 ,

Page 35: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

35

HFXC potentials in finite basis sets

Page 36: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

36

Hierarchy of approximations to the EXX potential

𝑣 X❑=𝑣S

HF + 1𝜌HF βˆ‘

𝑖=1

𝑁

(πœ– π‘–βˆ’πœ–π‘–HF )|πœ™π‘–

HF|2+ 𝜏

HF βˆ’πœπœŒHF

OAP

ECOAP

HFXC

Page 37: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

37

Summary

β€’ Orbital-averaged potentials (e.g., Slater) can be constructed by Kohnβ€’Sham inversion

β€’ Hierarchy or approximations to the OEP: OAP (Slater) < ECOAP < HFXC

β€’ ECOAP Slater potential KLI LHF

β€’ HFXC potential OEP

β€’ Same applies to all occupied-orbital functionals

Page 38: Efficient methods for computing exchange-correlation potentials for orbital-dependent functionals Viktor N. Staroverov Department of Chemistry, The University

38

Acknowledgments

β€’ Eberhard Engelβ€’ Leeor Kronik

for OEP benchmarks


Top Related