A simple birth process has rates qn,n+1=nλ. Its pgf solves Gt=λz(z−1)Gz, giving G(z,t)=e−λtz/[1−(1−e−λt)z]. Thus from one ancestor Xt is geometric with parameter e−λt and mean eλt. From i ancestors it is a sum of i such variables, hence negative binomial with mean ieλt.