Asymptotics of Quasi-Stationary Distributions of Small Noise Stochastic Dynamical Systems in Unbounded Domains

TL;DR

This paper investigates the asymptotic behavior of quasi-stationary distributions in small noise stochastic systems modeling biological populations, extending previous results to unbounded domains and continuous-time dynamics.

Contribution

It extends the analysis of QSD asymptotics to unbounded state spaces and continuous-time models, providing new support characterizations and extinction time bounds.

Findings

Limit points of QSD supported on interior attractors

Expected extinction times scale exponentially with system size

Established existence and tightness of QSD in unbounded domains

Abstract

We consider a collection of Markov chains that model the evolution of multitype biological populations. The state space of the chains is the positive orthant, and the boundary of the orthant is absorbing representing the extinction states of different population types. We are interested in the long-term behavior of the Markov chain away from extinction, under a small noise scaling. Under this scaling, the trajectory of the Markov process over any compact interval converges in distribution to the solution of an ordinary differential equation (ODE) evolving in the positive orthant. We study the asymptotic behavior of the quasi-stationary distributions (QSD) in this scaling regime. Our main result shows that the limit points of the QSD are supported on the union of interior attractors of the flow determined by the ODE. We also give lower bounds on expected extinction times which scale exponentially with the system size. Results of this type when the deterministic dynamical system obtained under the scaling limit is given by a discrete time evolution equation and the dynamics are essentially in a compact space (namely, the one step map is a bounded function) have been studied by Faure and Schreiber (2014). Our results extend these to a setting of an unbounded state space and continuous time dynamics. The proofs rely on uniform large deviation results for small noise stochastic dynamical systems and methods from the theory of dynamical systems. In general QSD for Markov chains with absorbing states and unbounded state spaces may not exist. We study one basic family of Binomial-Poisson models in the positive orthant where one can use Lyapunov function methods to establish existence of QSD and also to argue the tightness of the QSD of the scaled sequence of Markov chains. The results from the first part are then used to characterize the support of limit points of this sequence of QSD.