A dynamical low-rank approach to the chemical master equation

Stochastic reaction kinetics have increasingly been used to study cellular systems, with applications ranging from viral replication to gene regulatory networks and to signaling pathways. The underlying evolution equation, known as the chemical master equation (CME), can rarely be solved with traditional methods due to the huge number of degrees of freedom. We present a new approach to directly solve the CME by a dynamical low-rank approximation based on the Dirac-Frenkel-McLachlan variational principle. The new approach has the capability to substantially reduce the number of degrees of freedom, and to turn the CME into a computationally tractable problem. We illustrate the accuracy and efficiency of our methods in application to two examples of biological interest.