arXiv:2604.13213stat.MLcs.LG2026-04

用最优控制方法高效计算罕见事件的反应路径和速率。

Rare Event Analysis via Stochastic Optimal Control

论文配图:Rare Event Analysis via Stochastic Optimal Control
图 1 · 摘自论文原文
  • 将共谋函数估计转化为随机最优控制问题,通过反馈控制引导轨迹。
  • 在基准测试中,反应速率和平衡常数估计精度显著优于现有方法。
  • 适合研究生物分子构象变化、相变等罕见事件的科研人员。

诸如生物分子构象变化、相变和化学反应等罕见事件是许多物理系统行为的核心,但因其在无偏模拟中几乎不发生,计算研究极为困难。过渡路径理论(TPT)为这类事件提供了严格的统计框架:它刻画了两个指定稳态(反应物与产物)之间的反应轨迹集合,其中核心对象——共谋函数(committor function)——给出了系统接下来到达产物而非反应物的概率,编码了全部关键的动力学与热力学信息。本文提出一种新框架,将共谋函数估计建模为随机最优控制(SOC)问题。在此框架下,共谋函数定义了一个反馈控制策略(与其对数梯度成正比),可主动引导轨迹进入反应区域,从而实现反应路径的高效采样。为求解由此产生的击中时间控制问题,我们设计了两种互补目标:直接反向传播损失与具有首阶最优性保证的无偏值匹配损失。此外,针对受困于中间盆地的多重稳定性问题,引入一种保持反应流同时降低有效能垒的新采样过程。在多个基准系统上,该框架在共谋函数估计、反应速率和平衡常数方面均显著优于现有方法。

原文摘要 · Abstract (English)

Rare events such as conformational changes in biomolecules, phase transitions, and chemical reactions are central to the behavior of many physical systems, yet they are extremely difficult to study computationally because unbiased simulations seldom produce them. Transition Path Theory (TPT) provides a rigorous statistical framework for analyzing such events: it characterizes the ensemble of reactive trajectories between two designated metastable states (reactant and product), and its central object--the committor function, which gives the probability that the system will next reach the product rather than the reactant--encodes all essential kinetic and thermodynamic information. We introduce a framework that casts committor estimation as a stochastic optimal control (SOC) problem. In this formulation the committor defines a feedback control--proportional to the gradient of its logarithm--that actively steers trajectories toward the reactive region, thereby enabling efficient sampling of reactive paths. To solve the resulting hitting-time control problem we develop two complementary objectives: a direct backpropagation loss and a principled off-policy Value Matching loss, for which we establish first-order optimality guarantees. We further address metastability, which can trap controlled trajectories in intermediate basins, by introducing an alternative sampling process that preserves the reactive current while lowering effective energy barriers. On benchmark systems, the framework yields markedly more accurate committor estimates, reaction rates, and equilibrium constants than existing methods.

罕见事件最优控制反应路径共谋函数

Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。