This section is devoted to a numerical discussion of our MP2 correlation energy density implementation in the Python coding language . The resulting generator reads the relevant atomic orbital functions from a HF treatment using the PySCF library and performs the necessary MP2-based integrations and tensor contrac-tions. This procedure scales similar to other second-order perturbation theory derived models, and thus, we can employ some numerical approximation techniques to enhance the applicability and efficiency of our code. The effect on the accuracy can be discussed using the correlation energy result from the already implemented MP2 module in PySCF (refMP2).