I_{\vec{n},\vec{m}}=\frac{1}{N}\sum_{\vec{k}}J(\vec{k})e^{i\vec{k}(\vec{n}-\vec{m})