\hat{a}^+_{\vec{n}}=\frac{1}{\sqrt{N}}\sum_{\vec{k  }}\hat{a}^+_{\vec{k}}e^{-i\vec{k}\cdot\vec{n}