ORCA中解决SCF不收敛问题的关键字

413

大佬有没有时间输出一波这个table里面新版本的orca的命令？

chenjhust

请问在高斯09用b3lyp/6-311g(d,p)优化好的初始结构，经过multiwfn转化成ORCA（ORCA4.2.1）的输入文件计算更高精度的单点能时老是出现SCF不收敛的问题
，试过kdiis，调整积分精度都无法解决，请问有什么好的建议吗？

ionexchangeC

chenjhust 发表于 2024-3-7 17:15
请问在高斯09用b3lyp/6-311g(d,p)优化好的初始结构，经过multiwfn转化成ORCA（ORCA4.2.1）的输入文件计算更 ...

换用ORCA5，默认设置会使用TRAH辅助收敛
如果还不行并且TRAH步没有收敛趋势，那尝试使用noTRAH和slowconv关键词，并加大SCF迭代步数上限

zjxitcc

chenjhust 发表于 2024-3-7 17:15
请问在高斯09用b3lyp/6-311g(d,p)优化好的初始结构，经过multiwfn转化成ORCA（ORCA4.2.1）的输入文件计算更 ...

老铁不要短时间内立即重复提问http://bbs.keinsci.com/forum.php ... 283756&fromuid=2632

解决办法不唯一，你可以尝试他人的建议，我提供我认为最好的建议：不需要试kdiis、调整积分精度、slowconv、加大SCF迭代步数上限这些。写一个gjf文件，这里以水分子为例

1. %chk=h2o.chk
   
2. %mem=10GB
   
3. %nprocshared=4
   
4. #p UPW91B95/def2TZVPP nosymm int(nobasistransform) IOp(3/76=1000005000,3/77=0000005000,3/78=0731007310)
   
5. 
   
6. title
   
7. 
   
8. 0 3
   
9. O     0.000000    0.000000    0.062007
   
10. H     0.000000   -0.783976   -0.492052
    
11. H     0.000000    0.783976   -0.492052

复制代码

内存和核数根据自己机器情况修改。提交Gaussian任务，获得chk文件。运行

1. formchk h2o.chk h2o.fch
   
2. fch2mkl h2o.fch
   
3. orca_2mkl h2o_o -gbw

复制代码

获得h2o_o.inp和h2o_o.gbw文件。打开inp文件将前3行改成目标计算，例如

1. %pal nprocs 4 end
   
2. %maxcore 2500
   
3. ! UKS TightSCF RIJCOSX grid5 gridx5 def2/J def2-TZVPP/C RI-PWPB95 D3

复制代码

然后提交给ORCA算，会自动从gbw文件中读取轨道，SCF迅速收敛，进入双杂化计算部分。如果使用ORCA 5或更高版本，将grid5 gridx5换成defgrid3。

chenjhust

zjxitcc 发表于 2024-3-7 17:42
老铁不要短时间内立即重复提问http://bbs.keinsci.com/forum.php ... 283756&fromuid=2632

解决办法不 ...

谢谢老师，用了您的方法后确实很快就收敛了，现在有个问题
我很多结果是用这个指令算的

1. ! PWPB95 D3 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX grid4 gridx4 tightSCF noautostart miniprint nopop

复制代码

换成您推荐的这个，前面的需要重新算吗？

1. ! UKS TightSCF RIJCOSX grid5 gridx5 def2/J def2-TZVPP/C PWPB95 D3

复制代码

这二者有差异吗？

wzkchem5

chenjhust 发表于 2024-3-7 11:14
谢谢老师，用了您的方法后确实很快就收敛了，现在有个问题
我很多结果是用这个指令算的

后一种写法漏写了def2-TZVPP，得先把这个补上。
之后两种写法唯一的区别就是格点大小了，为了结果可比，格点必须统一。zjxitcc说的grid5 gridx5只是举个例子，如果你其他计算都是用grid4 gridx4算的，这里自然要改成grid4 gridx4。

zjxitcc

chenjhust 发表于 2024-3-7 18:14
谢谢老师，用了您的方法后确实很快就收敛了，现在有个问题
我很多结果是用这个指令算的

wzkchem5老师已经给你解释过了，主要是格点差异。我补充一条，fch2mkl小程序产生的inp文件中含有基组数据（数据量很大，没有在帖子里展示出来，你在自己电脑上打开inp文件便可看到）。而gjf文件用的基组名称就是def2TZVPP，传过来的基组数据就代表def2TZVPP，从而不需要在这一行写def2-TZVPP基组名称，这不同于你以前的计算。

chenjhust

zjxitcc 发表于 2024-3-7 19:07
wzkchem5老师已经给你解释过了，主要是格点差异。我补充一条，fch2mkl小程序产生的inp文件中含有基组数据 ...

好的，谢谢老师们

Keys

zjxitcc 发表于 2024-3-7 17:42
老铁不要短时间内立即重复提问http://bbs.keinsci.com/forum.php ... 283756&fromuid=2632

解决办法不 ...

老师好，我在使用molclus做构象搜索时，在使用orca计算其中一个构象的单点能的过程中，SCF无法收敛。
若使用您的方法，先在Gaussian中计算，使用同样计算级别加上

nosymm int(nobasistransform) IOp(3/76=1000005000,3/77=0000005000,3/78=0731007310)

关键词得出的轨道再给orca算，最终orca计算后得出的能量和其余orca计算正常收敛的构象的能量是否具有可比性？
或者说，在Gaussian计算中改变LDA和高阶交换泛函组合比例等相关参数是否影响转交的轨道和orca自己计算的轨道的可比性？
十分感谢解答

zjxitcc

Keys 发表于 2025-9-22 21:46
老师好，我在使用molclus做构象搜索时，在使用orca计算其中一个构象的单点能的过程中，SCF无法收敛。
若 ...

nosymm int(nobasistransform)要加。IOp(3/76=1000005000,3/77=0000005000,3/78=0731007310)是搭配PW91B95使用的，目的是凑出双杂化泛函PWPB95的自洽部分，如果你不用这个双杂化泛函，则不要写这个IOP。具有可比性。

Keys

zjxitcc 发表于 2025-9-22 21:55
nosymm int(nobasistransform)要加。IOp(3/76=1000005000,3/77=0000005000,3/78=0731007310)是搭配PW91B9 ...

谢谢老师，再请问一下，如果考虑溶剂模型是否也是二者输入文件中都加入如下SMD溶剂模型关键词即可？
Gaussian：

1. SCRF(SMD,solvent=THF)

复制代码

Orca:

1. smd true
   
2. SMDsolvent "THF"
   
3. end

复制代码

zjxitcc

Keys 发表于 2025-9-23 01:40
谢谢老师，再请问一下，如果考虑溶剂模型是否也是二者输入文件中都加入如下SMD溶剂模型关键词即可？
Gau ...

是的。也可以二者都用PCM。

Keys

zjxitcc 发表于 2025-9-23 08:43
是的。也可以二者都用PCM。

老师好，我模仿您的方法进行计算，在orca计算过程中SCF仍旧无法收敛，似乎并没有从gbw文件中读取轨道，不知道哪个步骤出了问题，请您再指导一下。
具体细节：
先在Gaussian中进行计算，关键词：

1. %chk=try.chk
   
2. %mem=96GB
   
3. %nprocs=48
   
4. #p UPW91B95/def2QZVPP nosymm int(nobasistransform) IOp(3/76=1000005000,3/77=0000005000,3/78=0731007310) SCRF(SMD,solvent=THF)

复制代码

计算结果正常，输出了chk文件，按照您的方法，使用mokit组件fch2mkl得到了try_o.inp和try_o.gbw

修改try_o.inp的前三行，并添加了溶剂模型：

1. ! UKS PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX tightSCF noautostart miniprint nopop
   
2. %maxcore 5500
   
3. %pal nprocs 48 end
   
4. %cpcm
   
5. smd true
   
6. SMDsolvent "THF"
   
7. end

复制代码

最后的结果仍旧是没有收敛，输出文件如下

zjxitcc

Keys 发表于 2025-9-24 20:25
老师好，我模仿您的方法进行计算，在orca计算过程中SCF仍旧无法收敛，似乎并没有从gbw文件中读取轨道，不 ...

（1）ORCA单点计算可以自动从gbw文件中读取轨道，但你写了noautostart，让它不自动读取了，没必要写这个。注意，如果你翻看其他帖子，有的帖子会告诉你写MOread和%moinp去读取轨道，这俩也是没必要写的。
（2）fch2mkl会产生try_o.inp文件，这里面已经写好了基组信息，不需要再写def2-QZVPP，但可以写辅助基。仔细看实例操作
http://bbs.keinsci.com/forum.php ... 349929&fromuid=2632
（3）ORCA DFT计算永远建议加上关键词defgrid3。除非你已经算了一批结果都没加defgrid3，那现在为了统一精度，不得已放弃加defgrid3。如果这批任务刚开始算，建议都加defgrid3。

Keys

zjxitcc 发表于 2025-9-24 20:51
（1）ORCA单点计算可以自动从gbw文件中读取轨道，但你写了noautostart，让它不自动读取了，没必要写这个 ...

没有系统学习研究过Orca，一般任务都在用Gaussian算。所以用Orca一直都是在照葫芦画瓢，理解非常粗浅。之前的计算都没有修改积分格点，这是最后一个构象，如果实在收敛不了，我都想如换个稍低一点的计算级别把所有构象都重新算一遍
刚刚修改关键词重新计算，但仍旧是SCF不收敛，请您再看一看

1. ! UKS PWPB95 D3 def2/J def2-QZVPP/C RIJCOSX tightSCF miniprint nopop
   
2. %maxcore 5500
   
3. %pal nprocs 48 end
   
4. %cpcm
   
5. smd true
   
6. SMDsolvent "THF"
   
7. end

复制代码