The Wigner function for a wave function $\Psi(\vec{r})$ is $$ W(\vec{r},\vec{k}) = \frac{1}{2\pi} \int dy e^{-i \vec{k} \cdot \vec{y}} \Psi^{*}(\vec{r}-\vec{y}/2) \Psi(\vec{r}+\vec{y}/2) . \tag{1} $$ Recently, I've come to learn of the generalized Cahill phase space representation $$ F^{s}(\alpha) = \frac{1}{\pi^{2}} \int d^{2}\xi e^{\alpha\xi^{*}-\alpha^{*}\xi+\frac{s}{2}|\xi|^{2}} \mathrm{Tr}D(\xi)\rho , \tag{2} $$ which is a continuous family of quasiprobability distribution functions with respect to the parameter $s\in\left[-1,1\right]$. Here, $D(\xi) = e^{\xi a^{\dagger} - \xi^{*} a}$ is the displacement operator, where $a$ and $a^{\dagger}$ are the annihilation and creation operators satisfying $[a, a^{\dagger}] = 1$. As $s$ interpolates from $-1$ to $0$ to $1$, it is able to take the form of the Glauber–Sudarshan $P$ function, the Wigner Function and the Husimi $Q$ function.
This means to me that the Wigner function can be represented in a coherent state expansion $W(\alpha)$, but to my knowledge not all states admit to a quantum optics description (like maybe a simple hydrogen atom), therefore can they be presented as such? Is it realizable only mathematically, or does it possess physical meaning in coherent state space?