From the explicit expressions for the spherical harmonics it is straightforward to derive that $$x = \rho \sqrt{\frac{2\pi}{3}}(Y_{1,-1}(\theta, \phi) - Y_{1,1}(\theta, \phi)),$$
$$y = i\rho \sqrt{\frac{2\pi}{3}}(Y_{1,-1}(\theta, \phi) + Y_{1,1}(\theta, \phi)),$$
$$z = \rho...