Double-Hybrid Density functionals for the condensed phase: gradients, stress tensor, and Auxiliary-Density Matrix Method acceleration

Due to their high accuracy, Double-Hybrid Density functionals emerged to important methods for molecular electronic-structure calculations. The high computational costs of double-hybrid calculations in condensed phase and the lack of efficient gradient implementations thereof inhibit a wide applicability for periodic systems. We present an implementation of gradients for Double-Hybrid functional theory into CP2K. The Auxiliary Density Matrix Method (ADMM) reduces the overhead from the Hartree-Fock calculations providing an efficient and accurate methodology to tackle condensed phase systems. First applications to water containing systems of different densities and molecular crystals pave the way for advanced studies. We present large benchmark systems to discuss the efficiency of our methodology on modern super computing hardware.