计算化学公社

标题: 求解环丙烷三个碳原子中有一个磁屏蔽张量不同 [打印本页]

作者
Author: chemicalchange 时间: 2023-9-9 14:38
标题: 求解环丙烷三个碳原子中有一个磁屏蔽张量不同
本帖最后由 chemicalchange 于 2023-9-9 15:07 编辑

各位老师好，我最近在学习核磁化学位移的计算，计算环丙烷的化学位移时注意到无论是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

作者
Author: zjxitcc 时间: 2023-9-9 14:53
能否在原帖中补充一下你在Gaussian和ORCA中用的关键词

作者
Author: chemicalchange 时间: 2023-9-9 15:04

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

已补充，谢谢

作者
Author: sobereva 时间: 2023-9-9 22:00
去掉溶剂模型再试、提升DFT积分格点精度再试

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

(, 下载次数 Times of downloads: 8)

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