第二章 搭台体系构建

第二章 · 搭台：体系构建

章节摘要：体系构建是把"我想研究的分子"变成"程序能推演的原子集合"的全过程：选力场、定拓扑、装盒子、填溶剂、加离子，五步环环相扣，任何一步的草率都会在第3章之后以崩溃或错误结论的形式讨债。本章以聚环氧乙烷熔体与 DPPC 脂膜为主案例、以经典蛋白流程为对照，把每一步的命令、输出与检查方法过一遍。读完你应当能独立搭出一个电荷中性、参数自洽、浓度明确的可运行体系。

学习目标

1. 面对一个新分子，能按"化学家族 → 现成力场 → 参数来源 → 水模型配套"的顺序做出选型，并说出理由与风险。
2. 用脚本核对手头力场的 LJ 参数、混合规则与 1-4 处理，确认拓扑文件引用的参数与力场目录里的表一致。
3. 跑通 pdb2gmx 的完整流程，读懂它的交互问答（端基、质子化、二硫键），并知道 pdb2gmx 认不了的分子（聚合物、脂质、药物小分子）该走哪条替代路线。
4. 用 editconf、solvate、genion 三条命令完成盒子、溶剂化与离子中和，并能在运行前后对上原子数的账。
5. 对一个 6 nm 立方盒子，口算出需要多少水分子、0.15 M 盐浓度对应多少对离子——这两个数字是体系构建的"九九表"。

核心概念速览

体系构建的本质是三本账：原子的账（拓扑里每个原子都有类型、电荷、键接）、电荷的账（净电荷必须中和或明确不中和）、浓度的账（溶剂与溶质的比例要能换算成实验可比的数）。这三本账对上，体系才算立住了。

子章节导航

• 2.1 力场选型与参数核对——选型决策表加上一段真正能跑的核对脚本，把"参数从哪来"变成"参数我验过"。PEO 的联合原子方案与 DPPC 的 CHARMM 方案都在本节定案。
• 2.2 pdb2gmx 拓扑生成实操——以多肽走一遍标准流程，把交互问答的每个选项讲透；再看聚合物走 CGenFF 与手写 itp 的替代路线，原子数怎么对账。
• 2.3 盒子溶剂化与加离子——盒子留多边、水填多少、离子加几对，全部用演算回答；顺带讲清周期性边界与最小镜像这两个贯穿全书的几何约定。

子章节之间的逻辑关系

三节是一条装配线：选型决定拓扑里写什么参数（2.1），拓扑决定装配时原子从哪来（2.2），装配完的体系决定后续一切模拟的物理含义（2.3）。倒着做会返工——先填水再想力场，换力场时连水模型都得跟着换，等于重搭。

[2.1 选型与核对] --参数表--> [2.2 拓扑落地] --原子清单--> [2.3 装配与对账]
       |                            |                          |
   决定"物理"                  决定"结构"                 决定"浓度与电荷"
       +------------ 装配完成的 gro + top 交棒 ------------> [第三章 grompp]
​

本章演算与脚本清单

关键数字速查

• 33.4：每立方纳米的水分子数，装配演算的通用常数。
• -d 1.0 底线：盒壁距离必须大于截断半径，柔性链建议 1.2。
• 0.15 M + 中和：离子数两件套，浓度按体积、中和按电荷，别混。
• r_c < 半盒宽：周期性自相互作用的红线。

带表进入实操，本章的三个节就是三张工序卡：选型（2.1）、办证（2.2）、装配（2.3），每张卡验收动作都在节末。

前置知识与后续延伸

前置：1.2 的六大势能项与混合规则——本章反复用到"σ、ε 从哪查、1-4 怎么缩放"这些概念；Linux 命令行基本操作。

延伸：2.3 装好的体系在 3.1 进入 mdp 参数逐项精读；本章建立的"对账"习惯在 4.3 排错时是第一反应——绝大多数崩溃的第一嫌疑人都是构建阶段留下的暗伤。脂膜装配的特殊性（双层定向、面积压强）在 7.2 案例里完整展开。