刘建川・实战 第 5 篇｜原子类型与力场参数：模拟的核心，看懂 [atomtypes] 块
上一篇我们搞定了非键作用的**参数开关——[nonbond_params]模块，而这一篇，我们聚焦它的参数源头**——[atomtypes] 原子类型块。
[atomtypes]是整个拓扑文件的**根基**，定义了所有原子的基础属性和力场参数，上一篇的非键参数、前几篇的原子信息，都依赖于这个模块。新手最容易忽略它的重要性，导致原子类型不匹配、参数报错，今天我们直击重点，快速吃透核心。
【14年实战经验】 [atomtypes] 是“参数母版”，看懂它，就能理解所有拓扑参数的由来，避开90%的参数不匹配报错， 新手不用死记参数，会核对、懂规则就够了。
一、核心认知：[atomtypes] 到底是什么？
核心作用：定义所有原子的**身份信息和力场参数，是TOP/ITP文件中所有模块（[atoms]、[nonbond_params]等）的参数基础，相当于原子参数字典**。
新手必记3点（不搞复杂理论）：
位置：通常位于TOP文件最开头（力场引用后），或ITP文件的[ moleculetype ]模块之前； 关联性：[atoms]模块的原子类型、[nonbond_params]的原子类型对，都必须来自这里的定义； 必要性：必需模块，缺失则模拟无法启动，参数错误则全体系报错。
二、重点参数拆解（新手只看核心，不记冗余）
[atomtypes]模块每一行对应一种原子类型，核心包含6列参数，重点掌握前5列，精准避坑，以amber力场为例：
1. 第一列：原子类型名称（type name）
原子的“身份标识”（如C3、O、N4），**大小写敏感**，必须与[atoms]、[nonbond_params]中一致（如[atoms]写C3，这里就不能写c3）。
2. 第二列：原子序号（at.num）
该原子在元素周期表的序号。
3. 第三列：原子质量（mass）
原子的相对原子质量（如C为12.01，O为16.00），力场固定，**新手严禁修改**，修改会导致模拟能量异常。
4. 第四列：原子电荷（charge）
原子的默认电荷（如H为0.418，O为-0.836），是静电作用的基础，与[atoms]模块的原子电荷需匹配。
5. 第五列：粒子类型（type）
把原子当作真实原子处理，此参数为：A，若是虚拟位点则为：D。
6. 第六列：范德华参数 σ（sigma）
与[nonbond_params]的σ一致，定义该原子范德华作用范围，力场固定，不修改。单位：nm。
7. 第七列：范德华参数 ε（epsilon）
与[nonbond_params]的ε一致，定义该原子范德华作用强度，力场固定，不修改。单位：kJ/mol。
重点提醒：所有参数均由力场固定，新手无需编写、无需修改。
三、实战衔接：与前后模块的核心关联（必看）
[atomtypes]是拓扑文件的“根基”，与前4篇模块紧密关联，理清关联才能避免报错：
**与[nonbond_params]（第4篇）**：[nonbond_params]的原子类型对（type1/type2），必须是这里定义的原子类型，否则提示“原子类型未找到”； **与[atoms]（第2篇）**：[atoms]模块中每个原子的“类型”，必须对应这里的原子类型，电荷需匹配，否则参数冲突； **与[bonds]/[angles]（第3篇）**：共价键、键角的参数，本质依赖这里的原子类型，原子类型错误会导致共价参数异常。
四、实战避坑：2个必核对+3个高频误区
必核对2点（核心避坑）
原子类型：[atoms]、[nonbond_params]中出现的原子类型，必须都在[atomtypes]中存在； 力场匹配：原子类型和参数，必须与所用力场的默认设置一致。
高频误区（新手必避）
随意修改原子类型名称：如把C3改成C4，导致与其他模块不匹配，模拟报错； 手动修改ε、σ或电荷：参数异常，导致非键作用、静电作用紊乱，体系崩溃； 忽略原子类型大小写：如C3和c3视为不同类型，提示“类型未找到”。
五、实战小结（4句口诀记牢）
atomtypes是根基，所有参数从这起； 原子类型要统一，大小写也需注意； 力场参数不修改，核对一致就可以； 关联模块无偏差，模拟顺利不费力。
