GROMACS 簇分析案例解读：如何联读最大簇、簇数与分布-模拟计算-科学指南针

【摘要】用纳米水滴演化案例理解 GROMACS 簇分析：先确认 cutoff，再联读最大簇、簇个数、平均簇大小和簇尺寸分布。

GROMACS 簇分析不是只画一条曲线，而是用统一的连通判据回答三个问题：母相还剩多少、体系分成了多少部分、离开母相的是单个分子还是小簇。对需要判断蒸发、团聚、解体或相分离过程的研究者而言，最大簇大小、簇个数、平均簇大小和簇尺寸分布应当放在同一时间窗口内解读。

下面以一个纳米级水滴在真空中受热的计算案例为边界，说明 gmx clustsize 的输出应怎样读。它提供的是一套阅读方法，不替代针对新体系的模型验证、参数选择或机理结论。

先给结论：四类输出各回答什么

簇指在选定距离阈值内彼此连通的一组粒子。判断簇之前，先明确 cutoff 的物理依据；对液体，常以径向分布函数第一峰后的第一谷作为候选。水的演示体系以 O-O 距离约 0.35 nm 为例，但这个数值不是其他分子或其他状态点的默认参数。

输出	直接回答的问题	不能单独证明什么
最大簇大小 Nmax(t)	主体连续相还剩多大	小簇的组成和数量
簇个数 Nc(t)	同一时刻有多少连通团	每个团的尺度分布
平均簇大小 <N>(t)	所有簇的典型平均尺度	是否仍存在一个主导母相
簇尺寸分布	单体、二聚体和更大簇各占多少	每次事件发生在何时

实用读法是：先看 Nmax 是否仍维持高平台，再用 Nc 和 <N> 判断体系是逐步散开还是同时形成很多碎片，最后用直方图确认碎片的尺寸组成。四项指向同一趋势时，解释更稳；若它们彼此矛盾，应回头核查 cutoff、周期性边界、分组对象和取样时间窗。

案例边界：一个纳米水滴的两种温度条件

演示模型由 1,000 个 SPC/E 水分子组成，水滴半径约 1.93 nm，置于 10 x 10 x 10 nm 的真空盒中，在 NVT 条件下比较 400 K 和 600 K 的 10 ns 演化。这里的盒子余量是分析前提之一：蒸气空间不足时，彼此靠近的气相分子可能被距离阈值重新连成表观大簇。

因此，这个案例能够支持的表述是“该设定下的簇输出呈现不同的演化节奏”。它不能直接推出所有水滴、所有力场、所有温控方式或任何其他材料体系都会以相同时间尺度解体。

先读最大簇：母相是否仍连续

最大簇曲线是最适合先看的总览图。高位平台通常意味着大部分分子仍属于同一个连通母相；缓慢下降说明主簇逐步损失成员；陡降到很小的数值，则提示原有连续母相已经不再占主导。曲线中的台阶高度还可帮助提出问题：一次减少 1 时，更接近单分子离开；一次减少多个成员时，需要结合轨迹快照和分布判断是否存在小簇脱离。

在这个演示中，400 K 条件下 Nmax 在 10 ns 内从 1,000 降到 961，仍以大簇为主体。600 K 条件下，主簇在约 0.1 ns 时间尺度内迅速下降到十几个成员的量级。两条曲线说明的是同一模型在这两组条件下的主簇保留差异；它们不是对“温度升高必然按某一固定速率蒸发”的普适定律。

怎样避免把曲线看得太满

1.确认比较对象使用同一个 cutoff。若每个条件各自改变阈值，曲线差异可能来自判据而非体系行为。

2.核对轨迹是否经过合理的 PBC 处理。跨边界的连续团被切开，会人为降低最大簇。

3.对多原子分子确认是否按整分子分组。把同一分子的原子拆开，可能让“粒子数”的含义发生变化。

4.用原子坐标或代表性帧检查异常转折。曲线负责提示位置，构型检查负责确认发生了什么。

再读簇个数与平均簇：体系是在变碎还是在重组

仅有 Nmax 下降，还不足以说明母相外的组成。若 Nc 上升，而 <N> 同时下降，常见的图像是一个大团逐步转化为更多、更小的连通团；若 Nmax 的变化不大但 Nc 波动，需要查看是否只是边缘出现短寿命小簇。

演示的高温条件中，簇个数一度达到 687 个，平均簇大小趋向 1 附近。这与最大簇的快速下降共同指向“主簇消失后，体系中存在大量小连通单元”的状态。低温条件中，主簇仍占主要部分，少量脱离粒子不会让 Nc、<N> 的变化具有同样的解体含义。

一个稳妥的分析顺序是把三条时间序列放在同一横轴上：标记 Nmax 的突变时间，再检查 Nc 是否同步增加、<N> 是否同步降低。若三个方向不同步，不要急于命名机理，而应检查是否有局部重聚、不同的取样间隔，或不合适的连通阈值。

再读分布：离开母相的是单体还是小簇

直方图补上了时间序列缺失的信息：它把选定时间窗内的簇按尺寸计数。单体峰高，说明独立分子占多数；二聚体和更大尺寸一侧的占比上升，说明不能只用“单分子逃逸”概括该窗口内的所有产物。分布必须与统计时间窗一起报告，因为早期、中期和后期的组成可能不同。

该案例统计后期窗口时，400 K 条件下单体约占 89.9%，二聚体约占 5.4%，更大簇约占 4.7%；600 K 条件下单体约占 71.3%，二聚体约占 15.9%，更大簇约占 12.9%。两组条件都以单体为主，但高温组的右尾更明显。这个比较支持“该案例高温条件下小簇比例更高”的描述；若要讨论具体脱离机制，仍需回到时间分辨轨迹和模型设定。

从输出到结论：一张联读检查表

观察到的组合	可提出的解释	下一步应核查
Nmax 高平台，Nc 接近 1	主体仍以一个连续大簇存在	边缘分子和短寿命小簇是否被漏掉
Nmax 下降，Nc 上升，<N> 下降	大簇正在转化为多个更小的连通团	cutoff 是否固定，PBC 是否正确
Nmax 突降，直方图单体/二聚体占比提高	主簇快速消失，小尺度产物增加	轨迹快照、取样时间窗和分组规则
Nmax 异常回升或长期偏大	可能存在重聚，也可能是判据或盒子问题	盒子余量、气相密度和 cutoff 敏感性

这张表的目的不是替代验证，而是把“看见一条曲线”转化为可复核的问题链。对于不同温度、不同液滴尺寸或不同材料体系，先固定比较口径，再讨论趋势，通常比直接比较单个数字更可靠。

报告结果时，把判据和窗口写在数字旁边

簇分析常见的误读来自“只报一个数字”。例如，报告某时刻 Nmax 为 50，却不说明这里的粒子单位、cutoff、PBC 处理和轨迹时间点，读者就无法判断这 50 代表连续母相、小团簇，还是某个边界处理下的表观连接。更可复核的写法是把每一项结果和分析口径放在一起：分析对象、距离定义、固定的 cutoff、是否按整分子计数、统计窗口与输出单位。

对于条件比较，建议同时保留三个层次的信息。第一层是时间序列，指出突变、平台或回升发生在哪里；第二层是该时间段的构型检查，确认曲线代表的结构变化；第三层是尺寸分布，说明变化后剩下的连通单元是什么。这样即使后续需要调整采样窗口或复算，也能明确哪些结论依赖于哪一层证据。