刘建川・GROMACS MDP 实战 第 6 篇｜盒子、边界与平衡模拟
大家好，我是刘建川（网名刘十三），14 年深耕 GROMACS 分子模拟领域，始终专注分子模拟实战教学。
上一篇我们掌握了压力控制参数，搞定了 NPT 平衡的核心设置；这一篇聚焦模拟的 “空间基础”—— 盒子、边界条件，再搭配平衡模拟的关键要点，内容适中、好记好用，衔接前几篇的知识点，帮你避开平衡阶段的常见坑。
【14 年实战经验】盒子和边界条件看似简单，却是模拟不出现 “边界效应” 的关键；平衡模拟的顺序和细节没做好，后续生产模拟必崩溃，新手重点记 “固定配置 + 正确顺序” 即可。
一、盒子与边界条件（核心配置）
盒子是分子模拟的 “容器”，边界条件用来消除容器边缘的不合理作用，两者配套设置，常规体系直接套用以下配置：
pbc = xyz ; 三维周期性边界条件，必开，消除边界效应模拟盒子类型，在建立模拟盒子 gro 文件时候，可指定，常用 cubic 立方体盒子；gro 文件最后一行三个数，是模拟盒子的 x、y、z 尺寸（单位 nm）。
✅ 关键说明：
pbc 参数：固定设为 xyz（三维启用），关闭（no）会导致分子在盒子边缘碰撞、体系崩溃，常规模拟绝不关闭；
盒子尺寸：建模时（如 gmx solvate）会自动根据体系大小设置，或者自己手动指定参数（gmx editconf -box）。对于蛋白溶液模拟来说确保蛋白分子与盒子边缘有足够距离（≥1.0 nm），避免周期镜像影响。
二、平衡模拟关键要点（衔接前几篇）
平衡模拟是模拟稳定的 “必经步骤”，核心是 “先稳温度、再稳密度”，顺序不能乱，参数沿用前几篇的通用配置，无需额外修改：
2.1 平衡模拟的正确顺序（必记）
能量最小化（emin）：消除原子不合理接触（第一篇、第二篇模板直接用）； NVT 平衡（nvt）：只控温、不控压，稳定体系温度（第四篇温控参数）； NPT 平衡（npt）：控温 + 控压，调整体系密度（第五篇温压组合模板）。
2.2 平衡模拟核心配置（直接抄）
; 沿用前几篇通用参数，无需修改integrator = md ; 动力学积分器dt = 0.002 ; 时间步2fs; 温控（第四篇）tcoupl = V-rescaletc-grps = Protein SOLtau_t = 0.1 0.1ref_t = 300 300; 控压（第五篇）pcoupl = Parrinello-Rahmanpcoupltype = isotropictau_p = 5.0ref_p = 1.0compressibility = 4.5e-5; 盒子与边界（本篇）pbc = xyz
✅ 关键说明：
顺序不能乱：跳过 NVT 直接跑 NPT，会导致体系崩溃；
参数不修改：新手遇到不明确的参数，不要放在 mdp 文件里，gromacs 大多有默认参数的；
验证标准：NVT 平衡后温度稳定在 300±5K，平衡后随模拟时间的温度曲线波动是非常小的；但是 NPT 平衡后密度平均值一般在设定值周围，随模拟时间的压力曲线波动是非常大的，这是正常现象。
三、新手常见避坑点
坑 1：关闭 pbc 或设为 no → 体系边缘分子碰撞，模拟崩溃，解决：固定设 pbc=xyz； 坑 2：平衡顺序颠倒（先 NPT 后 NVT） → 温压失衡，解决：严格按 “em→NVT→NPT” 顺序；
使用小结
盒子与边界：固定用 pbc=xyz； 平衡模拟：按 “能量最小化→NVT→NPT” 顺序，参数沿用前几篇通用配置； 核心原则：不跳过平衡、不乱改参数，就能避免大部分崩溃问题。
📌 下一篇预告：刘建川・GROMACS MDP 实战 第 7 篇｜输出设置与常见错误修正
盒子、边界和平衡是模拟的 “基础框架”，掌握这些能保证模拟稳定运行，但想要高效处理模拟输出文件、快速排查报错，还需要系统的实战指导。
近期我将开设 GROMACS 线下小班，14 年实战经验，手把手带你打通 GROMACS 全流程，从零基础入门到科研成果落地，让模拟真正成为你的科研核心利器。
