| 查看: 1570 | 回复: 0 | ||
[求助]
gromacs能量最小化之后，分子构型都变了！！
|
|
刚接触gromacs，想模拟一下小分子有机物。先拿苯练手！ 遇上一些问题，想请教一下各位！！谢啦 先讲一下我使用过程(避免是操作上的失误！！) 1、得到(PRODRG)苯的pdb文件(benzene.pdb)，建立一个空盒子（empty.gro）； 2、向空盒子中添加200个苯分子 genbox -cp empty.gro -ci benzene.pdb -nmol 200 -o benzene.gro ; 3、编写top文件和mdp文件，使用的是oplsaa力场，如附件！不知道写对没！！ ***这里遇到一个问题是苯它是平面分子，我的力场给了improper 二面角，但是还是观察苯上的原子不共面*** 4、grompp -f benzene.mdp -c benzene.gro -p system.top -o benzene.tpr 5、mdrun -v -deffnm benzene 6、 按理第2步得到的初始结构并不是合理的，需要先进行能量最小化步。但是得到的结果发现。能量最小过程后(超过迭代步数，能量并未收敛值规定值)，苯分子的构型完全变了，如图1,2 (前后比较！！) ###请各位看看，比较急！！金币不够再加！！！#### 谢谢啦！！  1.png  2.png |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : benzene.gro
- 附件 2 : 2.png
- 附件 3 : benzene.pdb
- 附件 4 : empty.gro
- 附件 5 : benzene.gro
- 附件 6 : system.top
- 附件 7 : md.mdp
- 附件 8 : 1.png
2015-08-23 16:58:33, 105.55 K
2015-08-23 17:03:04, 28.66 K
2015-08-23 17:05:56, 1.25 K
2015-08-23 17:06:07, 28 bytes
2015-08-23 17:06:18, 105.55 K
2015-08-23 17:06:39, 100 bytes
2015-08-23 17:06:50, 879 bytes
2015-08-23 17:08:07, 6.31 K
» 猜你喜欢
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
面上项目申报
已经有3人回复
-
遇见不省心的家人很难过
已经有22人回复
» 本主题相关商家推荐: (我也要在这里推广)
