NBO7 在orca4.2与G16A结果不同

lqhenwunai

我通过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.

请问这是怎么回事？