建川・GROMACS MDP 实战 第 2 篇｜通用 MDP 完整模板
大家好，我是刘建川（网名刘十三），14 年深耕 GROMACS 分子模拟领域，始终专注分子模拟实战教学。
上一篇我们理清了 MDP 文件的核心作用，这篇直接给大家可直接copy、稳定通用的 GROMACS 全套 MDP 模板，常规蛋白、小分子水溶液体系都能直接用。
【14 年实战经验】这套模板是我多年科研与教学总结的稳定版，收敛性好、报错率低，覆盖模拟全流程，入门到科研都能用。
一、能量最小化（emin.mdp）
integrator = steep ; 积分器：最陡下降法，能量最小化专用，不可替换nsteps = 50000 ; 最大迭代步数：50000步覆盖常规体系能量收敛需求emtol = 1000.0 ; 能量收敛阈值：最大力<1000 kJ/(mol·nm)时停止emstep = 0.01 ; 能量优化步长：0.01是通用安全值cutoff-scheme = Verlet ; 截断方案：GROMACS推荐的高效截断方式rcoulomb = 1.0 ; 库仑相互作用截断半径：1.0 nm（水溶液通用）rvdw = 1.0 ; 范德华相互作用截断半径：与库仑半径保持一致
✅ 关键说明：
核心作用：消除分子结构中原子间不合理的近距离接触，避免后续模拟崩溃；
修改建议：所有参数默认即可，仅当体系超大（如复杂复合物）时，可将 nsteps 调至 100000；
验证标准：运行完成后日志显示 “Energy minimization converged” 即达标。
integrator = md ; 积分器：分子动力学正则积分，平衡/生产模拟专用dt = 0.002 ; 时间步长：0.002 ps（2 fs），改大易导致模拟不稳定nsteps = 50000 ; 总步数：50000步×0.002 ps/步=100 ps，足够完成温度平衡tcoupl = V-rescale ; 温控算法：新手最稳的温控器，不易发散tc-grps = Protein SOL ; 控温分组：蛋白和溶剂分开控温，贴合实际体系tau_t = 0.1 0.1 ; 温控弛豫时间：0.1 ps快速稳定温度，两组需对应ref_t = 300 300 ; 参考温度：300 K（常温），可按需改为298/310 Kpbc = xyz ; 周期性边界条件：三维启用，消除边界效应cutoff-scheme = Verlet ; 截断方案：与能量最小化保持一致，保证计算逻辑统一rcoulomb = 1.0 ; 库仑相互作用截断半径：1.0 nm（水溶液通用）rvdw = 1.0 ; 范德华相互作用截断半径：与库仑半径保持一致nstxout-compressed = 500 ; 压缩轨迹输出频率：每500步输出一次，兼顾数据量和精度nstenergy = 100 ; 能量输出频率：每100步输出一次，便于监控温度变化
✅ 关键说明：
核心作用：固定体系体积，将温度稳定到目标值，为后续控压平衡打基础；
修改建议：仅 ref_t 可按需调整，其余参数严禁修改；
验证标准：查看能量文件（edr），温度曲线在目标值 ±5 K 内波动即达标。
三、NPT 平衡｜npt.mdp
integrator = md ; 积分器：分子动力学正则积分，平衡/生产模拟专用dt = 0.002 ; 时间步长：0.002 ps（2 fs），改大易导致模拟不稳定nsteps = 100000 ; 总步数：100000步×0.002 ps/步=200 ps，足够完成密度平衡tcoupl = V-rescale ; 温控算法：新手最稳的温控器，不易发散tc-grps = Protein SOL ; 控温分组：蛋白和溶剂分开控温，贴合实际体系tau_t = 0.1 0.1 ; 温控弛豫时间：0.1 ps快速稳定温度，两组需对应ref_t = 300 300 ; 参考温度：300 K（常温），可按需改为298/310 Kpcoupl = Parrinello-Rahman ; 控压算法：适用于各向同性体系，控压精度高pcoupltype = isotropic ; 控压类型：各向同性控压，水溶液/蛋白体系专用tau_p = 2.0 ; 控压弛豫时间：2.0 ps通用值，过短易导致体积波动ref_p = 1.0 ; 参考压力：1.0 bar（标准大气压），可改至1.01325 barcompressibility = 4.5e-5 ; 压缩系数：水的标准值，不可修改pbc = xyz ; 周期性边界条件：三维启用，消除边界效应cutoff-scheme = Verlet ; 截断方案：与能量最小化保持一致，保证计算逻辑统一rcoulomb = 1.0 ; 库仑相互作用截断半径：1.0 nm（水溶液通用）rvdw = 1.0 ; 范德华相互作用截断半径：与库仑半径保持一致nstxout-compressed = 500 ; 压缩轨迹输出频率：每500步输出一次，兼顾数据量和精度nstenergy = 100 ; 能量输出频率：每100步输出一次，便于监控温度变化
✅ 关键说明：
核心作用：在稳定温度的基础上，调整体系体积使压力达标，让体系密度接近真实值； 修改建议：仅 ref_p 可微调，其余控压参数是水溶液体系最优配置； 验证标准：密度曲线稳定在 1000 kg/m³ 左右，压力在 1.0 bar±0.5 bar 波动即达标。
四、生产模拟｜md.mdp
integrator = md ; 积分器：分子动力学正则积分，平衡/生产模拟专用dt = 0.002 ; 时间步长：0.002 ps（2 fs），改大易导致模拟不稳定nsteps = 5000000 ; 总步数：5000000步×0.002 ps/步=10 ns，可按需调整tcoupl = V-rescale ; 温控算法：新手最稳的温控器，不易发散tc-grps = Protein SOL ; 控温分组：蛋白和溶剂分开控温，贴合实际体系tau_t = 0.1 0.1 ; 温控弛豫时间：0.1 ps快速稳定温度，两组需对应ref_t = 300 300 ; 参考温度：300 K（常温），可按需改为298/310 Kpcoupl = Parrinello-Rahman ; 控压算法：适用于各向同性体系，控压精度高pcoupltype = isotropic ; 控压类型：各向同性控压，水溶液/蛋白体系专用tau_p = 2.0 ; 控压弛豫时间：2.0 ps通用值，过短易导致体积波动ref_p = 1.0 ; 参考压力：1.0 bar（标准大气压），可改至1.01325 barcompressibility = 4.5e-5 ; 压缩系数：水的标准值，不可修改pbc = xyz ; 周期性边界条件：三维启用，消除边界效应cutoff-scheme = Verlet ; 截断方案：与能量最小化保持一致，保证计算逻辑统一rcoulomb = 1.0 ; 库仑相互作用截断半径：1.0 nm（水溶液通用）rvdw = 1.0 ; 范德华相互作用截断半径：与库仑半径保持一致nstxout-compressed = 1000 ; 压缩轨迹输出频率：每1000步输出一次，减少文件体积nstenergy = 500 ; 能量输出频率：每500步输出一次，兼顾监控和文件大小
✅ 关键说明：
核心作用：正式采集模拟数据，用于后续的结构分析、动力学分析等科研场景；
修改建议：唯一必改项是 nsteps（控制模拟时长），输出频率可根据分析需求调整（如需要更高时间分辨率可改小）；
文件注意：生产模拟会生成大文件，建议提前规划存储空间。
使用小结
运行顺序：能量最小化 → NVT → NPT → 生产模拟，一步都不能跳过； 核心原则：新手 “只改必改项，不动默认项”，避免因参数混乱导致模拟失败； 体系适配：膜蛋白 / 核酸体系需调整控压类型 / 截断半径，后续文章会专项讲解。
📌 下一篇预告：刘建川・GROMACS MDP 实战 第 3 篇｜积分器与时间步参数
这套通用 MDP 模板能帮你快速上手基础模拟，但想要根据自身科研体系灵活调参、规避报错、做出高质量模拟数据，还需要系统掌握背后逻辑。
近期我将开设 GROMACS 线下小班，14 年实战经验，手把手带你打通 GROMACS 全流程，从零基础入门到科研成果落地，让模拟真正成为你的科研核心利器。
