Daubechies Wavelets for Linear Scaling Density Functional Theory

We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.

Many-body Greens function GW and Bethe-Salpeter study of the optical excitations in a paradigmatic model dipeptide

We study within the many-body Greens function GW and Bethe-Salpeter formalisms the excitation energies of a paradigmatic model dipeptide, focusing on the four lowest-lying local and charge-transfer excitations. Our GW calculations are performed at the self-consistent level, updating first the quasiparticle energies, and further the single-particle wavefunctions within the static Coulomb-hole plus screened-exchange approximation to the GW self-energy operator. Important level crossings, as compared to the starting Kohn-Sham LDA spectrum, are identified. Our final Bethe-Salpeter singlet excitation energies are found to agree, within 0.07 eV, with CASPT2 reference data, except for one charge-transfer state where the discrepancy can be as large as 0.5 eV. Our results agree best with LC-BLYP and CAM-B3LYP calculations with enhanced long-range exchange, with a 0.1 eV mean absolute error. This has been achieved employing a parameter-free formalism applicable to metallic or insulating extended or finite systems.