求助在gromacs模拟的过程中要限制某几组原子之间的距离

琦锅锅

wuzhiyi 发表于 2020-9-19 17:49
bond 是distance_restraints的特殊情况
bA 距离 kA 强度 后面两个是fep的

嗯嗯，好的，谢谢

琦锅锅

wuzhiyi 发表于 2020-9-19 17:49
bond 是distance_restraints的特殊情况
bA 距离 kA 强度 后面两个是fep的

好的，非常感谢

Jeffrey0521

新手求问，到底是disre还是dihre啊，我看手册都写的disre但是sob老师写的dihre。http://sobereva.com/10

喵星大佬

Jeffrey0521 发表于 2024-1-18 15:02
新手求问，到底是disre还是dihre啊，我看手册都写的disre但是sob老师写的dihre。http://sobereva.com/10

普通情况不要用[ distance_restrains ]去设置距离限制，直接用[ bonds ]里面函数类型6或者10去做谐振/分段的距离限制势。[ distance_restrains ]里面的距离限制是用来精修NMR解的结构用的，距离限制用的是时间平均的形式满足NOE信号的距离-6次方的期望，力常数也要在mdp里专门设置，不是正常理解的那种距离限制的形式

至于本文楼主说的那个目的，应该用pull去实现

Jeffrey0521

喵星大佬 发表于 2024-1-19 03:45
普通情况不要用[ distance_restrains ]去设置距离限制，直接用[ bonds ]里面函数类型6或者10去做谐振/分 ...

我想限制形成氢键的两个原子之间的距离也是这样设置吗？我在mdp里写入了力常数的，所以才问是写入的disre还是dihre

喵星大佬

Jeffrey0521 发表于 2024-1-19 15:45
我想限制形成氢键的两个原子之间的距离也是这样设置吗？我在mdp里写入了力常数的，所以才问是写入的disre ...

那个功能就不是用来做这个的，直接在[ bonds ]里面加就好

Newbee_Ccc

琦锅锅 发表于 2020-9-29 16:32
好的，非常感谢

楼主，请问现在是否解决了？我也按照这个方式，加上了限制但是MD的时候发现还是跑飞了

Newbee_Ccc

wuzhiyi 发表于 2020-9-19 17:49
bond 是distance_restraints的特殊情况
bA 距离 kA 强度 后面两个是fep的

您好！我参考了您那建议在top文件中加了这个命令，

[ intermolecular_interactions]
[ bonds ]
; ai     aj    type   bA      kA     bB      kB
4562  4576   6      2.536   1000.0
4514  4575   6      4.210   1000.0
随后又在mdp文件中加入了如下命令，
disre                    = Simple
; Force weighting of pairs in one distance restraint: Conservative or Equal

disre-weighting          = Conservative
; Use sqrt of the time averaged times the instantaneous violation
disre-mixed               = no
disre-fc                     = 500
disre-tau                   = 0
; Output frequency for pair distances to energy file
nstdisreout               = 100
但是当我做了短时间模拟后发现，并未限制住两个原子间的距离，能否麻烦您帮忙看看，可能是哪里出了问题呢？