计算化学公社

标题: 请问orca这种scf的情况有啥好办法吗？ [打印本页]

作者
Author: 413 时间: 2023-1-2 22:38
标题: 请问orca这种scf的情况有啥好办法吗？
本帖最后由 413 于 2023-1-2 23:11 编辑

大佬们，第一圈scf，这是不是只能调整初始结构了？可是这个结构是从文献中下载的想必是已经优化过的。! BP86 opt lanl2dz angs PRINTBASIS PRINTMOS def2/J  defgrid3 VERYTIGHTSCF
去掉verytightscf等没有改善

  --------------------------------------------------------------------------------------------
Iter.       energy          ||Error||_2       Shift    TRadius  Mac/Mic       Rej.
--------------------------------------------------------------------------------------------
WARNING: 11 diagonal Hessian elements are negative!
0    -19933.947993575650    1.851715e-01                   0.400 (TRAH MAcro) No
WARNING : negative HOMO - LUMO gap : -0.002214
0    dE -1.965799e-03    1.669000e-01 -3.9298e-03    0.021 (TRAH MIcro)
0    dE -2.306444e-03    8.106204e-02 -4.6078e-03    0.033 (TRAH MIcro)
0    dE -2.401580e-03    4.452627e-02 -4.7970e-03    0.036 (TRAH MIcro)
0    dE -2.430933e-03    1.904166e-02 -4.8563e-03    0.034 (TRAH MIcro)
0    dE -2.437123e-03    1.252820e-02 -4.8690e-03    0.033 (TRAH MIcro)
0    dE -2.442887e-03    7.360761e-03 -4.8808e-03    0.032 (TRAH MIcro)
1    -19933.950433610327    7.444445e-03                   0.480 (TRAH MAcro) No
WARNING : small HOMO - LUMO gap : 0.000426
1    dE -6.315391e-07    8.501462e-03 -1.2631e-06    0.000 (TRAH MIcro)
1    dE -3.063033e-06    4.336288e-03 -6.1260e-06    0.004 (TRAH MIcro)
1    dE -4.222101e-06    3.809155e-03 -8.4440e-06    0.005 (TRAH MIcro)
1    dE -6.699589e-06    4.046131e-03 -1.3398e-05    0.009 (TRAH MIcro)
1    dE -7.850745e-06    3.212134e-03 -1.5699e-05    0.012 (TRAH MIcro)
1    dE -8.490762e-06    1.903170e-03 -1.6978e-05    0.014 (TRAH MIcro)
1    dE -8.711531e-06    1.064994e-03 -1.7419e-05    0.015 (TRAH MIcro)
1    dE -8.778899e-06    6.406702e-04 -1.7554e-05    0.015 (TRAH MIcro)
2    -19933.950442507663    6.443776e-04                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000485
2    dE -1.918516e-08    3.887428e-04                         (NR MIcro)
2    dE -3.090945e-08    2.089152e-04                         (NR MIcro)
2    dE -3.261710e-08    1.646125e-04                         (NR MIcro)
2    dE -3.421582e-08    1.171871e-04                         (NR MIcro)
2    dE -3.492489e-08    8.970408e-05                         (NR MIcro)
2    dE -3.543493e-08    6.248423e-05                         (NR MIcro)
2    dE -3.577637e-08    2.566566e-05                         (NR MIcro)
2    dE -3.582480e-08    1.554480e-05                         (NR MIcro)
2    dE -3.583287e-08    1.338854e-05                         (NR MIcro)
2    dE -3.584741e-08    9.265650e-06                         (NR MIcro)
2    dE -3.585322e-08    5.858820e-06                         (NR MIcro)
2    dE -3.585581e-08    3.938254e-06                         (NR MIcro)
3    -19933.950442515230    2.106220e-05                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
3    dE -3.092034e-11    1.371533e-05                         (NR MIcro)
3    dE -4.345352e-11    5.359065e-06                         (NR MIcro)
3    dE -4.487958e-11    3.418172e-06                         (NR MIcro)
3    dE -4.529537e-11    2.506082e-06                         (NR MIcro)
3    dE -4.559169e-11    1.602808e-06                         (NR MIcro)
3    dE -4.574939e-11    7.567025e-07                         (NR MIcro)
3    dE -4.578447e-11    3.382672e-07                         (NR MIcro)
3    dE -4.579159e-11    1.881016e-07                         (NR MIcro)
4    -19933.950442505899    1.027344e-05                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
4    dE -9.512418e-12    2.299202e-06                         (NR MIcro)
4    dE -9.890366e-12    4.536606e-07                         (NR MIcro)
4    dE -9.899325e-12    3.091530e-07                         (NR MIcro)
4    dE -9.905387e-12    1.344416e-07                         (NR MIcro)
5    -19933.950442491743    9.976436e-06                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
5    dE -9.094870e-12    8.182287e-07                         (NR MIcro)
5    dE -9.132301e-12    4.268049e-07                         (NR MIcro)
5    dE -9.142867e-12    2.216063e-07                         (NR MIcro)
5    dE -9.145217e-12    1.448404e-07                         (NR MIcro)
6    -19933.950442530720    8.562137e-06                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
6    dE -6.826804e-12    1.677517e-06                         (NR MIcro)
6    dE -7.025496e-12    3.402775e-07                         (NR MIcro)
6    dE -7.030768e-12    2.315498e-07                         (NR MIcro)
6    dE -7.033766e-12    1.327943e-07                         (NR MIcro)
7    -19933.950442521949    8.439129e-06                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
7    dE -6.867907e-12    7.980775e-07                         (NR MIcro)
7    dE -6.907392e-12    3.436733e-07                         (NR MIcro)
7    dE -6.914173e-12    1.981457e-07                         (NR MIcro)
8    -19933.950442496975    9.491528e-06                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
8    dE -8.252073e-12    1.710341e-06                         (NR MIcro)
8    dE -8.457295e-12    6.994985e-07                         (NR MIcro)
8    dE -8.481940e-12    4.487147e-07                         (NR MIcro)
8    dE -8.493569e-12    2.215717e-07                         (NR MIcro)
8    dE -8.496988e-12    8.342446e-08                         (NR MIcro)
9    -19933.950442543046    1.095220e-05                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
9    dE -1.089342e-11    4.308515e-07                         (NR MIcro)
9    dE -1.090467e-11    2.136815e-07                         (NR MIcro)
9    dE -1.090670e-11    1.458736e-07                         (NR MIcro)
10    -19933.950442531743    1.147354e-05                         (NR MAcro)
WARNING : small HOMO - LUMO gap : 0.000486
10    dE -1.170280e-11    1.255824e-06                         (NR MIcro)
10    dE -1.181436e-11    5.678295e-07                         (NR MIcro)
10    dE -1.183452e-11    2.819044e-07                         (NR MIcro)
10    dE -1.183850e-11    1.724793e-07                         (NR MIcro)
11    -19933.950442537745    1.098830e-05                         (NR MAcro

作者
Author: zjxitcc 时间: 2023-1-2 23:25
本帖最后由 zjxitcc 于 2023-10-4 16:46 编辑

很简单，用高斯对当前结构算个单点，然后用fch2mkl小程序传轨道一下，ORCA极速收敛，然后你就可以继续做结构优化了。不过所有原子全用LANL2DZ未免有点粗暴

作者
Author: 413 时间: 2023-1-2 23:43

zjxitcc 发表于 2023-1-2 23:25
很简单，用高斯对当前结构算个单点，然后用fch2mkl小程序传轨道一下，ORCA极速收敛，然后你就可以继续做结 ...

感谢大师指点
我简直是猪脑，竟然忘了mokit这么好用的工具

作者
Author: sobereva 时间: 2023-1-3 00:56
没具体体系细节没法说，起码传个完整的输入文件

默认的SCF都不收敛，哪能再加verytightscf雪上加霜

杂化泛函往往比纯泛函SCF收敛容易，没绝对必要甭用纯泛函。而且就算用纯泛函，ORCA里用B97-3c也是比当前级别明显更好的选择。

trah经常表现得很拉胯。notrah再试，还不行加上KDIIS再试

作者
Author: 413 时间: 2023-1-3 02:04
本帖最后由 413 于 2023-1-3 02:05 编辑

sobereva 发表于 2023-1-3 00:56
没具体体系细节没法说，起码传个完整的输入文件

默认的SCF都不收敛，哪能再加verytightscf雪上加霜

谢谢老板回复

加上verytightscf是因为之前发现orca默认的收敛标准有时候会跟gaussian、以及实验结论相违背，有时候提高精度会有改善

换了pbe，用PBE或者BP86都是同流合污的选择

也试过notrah，能量崩了，还不如加上trah看起来有希望些，虽然trah也没能收敛

作者
Author: sobereva 时间: 2023-1-3 03:45

413 发表于 2023-1-3 02:04
谢谢老板回复

加上verytightscf是因为之前发现orca默认的收敛标准有时候会跟gaussian、以及实验结论相 ...

撑死了tightSCF也就足够了，verytightSCF完全没必要
换PBE没意义，都是纯泛函，SCF收敛性不会有多大差别
KDIIS比Trah见效几率大得多得多

作者
Author: 413 时间: 2023-1-3 08:21

sobereva 发表于 2023-1-3 03:45
撑死了tightSCF也就足够了，verytightSCF完全没必要
换PBE没意义，都是纯泛函，SCF收敛性不会有多大差别 ...

好的，社长。
我按照您的建议改下试试

作者
Author: 冰释之川 时间: 2023-1-3 08:41

zjxitcc 发表于 2023-1-2 23:25
很简单，用高斯对当前结构算个单点，然后用fch2mkl小程序传轨道一下，ORCA极速收敛，然后你就可以继续做结 ...

话说，需不需要高斯跑单点的时候计算级别保持一致？← v ←

作者
Author: zjxitcc 时间: 2023-1-3 13:21

冰释之川发表于 2023-1-3 08:41
话说，需不需要高斯跑单点的时候计算级别保持一致？← v ←

要的，比如这个例子里是BP86/LANL2DZ

作者
Author: 冰释之川 时间: 2023-1-3 15:17

zjxitcc 发表于 2023-1-3 13:21
要的，比如这个例子里是BP86/LANL2DZ

那如果有些泛函高斯里木有怎么办

作者
Author: zjxitcc 时间: 2023-1-3 16:37

冰释之川发表于 2023-1-3 15:17
那如果有些泛函高斯里木有怎么办

可以自定义的泛函，那就用自定义；连自定义都没办法的话，那就只能选一个最接近的了。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)