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 imple- mented 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 situ...

We extend our recently developed quantum-mechanical/molecular mechanics (QM/MM) approach [Dziedzic et al., J. Chem. Phys. 145, 124106 (2016)] to enable in situ optimization of the localized orbitals. The quantum subsystem is described with ONETEP linear-scaling density functional theory and the classical subsystem – with the AMOEBA polarizable force field. The two subsystems interact via multipolar electrostatics and are fully...

We present an overview of the ONETEP program for linear-scaling density functional theory (DFT) calculations with large basis set (planewave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the...

The correlated-band theory implemented as a combination of the local-density approximation with the dynamical mean-field theory is applied to PuO2. An insulating electronic structure, consistent with the experimental valence-band photoemission spectra, is obtained. The calculations yield a nonmagnetic ground state that is characterized by a noninteger filling of the plutonium 5f shell. The noninteger filling as well as the satellites...