求助高斯计算分子四态SCF不收敛问题

Vivar

各位老师好！我现在在计算一批分子的量化特征，先使用B3LYP-D3/6-31G*进行结构优化，然后再使用mulitwfn将结构优化后生成的fch文件（22、-1、1）转换成四个electrons states文件，得到四个gjf文件，有一些分子的部分gjf文件提交计算后是不收敛的，根据sob老师的教程，分别尝试了加关键词如SCF(novaracc,noincfock) int=acc2e=12、SCF=XQC、SCF=Fermi、SCF=noDIIS等方法，其中只有在加上SCF=XQC后才能达到收敛，但是在加了opt=stable后，又无法收敛，请各位老师给予指导意见，谢谢！

zjxitcc

（1）有stable=opt，没有opt=stable，这个在论坛上早有人提问，并且有准确回答，自己一搜就有。
（2）SCF(conver=6)一般用于粗计算，不应做为最终结果/报道结果，其收敛后应读取轨道 继续收敛至正常收敛限。
（3）int=acc2e=12对于g16是默认的，不需要写/不需要尝试。noDIIS完全没必要尝试。
（4）你上传的log文件中SCF在conver=6下已收敛，只不过g16认为此收敛限太松，不值得后续计算和分析，因此报错。解决办法就是去掉SCF(conver=6)。可以保留（建议保留）scf=xqc stable=opt。
（5）有一点提醒一下，如果你有整体电荷为2或-2，自旋多重度为1的体系要计算，这时候要注意该计算任务是否存在对称破缺单重态，不要只按闭壳层单重态计算，也不要只写UB3LYP就以为是计算对称破缺单重态了。如果不清楚对称破缺单重态，在此论坛上搜，过往有无数帖子。

Vivar

zjxitcc 发表于 2024-9-20 10:02
（1）有stable=opt，没有opt=stable，这个在论坛上早有人提问，并且有准确回答，自己一搜就有。
（2）SCF( ...

谢谢老师的回复，SCF(conver=6)这个关键词是multiwfn默认生成的，所以我就认为这个精度已经足够了，我之前的六千多个分子都是按照这样默认去计算的，是否都需要重新计算呢？我最终需要的参数放到这个附件里面了，请问这个精度是否已经满足我的需求了呢？

zjxitcc

Vivar 发表于 2024-9-20 10:53
谢谢老师的回复，SCF(conver=6)这个关键词是multiwfn默认生成的，所以我就认为这个精度已经足够了，我之 ...

不需要重新计算，但可以续算。我个人认为需要读取现有的轨道 继续收敛至默认收敛限。当然了，也会有人认为这样做有点僵化，这个问题没有标准答案，我只能说出我个人的观点。
有一点我想大多数人会同意：下次有新任务计算时，最好采用默认SCF收敛限。这批算完的任务你想怎么处理就怎么处理。