计算化学公社

标题: 采用B3LYP/6-31g（d.p）优化后出现问题 [打印本页]

作者
Author: xiaomidaxue 时间: 2019-7-29 12:34
标题: 采用B3LYP/6-31g（d.p）优化后出现问题
采用B3LYP/6-31g（d.p）优化后能够正常结束，但是打开log问题时提示如下错误：见附图1，并且HOMO/LUMO的数值存在问题。
请问这是什么原因产生的，我的优化结构opt+freq能够正常结束，但是出现上述错误后，导致我的HOMO/LUMO轨道出现较大问题。
谢谢！

作者
Author: zjxitcc 时间: 2019-7-29 13:27
“采用B3LYP/6-31g（d.p）优化后能够正常结束”，这句里没有一个名词与你的log文件能对的上。。。
（1）log文件里面显示你用的是3-21G，这样的计算没有意义；
（2）结尾明明是“Error termination”，你却说“正常结束”；
（3）看分子轨道要用.fchk文件看。
先把这些改了再看吧。

作者
Author: xiaomidaxue 时间: 2019-7-29 16:07

zjxitcc 发表于 2019-7-29 13:27
“采用B3LYP/6-31g（d.p）优化后能够正常结束”，这句里没有一个名词与你的log文件能对的上。。。
（1）lo ...

不好意思，上传的文件有问题。现在从新上传正确的文件。
谢谢！

作者
Author: zjxitcc 时间: 2019-7-29 16:56
本帖最后由 zjxitcc 于 2019-7-29 16:57 编辑

优化过程是正常的，没有问题。你用GV打开碰到的问题其实算是高斯本身的bug，极化率数值太大，导致显示不全变成星号****（在32519行）。
解决方法是把显示不全的这两行删去，或者自己修改为任一个格式正常的数。
如果想知道极化率数值，其实可以在文件末尾找到的，因此此处可以随便改。
再提醒2点：（1）这种体系结构优化用B3LYP是不够的，需要加色散；（2）不要再用log看轨道了。

作者
Author: xiaomidaxue 时间: 2019-7-29 17:45

zjxitcc 发表于 2019-7-29 16:56
优化过程是正常的，没有问题。你用GV打开碰到的问题其实算是高斯本身的bug，极化率数值太大，导致显示不全 ...

好的，非常感谢。这种体系论文直接写B3LYP进行优化，我是想在论文结构上就行修改，然后合成材料。
log主要是看HOMO/LUMO的数值。一般轨道电子云分布我也是通过fchk或者cube文件进行查看。

作者
Author: zjxitcc 时间: 2019-7-29 18:39

xiaomidaxue 发表于 2019-7-29 17:45
好的，非常感谢。这种体系论文直接写B3LYP进行优化，我是想在论文结构上就行修改，然后合成材料。
log主 ...

如果你所说的“论文”是别人已发表的文章，那是前人考虑不周（为啥没考虑，可能的因素很多），而你现在应该考虑色散，结果更有说服力；不加的话，碰到挑刺的审稿人就完了，看《坚持使用B3LYP算有机体系的人的下场》http://bbs.keinsci.com/thread-12773-1-1.html。

作者
Author: sobereva 时间: 2019-7-29 22:37
顺带提一句：用Multiwfn载入fch文件，进入主功能0，直接在文本窗口就能看见HOMO、LUMO的各种单位下的能量以及序号，比起直接看log/out文件省事得多

看电子密度分布的等值面图的话，用Multiwfn载入fch，用主功能5直接就能绘制（看手册4.5节的例子），完全不需要经历产生cube文件

看轨道的话，用Multiwfn载入fch文件直接就能看，比用gview省事得多而且速度快得多得多，效果还好得多，还免费，见
使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）
用VMD绘制艺术级轨道等值面图的方法（含演示视频）
http://sobereva.com/449（http://bbs.keinsci.com/thread-11550-1-1.html）

作者
Author: xiaomidaxue 时间: 2019-7-30 16:45

zjxitcc 发表于 2019-7-29 18:39
如果你所说的“论文”是别人已发表的文章，那是前人考虑不周（为啥没考虑，可能的因素很多），而你现在应 ...

好的，非常感谢！

作者
Author: xiaomidaxue 时间: 2019-7-30 16:46

sobereva 发表于 2019-7-29 22:37
顺带提一句：用Multiwfn载入fch文件，进入主功能0，直接在文本窗口就能看见HOMO、LUMO的各种单位下的能量以 ...

多谢sober老师，multiwfn真的很强大

欢迎光临计算化学公社 (http://bbs.keinsci.com/)