刘建川・实战 第6篇|体系定义：看懂 [system] 与 [molecules] 配置块

上一篇我们搞定了拓扑文件的**参数根基——[atomtypes]模块，明确了所有原子的身份和参数；这一篇，我们聚焦体系配置的总开关**——[system] 与 [molecules] 配置块。

这两个模块是TOP文件的**收尾核心**，负责定义整个模拟体系的基本信息和分子组成，前面所有模块（[atoms]、[nonbond_params]等）的参数，最终都要通过这两个模块整合，形成完整模拟体系。新手最容易在这里出错（如分子数量不匹配、名称不一致），导致模拟无法启动，今天直击重点、精准避坑。

【14年实战经验】 [system] 定体系，[molecules] 配分子，二者协同，才能让所有参数“落地生效”， 新手不用复杂操作，重点核对“名称一致、数量匹配”，就能避开核心报错。

一、核心认知：两个模块的核心作用（直白易懂）

[system] 和 [molecules] 仅存在于TOP主文件中，分工明确、缺一不可，核心作用一句话分清：

1. [system] 模块：定义模拟体系的“基本信息”

核心作用：给模拟体系“命名”、标注基本属性，无复杂参数，仅作标识和基础配置，新手无需修改核心内容。

新手必记2点：

1. 位置：通常位于TOP文件末尾，[molecules]模块之前；
2. 核心内容：仅1行，可任意书写，或者默认配置即可，无需改动。

2. [molecules] 模块：定义模拟体系的“分子组成”

核心作用：配置体系中所有分子的名称和数量，是连接TOP文件（参数）和结构文件（GRO/PDB）的关键，也是报错重灾区。

新手必记2点：

1. 位置：TOP文件最后一行，是整个拓扑文件的“收尾模块”；
2. 核心要求：分子名称、数量，必须与结构文件（GRO）、ITP文件的分子名称完全一致。

二、重点拆解：两个模块的实操要点（新手只看核心）

1. [system] 模块（极简实操，无需修改）

标准格式（1行搞定）： 如：Protein-Ligand-SOL

实战提醒：名称可自定义（如“DNA-Water”），但建议与体系实际组成一致，方便后续查看。

2. [molecules] 模块（核心重点，必核对）

标准格式（每一行对应一种分子）： Protein 1 LIG 1 SOL 1000 IONS 2

参数解读（核心必看）：

1. 第一列：分子名称（如Protein、LIG、SOL），必须与ITP文件的[ moleculetype ]模块中“molname”一致，与结构文件（GRO）中的分子名称一致；
2. 第二列：分子数量（如1、1000），必须与结构文件（GRO）中该分子的实际数量、顺序一致，不能多、不能少。

重点提醒：分子名称大小写敏感（如SOL和sol视为不同分子），数量不匹配会直接提示“分子数量不一致”，模拟无法启动。

三、实战衔接：与前后模块的核心关联（必看）

这两个模块是拓扑文件的“整合核心”，与前5篇模块紧密关联，理清关联才能避免报错：

1. **与[atomtypes]/[nonbond_params]（第4、5篇）**：通过分子名称关联ITP文件，间接调用原子类型和非键参数，确保参数统一；
2. **与[atoms]（第2篇）**：[molecules]的分子数量，需与[atoms]模块的原子总数、结构文件的原子总数匹配；
3. **与ITP文件**：[molecules]的分子名称，必须与对应ITP文件中[ moleculetype ]的“molname”完全一致，否则无法引用ITP参数。

四、实战避坑：3个必核对+3个高频误区

必核对3点（核心避坑）

1. 名称一致：[molecules]的分子名称，与ITP的molname、结构文件的分子名称完全一致；
2. 数量匹配：[molecules]的分子数量，与结构文件中该分子的实际数量、顺序一致；
3. 无多余分子：不添加体系中不存在的分子（如未加离子，却写IONS 2），否则报错。

高频误区（新手必避）

1. 分子名称大小写错误：如SOL写成sol，导致无法匹配ITP文件，提示“分子未定义”；
2. 分子数量不匹配：如结构文件有1000个水分子，[molecules]写SOL 999，提示“原子数量不一致”；

五、实战小结（4句口诀记牢）

1. system定体系名称，可随意更改；
2. molecules配分子，名称数量要匹配；
3. 关联ITP和结构，大小写别大意；
4. 核对无误再运行，模拟顺利不费力。