We present a method for the calculation of four-centre two-electron repulsion integrals in terms of localised non-orthogonal generalised Wannier functions (NGWFs). Our method has been implemented in the ONETEP program and is used to compute the Hartree-Fock exchange energy component of Hartree-Fock and Density Functional Theory (DFT) calculations with hybrid exchange-correlation functionals. As the NGWFs are optimised in situin terms of a systematically improvable basis set which is equivalent to plane waves, it is possible to achieve large basis set accuracy in routine calculations. The spatial localisation of the NGWFs allows us to exploit the exponential decay of the density matrix in systems with a band gap in order to compute the exchange energy with a computational effort that increases linearly with the number of atoms. We describe the implementation of this approach in the ONETEPprogram for linear-scaling first principles quantum mechanical calculations. We present extensive numerical validation of all the steps in our method. Furthermore, we find excellent agreement in energies and structures for a wide variety of molecules when comparing with other codes. We use our method to perform calculations with the B3LYP exchange-correlation functional for models of myoglobin systems bound with O2 and CO ligands and confirm that the same qualitative behaviour is obtained as when the same myoglobin models are studied with the DFT+U approach which is also available in ONETEP. Finally, we confirm the linear-scaling capability of our method by performing calculations on polyethylene and polyacetylene chains of increasing length.

1.
P.
Hohenberg
and
W.
Kohn
,
Phys. Rev.
136
,
B864
(
1964
).
2.
W.
Kohn
and
L. J.
Sham
,
Phys. Rev.
140
,
A1133
(
1965
).
3.
J.
Tao
,
J. P.
Perdew
,
V. N.
Staroverov
, and
G. E.
Scuseria
,
Phys. Rev. Lett.
91
,
146401
(
2003
).
4.
E.
Artacho
and
L.
Miláns del Bosch
,
Phys. Rev. A
43
,
5770
(
1991
).
5.
E.
Schwegler
,
M.
Challacombe
, and
M.
Head-Gordon
,
J. Chem. Phys.
106
,
9708
(
1997
).
6.
C.
Ochsenfeld
,
C. A.
White
, and
M.
Head-Gordon
,
J. Chem. Phys.
109
,
1663
(
1998
).
7.
S. A.
Maurer
,
D. S.
Lambrecht
,
D.
Flaig
, and
C.
Ochsenfeld
,
J. Chem. Phys.
136
,
144107
(
2012
).
8.
E.
Rudberg
,
E. H.
Rubensson
, and
P.
Salek
,
J. Chem. Theory Comput.
7
,
340
(
2011
).
9.
J.
Kussmann
and
C.
Ochsenfeld
,
J. Chem. Phys.
138
,
134114
(
2013
).
10.
M.
Guidon
,
J.
Hutter
, and
J.
VandeVondele
,
J. Chem. Theory Comput.
5
,
3010
(
2009
).
11.
M.
Del Ben
,
J.
Hutter
, and
J.
VandeVondele
,
J. Chem. Theory Comput.
8
,
4177
(
2012
).
12.
X.
Ren
,
P.
Rinke
,
V.
Blum
,
J.
Wieferink
,
A.
Tkatchenko
,
A.
Sanfilippo
,
K.
Reuter
, and
M.
Scheffler
,
New J. Phys.
14
,
053020
(
2012
).
13.
B. I.
Dunlap
,
J. W. D.
Connolly
, and
J. R.
Sabin
,
J. Chem. Phys.
71
,
3396
(
1979
).
14.
B. I.
Dunlap
,
J. W. D.
Connolly
, and
J. R.
Sabin
,
J. Chem. Phys.
71
,
4993
(
1979
).
15.
C.-K.
Skylaris
,
L.
Gagliardi
,
N. C.
Handy
,
A. G.
Ioannou
,
S.
Spencer
, and
A.
Willetts
,
J. Mol. Struct.: THEOCHEM
501–502
,
229
(
2002
).
16.
R.
Polly
,
H.-J.
Werner
,
F.
Manby
, and
P. J.
Knowles
,
Mol. Phys.
102
,
2311
(
2004
).
17.
P.
Merlot
,
T.
Kjaergaard
,
T.
Helgaker
,
R.
Lindh
,
F.
Aquilante
,
S.
Reine
, and
T. B.
Pedersen
,
J. Comput. Chem.
34
,
1486
(
2013
).
18.
B. I.
Dunlap
,
J. Mol. Struct.: THEOCHEM
529
,
37
(
2000
).
19.
M.
Lorenz
,
L.
Maschio
,
M.
Schütz
, and
D.
Usvyat
,
J. Chem. Phys.
137
,
204119
(
2012
).
20.
M.
Gibson
,
S.
Brandt
, and
S.
Clark
,
Phys. Rev. B
73
,
125120
(
2006
).
21.
M.
Marsman
,
J.
Paier
,
A.
Stroppa
, and
G.
Kresse
,
J. Phys.: Condens. Matter
20
,
064201
(
2008
).
22.
X.
Wu
,
A.
Selloni
, and
R.
Car
,
Phys. Rev. B
79
,
085102
(
2009
).
23.
C.-K.
Skylaris
,
A. A.
Mostofi
,
P. D.
Haynes
,
C. J.
Pickard
, and
M. C.
Payne
,
Comput. Phys. Commun.
140
,
315
(
2001
).
24.
P. J.
Stephens
,
F. J.
Devlin
,
C. F.
Chabalowski
, and
M. J.
Frisch
,
J. Phys. Chem.
98
,
11623
(
1994
).
25.
R.
Peverati
and
D. G.
Truhlar
,
Phys. Chem. Chem. Phys.
14
,
16187
(
2012
).
26.
D. D.
O'Regan
,
N. D. M.
Hine
,
M. C.
Payne
, and
A. A.
Mostofi
,
Phys. Rev. B
85
,
085107
(
2012
).
27.
D. D.
O'Regan
,
N. D. M.
Hine
,
M. C.
Payne
, and
A. A.
Mostofi
,
Phys. Rev. B
82
,
081102
(
2010
).
28.
C.-K.
Skylaris
,
P. D.
Haynes
,
A. A.
Mostofi
, and
M. C.
Payne
,
J. Chem. Phys.
122
,
084119
(
2005
).
29.
C.-K.
Skylaris
,
A. A.
Mostofi
,
P. D.
Haynes
,
O.
Diéguez
, and
M. C.
Payne
,
Phys. Rev. B
66
,
035119
(
2002
).
30.
A. A.
Mostofi
,
P. D.
Haynes
,
C.-K.
Skylaris
, and
M. C.
Payne
,
J. Chem. Phys.
119
,
8842
(
2003
).
31.
P. D.
Haynes
,
C.-K.
Skylaris
,
A. A.
Mostofi
, and
M. C.
Payne
,
J. Physics: Condens. Mat.
20
,
294207
(
2008
).
32.
J.
Almlöf
,
K.
Korsel
, and
K.
Faegri
 Jr.
,
J. Comput. Chem.
3
,
385
(
1982
).
33.
C.-K.
Skylaris
,
P. D.
Haynes
,
A. A.
Mostofi
, and
M. C.
Payne
,
Phys. Status Solidi B
243
,
973
(
2006
).
34.
M. R.
Jarvis
,
I. D.
White
,
R. W.
Godby
, and
M. C.
Payne
,
Phys. Rev. B
56
,
14972
(
1997
).
35.
N. D. M.
Hine
,
J.
Dziedzic
,
P. D.
Haynes
, and
C.-K.
Skylaris
,
J. Chem. Phys.
135
,
204103
(
2011
).
36.
O.
Vahtras
,
J.
Almlöf
, and
M. W.
Feyereisen
,
Chem. Phys. Lett.
213
,
514
(
1993
).
37.
P. D.
Haynes
and
M. C.
Payne
,
Comput. Phys. Commun.
102
,
17
(
1997
).
38.
Q.
Hill
, “
Development of more accurate computational methods within linear-scaling density functional theory
,” Ph.D. thesis (
University of Southampton
,
2010
).
39.
B.
Monserrat
and
P. D.
Haynes
,
J. Phys. A: Math. Theor.
43
,
465205
(
2010
).
40.
N. D. M.
Hine
,
P. D.
Haynes
,
A. A.
Mostofi
,
C.-K.
Skylaris
, and
M. C.
Payne
,
Comput. Phys. Commun.
180
,
1041
(
2009
).
41.
A. A.
Mostofi
,
C.-K.
Skylaris
,
P. D.
Haynes
, and
M. C.
Payne
,
Comput. Phys. Commun.
147
,
788
(
2002
).
42.
S.
Goedecker
and
G.
Scuseria
,
Comput. Sci. Eng.
5
,
14
(
2003
).
43.
V. N.
Staroverov
,
G. E.
Scuseria
,
J.
Tao
, and
J. P.
Perdew
,
J. Chem. Phys.
119
,
12129
(
2003
).
44.
CRC Handbook of Chemistry and Physics
, 3rd ed., edited by
D. R.
Lide
(
CRC
,
Boca Raton, FL
,
2002
).
45.
J.
Oláh
and
J.
Harvey
,
J. Phys. Chem. A
113
,
7338
(
2009
).
46.
D. J.
Cole
,
D. D.
O'Regan
, and
M. C.
Payne
,
J. Phys. Chem. Lett.
3
,
1448
(
2012
).
You do not currently have access to this content.