The long-range contribution 1 R AB is calculated using the smooth particle mesh Ewald method as implemented for tight binding approaches in CP2K, scaling as N log N with system size N . Exchange interactions based on γ EX are treated by implying the minimum image convention, thus restricting the sum over AB to neighbors within the unit cell.