Accurate, Large-Scale and Affordable Hybrid-PBE0 Calculations with GPU-Accelerated Supercomputers Laura E. Ratcliff Argonne Leadership Computing Facility, Argonne National Laboratory, Illinois 60439, USA Department of Materials, Imperial College London, London SW7 2AZ, UK E-mail:
[email protected] A. Degomme Department of Physics, Universität Basel, Klingelbergstr. 82, 4056 Basel, Switzerland José A. Flores-Livas Department of Physics, Universität Basel, Klingelbergstr. 82, 4056 Basel, Switzerland Stefan Goedecker Department of Physics, Universität Basel, Klingelbergstr. 82, 4056 Basel, Switzerland Luigi Genovese Univ. Grenoble Alpes, CEA, INAC-SP2M, L_Sim, F-38000, Grenoble, France E-mail:
[email protected] December 2017 Abstract. Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding – when not out of reach – if high quality basis sets are used. We present a highly efficient multiple GPU implementation of the exact exchange operator which allows hybrid functional density-functional theory calculations with systematic basis sets without additional approximations for up to a thousand atoms. This method is implemented in a portable real-space- based algorithm, released as an open-source package. With such a framework hybrid DFT calculations of high quality become accessible on state-of-the-art supercomputers within a time-to-solution of the same order of magnitude as traditional semilocal-GGA functionals. arXiv:1712.07973v1 [cond-mat.mtrl-sci] 21 Dec 2017 Accurate, Large-Scale and Affordable Hybrid-PBE0 Calculations with GPU-Accelerated Supercomputers 2 1. Introduction results coming from these codes employ hybrid XC functionals. The main reason is that, for typical Density-functional theory (DFT) in principle is not an systems, hybrid functional calculations are one to approximation, and should produce the exact ground- two orders of magnitude more expensive than DFT- state energy and density, but in practice the crucial GGA calculations.