I would rewrite the Hamiltonian using two independent creation and annihilation operators
$$
\begin{align}
\sqrt{\frac{\hbar}{2 m\omega_x}}(\hat{a}^{\dagger} + \hat{a}) & = x, \\
i\sqrt{\frac{\hbar m\omega_x}{2}}(\hat{a}^{\dagger} - \hat{a}) & = p_x, \\
\sqrt{\frac{\hbar}{2 m\omega_y}}(\hat{b}^{\dagger} + \hat{b}) & = y, \\
i\sqrt{\frac{\hbar m\omega_y}{2}}(\hat{b}^{\dagger} - \hat{b}) & = p_y.
\end{align}
$$
The Hamiltonian has the following form
$$\hat{\mathcal{H}} = \hbar \omega_x (\hat{a}^{\dagger}\hat{a} + 1/2) + \hbar \omega_y(\hat{b}^{\dagger}\hat{b} + 1/2).$$
Coherent state $|\alpha, \beta\rangle$ is defined as an eigenstate of annihilation operators
$$\hat{a}|\alpha, \beta\rangle = \alpha|\alpha,\beta\rangle \ \ \ \hat{b}|\alpha, \beta\rangle = \beta|\alpha,\beta\rangle.$$
In fact $|\alpha,\beta\rangle$ is a product of two independent coherent states $|\alpha,\beta\rangle = |\alpha\rangle \otimes |\beta\rangle$. When you initially start evolution with $|\alpha,\beta\rangle$ you will end up with $|\alpha e^{-i\omega_x t}, \beta e^{-i\omega_y t}\rangle$ at time $t$.
Coherent state
A collection of operators $\{\hat{a},\hat{b},\hat{a}^{\dagger},\hat{b}^{\dagger},1\}$ span the Heisenberg algebra. We follow the definition of Perelomov. General elements of the Heisenberg group can be written like this
$$G(\gamma,\alpha,\beta) = \exp\left(\gamma 1 + \alpha\hat{a}^{\dagger} - \alpha^{*}\hat{a} + \beta\hat{b}^{\dagger} - \beta^{*}\hat{b} \right) = \exp\left(\alpha\hat{a}^{\dagger} - \alpha^{*}\hat{a}\right) \exp\left(\beta\hat{b}^{\dagger} - \beta^{*}\hat{b} \right)\exp(1\gamma).$$
We choose a specific state in the Hilbert space, in our case it will be vacuum $|0\rangle$, and act with $G(\gamma,\alpha,\beta)$ on it
$$|\gamma,\alpha,\beta\rangle = G(\gamma,\alpha,\beta)|0\rangle.$$
Because $\exp(1\gamma)$ imprints just a global phase factor it is not a part of the definition of the coherent state. A coherent state is defined like this
$$|\alpha,\beta\rangle = \exp\left(\alpha\hat{a}^{\dagger} - \alpha^{*}\hat{a}\right) \exp\left(\beta\hat{b}^{\dagger} - \beta^{*}\hat{b} \right)|0\rangle$$
You can define similar coherent states for n-dimensional harmonic oscillator.