Selection rule on angular momentum quantum number for vector operators According to the Wigner-Eckart theorem, the matrix element of a vector operator in the basis of eigenstates of the square of the magnitude of the angular momentum ($J^2$) and the z-component of the angular momentum ($J_z$) is proportional to a Clebsch-Gordan coefficient (I am aware that this theorem applies to all spherical tensor operators, but seeing as I'm not very familiar with these other tensor operators, it's best to stick with vector operators, though I'll take an answer discussing more general tensor operators as well). From this, it's concluded that the selection rule on the angular momentum quantum number is
$$\Delta j=0, \pm1 \tag{1}$$
since the relevant Clebsch-Gordan coefficients vanish whenever (1) is violated.
On Wikipedia, the theorem is derived by showing that the matrix elements satisfy the same recursion relations as the Clebsch-Gordan coefficients. The problem is that I don't see how this leads to the selection rules (1). Why can't there also be some other sets of quantities (apart from something directly proportional to the Clebsch-Gordan coefficients) that satisfy the recursion relations, but that don't vanish whenever (1) is violated? The selection rules on $J_z$ can be found by directly applying the commutators between $J_z$ and the vector operator, without having to worry about the recursion relations. But I've not seen a similar thing being done to $J^2$. So how does one arrive at the selection rule (1)?
 A: The key point is that
$$
T^{1}_{m_1}\vert \ell_2,m_2\rangle \tag{1}
$$
transforms in exactly the same way as
\begin{align}
\vert 1,m_1\rangle\vert \ell_2,m_2\rangle \, .\tag{2}
\end{align}  You can see this because
\begin{align}
L_zT^{1}_{m_1}\vert \ell_2m_2\rangle &= 
[L_z,T^{1}_{m_1}]\vert \ell 2 m_2\rangle +
T^{1}_{m_1}L_z\vert \ell_2 m_2\rangle\, ,\\
&=(m_1+m_2)T^1_{m_1}\vert\ell_2 m_2\rangle
\end{align}
and likewise the actions
\begin{align}
L_\pm T^{1}_{m_1}\vert \ell_2 m_2\rangle 
&=[L_\pm,T^1_{m_1}]\vert \ell_2m_2\rangle + T^{1}_{m_1} L_\pm \vert \ell_2m_2\rangle\, ,\\
&\sim L_\pm \left[\vert 1m_1\rangle \vert \ell_2 m_2\rangle\right]
=\left[L_\pm\vert 1m_1\rangle\right]\vert\ell_2 m_2\rangle 
+\vert 1m_1\rangle\left[L_\pm \vert \ell_2m_2\rangle\right]
\end{align}
where by $\sim$ I mean the same factors with square roots come out, and
of course $[L_\pm,T^{1}_{m_1}]\propto T^1_{m\pm 1}$.
If you believe that (1) does indeed transform in the same way as (2) does, then it is clear that, if (2) yields states with $j$ values ranging from $\ell_2+1$ to $\vert \ell_2-1\vert$, then the combination in (1) will also yields states with values of $j$ in that range.  That’s basically the selection rule $\Delta j=0,\pm 1$.
Obviously if you do this not with a vector operator but a general tensor
$T^{\ell_1}$, you will find the range of $j$ is just the same as the range in $\ell_1\otimes\ell_2$.
