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

量化小菜鸡 · 发表于 Post on 2020-7-1 14:45:25

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

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

DLPNOCCSDTaugccpVTZ_06_36TS.out (30.8 KB, 下载次数 Times of downloads: 6)
您试过加aug的吗？我这个似乎还是不行，第一步总是基组不一样。

量化小菜鸡 · 发表于 Post on 2020-7-1 15:04:11

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

这个是fchk，mkl和out文件。您有空了帮我看看，我觉得是不是G16和ORCA对aug的弥散有区别啊。
链接：https://pan.baidu.com/s/1iEEZW6ck3nQzbFLsbmYHoQ
提取码：b5m5

zjxitcc · 发表于 Post on 2020-7-1 18:00:56

量化小菜鸡发表于 2020-7-1 15:04
这个是fchk，mkl和out文件。您有空了帮我看看，我觉得是不是G16和ORCA对aug的弥散有区别啊。
链接：http ...

网盘太慢了，我用你orca输出文件里的坐标直接建个gjf算了，比网盘下载快多了。
用fch2mkl生成的inp和mkl文件结果毫无问题，ORCA的RHF 1圈就收敛了。为节约空间，压缩包里只有gjf, log, inp, out四个文件

题外话：你2000+基函数，你确定你有相应的机器配置来算DLPNO-CCSD(T)？感觉够呛。。。

orca_examples.zip (110.63 KB, 下载次数 Times of downloads: 4)

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

zjxitcc 发表于 2020-7-1 18:00
网盘太慢了，我用你orca输出文件里的坐标直接建个gjf算了，比网盘下载快多了。
用fch2mkl生成的inp和mkl ...

谢谢了。我帮组里同学做的，他做benchmarking和方法的。他用超算能算的动，极限就到TZ了。ORCA就是吃内存太厉害了。

量化小菜鸡 · 发表于 Post on 2020-7-1 18:37:06

zjxitcc 发表于 2020-7-1 18:00
网盘太慢了，我用你orca输出文件里的坐标直接建个gjf算了，比网盘下载快多了。
用fch2mkl生成的inp和mkl ...

您有空可以试试，用默认的ORCA的aug-cc-pVTZ直接去读取算的波函数，似乎是不能一轮收敛。
我照着ORCA PrintBasis出来的基函数用G16算初猜，这种可以一轮直接收敛。非常感谢您帮助。

zjxitcc · 发表于 Post on 2020-7-1 19:44:35

本帖最后由 zjxitcc 于 2020-7-1 19:47 编辑

量化小菜鸡发表于 2020-7-1 18:37
您有空可以试试，用默认的ORCA的aug-cc-pVTZ直接去读取算的波函数，似乎是不能一轮收敛。
我照着ORCA Pr ...

我在#7楼就回复过了，最好不要用ORCA内置基组（如直接写名称），这只是有些时候能奏效，如果你经过几圈能收敛也行。但我仍强烈推荐用fch2mkl生成的inp文件，才能保证基组严格一致，也无需“ORCA PrintBasis出来的基函数用G16算初猜”这种操作。不然我让fch2mkl产生inp文件那不是多此一举了

sobereva · 发表于 Post on 2020-7-1 20:43:15

量化小菜鸡发表于 2020-7-1 14:45
您试过加aug的吗？我这个似乎还是不行，第一步总是基组不一样。

对于CH3NH2，使用下文的方法，用aug-cc-pVTZ的情况下（关键词#p aug-cc-pVTZ int=NoBasisTransform），ORCA照样一步就收敛
将Gaussian等程序收敛的波函数作为ORCA的初猜波函数的方法
http://sobereva.com/517（http://bbs.keinsci.com/thread-14905-1-1.html）

ITER    Energy       Delta-E       Max-DP    RMS-DP    [F,P]    Damp
            ***  Starting incremental Fock matrix formation  ***
                  *** Initiating the SOSCF procedure ***
                  *** Re-Reading the Fockian ***
                  *** Removing any level shift ***
ITER    Energy    Delta-E       Grad    Rot    Max-DP RMS-DP
  0 -95.25471790 -95.2547179018  0.000004  0.000004  0.000034  0.000000
            *** Restarting incremental Fock matrix formation ***
               **** Energy Check signals convergence ****
            ***Rediagonalizing the Fockian in SOSCF/NRSCF***

            *****************************************************
            *                   SUCCESS                      *
            *          SCF CONVERGED AFTER 1 CYCLES       *

结果和Gaussian给出的HF=-95.2547179完全一致，明显对aug-cc-pVTZ基组定义上并不存在区别。因此完全没必要在基组的差异性上纠结。

当体系较大，弥散函数很多的情况下，不同程序对于去除线性依赖基函数的做法不同。如果按我博文的做法在ORCA里没法一步收敛，线性依赖基函数的处理应当是导致问题的关键。
Gaussian里用IOp(3/32=2)可以避免自动去除线性依赖基函数，但对你那么大体系还用aug-cc-pVTZ，很可能导致数值不稳定。

你那个体系用aug-cc-pVTZ完全没必要，用个jun-cc-pVTZ足矣，还能显著减轻线性依赖问题，还能大幅节约耗时。

量化小菜鸡 · 发表于 Post on 2020-7-2 02:33:54

zjxitcc 发表于 2020-7-1 19:44
我在#7楼就回复过了，最好不要用ORCA内置基组（如直接写名称），这只是有些时候能奏效，如果你经过几圈能 ...

好的，谢谢。学习到新东西了

量化小菜鸡 · 发表于 Post on 2020-7-2 02:34:18

sobereva 发表于 2020-7-1 20:43
对于CH3NH2，使用下文的方法，用aug-cc-pVTZ的情况下（关键词#p aug-cc-pVTZ int=NoBasisTransform），OR ...

谢谢社长帮助

TDHFjiang · 发表于 Post on 2020-7-2 14:37:28

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

我使用的是aug-cc-pVTZ、aug-cc-pVTZ-PP基组生成的fchk, 用了MOKIT中生成的.gbw, 虽然在读取轨道信息时，出现如图1所示的问题，但是6圈后就收敛了(图2），（我没有完全采用MOKIT产生的inp, 所以会收敛较慢，通常情况下用MOKIT产生的inp应该会1-2圈收敛）C:\Users\yrjia\Desktop\Snipaste_2020-07-02_14-32-25.png C:\Users\yrjia\Desktop\Snipaste_2020-07-02_14-32-40.png

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

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

浏览过的版块