Minimum Tracking Linear Response Hubbard and Hund Corrected Density Functional Theory in CP2K

We present the implementation of the Hubbard (<i>U</i>) and Hund (<i>J</i>) corrected Density Functional Theory (DFT + <i>U</i> + <i>J</i>) functionality in the Quickstep program, which is part of the CP2K suite. The tensorial and Löwdin subspace representations are implemented and compared. Full analytical DFT + <i>U</i> + <i>J</i> forces are implemented and benchmarked for the tensorial and Löwdin representations. We also present the implementation of the recently proposed minimum-tracking linear-response method that enables the <i>U</i> and <i>J</i> parameters to be calculated on first-principles basis without reference to the Kohn-Sham eigensystem. These implementations are benchmarked against recent results for different materials properties including DFT + <i>U</i> band gap opening in NiO, the relative stability of various polaron distributions in TiO₂, the dependence of the calculated TiO₂ band gap on +<i>J</i> corrections, and, finally, the role of the +<i>U</i> and +<i>J</i> corrections for the computed properties of a series of the hexahydrated transition metals. Our implementation provides results consistent with those already reported in the literature from comparable methods. We conclude the contribution with tests on the influence of the Löwdin orthonormalization on the occupancies, calculated parameters, and derived properties.