[Kupershmidt](https://arxiv.org/abs/math/0001187) argues for $$\mu_{K}(j)=(pq)^j (1-pq)$$ as the "natural" $q$-geometric distribution (with a simple normalisation constant), because it produces the $q$-analog of the Pascal (negative binomial) distribution.