Shifted-exponential Poisson mixture ($G(z)=e^{\alpha (z-1)}\beta /(\beta +1-z)$)

When the random intensity is $\alpha +Y$ with $Y$ having an exponential distribution of rate $\beta$, the mixed count is the sum of independent Poisson distribution and support-zero geometric distribution counts. The Poisson parameter is $\alpha$ and the geometric success probability is $\beta /(\beta +1)$. At $\alpha =0$ the law is geometric.

For the shifted-exponential Poisson mixture, the displayed relation holds for $n\ge 2$, initialized by $p_0=e^{-\alpha }\beta /(\beta +1)$ and $p_1=(\alpha +1/(\beta +1))p_0$. It follows by differentiating the probability generating function and comparing coefficients. Nonnegativity follows from the independent Poisson-geometric convolution of independent random variables representation.