计算化学公社

标题: 双光子吸收截面值报错 [打印本页]

作者
Author: weixue 时间: 2021-11-5 09:33
标题: 双光子吸收截面值报错
sob老师，您好！我算了一个金属配合物的双光子吸收截面值，里面有Ru元素，用的SDD基组，其他原子用的6-31G*,泛函用的b3lyp。一直报下面的错误，请问应该如何解决

作者
Author: zjxitcc 时间: 2021-11-5 09:46
本帖最后由 zjxitcc 于 2021-11-5 10:56 编辑

从截图信息看是SCF不收敛，当然你可以在输出文件中再往上翻一翻确认是否SCF迭代超过最大圈数。如果是，有两种办法：
（1）查阅Dalton手册，寻找帮助SCF收敛的关键词，也可以阅读sob老师这篇博文《解决SCF不收敛问题的方法》http://sobereva.com/61寻找Dalton里有无相应功能；
（2）不查Dalton手册，在Gaussian中用相应的方法/基组算完了传轨道给Dalton，技巧在上次回帖时就介绍过了，这里再详细讲一下，先用高斯做计算

%chk=a.chk
#p B3LYP genecp nosymm int=nobasistransform
title
0 1
[坐标]
Ru 0
SDD
****
H C N O 0
6-31G(d)
****

复制代码

此时我们没加guess(only,save)，因为想让高斯算收敛。交给高斯算完后转化得到fch文件
formchk a.chk a.fch

从fch文件生成Dalton输入文件（.dal和.mol两个文件）
fch2dal a.fch

你可以在.dal文件中增加你要的关键词，高斯算好的轨道数据就在.dal文件里，基组信息在.mol文件里。fch2dal小程序免费开源，可在https://gitlab.com/jxzou/mokit下载代码自行编译。

作者
Author: weixue 时间: 2021-11-5 10:48

zjxitcc 发表于 2021-11-5 09:46
从截图信息看是SCF不收敛，当然你可以在输出文件中再往上翻一翻确认是否SCF迭代超过最大圈数。如果是，有 ...

您好！我想问一下，您遇到过这种不收敛的情况吗？我从dalton手册中没有找到具体的解决方案。方案2可以完全解决是吗？因为我们没有安装fch2dal小程序。

作者
Author: zjxitcc 时间: 2021-11-5 10:51
本帖最后由 zjxitcc 于 2021-11-5 10:55 编辑

weixue 发表于 2021-11-5 10:48
您好！我想问一下，您遇到过这种不收敛的情况吗？我从dalton手册中没有找到具体的解决方案。方案2可以完 ...

可以解决，因为我就是开发fch2dal小程序的，我能不知道么。。。高斯SCF收敛性本身很好，通常不加额外关键词就能顺利收敛；即使高斯要加关键词帮助收敛，也有很多经验和帖子供参考，还有sob老师这篇博文sob老师这篇博文《解决SCF不收敛问题的方法》http://sobereva.com/61，所以用起来得心应手。

看一例Dalton收敛到能量偏高的SCF解上去，然后从高斯传轨道过去1圈收敛得到低能解的例子《离线安装量子化学软件Dalton离线安装量子化学软件Dalton》此文第5部分

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