With the definition and derivation of the relevant quantities for e c (r) we can focus now on obtaining highly accurate energy densities. For this purpose, e c,ee (r) and other related quantities described in Sec. IV have been implemented using the PySCF 2.6.2 package . Crucial to this implementation is the full configuration interaction (FCI) to compute Ψ λ (called here λ-FCI), the minimizer of Ĥλ in Eq. 1.