Tags: blog, method
$$ \begin{equation} \hat{H} = \frac{\omega_{c}}{2}\left(\hat{A}^{\dagger}\hat{A} + \hat{A}\hat{A}^{\dagger}\right) \end{equation} $$
with $\hat{A} = \hat{a} f(\hat{n})$.