e scaling of correlated wave function methods lies, in their standard implementation, between O(N 5 ) and O(e N ), where N is proportional to the system size. erefore, performing accurate calculations on chemically meaningful systems requires i) approximations that can lower the computational scaling, and ii) e cient implementations that take advantage of modern massively parallel architectures.