The expensive cost of computing exact exchange in periodic systems limits the application range of density functional theory with hybrid functionals. To reduce the computational cost of exact change, we present a range-separated algorithm to compute electron repulsion integrals for Gaussian-type crystal basis. The algorithm splits the full-range Coulomb interactions into short-range and long-range parts, which are, respectively, computed in real and reciprocal space. This approach significantly reduces the overall computational cost, as integrals can be efficiently computed in both regions. The algorithm can efficiently handle large numbers of k points with limited central processing unit (CPU) and memory resources. As a demonstration, we performed an all-electron k-point Hartree–Fock calculation for LiH crystal with one million Gaussian basis functions, which was completed on a desktop computer in 1400 CPU hours.
REFERENCES
1.
D. R.
Bowler
and T.
Miyazaki
, “Calculations for millions of atoms with density functional theory: Linear scaling shows its potential
,” J. Phys.: Condens. Matter
22
, 074207
(2010
).2.
J.
VandeVondele
, U.
Borštnik
, and J.
Hutter
, “Linear scaling self-consistent field calculations with millions of atoms in the condensed phase
,” J. Chem. Theory Comput.
8
, 3565
–3573
(2012
).3.
C.
Pisani
and R.
Dovesi
, “Exact-exchange Hartree-Fock calculations for periodic systems. I. Illustration of the method
,” Int. J. Quantum Chem.
17
, 501
–516
(1980
).4.
R.
Dovesi
, C.
Pisani
, C.
Roetti
, and V. R.
Saunders
, “Treatment of Coulomb interactions in Hartree-Fock calculations of periodic systems
,” Phys. Rev. B
28
, 5781
–5792
(1983
).5.
F.
Gygi
and A.
Baldereschi
, “Self-consistent Hartree-Fock and screened-exchange calculations in solids: Application to silicon
,” Phys. Rev. B
34
, 4405
–4408
(1986
).6.
G.
Kresse
and D.
Joubert
, “From ultrasoft pseudopotentials to the projector augmented-wave method
,” Phys. Rev. B
59
, 1758
–1775
(1999
).7.
J.
Heyd
and G. E.
Scuseria
, “Efficient hybrid density functional calculations in solids: Assessment of the Heyd-Scuseria-Ernzerhof screened Coulomb hybrid functional
,” J. Chem. Phys.
121
, 1187
–1192
(2004
).8.
J.
Paier
, R.
Hirschl
, M.
Marsman
, and G.
Kresse
, “The Perdew–Burke–Ernzerhof exchange-correlation functional applied to the G2-1 test set using a plane-wave basis set
,” J. Chem. Phys.
122
, 234102
(2005
).9.
A. F.
Izmaylov
, G. E.
Scuseria
, and M. J.
Frisch
, “Efficient evaluation of short-range Hartree-Fock exchange in large molecules and periodic systems
,” J. Chem. Phys.
125
, 104103
(2006
).10.
J.
Spencer
and A.
Alavi
, “Efficient calculation of the exact exchange energy in periodic systems using a truncated Coulomb potential
,” Phys. Rev. B
77
, 193110
(2008
).11.
M.
Guidon
, J.
Hutter
, and J.
VandeVondele
, “Robust periodic Hartree-Fock exchange for large-scale simulations using Gaussian basis sets
,” J. Chem. Theory Comput.
5
, 3010
–3021
(2009
).12.
X.
Xu
and N. A. W.
Holzwarth
, “Projector augmented wave formulation of Hartree-Fock calculations of electronic structure
,” Phys. Rev. B
81
, 245105
(2010
).13.
R.
Sundararaman
and T. A.
Arias
, “Regularization of the Coulomb singularity in exact exchange by Wigner-Seitz truncated interactions: Towards chemical accuracy in nontrivial systems
,” Phys. Rev. B
87
, 165122
(2013
).14.
N.
Marzari
and D.
Vanderbilt
, “Maximally localized generalized Wannier functions for composite energy bands
,” Phys. Rev. B
56
, 12847
–12865
(1997
).15.
F.
Gygi
, “Compact representations of Kohn-Sham invariant subspaces
,” Phys. Rev. Lett.
102
, 166406
(2009
).16.
N.
Marzari
, A. A.
Mostofi
, J. R.
Yates
, I.
Souza
, and D.
Vanderbilt
, “Maximally localized Wannier functions: Theory and applications
,” Rev. Mod. Phys.
84
, 1419
–1475
(2012
).17.
F.
Gygi
and I.
Duchemin
, “Efficient computation of Hartree-Fock exchange using recursive subspace bisection
,” J. Chem. Theory Comput.
9
, 582
–587
(2013
).18.
R. A.
DiStasio
, B.
Santra
, Z.
Li
, X.
Wu
, and R.
Car
, “The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water
,” J. Chem. Phys.
141
, 084502
(2014
).19.
A.
Görling
, “Exact treatment of exchange in Kohn-Sham band-structure schemes
,” Phys. Rev. B
53
, 7024
–7029
(1996
).20.
S.
Chawla
and G. A.
Voth
, “Exact exchange in ab initio molecular dynamics: An efficient plane-wave based algorithm
,” J. Chem. Phys.
108
, 4697
–4700
(1998
).21.
X.
Wu
, A.
Selloni
, and R.
Car
, “Order-n implementation of exact exchange in extended insulating systems
,” Phys. Rev. B
79
, 085102
(2009
).22.
L.
Lin
, “Adaptively compressed exchange operator
,” J. Chem. Theory Comput.
12
, 2242
–2249
(2016
).23.
W.
Hu
, L.
Lin
, A. S.
Banerjee
, E.
Vecharynski
, and C.
Yang
, “Adaptively compressed exchange operator for large-scale hybrid density functional calculations with applications to the adsorption of water on silicene
,” J. Chem. Theory Comput.
13
, 1188
–1198
(2017
).24.
T. A.
Barnes
, T.
Kurth
, P.
Carrier
, N.
Wichmann
, D.
Prendergast
, P. R. C.
Kent
, and J.
Deslippe
, “Improved treatment of exact exchange in quantum espresso
,” Comput. Phys. Commun.
214
, 52
–58
(2017
).25.
T. D.
Kühne
, M.
Iannuzzi
, M.
Del Ben
, V. V.
Rybkin
, P.
Seewald
, F.
Stein
, T.
Laino
, R. Z.
Khaliullin
, O.
Schütt
, F.
Schiffmann
, D.
Golze
, J.
Wilhelm
, S.
Chulkov
, M. H.
Bani-Hashemian
, V.
Weber
, U.
Borštnik
, M.
Taillefumier
, A. S.
Jakobovits
, A.
Lazzaro
, H.
Pabst
, T.
Müller
, R.
Schade
, M.
Guidon
, S.
Andermatt
, N.
Holmberg
, G. K.
Schenter
, A.
Hehn
, A.
Bussy
, F.
Belleflamme
, G.
Tabacchi
, A.
Glöß
, M.
Lass
, I.
Bethune
, C. J.
Mundy
, C.
Plessl
, M.
Watkins
, J.
VandeVondele
, M.
Krack
, and J.
Hutter
, “CP2K: An electronic structure and molecular dynamics software package—Quickstep: Efficient and accurate electronic structure calculations
,” J. Chem. Phys.
152
, 194103
(2020
).26.
R.
Dovesi
, F.
Pascale
, B.
Civalleri
, K.
Doll
, N. M.
Harrison
, I.
Bush
, P.
D’Arco
, Y.
Noël
, M.
Rérat
, P.
Carbonnière
, M.
Causà
, S.
Salustro
, V.
Lacivita
, B.
Kirtman
, A. M.
Ferrari
, F. S.
Gentile
, J.
Baima
, M.
Ferrero
, R.
Demichelis
, and M.
De La Pierre
, “The CRYSTAL code, 1976-2020 and beyond, a long story
,” J. Chem. Phys.
152
, 204111
(2020
).27.
Q.
Sun
, X.
Zhang
, S.
Banerjee
, P.
Bao
, M.
Barbry
, N. S.
Blunt
, N. A.
Bogdanov
, G. H.
Booth
, J.
Chen
, Z.-H.
Cui
, J. J.
Eriksen
, Y.
Gao
, S.
Guo
, J.
Hermann
, M. R.
Hermes
, K.
Koh
, P.
Koval
, S.
Lehtola
, Z.
Li
, J.
Liu
, N.
Mardirossian
, J. D.
McClain
, M.
Motta
, B.
Mussard
, H. Q.
Pham
, A.
Pulkin
, W.
Purwanto
, P. J.
Robinson
, E.
Ronca
, E. R.
Sayfutyarova
, M.
Scheurer
, H. F.
Schurkus
, J. E. T.
Smith
, C.
Sun
, S.-N.
Sun
, S.
Upadhyay
, L. K.
Wagner
, X.
Wang
, A.
White
, J. D.
Whitfield
, M. J.
Williamson
, S.
Wouters
, J.
Yang
, J. M.
Yu
, T.
Zhu
, T. C.
Berkelbach
, S.
Sharma
, A. Y.
Sokolov
, and G. K.-L.
Chan
, “Recent developments in the PySCF program package
,” J. Chem. Phys.
153
, 024109
(2020
).28.
V.
Blum
, R.
Gehrke
, F.
Hanke
, P.
Havu
, V.
Havu
, X.
Ren
, K.
Reuter
, and M.
Scheffler
, “Ab initio molecular simulations with numeric atom-centered orbitals
,” Comput. Phys. Commun.
180
, 2175
–2196
(2009
).29.
S. G.
Balasubramani
, G. P.
Chen
, S.
Coriani
, M.
Diedenhofen
, M. S.
Frank
, Y. J.
Franzke
, F.
Furche
, R.
Grotjahn
, M. E.
Harding
, C.
Hättig
, A.
Hellweg
, B.
Helmich-Paris
, C.
Holzer
, U.
Huniar
, M.
Kaupp
, A.
Marefat Khah
, S.
Karbalaei Khani
, T.
Müller
, F.
Mack
, B. D.
Nguyen
, S. M.
Parker
, E.
Perlt
, D.
Rappoport
, K.
Reiter
, S.
Roy
, M.
Rückert
, G.
Schmitz
, M.
Sierka
, E.
Tapavicza
, D. P.
Tew
, C.
van Wüllen
, V. K.
Voora
, F.
Weigend
, A.
Wodyński
, and J. M.
Yu
, “TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations
,” J. Chem. Phys.
152
, 184107
(2020
).30.
J.
VandeVondele
, M.
Krack
, F.
Mohamed
, M.
Parrinello
, T.
Chassaing
, and J.
Hutter
, “Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach
,” Comput. Phys. Commun.
167
, 103
–128
(2005
).31.
G.
Lippert
, J.
Hutter
, and M.
Parrinello
, “A hybrid Gaussian and plane wave density functional scheme
,” Mol. Phys.
92
, 477
–488
(1997
).32.
D.
Marx
and J.
Hutter
, “Ab initio molecular dynamics: Basic theory and advanced methods
,” Angew. Chem., Int. Ed.
48
, 9404
–9405
(2009
).33.
C.
Hartwigsen
, S.
Goedecker
, and J.
Hutter
, “Relativistic separable dual-space Gaussian pseudopotentials from H to Rn
,” Phys. Rev. B
58
, 3641
–3662
(1998
).34.
A.
Willand
, Y. O.
Kvashnin
, L.
Genovese
, Á.
Vázquez-Mayagoitia
, A. K.
Deb
, A.
Sadeghi
, T.
Deutsch
, and S.
Goedecker
, “Norm-conserving pseudopotentials with chemical accuracy compared to all-electron calculations
,” J. Chem. Phys.
138
, 104109
(2013
).35.
W.
Gao
and J. R.
Chelikowsky
, “Accuracy of partial core corrections using fourier transforms in pseudopotential-density functional theory
,” J. Chem. Theory Comput.
14
, 6515
–6520
(2018
).36.
G.
Prandini
, A.
Marrazzo
, I. E.
Castelli
, N.
Mounet
, and N.
Marzari
, “Precision and efficiency in solid-state pseudopotential calculations
,” npj Comput. Mater.
4
, 72
(2018
).37.
J.-B.
Lu
, D. C.
Cantu
, M.-T.
Nguyen
, J.
Li
, V.-A.
Glezakou
, and R.
Rousseau
, “Norm-conserving pseudopotentials and basis sets to explore lanthanide chemistry in complex environments
,” J. Chem. Theory Comput.
15
, 5987
–5997
(2019
).38.
K.
Lejaeghere
, G.
Bihlmayer
, T.
Björkman
, P.
Blaha
, S.
Blügel
, V.
Blum
, D.
Caliste
, I. E.
Castelli
, S. J.
Clark
, A.
Dal Corso
, S.
de Gironcoli
, T.
Deutsch
, J. K.
Dewhurst
, I.
Di Marco
, C.
Draxl
, M.
Dułak
, O.
Eriksson
, J. A.
Flores-Livas
, K. F.
Garrity
, L.
Genovese
, P.
Giannozzi
, M.
Giantomassi
, S.
Goedecker
, X.
Gonze
, O.
Grånäs
, E. K. U.
Gross
, A.
Gulans
, F.
Gygi
, D. R.
Hamann
, P. J.
Hasnip
, N. A. W.
Holzwarth
, D.
Iuşan
, D. B.
Jochym
, F.
Jollet
, D.
Jones
, G.
Kresse
, K.
Koepernik
, E.
Küçükbenli
, Y. O.
Kvashnin
, I. L. M.
Locht
, S.
Lubeck
, M.
Marsman
, N.
Marzari
, U.
Nitzsche
, L.
Nordström
, T.
Ozaki
, L.
Paulatto
, C. J.
Pickard
, W.
Poelmans
, M. I. J.
Probert
, K.
Refson
, M.
Richter
, G.-M.
Rignanese
, S.
Saha
, M.
Scheffler
, M.
Schlipf
, K.
Schwarz
, S.
Sharma
, F.
Tavazza
, P.
Thunström
, A.
Tkatchenko
, M.
Torrent
, D.
Vanderbilt
, M. J.
van Setten
, V.
Van Speybroeck
, J. M.
Wills
, J. R.
Yates
, G.-X.
Zhang
, and S.
Cottenier
, “Reproducibility in density functional theory calculations of solids
,” Science
351
, aad3000
(2016
).39.
M.
Krack
and M.
Parrinello
, “All-electron ab-initio molecular dynamics
,” Phys. Chem. Chem. Phys.
2
, 2105
–2112
(2000
).40.
E.
Kucukbenli
, M.
Monni
, B. I.
Adetunji
, X.
Ge
, G. A.
Adebayo
, N.
Marzari
, S.
de Gironcoli
, and A.
Dal Corso
, “Projector augmented-wave and all-electron calculations across the periodic table: A comparison of structural and energetic properties
,” arXiv:1404.3015 [cond-mat.mtrl-sci]
(2014
).41.
S. V.
Levchenko
, X.
Ren
, J.
Wieferink
, R.
Johanni
, P.
Rinke
, V.
Blum
, and M.
Scheffler
, “Hybrid functionals for large periodic systems in an all-electron, numeric atom-centered basis framework
,” Comput. Phys. Commun.
192
, 60
–69
(2015
).42.
T.
Biktagirov
, W. G.
Schmidt
, and U.
Gerstmann
, “Calculation of spin-spin zero-field splitting within periodic boundary conditions: Towards all-electron accuracy
,” Phys. Rev. B
97
, 115135
(2018
).43.
C.-N.
Yeh
, A.
Shee
, Q.
Sun
, E.
Gull
, and D.
Zgid
, “Relativistic self-consistent GW: Exact two-component formalism with one-electron approximation for solids
,” Phys. Rev. B
106
, 085121
(2022
).44.
F.
Gygi
, “All-electron plane-wave electronic structure calculations
,” J. Chem. Theory Comput.
19
, 1300
–1309
(2023
).45.
G.
Lippert
, J.
Hutter
, and M.
Parrinello
, “The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations
,” Theor. Chem. Acc.
103
, 124
–140
(1999
).46.
M. D.
Towler
, A.
Zupan
, and M.
Causà
, “Density functional theory in periodic systems using local Gaussian basis sets
,” Comput. Phys. Commun.
98
, 181
–205
(1996
).47.
E.
Schwegler
and M.
Challacombe
, “Linear scaling computation of the Fock matrix. IV. Multipole accelerated formation of the exchange matrix
,” J. Chem. Phys.
111
, 6223
–6229
(1999
).48.
J.
Heyd
, G. E.
Scuseria
, and M.
Ernzerhof
, “Hybrid functionals based on a screened Coulomb potential
,” J. Chem. Phys.
118
, 8207
–8215
(2003
).49.
C. J.
Tymczak
, V. T.
Weber
, E.
Schwegler
, and M.
Challacombe
, “Linear scaling computation of the Fock matrix. VIII. Periodic boundaries for exact exchange at the γ point
,” J. Chem. Phys.
122
, 124105
(2005
).50.
L.
Maschio
, D.
Usvyat
, F. R.
Manby
, S.
Casassa
, C.
Pisani
, and M.
Schütz
, “Fast local-MP2 method with density-fitting for crystals. I. Theory and algorithms
,” Phys. Rev. B
76
, 075101
(2007
).51.
Š.
Varga
, “Long-range analysis of density fitting in extended systems
,” Int. J. Quantum Chem.
108
, 1518
–1527
(2008
).52.
L.
Maschio
and D.
Usvyat
, “Fitting of local densities in periodic systems
,” Phys. Rev. B
78
, 073102
(2008
).53.
M.
Del Ben
, J.
Hutter
, and J.
VandeVondele
, “Electron correlation in the condensed phase from a resolution of identity approach based on the Gaussian and plane waves scheme
,” J. Chem. Theory Comput.
9
, 2654
–2671
(2013
).54.
T.
Yamada
, R. P.
Brewster
, and S.
Hirata
, “Asymptotic expansion of two-electron integrals and its application to Coulomb and exchange lattice sums in metallic, semimetallic, and nonmetallic crystals
,” J. Chem. Phys.
139
, 184107
(2013
).55.
Q.
Sun
, T. C.
Berkelbach
, J. D.
McClain
, and G. K.
Chan
, “Gaussian and plane-wave mixed density fitting for periodic systems
,” J. Chem. Phys.
147
, 164119
(2017
).56.
A.
Irmler
, A. M.
Burow
, and F.
Pauly
, “Robust periodic Fock exchange with atom-centered Gaussian basis sets
,” J. Chem. Theory Comput.
14
, 4567
–4580
(2018
).57.
X.
Wang
, C. A.
Lewis
, and E. F.
Valeev
, “Efficient evaluation of exact exchange for periodic systems via concentric atomic density fitting
,” J. Chem. Phys.
153
, 124116
(2020
).58.
J.
Lee
, A.
Rettig
, X.
Feng
, E.
Epifanovsky
, and M.
Head-Gordon
, “Faster exact exchange for solids via occ-RI-K: Application to combinatorially optimized range-separated hybrid functionals for simple solids with pseudopotentials near the basis set limit
,” J. Chem. Theory Comput.
18
, 7336
–7349
(2022
).59.
S.
Sharma
, A. F.
White
, and G.
Beylkin
, “Fast exchange with Gaussian basis set using robust pseudospectral method
,” J. Chem. Theory Comput.
18
, 7306
–7320
(2022
).60.
H.-Z.
Ye
and T. C.
Berkelbach
, “Fast periodic Gaussian density fitting by range separation
,” J. Chem. Phys.
154
, 131104
(2021
).61.
H.-Z.
Ye
and T. C.
Berkelbach
, “Tight distance-dependent estimators for screening two-center and three-center short-range Coulomb integrals over Gaussian basis functions
,” J. Chem. Phys.
155
, 124106
(2021
).62.
S. J.
Bintrim
, T. C.
Berkelbach
, and H.-Z.
Ye
, “Integral-direct Hartree-Fock and Møller-Plesset perturbation theory for periodic systems with density fitting: Application to the benzene crystal
,” J. Chem. Theory Comput.
18
, 5374
–5381
(2022
).63.
Q.
Sun
, “Various integral estimations and screening schemes for extended systems in PySCF
,” arXiv:2302.11307 [physics.chem-ph]
(2023
).64.
L.
Füsti-Molnár
and P.
Pulay
, “The Fourier transform Coulomb method: Efficient and accurate calculation of the Coulomb operator in a Gaussian basis
,” J. Chem. Phys.
117
, 7827
–7835
(2002
).65.
J.
Paier
, C. V.
Diaconu
, G. E.
Scuseria
, M.
Guidon
, J.
VandeVondele
, and J.
Hutter
, “Accurate Hartree-Fock energy of extended systems using large Gaussian basis sets
,” Phys. Rev. B
80
, 174114
(2009
).© 2023 Author(s). Published under an exclusive license by AIP Publishing.
2023
Author(s)
You do not currently have access to this content.
