【摘要】 系统讲清 GROMACS 中冻结与位置限制的区别、参数写法、方向控制、索引组设置与 posre 文件引用方法。

科研人必备干货｜全套科研知识库

常用科研资料汇总，点击领取即可获取。

 

在 GROMACS 中，冻结（Freeze）和位置限制（Position Restraints）都能降低部分原子的运动自由度，但两者的物理含义和实现方式并不相同。冻结是直接阻止原子在指定方向上的坐标更新，位置限制则是在原子偏离参考位置时施加恢复力，让原子在参考位置附近振动。

 

冻结和位置限制的核心区别

如果目标是让某一组原子在模拟过程中保持绝对不动，通常使用冻结；如果目标是保留局部热振动，但限制其偏离参考结构的幅度，则更适合使用位置限制。

方法

作用机制

典型用途

Freeze

禁止指定方向上的坐标更新

固定基底、限制墙面、只允许单方向运动

Position Restraints

偏离参考位置时施加恢复力

平衡阶段稳定骨架、保持结构参考位形

 

冻结参数怎么写

GROMACS 中冻结主要通过 .mdp 文件里的两个参数控制：

  • freezegrps：指定需要冻结的原子组

  • freezedim：指定该组在 X/Y/Z 三个方向上的冻结状态

 

完全冻结一个原子组

freezegrps = MOL freezedim = Y Y Y

这表示 MOL 组在 X、Y、Z 三个方向全部冻结，模拟过程中保持固定。

 

只冻结部分方向

freezegrps = MOL freezedim = Y Y N

这表示 X、Y 方向冻结，Z 方向仍可运动，适合只允许分子沿某一方向滑移或扩散的设置。

 

多组冻结时要注意什么

如果要同时冻结多个原子组，每个原子组都必须对应三个方向参数。例如：

freezegrps = Group_A Group_B freezedim = Y Y Y Y Y N

这里 Group_A 对应完全冻结，Group_B 只允许沿 Z 方向运动。组数量和方向参数数量必须严格匹配，否则容易在预处理阶段出现设置错误。

 

位置限制如何生成

位置限制通过拓扑中的 [ position_restraints ] 定义，通常最方便的做法是直接用 gmx genrestr 自动生成：

gmx genrestr -f conf.gro -o posre.itp -fc 1000 1000 1000

生成后的 posre.itp 会为所选原子组写入三方向的力常数。常见字段包括：

  • atom：原子编号

  • funct：位置限制势函数类型，常规为 1

  • fcx/fcy/fcz：X、Y、Z 三方向的限制力常数

 

不同方向的位置限制怎么控制

位置限制不要求三个方向使用相同力常数。例如：

gmx genrestr -f conf.gro -o posre.itp -fc 1000 1000 0

这表示只在 X、Y 方向施加限制，而 Z 方向自由。对于界面模型、膜体系和表面吸附体系，这种各向异性限制非常常见。

 

posre 文件怎么启用

生成 posre.itp 后，还需要在对应分子拓扑中通过宏引用：

#ifdef POSRES #include "posre.itp" #endif

然后在 .mdp 中加入：

define = -DPOSRES

如果不启用这个宏，对应的位置限制不会被加载。

 

grompp 为什么还要加 -r

使用位置限制时，gmx grompp 需要明确参考坐标文件：

gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr

其中 -c 是模拟起始结构，-r 是位置限制参考坐标。若缺少 -r 或参考文件不匹配，通常会导致位置限制设置无法按预期生效。

 

FAQ

冻结和位置限制可以一起用吗？

可以，但需要明确它们分别作用在哪些原子组，以及是否会造成不必要的重复约束。对于完全固定的基底，通常只用冻结即可；对于想保持构型但保留热振动的结构，更适合用位置限制。

冻结是不是比位置限制更稳定？

不一定。冻结更“强”，但并不总是更合理。若体系本应保留局部热运动，完全冻结反而可能引入不自然的边界效应。选择哪一种方式，关键取决于模型物理目标。

 

结语

理解 Freeze 与 Position Restraints 的差异，远比记住单个命令更重要。只要把目标原子组、方向控制、参考结构和拓扑引用关系理清，GROMACS 中的大多数约束设置问题都能更快定位和修正。