\hat{H}=\int \mathrm{d^3 k} \hbar \omega_{\vec{k}} \left(\hat{a}^\dagger_{\vec{k}}\hat{a}_{\vec{k}}+\frac{1}{2}\right)