Bosonic statistics from commuting creation operators ($[a_p^{†},a_q^{†}]=0$)

Commuting creation operators make multiparticle states invariant under exchange of their labels. With $[a,a^{†}]=1$, the normalized number states $(a^{†}{)}^r|0⟩/\sqrt{r!}$ exist for all $r\ge 0$, giving unrestricted bosonic occupation numbers. The single-mode thermal sum is $Z={\sum }_{r\ge 0}{e}^{-\beta Er}=(1-e^{-\beta E}{)}^{-1}$ and the mean occupation is $(e^{\beta E}-1{)}^{-1}$. Thus the exchange symmetry and the Bose-Einstein distribution are two expressions of the same bosonic algebra.