双光子吸收截面值报错

weixue

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

zjxitcc

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

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

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

weixue

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

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

zjxitcc

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

1. %chk=a.chk
   
2. #p B3LYP genecp nosymm int=nobasistransform
   
3. 
   
4. title
   
5. 
   
6. 0 1
   
7. [坐标]
   
8. 
   
9. Ru 0
   
10. SDD
    
11. ****
    
12. H C N O 0
    
13. 6-31G(d)
    
14. ****
    
15. 
    

复制代码

此时我们没加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下载代码自行编译。