Skip to content

配对持久性与驻留时间

本文是一份面向初学者的 配对持久性 分析指南——MolPy 怎么衡量两个粒子维持"成键"（在距离截断内）多长时间，又怎么把这个时长转成驻留时间相关函数。常见用途包括水中的氢键动力学和电解质里的离子对寿命。除了时间相关函数的基本概念之外不需要别的背景；如果已经看过扩散与离子输运指南，那基础知识已经够了。

底层靠 Rust（molrs）完成每对粒子逐帧的簿记；MolPy 层负责提取坐标和盒子尺寸，返回类型化结果。

本文使用的约定

  • 配对指一个参考粒子 \(i\) 和一个伙伴粒子 \(j\)。
  • 如果一对粒子在某帧的最小镜像距离在截断范围内，就认为它们处于成键状态。可以使用两个截断：内截断 \(r_0\)（形成）和外截断 \(r_1\ge r_0\)（断裂）。
  • \(\langle\cdots\rangle_t\) 表示对时间原点取平均；\(\tau\) 是滞后时间。
  • 单位：长度 Å，时间 ps。相关函数 \(C(\tau)\) 无量纲。

1. 寿命为什么需要专门的工具

扩散系数和电导率（见输运指南）能告诉你物质移动了多远，但没法告诉你某种关联维持了多久。然而很多性质——质子转移速率、离子对稳定性、溶剂化壳层交换——取决于接触的寿命，而非体相迁移率。

处理这类问题的难点在于，接触很少是干脆利落的：两个粒子到了截断附近，可能每皮秒里反复进出好多次，但配对并没有真正断裂。一个合格的寿命度量必须明确什么才算"同一键持续存在"。这个选择就是本文要讨论的核心。

2. 存活相关函数

定义配对 \((i,j)\) 在时间 \(t\) 产生、在时间 \(t+\tau\) 被观测的存活指示函数：

\[ S_{ij}(t, t+\tau) = \begin{cases} 1 & \text{配对在 } t+\tau \text{ 时被视为仍存活},\\ 0 & \text{否则}. \end{cases} \]

持久性（驻留时间）相关函数对所有参考粒子和时间原点求平均，并按参考粒子归一化：

\[ \boxed{\;C(\tau) = \Big\langle\,\frac{1}{N_i}\sum_i\sum_j S_{ij}(t,\,t+\tau)\,\Big\rangle_t\;} \]

两个极限值可以帮助理解：

  • \(C(0)\) 是平均配位数——落在形成截断内的伙伴粒子平均数目。（\(\tau=0\) 的时候每个刚产生的配对自然都是存活的。）
  • \(C(\tau)\) 随键断裂而衰减；衰减时间就是驻留时间。用 \(C(\tau)\approx C(0)\,e^{-\tau/\tau_\text{res}}\) 拟合（或者直接对曲线积分）就能得到一个单一的寿命值。

接下来要确定的只有 \(S_{ij}\) 内部怎么判存活。

3. 存活的三种定义

文献里有多种存活规则；MolPy 实现了三种定义明确且广泛使用的。它们共享相同的产生条件——配对只在 \(t\) 时刻处于内截断 \(r_0\) 范围内时才算产生——区别在于如何判断已产生的配对到 \(t+\tau\) 时是否还存活。

3.1 连续存活

最严格的规则（Rapaport, 1983）：配对必须从 \(t\) 到 \(t+\tau\) 的每一帧都保持在存活截断 \(r_1\) 之内。一旦离开，对于该原点而言就永远死亡：

\[ S^\text{cont}_{ij}(t,t+\tau) = \prod_{s=0}^{\tau}\Theta\!\big(r_1 - r_{ij}(t+s)\big), \]

其中 \(\Theta\) 为阶跃函数。结果称为连续驻留时间——通常很短，因为任何一次瞬时离轨都会终结键合。使用单一截断（\(r_1=r_0\)）即可得到经典定义。

3.2 间歇存活

最宽松的规则（Luzar & Chandler, 1996）：配对在 \(t+\tau\) 时刻只要处于 \(r_1\) 范围内就算存活，中间是否离开过无所谓。允许重新形成：

\[ S^\text{int}_{ij}(t,t+\tau) = \Theta\!\big(r_1 - r_{ij}(t+\tau)\big). \]

结果称为间歇相关函数，衰减速度体现的是包含重新穿越的结构寿命。连续和间歇两种相关函数给出了真实值的上下界，它们的比值则反映了键的振荡程度。

3.3 稳态图像（SSP）

一种折中方案（Laage & Hynes, 2008），在不丢弃真正存活事件的前提下抑制振荡效应。使用两个截断：配对在内截断 \(r_0\) 范围内产生（"稳态反应物"），之后只要保持在外截断 \(r_1\) 范围内就算存活（"稳态产物"）；必须穿越外截断才被判定为断裂：

\[ S^\text{SSP}_{ij}(t,t+\tau) = \Theta\!\big(r_0 - r_{ij}(t)\big)\prod_{s=0}^{\tau}\Theta\!\big(r_1 - r_{ij}(t+s)\big), \qquad r_1\ge r_0. \]

\(r_0\) 和 \(r_1\) 之间的空隙相当于一个缓冲带：粒子游荡到刚好超出接触距离时不会立即被视为离开。对于离子对分析，SSP 是推荐的默认方法。

与 tame 定义的关系

这与 tame 的 persist 配方对应。tame 命名了两种定义：IMM（Impey–Madden–McDonald, 1983），使用单一截断加上一个忽略短于 \(t^*\) 离轨的容差时间；以及 SSP（Laage–Hynes）。MolPy 的 continuous/intermittent/ssp 使用明确的、不带时间容差的判据覆盖了相同的物理含义；IMM 容差时间变体有意未予复现（其已发布的实现中存在 timestep 越界问题）。

4. 使用 Persist

Persist 接受 "t1,t2:method:r0[,r1]" 格式的标签：

  • t1,t2——参考粒子和伙伴粒子的原子类型（例如阳离子和阴离子）。
  • method——continuous、intermittent 或 ssp。
  • r0[,r1]——内截断（以及可选的外截断），单位 Å。只给一个值时设置 \(r_1=r_0\)。
from molpy.compute import Persist

p = Persist(
    tags=[
        "3,4:ssp:3.0,4.0",        # 阳离子-阴离子，稳态图像
        "1,1:continuous:3.5",      # 同种粒子，单一截断连续模式
    ],
    max_dt=30.0,    # 最长滞后时间，ps
    dt=0.01,        # 时间步长，ps
)
result = p(trajectory)

C = result.correlations["3,4:ssp:3.0,4.0"]
C[0]            # 平均配位数（每个参考离子的伙伴数）
result.time     # 滞后时间 tau，ps

当 t1 == t2 时，同种粒子的自配对（\(i=j\)）会被自动排除。

5. 从持久性到配对扩散

持久性不止是描述性工具——它还允许按键合状态分解输运系数。将输运指南 §2 中的区分扩散相关函数与存活函数 \(f(r_{ij};s)\) 结合，即可得到扩散系数的配对贡献：1

\[ D^\text{d,pairing}_{\alpha\beta} = \lim_{\tau\to\infty}\frac{1}{6\tau N} \sum_i\sum_{j\ne i}\big\langle\Delta\mathbf{r}_{i,\alpha}(\tau)\cdot\Delta\mathbf{r}_{j,\beta}(\tau)\,f(r_{ij};s)\big\rangle, \]

这个式子只计入仍保持关联的伙伴的位移（"对内"通道）。它区分了配对离子和自由离子的扩散方式——对于离子对是否携带净电荷这个问题直接相关。仅在持久性计数已收敛并且 MDC 随时间呈线性时才能解读结果。

6. 陷阱检查清单

  1. 用单一截断而振荡又严重 → continuous 寿命会坍缩到帧间隔级别。改用带缓冲的 ssp（\(r_1 > r_0\)），或者用 intermittent。
  2. 截断不参考 RDF，随意取值 → \(r_0\) 应该在 \(g(r)\) 的第一配位壳层极小值处选取；随意选截断会让寿命失去意义。
  3. max_dt 小于实际寿命 → \(C(\tau)\) 在窗口内永远不会衰减，读不出驻留时间。
  4. 采样太稀疏 → 两种存活乘积都会遗漏快速的重新穿越；采样密度要足以解析接触动力学。
  5. 跨定义比较 → 连续、间歇和 SSP 这三种定义本身就给出不同的数值；始终要报告用的是哪一种。

7. 参考文献

  • D. C. Rapaport, Mol. Phys. 50, 1151 (1983) — 连续与间歇氢键相关函数。
  • A. Luzar, D. Chandler, Nature 379, 55 (1996); Phys. Rev. Lett. 76, 928 (1996) — 间歇相关与氢键动力学。
  • A. Luzar, J. Chem. Phys. 113, 10663 (2000) — 持久性/驻留时间定义。
  • R. W. Impey, P. A. Madden, I. R. McDonald, J. Phys. Chem. 87, 5071 (1983) — 带容差时间的驻留时间（IMM）。
  • D. Laage, J. T. Hynes, J. Phys. Chem. B 112, 14230 (2008) — 稳态图像（SSP）。

参见

  1. H. Gudla, Y. Shao et al., J. Phys. Chem. Lett. 12, 8460 (2021) — 通过持久性加权的区分相关函数计算扩散的配对贡献。 ↩