GROMACS体系能量最小化时加入的聚合物分子构型一直出错求助

查看数: 3882 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2021-12-7 15:23

本帖最后由牛文文于 2021-12-7 15:27 编辑各位老师同学大家好，我目前在尝试用质子化后的OPBI十二聚体掺杂磷酸进行模拟。力场用的是GAFF力场，用RESP电荷进行拟合。质子化后的OPBI十二聚体（带24个正电荷）的 ...

回复 Reply

sobereva 发表于 Post on 2021-12-10 04:29:19

牛文文发表于 2021-12-9 18:47
老师我修改了键角之后发现还是不行，看了看二面角应该也存在问题。但是二面角太复杂了，我打算试试楼里提 ...

跟优化没直接关系

牛文文 发表于 Post on 2021-12-9 18:47:47

sobereva 发表于 2021-12-9 06:25
先把能明显看出不合理的地方改了，如果还不行再说别的

老师我修改了键角之后发现还是不行，看了看二面角应该也存在问题。但是二面角太复杂了，我打算试试楼里提到的ztop了。
不过我有一个疑问，就是我这个问题的根本原因还是构型优化不合理导致的吗？
谢谢老师回复

牛文文 发表于 Post on 2021-12-9 09:01:32

心向暖阳发表于 2021-12-8 18:01
试试ztop

我先试试手动修改键角，之后可以研究研究ztop

sobereva 发表于 Post on 2021-12-9 06:25:09

牛文文发表于 2021-12-8 16:56
谢谢sob老师，仔细研究了下itp文件，确实出错的地方键角会偏小。请问我想手动修改的话，除了键角这里修改 ...

先把能明显看出不合理的地方改了，如果还不行再说别的

心向暖阳 发表于 Post on 2021-12-8 18:01:50

试试ztop

牛文文 发表于 Post on 2021-12-8 16:56:17

sobereva 发表于 2021-12-8 08:28
你先弄一个简化模型，只保留其中一个单元并且边缘用氢封端，然后跑跑，如果也有这个问题，仔细检查itp里的a ...

谢谢sob老师，仔细研究了下itp文件，确实出错的地方键角会偏小。请问我想手动修改的话，除了键角这里修改，还有其他地方需要修改吗？
gromacs初学者在这方面还没有经验，谢谢您的回答

sobereva 发表于 Post on 2021-12-8 08:28:26

你先弄一个简化模型，只保留其中一个单元并且边缘用氢封端，然后跑跑，如果也有这个问题，仔细检查itp里的angle项参数的合理性

成键方式不影响量化计算，但对于acpype能正确指认参数有至关重要的影响，给acpype提供的mol2文件里的形式键级一定要保证正确

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register