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

好的，谢谢老师们

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

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

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

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

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

复制代码

这二者有差异吗？

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

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

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

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

zjxitcc 发表于 2022-12-17 14:13
在高斯中用同样的泛函和基组算单点（加关键词nosymm int=nobasistransform），算完获得chk或fch文件，例 ...

非常感谢您详细的解答！

youxx189 发表于 2022-12-17 13:22
您好 请教一下   怎样实现高斯算好的轨道传到orca啊

在高斯中用同样的泛函和基组算单点（加关键词nosymm int=nobasistransform），算完获得chk或fch文件，例如叫h2o.fch，执行

1. fch2mkl h2o.fch

复制代码

生成两个文件h2o_o.mkl和h2o_o.inp，内含坐标、基组和轨道信息。你只需打开h2o_o.inp文件加上你要的泛函、其他你想要的关键词。执行

1. orca_2mkl h2o_o -gbw

复制代码

生成ORCA波函数文件h2o_o.gbw。提交h2o_o.inp给ORCA，SCF 1圈收敛。fch2mkl小程序是开源的，见https://gitlab.com/jxzou/mokit

zjxitcc 发表于 2021-9-30 16:04
如果你收敛有困难（包括 成功收敛需要很多圈），一定要检验波函数稳定性，此时电子结构复杂，可能辛辛苦 ...

您好 请教一下   怎样实现高斯算好的轨道传到orca啊  

Freeman 发表于 2021-9-30 09:00
现在ORCA5推出了TRAH-SCF，声称绝对（always）可以收敛。不过亲测特别特别慢，看图

准确来说，是在没有数值噪音问题的情况下，数学上保证可以收敛。如果格点不够大或者积分精度不够高，不排除即便TRAH也不收敛。可以类比高斯的scf=qc，也是牺牲了迭代次数换取收敛概率。但和scf=qc容易收敛到不稳定解不同，用TRAH收敛到不稳定解的概率反而比不用TRAH低。

Freeman 发表于 2021-9-30 16:00
现在ORCA5推出了TRAH-SCF，声称绝对（always）可以收敛。不过亲测特别特别慢，看图

如果你收敛有困难（包括 成功收敛需要很多圈），一定要检验波函数稳定性，此时电子结构复杂，可能辛辛苦苦得到的也未必是稳定波函数，还要再优化波函数。

PS：你图中的情况，还不如直接用高斯算完了传轨道给ORCA来得快。

现在ORCA5推出了TRAH-SCF，声称绝对（always）可以收敛。不过亲测特别特别慢，看图