24小时热门版块排行榜    

查看: 716  |  回复: 5

fight1man

金虫 (小有名气)

[求助] Gaussion 加水进行优化计算出错,求高手进来看看啊,感激不尽,焦急阿

一样的输入文件算法,但是对单体A和B分别计算错误不成功,奇怪的是,而对通过AUTODOCK对接得到的A和B的结合物C能计算成功,非常纳闷之中,请高手指点谢谢
附上输入文件
%chk=D:\xy\lv\danti\DMCDPM3h2o.chk
%mem=500MB
%rwf=d:\1,500mb,d:\2,500mb,d:\3,500mb,d:\4,500mb,d:\5,500mb
# hf scrf=(ipcm,solvent=water) geom=connectivity
和出错的输出文件[ Last edited by fight1man on 2012-10-25 at 14:51 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : BCDh2o.out
  • 2012-10-25 13:12:33, 60.26 K
  • 附件 2 : BCDh2o.gjf
  • 2012-10-25 13:21:11, 11.51 K
  • 附件 3 : DM1010g2.out
  • 2012-10-25 14:51:02, 64.05 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxyuan09

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2012-10-26 12:20:38
fight1man: 金币+25, 有帮助 2012-10-28 19:32:57
1、出错原因
(1)、BCDho.out
这个是自洽场(SCF)迭代不收敛,很多时候不太容易解决,可参照本版置顶帖Gaussian FAQ的2.3节的相关内容来解决。
http://muchong.com/bbs/viewthread.php?tid=1566146
(2)、DM1010g2.out
这个是IPCM溶剂模型计算中迭代不收敛。
可以在关键词那一行(#打头)添加IOP(1/6=N)来解决,默认N=200,你的输出文件就是200循环还没有收敛出错。你可以设置的大一点。

2、建议。
你使用的计算方法(hf)太粗错,使得计算容易出各种问题。
你的体系应该是一个客体分子进入到环糊精的腔体中,属于弱相互作用,hf方法处理的实在是太差了,可以使用专用的处理弱相互作用的方法,比如B97D。
另外基组没有指定的话是使用sto-3g,也太粗糙,应该弱相互作用应该有极化函数(比如6-31G*),最好也有弥散函数(比如6-31+G*)。
即:建议使用B97D/6-31G*或者B97D/6-31+G*。不过你的体系有点大,带了加号算不动,先用前一种方法试试看吧。
2楼2012-10-25 20:16:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-10-26 12:20:43
fight1man: 金币+25, 有帮助, 谢谢啊,看来我得换理论方法基组看看了。可是我的体系要加水溶剂模型,您觉得我用哪个理论方法算比较合适啊。谢谢 2012-10-28 19:32:42
你的问题出在IPCM溶剂模型计算中迭代不收敛。而且从 Iteration number XXX MaxDiff 后面的数值可以看出,越迭代距离收敛限越远,所以增加迭代次数是徒劳。
如果没有特别要求必须使用IPCM模型的话,建议将IPCM关键词改为PCM。

你计算时使用的理论方法和基组明显不合适,这在LS已经提到了。但这并非是导致出错的原因。
3楼2012-10-25 23:41:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fight1man

金虫 (小有名气)

引用回帖:
2楼: Originally posted by cxyuan09 at 2012-10-25 20:16:28
1、出错原因
(1)、BCDho.out
这个是自洽场(SCF)迭代不收敛,很多时候不太容易解决,可参照本版置顶帖Gaussian FAQ的2.3节的相关内容来解决。
http://muchong.com/bbs/viewthread.php?tid=1566146
(2)、DM1010g2.o ...

嗯,谢谢您的回复,因为我想用水溶剂模型,它的算法不是有限制吗,所以我先选了HF 本来想用DFT的 但是太耗时 时间不是很宽裕
我想问下您说的B97D可以加溶剂水吗,我这个体系有这个算法 一般大概需要多少时间啊 谢谢
4楼2012-10-28 19:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fight1man

金虫 (小有名气)

引用回帖:
2楼: Originally posted by cxyuan09 at 2012-10-25 20:16:28
1、出错原因
(1)、BCDho.out
这个是自洽场(SCF)迭代不收敛,很多时候不太容易解决,可参照本版置顶帖Gaussian FAQ的2.3节的相关内容来解决。
http://muchong.com/bbs/viewthread.php?tid=1566146
(2)、DM1010g2.o ...

词那一行(#打头)添加IOP(1/6=N)来解决,默认N=200,尝试了不行。 理论方法换成B97D 加水溶剂模型 还是出错
5楼2012-10-28 19:57:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fight1man

金虫 (小有名气)

引用回帖:
3楼: Originally posted by sobereva at 2012-10-25 23:41:40
你的问题出在IPCM溶剂模型计算中迭代不收敛。而且从 Iteration number XXX MaxDiff 后面的数值可以看出,越迭代距离收敛限越远,所以增加迭代次数是徒劳。
如果没有特别要求必须使用IPCM模型的话,建议将IPCM关键词 ...

换成了PCM模型 还是出错哈
6楼2012-10-28 19:58:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fight1man 的主题更新
信息提示
请填处理意见