Problem 1: Death Process¶
In this problem, we study a death process, which has the unique reaction $A \to \emptyset$ with rate constant $\delta$.
a) Write the ODE for the mean population.
Show hint
b) Write down the Chemical Master Equation for the stochastic kinetics. Qualitatively compare the two.
Show hint
- From state $n$, any individual may die, so the transition rate $n \mapsto n-1$ is $n \cdot \delta$.
- From state $n+1$, we can transition to $n$ at rate $\delta \cdot (n+1)$.
CME:
$$
\frac{dP_n}{dt} = \delta \cdot (n+1) \cdot P_{n+1} - \delta \cdot n \cdot P_n
$$
c) Find the expected value of the extinction time $T_0$ when starting from an initial population of $n$.
Show hint
- First event is a death with exponential waiting time mean $1/\delta n$.
- After the first death, the population is $n-1$.
Recurrence:
$$
T_0(n) = \frac{1}{\delta n} + T_0(n-1), \quad T(0) = 0
$$
d) Find the variance of the extinction time $T_0$.
Show hint
- Variance of an exponential with rate $\lambda$ is $1/\lambda^2$.
- Then sum the variances of each waiting time.
e) Argue that $\mathbb{P}(T_0(n) \leq t) = (1 - e^{-\delta t})^n$.
Show hint
f) Find the derivative $\frac{d}{dt}\mathbb{P}(T_0(n) \leq t)$.
Show hint
Use the chain rule.
g) Consider two independent identical death processes, one for species $A$ starting at an initial population of $n_A$ and one for species $B$ starting at an initial population of $n_B$. Denoting by $T_A$ the extinction time of $A$ and by $T_B$ that of $B$, find (an estimation of) the probability that $T_A \leq T_B$.
Show hint
Use the law of total probability and the fact that the derivative from part f) is the probability density function of $T_0(n)$.
Problem 2: Birth and Death Process¶
In this problem, we add a birth reaction to the death process from the previous problem, i.e., the reaction $A \to A + A$ with rate constant $\gamma$.
a) Write down an ODE describing the deterministic kinetics.
Show hint
b) Write down the Chemical Master Equation for the stochastic kinetics. Qualitatively compare the two.
Show hint
- Transitions: $n \mapsto n+1$ at rate $\gamma n$ and $n \mapsto n-1$ at rate $\delta n$.
CME for $n \geq 1$:
$$
\frac{dP_n}{dt} = \gamma (n-1) P_{n-1} + \delta (n+1) P_{n+1} - (\gamma+\delta) n P_n
$$
c) Let $q_n$ be the probability of extinction when starting from an initial population of $n$. Find a recursion formula for $q_n$.
Show hint
d) Find a solution of the recusion formula for $q_n$.
Show hint
Make the ansatz $q_n = r^n$ for some $r\in [0,1]$ and solve the resulting quadratic equation. Then distinguish the cases $\gamma > \delta$ and $\gamma \leq \delta$.Intuition: when births are stronger than deaths, the population can "escape" extinction.
e) For the case of almost-sure extinction, find the expected extinction time. Evaluate this formula for $n=1,2,3$.