[论文解读] Combinatorics of chemical reaction systems

Nicolas Behr, Gérard Duchamp+1 more|arXiv (Cornell University)|Dec 18, 2017

Protein Structure and Dynamics参考文献 26被引用 6

一句话总结

本文提出了一种化学反应系统的随机力学框架，将组合方法与统计物理统一，推导出概率分布时间演化的精确解析解。通过用Sobolev-Jacobi正交多项式替代McQuarrie的有缺陷的特征基，解决了二元反应系统中长期存在的问题，实现了对全部六种基本单物种反应类型的闭式解。

ABSTRACT

We propose a concise stochastic mechanics framework for chemical reaction systems that allows to formulate evolution equations for three general types of data: the probability generating functions, the exponential moment generating functions and the factorial moment generating functions. This formulation constitutes an intimate synergy between techniques of statistical physics and of combinatorics. We demonstrate how to analytically solve the evolution equations for all six elementary types of single-species chemical reactions by either combinatorial normal-ordering techniques, or, for the binary reactions, by means of Sobolev-Jacobi orthogonal polynomials. The former set of results in particular highlights the relationship between infinitesimal generators of stochastic evolution and parametric transformations of probability distributions.

研究动机与目标

• 开发一种统一的随机力学形式化方法，将组合生成函数与连续时间马尔可夫链整合于化学反应系统中。
• 解决二元化学反应缺乏精确解析解的问题，特别是涉及自催化或双分子过程的反应。
• 通过引入Sobolev-Jacobi正交多项式作为正确的特征基，解决McQuarrie特征基方法在二元反应中的不一致性。
• 证明半线性正规排序技术和生成函数方法可在全部六种基本单物种反应类型中获得精确解。
• 为将该框架推广至多物种和图重写系统奠定基础，表明其在单物种系统之外也具有更广泛的应用潜力。

提出的方法

• 使用纯态的向量空间及其上的概率分布来构建随机力学框架，重点关注生成函数。
• 应用概率生成函数、指数矩生成函数和阶乘矩生成函数作为演化变量。
• 利用随机过程的无穷小生成元定义时间演化，明确区分状态空间与演化算子。
• 采用组合正规排序技术求解非二元反应的演化方程，利用已知的生成函数恒等式。
• 对于二元反应，用Sobolev-Jacobi正交多项式替代McQuarrie的特征基，证明其构成生成元的一致特征基。
• 采用Gurappa-Panigrahi框架提出Sobolev-Jacobi多项式的新型定义，为推广至多物种系统提供可能。

实验结果

研究问题

• RQ1 如何构建一个随机力学框架，以统一化学反应系统中组合生成函数与连续时间马尔可夫链？
• RQ2 为何McQuarrie的特征基方法在二元化学反应中失效，如何修正？
• RQ3 Sobolev-Jacobi正交多项式能否作为二元反应系统无穷小生成元的有效且完备的特征基？
• RQ4 在演化半群和底层生成元方面，非二元与二元反应系统在结构上存在何种差异？
• RQ5 所提出的形式化方法在多物种或图重写随机转移系统中可推广到何种程度？

主要发现

• 本文利用组合正规排序技术，为所有六种基本单物种化学反应的概率分布时间演化提供了精确的闭式解。
• 对于二元反应，该框架通过用Sobolev-Jacobi正交多项式替代McQuarrie的特征基，解决了其方法中的根本性不一致问题，证明后者构成有效且完备的特征基。
• 二元反应系统的解被解析推导并明确呈现于命题6中，图2e和2f提供了说明性示例。
• Sobolev-Jacobi多项式的应用在正交多项式与化学反应动力学之间建立了新联系，暗示反应生成元中存在更深层次的代数结构。
• 引入了基于Gurappa-Panigrahi框架的Sobolev-Jacobi多项式的新定义，为将结果推广至多物种系统开辟了新途径。
• 该框架表明，矩生成函数和阶乘矩按明确定义的方程演化，其解可通过所提出的形式化方法获得，即使对于复杂反应类型亦然。