为何用ORCA读取G16算出的初猜第一步能量相差极大？

量化小菜鸡 · 发表于 Post on 2020-6-30 14:18:35

我使用社长的方法用G16产生aug-cc-pVTZ的波函数去作为ORCA的初猜，目的是更好的进行下一步的DLPNO-CCSD(T)/aug-cc-pVTZ的计算。因为是aug-cc系列，我使用了 int=NoBasisTransform nosymm的关键词。最终能量为HF=-1344.3638469
但用ORCA做接下来的计算的话，可见第一步能量相差较大，虽然基组验证一样。此外，因为MAX-DP极大，第二步能量巨变。

这个结果让我觉得十分的迷惑，感觉问题出在第一步的算法上。谢谢大家帮忙。另外我发觉直接用G16的aug-cc-pVTZ的话，就算加NoBasisTransform，跟ORCA用的也有些许不同。

量化小菜鸡 · 发表于 Post on 2020-6-30 14:19:22

---------------------
INITIAL GUESS: MOREAD
---------------------
Guess MOs are being read from file: 06_36TS-guess.ges
Input Geometry matches current geometry (good)
Input basis set matches current basis set (good)
MOs were renormalized
MOs were reorthogonalized (Cholesky)
                  ------------------
                  INITIAL GUESS DONE ( 4.1 sec)
                  ------------------
--------------
SCF ITERATIONS
--------------
ITER    Energy       Delta-E       Max-DP    RMS-DP    [F,P]    Damp
            ***  Starting incremental Fock matrix formation  ***
  0  -1343.8608293347 0.0000000000001252.06273066  3.82819878  0.3028578 0.7000
  1  29624.205102595030968.065931929730371.14597152  1.13787707 1692.9559380 0.7000

hebrewsnabla · 发表于 Post on 2020-6-30 14:32:25

本帖最后由 hebrewsnabla 于 2020-6-30 14:37 编辑

首先基函数数目一样吗？高斯算的时候要指定5D 7F

你是用Multiwfn把fchk转成mkl么？Multiwfn的这个功能有一点问题（一部分f,g基函数对应的MO系数要乘-1），你可以用MOKIT试试（https://gitlab.com/jxzou/mokit）

搞不定就上传fchk和mkl文件

zjxitcc · 发表于 Post on 2020-6-30 14:38:47

本帖最后由 zjxitcc 于 2020-6-30 14:41 编辑

mkl文件里的有些G, F, H函数的系数需要乘以-1，这点在MOKIT里的fch2mkl小程序考虑了。
你这基组很大，光写5D是不够的，还要加7F。
你可以传压缩后的fchk文件，或者传百毒网盘。

量化小菜鸡 · 发表于 Post on 2020-6-30 14:38:55

hebrewsnabla 发表于 2020-6-30 14:32
首先基函数数目一样吗？高斯算的时候要指定5D 7F

你是用Multiwfn把fchk转成mkl么？Multiwfn的这个功能有 ...

谢谢，我用的是5D（查高斯输出了），否则ORCA第一步应该不会说match the basis set。 mkl比较大，直接上传不上去。
我正在读manual，感觉问题像出在Cholesky inversion。他讲有linear dependence这个算法容易出问题，然后能量会大跳。

量化小菜鸡 · 发表于 Post on 2020-6-30 14:40:43

zjxitcc 发表于 2020-6-30 14:38
mkl文件里的有些G, F, H函数的系数需要乘以-1，这点在MOKIT里的fch2mkl小程序考虑了。
光5D是不够的，还要 ...

好的，谢谢，那我用printbasis一个一个比比看看。但ORCA第一步说基组符合应该就是对的上吧。

zjxitcc · 发表于 Post on 2020-6-30 14:42:13

量化小菜鸡发表于 2020-6-30 14:40
好的，谢谢，那我用printbasis一个一个比比看看。但ORCA第一步说基组符合应该就是对的上吧。

不一定，这并没有必然关系。而且传轨道最好不要用ORCA内置基组（如直接写名称），而用fch2mkl生成的inp文件，才能保证基组严格一致。

量化小菜鸡 · 发表于 Post on 2020-6-30 14:46:05

zjxitcc 发表于 2020-6-30 14:42
不一定，这并没有必然关系。而且传轨道最好不要用ORCA内置基组（如直接写名称），而用fch2mkl生成的inp文 ...

好的，学习了。您不讲我都不知道

量化小菜鸡 · 发表于 Post on 2020-6-30 14:48:13

zjxitcc 发表于 2020-6-30 14:38
mkl文件里的有些G, F, H函数的系数需要乘以-1，这点在MOKIT里的fch2mkl小程序考虑了。
你这基组很大，光写 ...

7F也有的，忘记说了。我用过默认aug-cc-pVTZ，还用过gen从EMSL的搞下来的。我试试怎么分享啊，网盘在外面可慢

zjxitcc · 发表于 Post on 2020-6-30 15:19:36

量化小菜鸡发表于 2020-6-30 14:48
7F也有的，忘记说了。我用过默认aug-cc-pVTZ，还用过gen从EMSL的搞下来的。我试试怎么分享啊，网盘在外面 ...

暂时不用分享fch文件。也不用gen，你就用高斯里的aug-cc-pVTZ就行。用MOKIT里的fch2mkl，仅编译组件fch2mkl即可。程序都是开源的，自己弄一下几分钟就行了。

sobereva · 发表于 Post on 2020-6-30 19:00:53

Gaussian里用Dunning系列基组默认就用5d 7f，完全没必要刻意写，你直接在输入文件里就会看到以下信息
Standard basis: Aug-CC-pVTZ (5D, 7F)
明显可见默认就是用的5d 7f

量化小菜鸡 · 发表于 Post on 2020-7-1 03:09:34

sobereva 发表于 2020-6-30 19:00
Gaussian里用Dunning系列基组默认就用5d 7f，完全没必要刻意写，你直接在输入文件里就会看到以下信息
Stan ...

对的，默认是5D 7F。

量化小菜鸡 · 发表于 Post on 2020-7-1 03:12:20

zjxitcc 发表于 2020-6-30 15:19
暂时不用分享fch文件。也不用gen，你就用高斯里的aug-cc-pVTZ就行。用MOKIT里的fch2mkl，仅编译组件fch2m ...

谢谢，我尝试了用MOKIT，第二轮能量不跳了，但还是不能一轮两轮就收敛。至少趋势是对的 :)
ORCA有弥散的收敛好一般，高斯直接上去算就收敛了。

sobereva · 发表于 Post on 2020-7-1 06:47:25

量化小菜鸡发表于 2020-6-30 14:48
7F也有的，忘记说了。我用过默认aug-cc-pVTZ，还用过gen从EMSL的搞下来的。我试试怎么分享啊，网盘在外面 ...

Multiwfn的那个问题已经解决，见

解决了一个Multiwfn产生mkl文件作为ORCA当初猜波函数的问题
http://bbs.keinsci.com/thread-18220-1-1.html

现在对于aug-cc-pVTZ这样含d以上角动量的情况也可以非常容易地通过Multiwfn来把Gaussian收敛的波函数当ORCA的初猜波函数。如果有问题，欢迎反馈。

量化小菜鸡 · 发表于 Post on 2020-7-1 13:12:17

sobereva 发表于 2020-7-1 06:47
Multiwfn的那个问题已经解决，见

解决了一个Multiwfn产生mkl文件作为ORCA当初猜波函数的问题

谢谢社长，这会就去试试。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] 为何用ORCA读取G16算出的初猜第一步能量相差极大？

评分 Rate

浏览过的版块