24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1750  |  回复: 13

亮剑随风

金虫 (小有名气)

[求助] 高斯opt+fre运算老出错,求高人指点已有1人参与

现把我的gjf和log文件附上,恳请高斯方面的高人帮忙给看看。。。到底是构型方面的问题还是结构太大的原因?新人刚接触高斯计算,本科毕业设计到现在还没啥进展。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : BL-03.log
  • 2014-04-25 14:58:47, 126.17 K
  • 附件 2 : BL-03.gjf
  • 2014-04-25 15:02:27, 3.85 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

庸者常立志,智者立长志!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

Convergence failure的信息。
看样子scf没收敛,我还没用gview看一下,造成这些问题的有可能是结构不合理或是其他原因暂时还不好判断。

另外建议做计算过程在计算参数这一行采用#p来显示详细的迭代信息,这样便于我们进行诊断,另外建议保存chk文件,它储存了很多波函数的信息,一来便于分析,二来就是read check是一种常用的scf收敛技巧。

另外站内有g09的中文手册最好看一下,借助毕业论文的机会了解一下量子化学这个领域开开眼界也不错,否则这个毕业论文就是浪费时间了。

结构似乎比较大,我可以尝试帮你做个预优化,达到6-31+g*的水平,其他的你自己来。我认为这么大的体系做频率分析4cores估计有些困难。

本科毕业论文要求不能太高。。。。
不一定挂在论坛,计算问题问题欢迎留言。
3楼2014-04-25 15:49:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

FoggyEM

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
最后一个原子坐标
H                  0.93815032    2.73869541    1.45776498
两个回车,然后把后面的数据都删掉。应该就可以了。

也就是这样
……(原子坐标)
H                  0.93815032    2.73869541    1.45776498
回车
回车
(剩下的都删了)
Everysinglestepoftheway,Ipaid.
2楼2014-04-25 15:39:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亮剑随风

金虫 (小有名气)

引用回帖:
2楼: Originally posted by FoggyEM at 2014-04-25 15:39:08
最后一个原子坐标
H                  0.93815032    2.73869541    1.45776498
两个回车,然后把后面的数据都删掉。应该就可以了。

也就是这样
……(原子坐标)
H                  0.93815032    2.73869 ...

好的!我试一下
庸者常立志,智者立长志!
4楼2014-04-25 16:14:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
4楼: Originally posted by 亮剑随风 at 2014-04-25 16:14:59
好的!我试一下...

和输入文件没关系,不用去修改格式,log文件的错误已经看到了。另外,有一点希望你能够明白,你的体系含有氢键,用 B3LYP结果会比较糟糕。这些内容是你应该查阅文献或资料得到的,目前帮你做预优化用的是wB97xd,文件稍微晚些传上来。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2014-04-25 16:28:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oyezzy

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by FoggyEM at 2014-04-25 15:39:08
最后一个原子坐标
H                  0.93815032    2.73869541    1.45776498
两个回车,然后把后面的数据都删掉。应该就可以了。

也就是这样
……(原子坐标)
H                  0.93815032    2.73869 ...

如果如你所讲的话,#行的 geom=connectivity 也得删除,否则直接回报错。
初来乍到多多关照
6楼2014-04-25 16:56:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亮剑随风

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 亮剑随风 at 2014-04-25 16:14:59
好的!我试一下...

果断不行!

[ 发自小木虫客户端 ]
庸者常立志,智者立长志!
7楼2014-04-26 07:14:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亮剑随风

金虫 (小有名气)

引用回帖:
5楼: Originally posted by 卡开发发 at 2014-04-25 16:28:22
和输入文件没关系,不用去修改格式,log文件的错误已经看到了。另外,有一点希望你能够明白,你的体系含有氢键,用 B3LYP结果会比较糟糕。这些内容是你应该查阅文献或资料得到的,目前帮你做预优化用的是wB97xd,文 ...

我的毕业设计主要是寻找两种分子的相互作用!即寻找分子间氢键。另外请教下(b3lyp不是最小机组么?)log文件中报的是l502错误,我现在将scf=tight改为scf=qc正在尝试。

[ 发自小木虫客户端 ]
庸者常立志,智者立长志!
8楼2014-04-26 07:23:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xrmc

银虫 (初入文坛)

引用回帖:
8楼: Originally posted by 亮剑随风 at 2014-04-26 07:23:04
我的毕业设计主要是寻找两种分子的相互作用!即寻找分子间氢键。另外请教下(b3lyp不是最小机组么?)log文件中报的是l502错误,我现在将scf=tight改为scf=qc正在尝试。
...

b3lyp是方法吧,不是基组
9楼2014-04-26 08:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亮剑随风

金虫 (小有名气)

引用回帖:
9楼: Originally posted by xrmc at 2014-04-26 08:04:44
b3lyp是方法吧,不是基组...

抱歉,是6-31G(p)基组。

[ 发自小木虫客户端 ]
庸者常立志,智者立长志!
10楼2014-04-26 08:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 亮剑随风 的主题更新
信息提示
请填处理意见