Yau Awards Archive 2020 — 2025

Topics选题方案 ／ 化学奖 ／ C07

C07

笔记本算力预算下水中小分子扩散系数的力场可分辨性：有限尺寸校正后两个力场差异何时才在统计上成立

推荐优先级 中分族 分子动力学资源需求 纯 CPU，多核后台连跑数周技能取向 命令行 MD + 时间序列统计失败风险 中高

1 · 研究问题

在纯 CPU 笔记本可承受的模拟预算内（盒子 ≤ 3,000 原子、累计轨迹 ≤ 300 ns），经 Yeh–Hummer 有限尺寸校正后，两个常用力场组合（如 GAFF2/TIP3P 与 OPLS-AA/TIP4P-2005）给出的小分子水中自扩散系数之差，需要多长的累计轨迹与多少条独立重复才能在 95% 置信区间上与零区分开？该阈值与分子大小、盒子尺寸如何标度？

2 · 研究背景与空白

自扩散系数 D 由均方位移（mean squared displacement, MSD）长时线性段斜率给出：D = lim MSD/(6t)，是最容易与实验对照的液相输运性质，也是力场验证的常规靶标。它有两个众所周知的陷阱：其一，周期边界下的流体力学自相互作用使 D 系统性偏小，需 Yeh–Hummer 校正 D∞ = D_PBC + kB T ξ/(6πηL)（ξ ≈ 2.837297，L 为盒长，η 为剪切黏度），小盒子里可达 20–40%；其二，MSD 斜率的统计误差只以 √t 收敛，故"两个力场差 10%"在短轨迹上常常是噪声。

已有工作走到哪一步：Yeh 与 Hummer（2004）确立有限尺寸校正公式；大量力场评测论文报告 GAFF/OPLS/CHARMM 配不同水模型下小分子 D 与实验的偏差（典型 20–50%）；TIP3P 水自扩散系数比实验高约 2 倍是公认缺陷。

空白在：这些工作几乎都在集群上跑几百纳秒到微秒，报告的是"力场 A 比力场 B 更接近实验"这类结论，没有人从相反方向问：在一台笔记本的算力上限内，这类结论到底还能不能成立。这是把算力预算显式写进研究设计的"可分辨性/统计功效"问题，属于换问题方向：资源约束不是妥协而恰是研究对象，正负结论都对同类小资源使用者有直接指导意义。

3 · 可检验假设

H1 在 2,500 原子量级的盒子中，单条 20 ns 轨迹给出的溶质 D 的相对标准误差 ≥ 15%；要把两个力场组合 20% 的 D 差异在 95% 置信区间上分辨开，需要累计 ≥ 100 ns（≥ 5 条独立重复）。

H2 有限尺寸校正项在本课题可及的盒子尺寸（L ≈ 3–4 nm）上贡献 ≥ 15% 的 D，其量级大于同一体系两个力场之间的差异；即不做有限尺寸校正的力场对比结论会被系统性污染。若 H2 被否证（校正项远小于力场差异），说明小盒子对力场排序无害，这对小资源使用者是一个正面结论。

4 · 量化验收标准

  1. 方法学校验（硬门槛）：用自建管线模拟纯 TIP3P 水（约 2,000–3,000 原子），计算其自扩散系数并与该水模型的公认文献值对比（TIP3P 在 298 K 的 D 约为 5.5–6.0 × 10⁻⁵ cm²/s 量级，须以原文核实），经有限尺寸校正后偏差 ≤ 10%；同时复现同一模型的密度（偏差 ≤ 1%）。此条不过关意味着系综设置、时间步长、静电处理或 MSD 拟合区间有误，后续全部力场对比结论无效。
  2. 统计口径预先写死：D 由 MSD 在明确指定的线性拟合区间（例如 20–80% 轨迹长度，区间选择规则事先冻结并公布）上最小二乘得到；误差棒来自 ≥ 5 条独立初速度重复的样本标准误，并额外用块平均（block averaging）交叉验证误差估计；所有力场差异的显著性用自助重抽样给 95% 置信区间，不用单条轨迹的拟合残差冒充误差棒（这是该领域最常见的错误）。
  3. 收敛性检查（必需）：报告能量守恒/温度控制的漂移量、NPT 平衡阶段的密度时间序列、以及 D 对拟合区间选择的敏感性（至少 3 种区间）。报告 MSD 的对数斜率是否达到 1.0（确认进入扩散区而非亚扩散区）。
  4. 有限尺寸校正必须自行实现并验证：至少用 2 个不同盒长（如 L ≈ 3.0 nm 与 3.8 nm）做外推对照，检验 D_PBC 对 1/L 的线性关系斜率与 Yeh–Hummer 公式预测值一致（偏差 ≤ 25%）；黏度 η 由同一模拟的 Green–Kubo 或文献值给出（取法事先写定）。
  5. 交付一张"可分辨的 D 相对差异 vs 累计轨迹长度"的功效曲线，并在图上标出本机 1 周 / 4 周预算线。
  6. 全部拓扑文件、mdp/输入文件、分析脚本与原始 MSD 数据开源，一条命令可重跑。

5 · 数据与工具

用途 来源 / 工具
MD 引擎（主） GROMACS（LGPL，conda-forge 或源码编译）。CPU 并行成熟，是本课题唯一推荐的主引擎
MD 引擎（不推荐作主力） OpenMM。官方与社区反复指出其为 GPU 设计、CPU 平台性能不佳（有报告：48 核节点上 95,000 原子体系仅约 0.2 ns/day）。本课题仅在体系搭建/力场转换环节使用其生态（如 openff-toolkit），不用它跑生产轨迹
力场与拓扑 GAFF2（经 acpype/AmberTools 或 OpenFF 工具链生成）、OPLS-AA（LigParGen 在线服务，可用性与许可须核实）、OpenFF Sage（开源，openff-toolkit 直接生成）。水模型 TIP3P / TIP4P-2005 均随 GROMACS 分发
实验对照 公开的小分子水中自扩散系数实验值（如甲醇、乙腈、尿素等，见物性手册与文献表格）。仅用于校验与对比，不计入本项目的数据贡献
分析 GROMACS gmx msd（内置 MSD 与 Einstein 拟合）、MDAnalysis 或 MDTraj（Python 侧独立重算，用于交叉校验）、自写 Yeh–Hummer 校正与块平均
算力 纯 CPU，这是全部路线中算力最紧的一条，必须在第 4 周实测本机数字。可用的公开标尺：8 核 CPU 上 41,175 原子体系约 0.661 ns/day。按原子数近似线性外推，2,500–3,000 原子体系在 8 核上约 9–11 ns/day（此外推必须由本机基准测试确认，不得直接采信）。据此，4 周后台连跑可得约 250–300 ns 累计轨迹，恰好落在"5 条 × 20 ns × 两个力场 + 校验"这一预算内。超出范围的方案（明确列出）：任何 > 10,000 原子的盒子、微秒级轨迹、自由能微扰/伞形采样、蛋白-配体体系——一律不做

6 · 方法路径

  1. 装 GROMACS（优先用 conda-forge 预编译包，避免编译链问题），跑通官方水盒子教程，在本机实测 ns/day 并据此重新确定盒子大小与轨迹长度（第 1–4 周，此步的输出直接决定后续方案规模）。
  2. 做纯 TIP3P 水的方法学校验硬门槛：密度 + 自扩散系数（含有限尺寸校正）对照文献值。
  3. 核实工具能力边界：确认 GROMACS 的 gmx msd 拟合区间默认取法、是否内置有限尺寸校正（几乎肯定不内置，须自行实现）、Green–Kubo 黏度是否可直接输出。这三项须在第 8 周前查清。自行实现的校正必须用双盒长外推验证——这一步本身可作为项目的方法学贡献写入。
  4. 搭建溶质体系（1 个小分子 + 水，两套力场组合），做 NPT 平衡与生产轨迹，每套 ≥ 5 条独立重复。
  5. 计算 D_PBC、施加有限尺寸校正、给出块平均与重复间误差棒；做拟合区间敏感性分析。
  6. 构造功效曲线：对轨迹长度做子采样（取前 2/5/10/20 ns），画"可分辨的相对差异 vs 累计轨迹长度"。
  7. 独立交叉校验：用 MDAnalysis 在 Python 侧完全独立地重算 MSD 与 D，确认与 gmx msd 结果一致（偏差 ≤ 3%）；再用一个盒长完全不同的体系确认标度关系。

7 · 新颖性边界

  • 本课题不声称提出新力场或新校正方法，不声称发现 TIP3P 扩散过快（这是公认缺陷），不声称给出比文献更准的 D 值——本课题的算力预算注定做不到这一点，任何相反的声称都会被评委在三十秒内问倒。
  • 已有工作具体完成了什么：Yeh 与 Hummer（J. Phys. Chem. B, 2004）给出周期边界下的有限尺寸校正公式；大量力场评测论文在集群上以几百纳秒到微秒的轨迹报告小分子 D 与实验的偏差（典型 20–50%）并对力场排序。这些工作的算力前提与本课题完全不同。
  • 本项目的贡献是换问题方向且为主结论：把算力预算显式写进研究设计，回答"在纯 CPU 笔记本上，力场对比结论的可分辨性下限在哪里"，输出一条功效曲线与一条可操作判据（"累计轨迹短于 X ns 时，不要报告力场之间 Y% 以内的 D 差异"）。附带贡献是自行实现并双盒长验证的有限尺寸校正流程。
  • 为什么有价值：小资源使用者（中学、本科课程、小课题组）大量在做这类对比，却普遍用单条短轨迹的拟合残差当误差棒，从而报出不存在的差异。本课题给出的是他们可以直接照用的门槛。
  • "差异不显著"同样是有效结论——而且在本课题中是预期的主结论之一：明确写出"在笔记本预算内两个力场的 D 差异不可分辨"，并用重复间误差棒证明这不是因为做得不够仔细，而是统计精度的物理上限。

8 · 决策门槛（go / no-go）

  • 第 4 周末（最关键的一道门）：本机 GROMACS 基准测试必须给出 ≥ 5 ns/day @ 2,500 原子。若低于此值，立即降级：其一，把体系缩到 1,200–1,500 原子（约 400 个水分子 + 溶质），代价是有限尺寸校正项更大——但这反而强化了 H2 的检验，主结论框架不受损；其二，把研究对象从"溶质扩散"换成"纯水自扩散"（水分子数多，统计量天然充足，同样长度的轨迹误差棒小得多），功效曲线的主结论完全保留。
  • 第 8 周末：纯 TIP3P 水的方法学校验必须通过（D 偏差 ≤ 10%、密度偏差 ≤ 1%）。若失败，先查静电处理（PME 参数）、时间步长与约束设置；仍失败则降级为只做"同一力场下不同盒长的有限尺寸标度检验"（这部分不需要与文献 D 值对齐，只需内部自洽），H2 完全保留，H1 改为"同一力场两个独立实现之间的可分辨性"。
  • 第 12 周前必须查清（不要边做边发现）：OPLS-AA 拓扑的获取途径（LigParGen 的可用性与使用条款）。若不可用，改用 OpenFF Sage 与 GAFF2 作为两个对比力场——方案的主结论不依赖于具体是哪两个力场。
  • 已知困难即研究内容：MSD 线性拟合区间的选择带有主观性。把区间选择的敏感性分析与"如何给出可复现的区间判据"写成一节，本身就是交付物。
  • 预算裁剪顺序：先砍第二个盒长（放弃双盒长外推验证，改引用 Yeh–Hummer 原文系数），再砍重复次数（5 → 3，误差棒变宽但功效曲线仍可画），再砍第二个力场（只做单力场的功效分析），最后砍 Green–Kubo 黏度（改用文献 η 值）。砍到只剩"纯水 + 单力场 + 3 重复 + 功效曲线"时主结论仍成立。
  • 选择前提：仅在学生有一台可连续后台运行数周的机器时启动。若机器需频繁携带、休眠、断电，改选可分批可中断的 C05 或 C06。