In our past work, we described the advantages of pursuing topological partitioning of the electron energy density, 𝛸(𝒓) (c.f., Eq. ( )). The topological features of 𝛸(𝒓) resembles those of the electron density -although effectively scaled by the value of the local electronegativity, 𝜒̅ (𝒓) -with sharp cusps at nuclear coordinates and with the function decaying to zero at large distances from the nuclei. Whereas our initial study relied on an orbital averaging approach (viz. Eq. ( ) to compute 𝛸(𝒓) in molecules, we here take on the challenge of extending our approach to crystals.