A general method for the computation of probabilities in systems of first order chemical reactions

We present a general method for the computation of molecular population distributions in a system of first-order chemical reactions. The method is based on the observation that the molecules in first-order reactions do not interact with each other. Therefore, the population distributions are a convolution of densities for one molecule. With this method one can study reactions involving any number of molecules. Such analysis is demonstrated on several examples, including an enzyme catalyst reaction and a first-order reaction chain.