计算化学公社

 找回密码 Forget password
 注册 Register

关于在gromacs中利用UFF力场模拟的问题

查看数: 25811 | 评论数: 9 | 收藏 Add to favorites 13
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2015-1-27 17:03

正文摘要:

最近在做一个包含Ir原子的有机配合物,分子式如图一。想利用gromacs模拟一下,但是一直没有做出来,还请各位帮我看看!先谢过了! 主要问题是Ir只能使用UFF力场,我利用OBGMX(http://software-lisc.fbk.eu/obgmx ...

回复 Reply

sun666 发表于 Post on 2021-6-2 09:40:01
请教楼主,金属配合物的问题解决了吗。
tjuptz 发表于 Post on 2020-12-23 22:29:47
youxx189 发表于 2020-8-28 13:19
您好请问一下,Ir原子在UFF力场下能量最小化后结构发生了变形,这是什么导致的啊,谢谢老师!

力场参数不合适呗
youxx189 发表于 Post on 2020-8-28 13:19:41
您好请问一下,Ir原子在UFF力场下能量最小化后结构发生了变形,这是什么导致的啊,谢谢老师!
yaochuang 发表于 Post on 2015-1-29 08:31:45
sobereva 发表于 2015-1-28 23:32
grompp加上-maxwarn xxx就可以设定warning达到xxx时才不允许执行

啊~  原来是这样,谢谢sob老师! 哈哈~
sobereva 发表于 Post on 2015-1-28 23:32:20
yaochuang 发表于 2015-1-28 15:45
谢谢,sob老师,我把原子名改了就可以了!  现在这个版本的gmx中超过20个不匹配就不能继续进行了。 以前 ...

grompp加上-maxwarn xxx就可以设定warning达到xxx时才不允许执行
yaochuang 发表于 Post on 2015-1-28 15:45:22
sobereva 发表于 2015-1-28 12:15
这只是警告,不是报错,可以无视。

.itp中有这样的部分

谢谢,sob老师,我把原子名改了就可以了!  现在这个版本的gmx中超过20个不匹配就不能继续进行了。 以前我总是在pdb中改原子名,看来是不行的。
sobereva 发表于 Post on 2015-1-28 12:15:27
这只是警告,不是报错,可以无视。

.itp中有这样的部分
[ atoms ]
; nr type  resnr    residue    atom     cgnr    charge       mass
1        C_R        1        UFF        C_R        1        0.025723        12.0107                       
2        C_R        1        UFF        C_R        2        -0.097501        12.0107                       
3        C_R        1        UFF        C_R        3        -0.119779        12.0107       
...
第2列是原子类型,第5列是原子名。实际上,只要你的pdb文件里的原子顺序和这部分的原子顺序一致就行了,就可以正常计算,名字不需要非得一致。但是grompp为了谨慎起见,发现第5列的原子名和结构文件里的原子名对不上,就会警告一下。要想避免警告,就把上面的第5列和pdb文件里的原子名弄一致。
sobereva 发表于 Post on 2015-1-28 01:39:47
本身gromacs就不支持UFF力场,而且UFF力场形式本身不能直接被gromacs所支持。OBGMX的作者把UFF力场形式做了一些近似,才勉强使之能用在gromacs里面。

最好把你详细的操作步骤和相关的文件都贴出来,否则说不清楚。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 16:53 , Processed in 0.166998 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list