I would like to convince myself of the following relationship in an astrophysical context:
\begin{aligned} & \sum_{m}\sum_{m^{\prime}}\left\langle a_{\ell m} a_{\ell m}^* a_{\ell m^{\prime}} a_{\ell m^{\prime}}^*\right\rangle \\ & =\left[\sum_m\left\langle a_{\ell m} a_{\ell m}^* a_{\ell m} a_{\ell m}^*\right\rangle+\sum_{m^{\prime}\neq m} \left\langle a_{\ell m} a_{\ell m}^* a_{\ell m^{\prime}} a_{\ell m^{\prime}}^*\right\rangle\right] \\ & =\left(3 C_{\ell}^2(2 \ell+1)+2 \ell C_{\ell}^2(2 \ell+1)\right). \end{aligned}
I apply Wick's theorem for computing the expectation of $a_{\ell m} $:
EDIT 1: I think that my issue of understanding is about the second sum when $m\neq m'$. I can't represent myself the number of terms in this sum ( equal to $2\ell\ \times(2\ell +1)$)
EDIT 2: I think the keypoint is the number of different pairs ($m,m^{\prime}$).
$A_{n}^{p}= \dfrac{(2\ell+1)!}{(2\ell-1)!}= (2\ell+1)\times 2\ell$
Is this interpretation correct?