$$ \newcommand{\op}[1]{\hat{#1}} \renewcommand{\vec}[1]{\mathbf{#1}} \newcommand{\ket}[1]{|{#1}\rangle} \newcommand{\bra}[1]{\langle {#1}|} $$ Let's say we have a spherically symmetric system, like a hydrogen atom, in which we label the eigenstates of the three (commuting) observables $\op H$, $\op{\vec L}^2$ and $\op L_3$, the hamiltonian, the squared angular momentum and its third component, such that $$ \op H\ket{n,l,m}=E_n\ket{n,l,m}\\ \op{\vec L}^2\ket{n,l,m}=l(l+1)\ket{n,l,m}\\ \op L_3\ket{n,l,m}=m\ket{n,l,m} $$ A state with $l=0$, which is of the type $$ \ket{\psi}=\sum_{n\in\mathbb N}a_n\ket{n,0,0} $$ has a radially symmetric wave function $\psi(\vec x)$ so the expected value $\bra{\psi}\op x_i\ket{\psi}$ has to be $0$ for each $i=1,2,3$. This tells us that $\ket{n,0,0}$ and $\op x_i\ket{n,0,0}$ are orthogonal states, so $\op x_i\ket{n,0,0}$ is in a subspace with $l\ge 1$.
Now, the spherical harmonics with $l=1$ are three linear combinations of the spatial coordinates $x_1,x_2,x_3$ and are linearly independent, so we could write $x,y,z$ as linear combinations of $Y_{1,-1}$, $Y_{1,0}$ and $Y_{1,1}$, e.g. $x=\frac12(Y_{1,1}+Y_{1,-1})$. Can this be seen as a "proof" that $\op x_i\ket{n,0,0}$ is a state with $l=1$? I get the broader picture, but I'm stuck mathematically: are there operators associated to the spherical harmonics which we can use in this case? Could they be just, for example, $\op Y_{1,1}=\op x+i\op y$? That seems just circular reasoning to me. And, having "defined" these spherical harmonic operators, how does applying the operator $\frac12(\op Y_{1,1}+\op Y_{1,-1})$ to $\ket{n,0,0}$ bring it to a state with $l=1$?