one may encounter points in space where 𝜌𝜌 � 𝑏𝑏 (𝐫𝐫) < 0 if some of the natural orbitals have occupation numbers 0 < 𝑛𝑛 𝑖𝑖 < 1. A negative effective density is not physical, and it can produce spurious results. To tackle this issue in our earlier versions of DC functionals, we set 𝜌𝜌 � 𝑏𝑏 (𝐫𝐫) and ∇ 𝐫𝐫 𝜌𝜌 � 𝑏𝑏 (𝐫𝐫) to zero for quadrature points r where 𝜌𝜌 � 𝑏𝑏 (𝐫𝐫) < 0: