An important cornerstone in the development of modern DFT approximations has been the introduction of spin-resolved same-spin (SS) and opposite-spin (OS) components of the correlation energy . Thus, it would also be useful to spin-resolve each term of e c (r). To do that, we can decompose the underlying ∆γ λ (r) = γ λ -γ 0 and ρ λ (r) into spin-up and spin-down components. This yields spin-up and spin-down contributions for e J,λ c (r) and e K,λ c (r), along with spin-resolution of their corresponding λ-averages (Eqs. 35b and 35c). For closedshell systems, each spin component contributes half of the total e J c (r) and e K c (r) values. More interesting is the spin-resolved e ee c (r), which requires decomposing the two-body pair density P λ 2 (r, r 2 ) into SS and OS contributions,