eQE 2.0: Subsystem DFT beyond GGA functionals

By adopting a divide-and-conquer strategy, subsystem-DFT (sDFT) can dramatically reduce the computational cost of large-scale electronic structure calculations. The key ingredients of sDFT are the nonadditive kinetic energy and exchange-correlation functionals which dominate it's accuracy. Even though, available semilocal nonadditive functionals find a broad range of applications, their accuracy is somewhat limited especially for those systems where achieving balance between exchange-correlation interactions on one side and nonadditive kinetic energy on the other is crucial. In eQE 2.0, we improve dramatically the accuracy of sDFT simulations by (1) implementing nonlocal nonadditive kinetic energy functionals based on the LMGP family of functionals; (2) adapting Quantum ESPRESSO's implementation of rVV10 and vdW-DF nonlocal exchange-correlation functionals to be employed in sDFT simulations; (3) implementing "deorbitalized" meta GGA functionals (e.g., SCAN-L). We carefully assess the performance of the newly implemented tools on the S22-5 test set. eQE 2.0 delivers excellent interaction energies compared to conventional Kohn-Sham DFT and CCSD(T). The improved performance does not come at a loss of computational efficiency. We show that eQE 2.0 with nonlocal nonadditive functionals retains the same linear scaling behavior achieved previously in eQE 1.0 with semilocal nonadditive functionals. Program summary: Program title: eQE CPC Library link to program files: https://doi.org/10.17632/g55n9xmnt2.1 Developer's repository link: https://gitlab.com/Pavanello/eqe Licensing provisions: GPLv2 Programming language: Fortran 90 External routines/libraries: BLACS, MPI Nature of problem: Solving the electronic structure of molecules and materials with subsystem density functional theory Solution method: This version of eQE uses subsystem Density-Functional Theory (DFT) to compute the electronic structure of molecules and materials. Subsystem DFT is a divide-and-conquer version of DFT that can be massively parallelized and is linear scaling in work and data. However, it requires the use of density functionals for kinetic and exchange-correlation energies. eQE can now employ nonlocal functionals as well as meta-GGA functionals for these energy terms. Additional comments: Url of stable release http://eqe.rutgers.edu