计算化学公社

标题: ORCA计算结构不收敛 [打印本页]

作者
Author: dcckin 时间: 2016-9-18 10:51
标题: ORCA计算结构不收敛
各位老师，我想问一下，我用ORCA 算个结构，结果在100次迭代后，SCF一直相差在10-4，不能收敛，该怎么办呢，
如果不直接算结构，可以不可以先算其他参量比如能量，将能量最低点的结构作为出猜，还有没有其他方法，谢谢！

作者
Author: Sevenyhlee 时间: 2016-9-18 12:02
不说出具体的体系和计算方法，或直接贴出inp文件，不太好说针对的办法。
常规思路应该是 1 检查结构是否合理；2 看看计算方法是否适合当前体系；3 关键词行加入SlowConv，同时加大scf的maxiter试一下？
ORCA的Input Library文档中有收敛问题的解决方案，可详细对照。

https://sites.google.com/site/orcainputlibrary/

作者
Author: dcckin 时间: 2016-9-18 15:14

Sevenyhlee 发表于 2016-9-18 12:02
不说出具体的体系和计算方法，或直接贴出inp文件，不太好说针对的办法。
常规思路应该是 1 检查结构是否合 ...

这是我的输入文件，不吝赐教
另外，我的网址打不开，能不能传个PDF文件给我，谢谢

!UKS B3LYP SV(P) SV/J OPT TightSCF grid5 VerySlowConv LSHIFT KDIIS SOSCF scfconv6

%basis newgto Pr
  0 6
1    961538.755282000056       0.006478572811
2    427350.557902999979       -0.003087733061
3    189933.581290000002       0.012936227713
4    84414.925017999994       0.005951443460
5    37517.744451999999       0.027640099371
6    16674.553090000001       0.040309057385
  0 1
7    7410.912484000000       1.000000000000
  0 1
8    3293.738882000000       1.000000000000
  0 1
9    1463.883948000000       1.000000000000
  0 1
10       650.615088000000       1.000000000000
  0 1
11       289.162261000000       1.000000000000
  0 1
12       128.516561000000       1.000000000000
  0 1
13       57.118471000000       1.000000000000
  0 1
14       25.385987000000       1.000000000000
  0 1
15       11.282661000000       1.000000000000
  0 1
16       5.014516000000       1.000000000000
  0 1
17       2.228674000000       1.000000000000
  0 1
18       0.990522000000       1.000000000000
  0 1
19       0.440232000000       1.000000000000
  0 1
20       0.195659000000       1.000000000000
  0 1
21       0.086959000000       1.000000000000
  0 1
22       0.038649000000       1.000000000000
  0 1
23       0.017177000000       1.000000000000
  1 5
1    12926.877106000000       0.003978702409
2    5170.750842000000       0.005409657504
3    2068.300337000000       0.027007359808
4       827.320135000000       0.086094511882
5       330.928054000000       0.260534549042
  1 1
6       132.371222000000       1.000000000000
  1 1
7       52.948489000000       1.000000000000
  1 1
8       21.179395000000       1.000000000000
  1 1
9       8.471758000000       1.000000000000
  1 1
10       3.388703000000       1.000000000000
  1 1
11       1.355481000000       1.000000000000
  1 1
12       0.542193000000       1.000000000000
  1 1
13       0.216877000000       1.000000000000
  1 1
14       0.086751000000       1.000000000000
  1 1
15       0.034700000000       1.000000000000
  1 1
16       0.013880000000       1.000000000000
  2 4
1       877.609164000000       0.004777805591
2       319.130605000000       0.027279230492
3       116.047493000000       0.148384186060
4       42.199088000000       0.459132970274
  2 1
5       15.345123000000       1.000000000000
  2 1
6       5.580045000000       1.000000000000
  2 1
7       2.029107000000       1.000000000000
  2 1
8       0.737857000000       1.000000000000
2 1
9       0.268312000000       1.000000000000
  2 1
10       0.097568000000       1.000000000000
  2 1
11       0.035479000000       1.000000000000
  2 1
12       0.012902000000       1.000000000000
  3 4
1       20.929801000000       0.103680288430
2       6.976600000000       0.316043027300
3       2.325533000000       0.483844050659
4       0.775178000000       0.364645629683
  3 1
5       0.258393000000       0.158110748793
  3 1
6       0.086131000000       1.000000000000
  end
end
%scf
         maxiter 500
         SOSCFStart 0.00033
end
* xyz 0 2

Pr                0.00000000 0.00000000 4.49209540
Ni                0.00000000 0.00000000 0.00000000
O                0.00000000 1.86800000 0.00000000
O                0.00000000 0.00000000 2.05524800
Pr                0.00000000 0.00000000 8.03990460
Sr                1.86800000 1.86800000 10.75809540
Sr                1.86800000 1.86800000 1.77390460
Al                1.86800000 1.86800000 6.26600000
O                1.86800000 0.00000000 0.00000000
O                1.86800000 0.00000000 6.26600000
O                0.00000000 1.86800000 6.26600000
O                0.00000000 0.00000000 10.47675200
O                1.86800000 1.86800000 8.32124800
O                1.86800000 1.86800000 4.21075200
*

作者
Author: sobereva 时间: 2016-9-18 16:30
到底是几何优化不收敛还是SCF不收敛？
“不直接算结构”是指不优化结构？什么叫将能量最低点的结构作为出猜？不优化怎么得到能量最低点？

作者
Author: dcckin 时间: 2016-9-19 08:55
谢谢你的回答。我是SCF不收敛，我想做的是优化结构，但输出文件出错，提示SCF不收敛。我意思是有没有其他办法，得到一个较合适的初步结构，再进行结构优化。

作者
Author: fwan1006 时间: 2022-11-29 13:32

sobereva 发表于 2016-9-18 16:30
到底是几何优化不收敛还是SCF不收敛？
“不直接算结构”是指不优化结构？什么叫将能量最低点的结构作为出 ...

sob老师，几何优化不收敛的话有什么通用的解决方法吗

作者
Author: wzkchem5 时间: 2022-11-29 15:12

fwan1006 发表于 2022-11-29 06:32
sob老师，几何优化不收敛的话有什么通用的解决方法吗

http://sobereva.com/164

作者
Author: sobereva 时间: 2022-11-30 09:06

fwan1006 发表于 2022-11-29 13:32
sob老师，几何优化不收敛的话有什么通用的解决方法吗

7L贴的我的博文里，使用较小步长上限、尝试不同的优化算法、用更好的Hessian等等，对于各种量化程序都是共通的

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