The angular momentum operator $L_z$ can be expressed in terms of $L_z = xp_y - yp_x$ where $ p = \hat{p} = -i\sqrt{hwm/2}(a^{\dagger} - a)$ and $x = \hat{x}\sqrt{\hbar/2mw}(a^{\dagger} + a)$.
We want to express $L_z$ in terms of the ladder operators $a$ and $a^{\dagger}$. By direct insertion we then have
$L_z = i \frac{\hbar}{2}((a_x^{\dagger} + a_x)(a_y^{\dagger} - a_y) - (a_y^{\dagger} + a_y)(a_x^{\dagger} - a_x))$
Here's where my problems begin. I've multiplied it out, and tried to invoke both some commutation relations and using the fact that $[a, a^{\dagger}] = 1$ to rewrite the expression, but I cant seem to get the correct answer. Im not that secure on the commutations and relations that can be used the rewrite these operators, obviously. Specifically, I'm not sure what commutations are valid when we're concerned with two directions.
The correct answer is $L_z = i\hbar(a_xa_y^{\dagger} - a_x^{\dagger}a_y)$. Would someone be kind enough to give me some hints or tricks to better manipulate these?