Tensor approximation of stationary distributions of chemical reaction networks

Abstract

We prove that the stationary distribution of a system of reacting species with a weakly- reversible reaction network of zero deficiency in the sense of Feinberg admits tensor- structured approximation of complexity which scales linearly with respect to the number of species and logarithmically in the maximum copy numbers as well as in the desired accuracy. Our results cover the classical mass-action and also Michaelis-Menten kinetics which correspond to two widely used classes of propensity functions, and also to arbitrary combinations of those. New rank bounds for tensor-structured approximations of the PDF of a truncated one-dimensional Poisson distribution are an auxiliary result of the present paper, which might be of independent interest. The present work complements recent results obtained by the authors jointly with M. Khammash and M. Nip on the tensor-structured numerical simulation of the evolution of system states distributions, driven by the Kolmogorov forward equation of the system, known also as the chemical master equation, or CME for short. For the two kinetics mentioned above we also analyze the low-rank tensor structure of the CME operator.