求解环丙烷三个碳原子中有一个磁屏蔽张量不同

chemicalchange

各位老师好，我最近在学习核磁化学位移的计算，计算环丙烷的化学位移时注意到无论是gaussian16A.03还是orca5.0.4输出的三个碳原子中的一个碳原子（以及上面的氢）的磁屏蔽张量总是有一定差异，输入的结构确认是严格的D3h点群。gaussian计算的结果如下（NMR计算关键词见帖末）：
SCF GIAO Magnetic shielding tensor (ppm):
      1  C    Isotropic =   195.4220   Anisotropy =    49.8843
   XX=   189.1310   YX=     0.0000   ZX=    -0.0000
   XY=    -0.0000   YY=   168.4568   ZY=     0.0000
   XZ=    -0.0000   YZ=    -0.0000   ZZ=   228.6781
   Eigenvalues:   168.4568   189.1310   228.6781
      2  C    Isotropic =   195.4209   Anisotropy =    49.8694
   XX=   173.6297   YX=     8.9531   ZX=    -0.0000
   XY=     8.9460   YY=   183.9659   ZY=    -0.0000
   XZ=     0.0000   YZ=     0.0000   ZZ=   228.6672
   Eigenvalues:   168.4632   189.1324   228.6672
      3  C    Isotropic =   195.4209   Anisotropy =    49.8694
   XX=   173.6297   YX=    -8.9531   ZX=     0.0000
   XY=    -8.9460   YY=   183.9659   ZY=     0.0000
   XZ=     0.0000   YZ=     0.0000   ZZ=   228.6672
   Eigenvalues:   168.4632   189.1324   228.6672
      4  H    Isotropic =    31.2902   Anisotropy =    14.3596
   XX=    25.6786   YX=    -0.0000   ZX=     0.0000
   XY=    -0.0000   YY=    32.5671   ZY=     3.8732
   XZ=    -0.0000   YZ=     9.3114   ZZ=    35.6249
   Eigenvalues:    25.6786    27.3287    40.8633
      5  H    Isotropic =    31.2902   Anisotropy =    14.3596
   XX=    25.6786   YX=    -0.0000   ZX=    -0.0000
   XY=    -0.0000   YY=    32.5671   ZY=    -3.8732
   XZ=     0.0000   YZ=    -9.3114   ZZ=    35.6249
   Eigenvalues:    25.6786    27.3287    40.8633
      6  H    Isotropic =    31.2888   Anisotropy =    14.3552
   XX=    30.8456   YX=     2.9824   ZX=     3.3536
   XY=     2.9803   YY=    27.4002   ZY=     1.9385
   XZ=     8.0593   YZ=     4.6575   ZZ=    35.6206
   Eigenvalues:    25.6796    27.3278    40.8590
      7  H    Isotropic =    31.2888   Anisotropy =    14.3552
   XX=    30.8456   YX=     2.9824   ZX=    -3.3536
   XY=     2.9803   YY=    27.4002   ZY=    -1.9385
   XZ=    -8.0593   YZ=    -4.6575   ZZ=    35.6206
   Eigenvalues:    25.6796    27.3278    40.8590
      8  H    Isotropic =    31.2888   Anisotropy =    14.3552
   XX=    30.8456   YX=    -2.9824   ZX=     3.3536
   XY=    -2.9803   YY=    27.4002   ZY=    -1.9385
   XZ=     8.0593   YZ=    -4.6575   ZZ=    35.6206
   Eigenvalues:    25.6796    27.3278    40.8590
      9  H    Isotropic =    31.2888   Anisotropy =    14.3552
   XX=    30.8456   YX=    -2.9824   ZX=    -3.3536
   XY=    -2.9803   YY=    27.4002   ZY=     1.9385
   XZ=    -8.0593   YZ=     4.6575   ZZ=    35.6206
   Eigenvalues:    25.6796    27.3278    40.8590
请问这是数值误差导致的吗？结构是对称的为什么有一个碳的磁屏蔽张量会有差异呢？
gaussian
关键词
nmr wB97XD/def2SVP scrf(solvent=chloroform)
orca
关键词
! wB97X-V def2-SVP def2/J RIJCOSX nmr tightSCF noautostart miniprint nopop
%maxcore  2800
%pal nprocs  36 end
! CPCM(Chloroform)
坐标
C                  0.00000000    0.86830100    0.00000000
C                  0.75197000   -0.43415000    0.00000000
C                 -0.75197000   -0.43415000    0.00000000
H                  0.00000000    1.46165100    0.91788000
H                  0.00000000    1.46165100   -0.91788000
H                 -1.26582600   -0.73082500   -0.91788000
H                 -1.26582600   -0.73082500    0.91788000
H                  1.26582600   -0.73082500    0.91788000
H                  1.26582600   -0.73082500   -0.91788000

zjxitcc

能否在原帖中补充一下你在Gaussian和ORCA中用的关键词

chemicalchange

zjxitcc 发表于 2023-9-9 14:53
能否在原帖中补充一下你在Gaussian和ORCA中用的关键词

已补充，谢谢

sobereva

去掉溶剂模型再试、提升DFT积分格点精度再试

北京科音基础量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）里也专门提到过这个问题