Similarly, one can investigate the accuracy of the structures of such non-covalently bound complexes obtained with sDFT geometry optimizations. Such studies were enabled by several new analytical gradient implementations for FDE/sDFT have appeared in the past years, including a Slater-function-based code in the Amsterdam Density Functional FDE-module [141], and an implementation in the subsystem quantum chemistry