The branching ratio of unimolecular decomposition can be evaluated by solving the rate equations. Recent advances in automated reaction path search methods have enabled efficient construction of the rate equations based on quantum chemical calculations. However, it is still difficult to solve the rate equations composed of hundreds or more elementary steps. This problem is especially serious when elementary steps that occur in highly different timescales coexist. In this article, we introduce an efficient approach to obtain the branching ratio from a given set of rate equations. It has been derived from a recently proposed rate constant matrix contraction (RCMC) method, and termed full-RCMC (f-RCMC). The f-RCMC gives the branching ratio without solving the rate equations. Its performance was tested numerically for unimolecular decomposition of C3H5 and C4H5. Branching ratios obtained by the f-RCMC precisely reproduced the values obtained by numerically solving the rate equations. It took about 95 h to solve the rate equations of C4H5 consisting of 234 elementary steps. In contrast, the f-RCMC gave the branching ratio in less than 1 s. The f-RCMC would thus be an efficient alternative of the conventional kinetic simulation approach.
All Science Journal Classification (ASJC) codes
- Computational Mathematics