刘建川·GROMACS TOP实战 第2篇｜原子与电荷：模拟的起点，看懂[atoms]块

大家好，我是刘建川（网名刘十三），14年深耕GROMACS分子模拟领域，始终专注分子模拟实战教学。

在TOP/ITP文件所有模块里，**[atoms]块是真正的“核心中的核心”——它定义分子里每一个原子的身份、类型、电荷、质量，是GROMACS认识分子的第一步，也是90%拓扑报错的根源**。

这一篇我们不讲虚的，直接带你看懂[atoms]每一列含义、怎么核对、怎么避错，让你从此不再怕拓扑报错。

【14年实战经验】 [atoms]对得上，模拟就稳一半；[atoms]对不上，一切都白搭。 只要原子数、电荷、名称和结构文件一致，大部分拓扑崩溃直接解决。

一、先搞懂：[atoms]块到底是什么？

**[atoms]是分子的“身份证表”**，每一行代表一个原子，每一列代表一项属性。GROMACS靠这张表知道：这是什么原子、带多少电荷、属于哪个残基、质量多少。

1. 位置：一定出现在ITP文件里，紧跟**[moleculetype]**之后
2. 地位：**必须有，不能缺，不能乱**
3. 规则：**行数 = 原子总数；列数固定，顺序固定**

二、[atoms]标准列详解（逐列看懂）

标准[atoms]一共有7列，我给你用最直白的语言讲清楚：

1. **原子序号**：从1开始递增，必须连续、不重复、不颠倒
2. **原子类型**：力场定义的原子种类（如C1、CA、OH、N1等）
3. **残基序号**：同一个分子一般填1即可
4. **残基名**：如UNK、LIG、SOL、ALA
5. **原子名**：必须和GRO/PDB完全一致（如C1、O2、H3），如果不一致，会提示警告，然后丢弃GRO/PDB名称，采用ITP里面的名称。
6. **电荷分组**：一般填0即可
7. **电荷**：原子的局部电荷，决定静电相互作用
8. **质量**：各个原子的相对原子质量

简单记：**序号、类型、残基号、残基名、原子名、分组、电荷、质量**（有些会省略电荷和质量，则会自动去force field.itp文件里面去找）。

三、实战最关键：3个必须核对的点

这3点对了，模拟基本不会因为[atoms]崩溃。

1. 原子行数 = 结构文件（GRO/PDB）原子总数
2. 原子名必须完全一样（大小写、数字、符号都要一致）
3. 分子总电荷尽量为0（带电体系除外）

不一致就一定会报错：原子不匹配、分子找不到、体系爆炸。

四、新手最常踩的3个坑

误区1：原子type大小写不一致

GROMACS大小写敏感，C1 ≠ c1，必须完全一致。

误区2：原子行数对不上

多一行、少一行，直接启动失败。

误区3：电荷乱改

电荷不匹配会导致体系瞬间崩溃、温度飙升、原子飞散。

五、实战小结：[atoms]块4句口诀

1. [atoms]是分子身份证
2. 行数=原子数，一列不能少
3. 原子名必须和结构文件完全一致

📌 下一篇预告： 刘建川·GROMACS TOP实战 第3篇｜键、角、二面角：分子骨架怎么看懂

我开设的GROMACS 线下小班将于4.10-12号正式开课，14年实战经验手把手教学，打通GROMACS全流程实操，具体课程详情、学习大纲及报名方式，点击下方进一步了解，或者wx:simuhomesci

报名详情：https://www.fangzhenxiu.com/train/1275461229226983424/

注：后续想学习GROMACS哪方面的内容，欢迎大家在告诉我！