System size expansion using Feynman rules and diagrams

TL;DR

This paper introduces a diagrammatic Feynman rule approach to the Chemical Master Equation, enabling systematic calculation of fluctuations and corrections beyond the linear noise approximation in stochastic biochemical systems.

Contribution

It develops a novel diagrammatic perturbation method for the Chemical Master Equation, allowing higher-order correlation calculations and improved noise spectrum analysis.

Findings

Derived closed-form leading order corrections to the noise spectrum

Applied method to noise-induced oscillations in the Brusselator

Demonstrated accuracy of the approach for large fluctuations

Abstract

Few analytical methods exist for quantitative studies of large fluctuations in stochastic systems. In this article, we develop a simple diagrammatic approach to the Chemical Master Equation that allows us to calculate multi-time correlation functions which are accurate to a any desired order in van Kampen's system size expansion. Specifically, we present a set of Feynman rules from which this diagrammatic perturbation expansion can be constructed algorithmically. We then apply the methodology to derive in closed form the leading order corrections to the linear noise approximation of the intrinsic noise power spectrum for general biochemical reaction networks. Finally, we illustrate our results by describing noise-induced oscillations in the Brusselator reaction scheme which are not captured by the common linear noise approximation.