By adopting a divide-and-conquer strategy, subsystem-DFT (sDFT) can dramatically reduce the computational cost of large-scale electronic structure calculations. The key ingredients sDFT are nonadditive kinetic energy and exchange-correlation functionals which dominate it's accuracy. Even though, semilocal find broad range applications, their accuracy is somewhat limited especially for those sys...