 -\sum_{\vec{n},\vec{m}}I_{\vec{n},\vec{m}}\hat{a}^+  _{\vec{n}}\hat{a}_{\vec{n}}\langle \hat{a}^+_{\vec{n}}\hat{b}^+_{\vec{m}}\rangle