[论文解读] The Berry phase and the pump flux in stochastic chemical kinetics

Nikolai A. Sinitsyn, Ilya Nemenman|arXiv (Cornell University)|Dec 11, 2006

Spectroscopy and Quantum Chemical Studies被引用 7

一句话总结

本文提出了一套几何贝里相位框架，用于描述经典随机化学动力学中的泵电流，特别针对酶周转的两态马尔可夫模型。通过推导在周期性速率调制下粒子跃迁的完整计数统计，识别出一种纯粹几何的通量贡献——与经典电流不同，该贡献源于参数空间中的绝热演化，其预测结果可被单分子实验测量。

ABSTRACT

We study a classical two-state stochastic system in a sea of substrates and products (absorbing states), which can be interpreted as a single Michaelis-Menten catalyzing enzyme or as a channel on a cell surface. We introduce a novel general method and use it to derive the expression for the full counting statistics of transitions among the absorbing states. For the evolution of the system under a periodic perturbation of the kinetic rates, the latter contains a term with a purely geometrical (the Berry phase) interpretation. This term gives rise to a pump current between the absorbing states, which is due entirely to the stochastic nature of the system. We calculate the first two cumulants of this current, and we argue that it is observable experimentally.

研究动机与目标

• 开发一个通用的理论框架，用于分析在周期性驱动下的经典随机马尔可夫过程中的粒子通量。
• 在最小的米氏- Menten 酶模型中，识别并量化一个几何（贝里相位）贡献到净通量。
• 提供一种方法，不仅计算平均通量，还计算完整计数统计，包括方差和高阶累积量。
• 证明泵电流的实验可观测性，特别是在单分子或低浓度条件下。
• 将量子节拍器中的几何相位概念拓展至经典随机化学动力学，扩展其适用范围。

提出的方法

• 将随机化学动力学形式化为具有吸收态（底物和产物）和中心态（酶-底物复合物）的连续时间马尔可夫链。
• 在 master 方程中引入计数场，以推导粒子跃迁的完整计数统计的生成函数。
• 对时变跃迁率应用绝热近似，假设速率常数的缓慢周期性调制。
• 通过提取修正生成矩阵的最低实部特征值，获得累积量生成函数。
• 从时变生成器的基态本征态计算贝里曲率，从而实现对泵电流的几何解释。
• 利用驱动频率的微扰理论，推导泵电流前两个累积量（平均值和方差）的显式表达式。

实验结果

研究问题

• RQ1 在具有周期性驱动的经典随机系统中，能否识别出纯粹几何（贝里相位）的粒子通量贡献？
• RQ2 在最小的酶促模型中，泵电流如何依赖于周期性速率调制的相位和振幅？
• RQ3 在大小和可探测性方面，经典电流与几何泵电流之间存在何种关系？
• RQ4 泵电流及其涨落能否在单分子或低酶浓度条件下被实验分辨？
• RQ5 贝里相位形式化在多大程度上可推广至其他马尔可夫化学动力学系统？

主要发现

• 泵电流完全源于完整计数统计中的几何（贝里相位）项，与经典电流无关，且与一个周期内参数空间中封闭区域的面积成正比。
• 泵电流按 $ J_{\text{pump}} \sim \omega / \langle K \rangle $ 缩放，其中 $ \omega $ 为驱动频率，$ \langle K \rangle $ 为平均总速率。
• 泵电流的方差因额外因子 $ \delta k_{1} \delta k_{-2} / (k_{1} k_{-2}) $ 而被抑制，使其比经典电流更难探测。
• 在典型酶参数下（$ k_2 \sim 10 \, \text{s}^{-1} $，$ K_M = 1 \, \text{mM} $，$ \delta k / k \sim 10\% $），泵电流达到 $ \sim 10^{-2}/T_0 $，其中 $ T_0 $ 为驱动周期。
• 对于持续数天的实验，泵电流与经典电流的标准差相当，表明在长期测量中具有可探测性。
• 在约 $ 10^6 $ 个酶分子下，泵电流及其方差可在一分钟内被观测到，表明在实际条件下具备实验可行性。