求助Dalton报错QC-SCF not implemented/tested for DFT

Voland

如题，最近我在用Dalton进行磷光光谱计算，这里根据Sob老师的意见，先用Gaussian进行Tn的几何优化，然后再通过Dalton计算振子强度等数据。但是过程中遇到了QC-SCF not implemented/tested for DFT的报错，看起来是QC-SCF目前并不支持DFT方法。那么，在此之前首先是SCF的DIIS不收敛的问题，我看了下手册，结合Sob老师在Gaussian中提到的能够解决SCF不收敛的建议，调整了.SHIFT（目前设置为-4.0），初猜结构也是Gaussian的DFT方法激发态几何优化收敛所得（泛函B3LYP，我注意到了Dalton这里的差别，不过这里不是重点），但是仍然感觉收敛非常缓慢，无意义地增加SCF ITER也不是个好办法（默认60）。请问一下还有什么办法可以加快收敛吗？
实在没辙那我只能用TDHF之类的方法做了，也不知道能不能行。
十分希望收到大家的建议，非常感谢！

Voland

噢，不好意思忘了说了，现在60圈大概收敛到D-2的水平吧，距离D-5还有相当的距离

zjxitcc

这还不简单用高斯算完之后，用fch2dal小程序生成.dal和.mol文件，Dalton飞速SCF收敛。举个简单的例子，高斯输入文件

1. %chk=h2o.chk
   
2. %mem=4GB
   
3. %nprocshared=2
   
4. #p RHF/def2TZVP nosymm int=nobasistransform
   
5. 
   
6. Title Card Required
   
7. 
   
8. 0 1
   
9. O       -0.31210987    0.61797752    0.0
   
10. H        0.64789013    0.61797752    0.0
    
11. H       -0.63256446    1.52291335    0.0

复制代码

算完之后执行

1. fch2dal h2o.chk

复制代码

生成h2o.dal和h2o.mol。交给Dalton做计算，SCF 1圈收敛。

如果你要算B3LYP，gjf文件里就把RHF改成B3LYP，然后打开生成的.dal文件，把其中的方法部分改成

1. **WAVE FUNCTIONS
   
2. .DFT
   
3. B3LYP

复制代码

其他你认为必要的关键词也可以自行添加。B3LYP的话，可能要花个5圈收敛（这是由于各个程序的默认积分格点不一样，B3LYP定义也未必一样）。

fch2dal小程序免费、开源，见https://gitlab.com/jxzou/mokit。也可以到https://gitlab.com/jxzou/mokit/-/releases下载Win预编译版本，直接使用。

Voland

zjxitcc 发表于 2022-8-20 13:46
这还不简单用高斯算完之后，用fch2dal小程序生成.dal和.mol文件，Dalton飞速SCF收敛。举个简单的例子 ...

非常感谢，然而我觉得你应该没明白我的意思，我的重点不是说输入文件，而是说现在面临Dalton的DIIS不收敛的问题，您拿一个水来算我自然也能正常运行并得到预期的结果，但是我的体系并不是那么简单的体系，而这个问题您是否有什么好办法呢？

zjxitcc

Voland 发表于 2022-8-20 14:13
非常感谢，然而我觉得你应该没明白我的意思，我的重点不是说输入文件，而是说现在面临Dalton的DIIS不收敛 ...

上面只是举个例子，我要是贴个100原子坐标，来几个过渡金属，大家鼠标上下滚动老久、不方便看帖啊。你体系就算再复杂，你用高斯能收敛的话，传给Dalton都是一样的操作步骤啊，没什么不同。
fch2dal不仅产生Dalton输入文件，文件里面还包含高斯算好的轨道（经过转换的）的，所以才能让Dalton极速收敛。

我算过很多复杂的例子，传轨道的方法也包括CASSCF，操作上没什么不同。

Voland

zjxitcc 发表于 2022-8-20 14:18
我上面只是举个例子啊，我要是贴个100原子坐标，来几个过渡金属，大家鼠标上下滚动老久、不方便看帖啊 ...

行，了解了，我试试看，多谢

zjxitcc

Voland 发表于 2022-8-20 14:23
行，了解了，我试试看，多谢

注意高斯gjf提交计算之前，输入文件关键词最好加上nosymm int=nobasistransform哈，这样传轨道更保险，Dalton收敛更快。

Voland

zjxitcc 发表于 2022-8-20 14:27
注意高斯gjf提交计算之前，输入文件关键词最好加上nosymm int=nobasistransform哈，这样传轨道更保险，Da ...

非常感谢，已经拿到能用的结果了。
祝您的MOKIT开发一切顺利！

zhouyulin

Voland 发表于 2022-8-21 22:45
非常感谢，已经拿到能用的结果了。
祝您的MOKIT开发一切顺利！

楼主，想请教一下你是怎么解决的，我这边也遇到了相同的问题

wzkchem5

zhouyulin 发表于 2025-6-27 13:33
楼主，想请教一下你是怎么解决的，我这边也遇到了相同的问题

楼上几个帖子已经把解决方案写得清清楚楚的了

zhouyulin

wzkchem5 发表于 2025-6-27 13:58
楼上几个帖子已经把解决方案写得清清楚楚的了

谢谢老师，只不过我对原先gaussian生成的T1的chk文件已经做了nosymm int=nobasistransform guess=read处理得到fchk文件，最后通过multiwfn程序得到mol文件，但是还是出现了相同的报错，所以想问一下还有什么解决方法吗

zjxitcc

zhouyulin 发表于 2025-6-27 15:39
谢谢老师，只不过我对原先gaussian生成的T1的chk文件已经做了nosymm int=nobasistransform guess=read处 ...

3L的回答写得非常详细，是用MOKIT的fch2dal小程序 从Gaussian传轨道给Dalton，建议您再阅读几遍。如果你觉得仍不过瘾，可以阅读
《利用MOKIT从Gaussian向其他量化程序传轨道》
fch2dal小程序会同时产生.dal和.mol文件，包含体系坐标，基组和赝势（如果有用到赝势的话），转化好的轨道系数。如果觉得哪里看不懂，可以具体问。

zhouyulin

zjxitcc 发表于 2025-6-27 15:50
3L的回答写得非常详细，是用MOKIT的fch2dal小程序 从Gaussian传轨道给Dalton，建议您再阅读几遍。如果你 ...

谢谢老师，我是第一次接触MOKIT程序，所以并不清楚。这边试着使用windows的fch2dal程序，将fchk文件拖入fch2dal.exe中，但是不知道为什么总是无法运行，这种情况是我的操作有问题吗

zhouyulin

zjxitcc 发表于 2025-6-27 15:50
3L的回答写得非常详细，是用MOKIT的fch2dal小程序 从Gaussian传轨道给Dalton，建议您再阅读几遍。如果你 ...

而且不知道为什么fch2inp.exe文件能够运行，并且能够生成相应的inp文件

zhouyulin

zjxitcc 发表于 2025-6-27 15:50
3L的回答写得非常详细，是用MOKIT的fch2dal小程序 从Gaussian传轨道给Dalton，建议您再阅读几遍。如果你 ...

老师，抱歉打扰一下，有个关于MOKIT的问题我还是不太理解，fch2inp.exe和fch2dal.exe这两个程序应该都属于能够直接在window下使用的小程序吧。为什么将fchk文件拖入fch2inp.exe后，能够生成相应的inp文件，而将fchk文件拖入fch2dal.exe后，则会直接闪退呢，是我的使用方法存在问题。而且不足fch2inp.exe，其它的几个exe程序也出现了类似的情况