刘建川・GROMACS TOP 实战 第 3 篇｜键、角、二面角：分子骨架怎么看懂
大家好，我是刘建川（网名刘十三），14年深耕GROMACS分子模拟领域，始终专注分子模拟实战教学。
上一篇我们吃透了[atoms]块，搞懂了原子的“身份信息”；这一篇我们聚焦分子的**骨架结构——[bonds]、[angles]、[dihedrals]三大模块，它们共同决定分子的空间形态，也是拓扑报错的高频重灾区**。
新手常忽略这三个模块，觉得“参数默认就行”，结果导致分子结构扭曲、原子重叠、模拟崩溃。这一篇我们直击核心，讲懂每个模块的作用、参数含义和避错技巧，帮你稳住分子骨架。
【14年实战经验】 原子是“零件”，键、角、二面角是“连接件”，零件对了，连接件出错，分子散架、异形； 三大模块的参数的核心是“匹配力场、贴合结构”，无需乱改，但必须会核对。
一、先搞懂：三大模块的核心作用（一句话分清）
[bonds]、[angles]、[dihedrals]均位于ITP文件中，紧跟[atoms]模块之后，三者分工明确、缺一不可，共同维持分子的稳定空间结构：
**[bonds]（键模块）**：定义原子间的共价键，控制键长，相当于“分子的骨架连杆”，决定原子间的距离； **[angles]（角模块）**：定义三个相连原子的夹角，控制键角，相当于“骨架的弯折角度”，决定分子的平面形态； **[dihedrals]（二面角模块）**：定义四个相连原子的扭转角度，控制分子的扭转形态，相当于“骨架的旋转关节”，决定分子的空间折叠方式。
简单记：**键控长、角控弯、二面角控扭转**，三者协同，分子才不会“散架”“扭曲”。
二、逐模块详解：核心参数+实战重点
三大模块的格式相似，均为“原子序号+参数”，无需新手手动编写（力场自动生成），但必须会看懂、会核对，避免参数 mis match 报错。
1. [bonds]模块（键长控制核心）
核心作用：规定两个相连原子的键长范围，防止键长过长或过短，导致分子结构异常。
关键参数（2个核心）：
**原子对**：必须是[atoms]中存在的、实际相连的原子，序号顺序可互换； **键长参数**：力场默认值（如C-C键约0.154nm），新手无需修改，修改会导致键断裂或原子重叠； **键常数参数**：该参数控制键的振动，可以等效理解为其参数不同是区分单键、双键。
实战注意：原子对必须真实存在，不能出现“未相连的原子成对”。
2. [angles]模块（键角控制核心）
核心作用：规定三个相连原子（A-B-C）的夹角范围，维持分子的平面或空间角度。
关键参数（2个核心）：
**原子序号**：三个原子必须连续相连（中间原子为B，是夹角的顶点），序号顺序不能乱（如C1-C2-O3，顶点为C2）； **键角参数**：力场默认值（如C-C-O键角约120°），新手禁止随意修改，否则会导致分子骨架弯折异常。
实战注意：原子序号的顺序是关键，颠倒顺序会导致键角识别错误，模拟报错。
3. [dihedrals]模块（扭转角控制核心）
核心作用：规定四个相连原子（A-B-C-D）的扭转角度，控制分子的折叠、旋转，是影响分子构象的关键。
关键参数（2个核心）：
**原子序号**：四个原子必须连续相连（中间两个原子B、C是扭转轴），序号顺序固定； **扭转角参数**：力场默认值，不同构象对应不同扭转角，新手无需修改，仅需核对原子序号正确。
实战注意：二面角模块缺失不会直接导致模拟崩溃，但会导致分子构象混乱，模拟数据无效。
三、实战关键：3个必须核对的点（避坑核心）
三大模块的报错，90%源于“参数与结构不匹配”，核对以下3点，即可避开大部分问题：
原子序号：[bonds]/[angles]/[dihedrals]中的所有原子序号，必须在[atoms]模块中存在，无多余、无缺失； 连接关系：原子序号的连接关系，必须与分子结构完全一致； 参数一致性：同一类键/角（如所有C-C键）的参数一致，不随意修改。
重点提醒：若手动修改过键长、键角参数，模拟时会出现“键断裂”“原子飞散”“温度异常飙升”等问题，新手尽量保持力场默认参数。
四、新手高频误区：3个最容易踩的坑
误区1：随意删除二面角模块
错误认知：“二面角不重要，删了也能跑”；后果：分子构象混乱，模拟结果无意义。
误区2：颠倒原子序号顺序
错误操作：[angles]中写反原子序号顺序；后果：GROMACS识别错误，模拟崩溃。
误区3：手动修改键长/键角参数
错误操作：觉得“键长不对”，手动修改参数；后果：键的作用力异常，分子结构瞬间崩溃。
五、实战小结：三大模块4句口诀
键控长、角控弯、二面角控扭转； 原子序号要存在，连接关系要对得上； 力场参数不乱改，默认参数最稳妥； 三者齐全无缺失，分子骨架才稳定。
