 -\sum_{\vec{k}}\hat{a}^+_{\vec{k}}\hat{a}_{\vec{k}}  \frac{1}{N}\sum_{\vec{q}}J(\vec{q})\langle \hat{a}^+_{\vec{q}}\hat{b}^+_{-\vec{q}}\rangle 