Spectral theory of metastability and extinction in a branching-annihilation reaction

TL;DR

This paper uses spectral methods to analytically study the metastability and extinction times in a reaction process involving branching and annihilation, providing accurate results that match numerical simulations.

Contribution

It introduces a spectral approach combining generating functions and Sturm-Liouville theory to analyze complex birth-death processes with new analytical insights.

Findings

Analytical probability distribution for metastable states

Accurate extinction time statistics

Resolution of boundary condition issues in spectral methods

Abstract

We apply the spectral method, recently developed by the authors, to calculate the statistics of a reaction-limited multi-step birth-death process, or chemical reaction, that includes as elementary steps branching A->2A and annihilation 2A->0. The spectral method employs the generating function technique in conjunction with the Sturm-Liouville theory of linear differential operators. We focus on the limit when the branching rate is much higher than the annihilation rate, and obtain accurate analytical results for the complete probability distribution (including large deviations) of the metastable long-lived state, and for the extinction time statistics. The analytical results are in very good agreement with numerical calculations. Furthermore, we use this example to settle the issue of the "lacking" boundary condition in the spectral formulation.