“…Several schemes have been proposed in the last decades to extract from ρ(r) as much chemical information as possible, which involve its visualization (profiles, 2D maps, 3D graphical representations), partition according to different schemes, topological analysis, etc. [4][5][6][7][8][9][10][11][12][13][14][15] In this respect, one of the most powerful and popular techniques is represented by Bader's quantum theory of atoms in molecules (QTAIM), which relies on the topological analysis of ρ(r). 16 If a routine analysis of the wave-function |Ψ of a small system is a non particularly demanding task from a computational point of view (when compared to the convergence of the self-consistent field, SCF, procedure or to the analytical evaluation of energy gradients), this is clearly no more the case either when rather sophisticated techniques are adopted (such as QTAIM) or when large systems are studied.…”