Quasistationarity and extinction for population processes

Damian Clancy
2024-12-10
Abstract:We consider stochastic population processes that are almost surely absorbed at the origin within finite time. Our interest is in the quasistationary distribution, $\boldsymbol{u}$, and the expected time, $\tau$, from quasistationarity to extinction, both of which we study via WKB approximation. This approach involves solving a Hamilton-Jacobi partial differential equation specific to the model. We provide conditions under which analytical solution of the Hamilton-Jacobi equation is possible, and give the solution. This provides a first approximation to both $\boldsymbol{u}$ and $\tau$. We provide further conditions under which a corresponding `transport equation' may be solved, leading to an improved approximation of $\boldsymbol{u}$. For multitype birth and death processes, we then consider an alternative approximation for $\boldsymbol{u}$ that is valid close to the origin, provide conditions under which the elements of this alternative approximation may be found explicitly, and hence derive an improved approximation for $\tau$. We illustrate our results in a number of applications.
Probability
What problem does this paper attempt to address?
The core problems that this paper attempts to solve are two key quantities in the stochastic population process that is almost surely absorbed to zero within a finite time: the quasi - stationary distribution \(u\) and the expected time \(\tau\) from the quasi - stationary state to extinction. Specifically, the author uses the WKB (Wentzel, Kramers, Brillouin) approximation method to study these quantities and provides solutions to these problems by solving the Hamilton - Jacobi partial differential equation specific to the model. ### Specific description of the problem 1. **Quasi - stationary distribution \(u\)**: - Before the population finally goes extinct, the system will reach an apparent equilibrium state, that is, the quasi - stationary distribution \(u\). - The quasi - stationary distribution describes the state distribution of the population over a long period before extinction occurs. 2. **Expected time \(\tau\) from the quasi - stationary state to extinction**: - This time is from when the system enters the quasi - stationary distribution until the final extinction. - This time is usually a random variable with an exponential distribution, and its distribution is completely determined by its expected value \(\tau\). ### Method of the solution To study these two key quantities, the author adopts the following methods: - **WKB approximation**: - The WKB approximation is an asymptotic analysis method, especially suitable for the case of a large parameter \(N\). - The author assumes that the quasi - stationary distribution \(u\) can be expressed as \(u_x = M_N\exp(-N V(y)-V_0(y)+O(1/N))\), where \(y = x/N\), \(V(y)\) and \(V_0(y)\) are functions independent of \(N\), and \(M_N\) is a constant independent of \(y\). - **Hamilton - Jacobi partial differential equation**: - Through the WKB approximation, the author derives a Hamilton - Jacobi partial differential equation to determine \(V(y)\). - Further, the author also derives a "transport equation" to improve the estimate of \(V_0(y)\). - **Conditions and results**: - The author provides a method for analytically solving the Hamilton - Jacobi equation under certain conditions and gives specific solutions. - For the multi - type birth - death process, the author also provides an alternative approximation method that is valid near the origin, thus obtaining a more accurate estimate of \(\tau\). ### Application examples The paper provides multiple application examples, including but not limited to: - **SIS infection model**: The SIS model with heterogeneous susceptibility and infection period. - **Multi - type birth - death process**: The process with linear birth rate and quadratic death rate. - **Competition process**: The model involving competition among multiple species. Through these methods, the author can better understand and predict the behavior of the population before extinction, which is of great significance for fields such as ecology, epidemiology, and immunology.