计算化学公社

标题: 求助：如何在x，y方向上限制几个原子的质心 [打印本页]

作者
Author: 1900 时间: 2021-8-27 16:32
标题: 求助：如何在x，y方向上限制几个原子的质心
我想要在z方向对分子一端的几个原子进行拉伸模拟，但是不希望这些原子在xy平面上平移太多。所以想只在xy平面上限制一下这几个原子的质心。

作者
Author: 园园乖 时间: 2021-8-27 22:03
可以在mdp文件中加入define=-DPOSRES，然后在itp文件中[ position_restraints ]字段中加入#ifdef POSRES，再在相应的原子后面分别设置限制力，（只xy平面限制的话，1000 1000 0，应该就可以）

作者
Author: 1900 时间: 2021-9-1 09:57

园园乖发表于 2021-8-27 22:03
可以在mdp文件中加入define=-DPOSRES，然后在itp文件中[ position_restraints ]字段中加入#ifdef POSRES， ...

多谢回复。但是这样只能限制原子位置，无法做到只限制质心。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)