In this paper, we introduce and examine a fractional linear birth–death process Nν(t), t>0, whose fractionality is obtained by replacing the time derivative with a fractional derivative in the system of difference-differential equations governing the state probabilities pkν(t), t>0, k≥0. We present a subordination relationship connecting Nν(t), t>0, with the classical birth–death process N(t), t>0, by means of the time process T2ν(t), t>0, whose distribution is related to a time-fractional diffusion equation.
¶
We obtain explicit formulas for the extinction probability p0ν(t) and the state probabilities pkν(t), t>0, k≥1, in the three relevant cases λ>μ, λ<μ, λ=μ (where λ and μ are, respectively, the birth and death rates) and discuss their behaviour in specific situations. We highlight the connection of the fractional linear birth–death process with the fractional pure birth process. Finally, the mean values $\mathbb{E}N_{\nu}(t)$ and $\operatorname{\mathbb{V}ar}N_{\nu}(t)$ are derived and analyzed.