计算化学公社

标题: NBO7 在orca4.2与G16A结果不同 [打印本页]

作者
Author: lqhenwunai 时间: 2020-5-28 17:22
标题: NBO7 在orca4.2与G16A结果不同
本帖最后由 lqhenwunai 于 2020-5-28 17:23 编辑

我通过G16A和orca4.2分别调用NBO7进行计算。

ORCA4.2的关键词为：
!pal BP86 NBO  cc-pvtz nori nopop
最后的能量为FINAL SINGLE POINT ENERGY    -79.853600707855

G16A的关键词为
#p bp86/cc-pvtz pop=nbo6
最后的能量为SCF Done:  E(RB-P86) =  -79.8534982765

能量的不同可能来自格点精度的的不同，这还好。但是如果比较NBO的E2能量的话,二者就有很大不同了。对于G16A，有

SECOND ORDER PERTURBATION THEORY ANALYSIS OF FOCK MATRIX IN NBO BASIS

   Threshold for printing: 0.50 kcal/mol
                                                      E(2) E(NL)-E(L) F(L,NL)
   Donor (L) NBO             Acceptor (NL) NBO    kcal/mol a.u. a.u.
===============================================================================

within unit  1
4. BD ( 1) C  1- H  6    14. BD*( 1) C  2- H  3    2.48 0.74 0.038
5. BD ( 1) C  1- H  7    16. BD*( 1) C  2- H  5    2.48 0.74 0.038
6. BD ( 1) C  1- H  8    15. BD*( 1) C  2- H  4    2.48 0.74 0.038
6. BD ( 1) C  1- H  8    43. RY ( 2) C  2          0.50 1.52 0.025
7. BD ( 1) C  2- H  3    11. BD*( 1) C  1- H  6    2.48 0.74 0.038
8. BD ( 1) C  2- H  4    13. BD*( 1) C  1- H  8    2.48 0.74 0.038
8. BD ( 1) C  2- H  4    18. RY ( 2) C  1          0.50 1.52 0.025
9. BD ( 1) C  2- H  5    12. BD*( 1) C  1- H  7    2.48 0.74 0.038

但是对于orca4.2, 打印出的E2还包括CR等轨道：
SECOND ORDER PERTURBATION THEORY ANALYSIS OF FOCK MATRIX IN NBO BASIS

   Threshold for printing: 0.50 kcal/mol
                                                      E(2) E(NL)-E(L) F(L,NL)
   Donor (L) NBO             Acceptor (NL) NBO    kcal/mol a.u. a.u.
===============================================================================

within unit  1
1. CR ( 1) C  1          10. BD*( 1) C  1- C  2    2.29 11.83 0.147
1. CR ( 1) C  1          11. BD*( 1) C  1- H  6    1.12 11.97 0.103
1. CR ( 1) C  1          12. BD*( 1) C  1- H  7    1.12 11.97 0.103
1. CR ( 1) C  1          13. BD*( 1) C  1- H  8    1.12 11.97 0.103
1. CR ( 1) C  1          19. RY ( 3) C  1          5.58 13.31 0.243
1. CR ( 1) C  1          20. RY ( 4) C  1          1.24 14.09 0.118
1. CR ( 1) C  1          21. RY ( 5) C  1          19.96 16.93 0.519
1. CR ( 1) C  1          22. RY ( 6) C  1          9.56 15.94 0.349
1. CR ( 1) C  1          24. RY ( 8) C  1          1.47 13.56 0.126
1. CR ( 1) C  1          25. RY ( 9) C  1          1.28 13.87 0.119
1. CR ( 1) C  1          28. RY (12) C  1          2.58 13.29 0.165

特别的，由ORCA4.2/nbo7得到了
8. BD ( 1) C  2- H  4    15. BD*( 1) C  2- H  4    6.25 0.75 0.061
即从BD8 到BD*15的E2能量为6.25kcal/mol, 而在G16A中，所有的pair的E2能量都小于3kcal/mol.

请问这是怎么回事？

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