计算化学公社

标题: Gaussian 计算Rotatory Strengths(length)符号问题 [打印本页]

作者
Author: ice 时间: 2022-3-22 11:09
标题: Gaussian 计算Rotatory Strengths(length)符号问题

师兄您好！

最近我使用Gaussian进行tddft对某两个镜像对称物质计算,

#p td=(full, nstates=10) cam-b3lyp/genecp em=gd3bj geom=checkpoint maxdisk=100GB density Integral=(FineGrid,Acc2E=10) nosymm

发现输出的结果：（以R-MeN8-EuI2tddft.log为例）

Ground to excited state transition electric dipole moments (Au):
   state       X          Y          Z       Dip. S.    Osc.
      1       0.0084    0.0002    0.0132    0.0002    0.0000
      2       -0.0308    0.0178    0.0171    0.0016    0.0001
      3       0.0044    0.1280    -0.1103    0.0286    0.0015
      4       0.0723    0.0633    0.0676    0.0138    0.0007
      5       -0.0653    0.0355    0.0398    0.0071    0.0004
      6       -0.0720    0.0362    0.0359    0.0078    0.0004
      7       0.0708    -0.0099    0.1308    0.0222    0.0012
      8       0.0010    0.0001    -0.0001    0.0000    0.0000
      9       -0.1323    0.0780    0.0760    0.0294    0.0018
      10       -0.0659    -0.1654    0.0571    0.0350    0.0021

Ground to excited state transition magnetic dipole moments (Au):
   state       X          Y          Z
      1       -0.1599    0.1211    0.0549
      2       -0.3497    0.0033    -0.6582
      3       0.8659    -1.0330    -1.2110
      4       -1.4316    1.0833    0.3504
      5       -1.0360    0.0662    -1.7479
      6       -0.9879    -0.0220    -1.8668
      7       -2.2831    2.0937    1.3268
      8       0.0708    0.0036    0.1467
      9       -1.7692    0.0421    -3.2163
      10       0.2484    0.1698    0.7498

Rotatory Strengths (R) in cgs (10**-40 erg-esu-cm/Gauss)
   state       XX       YY       ZZ    R(length)
      1       0.9472    -0.0140    -0.5120    0.1404
      2       -7.6062    -0.0420    7.9540    0.1019
      3       -2.7092    93.5297 -94.4297    -1.2031
      4       73.2168 -48.4597 -16.7479    2.6697
      5    -47.8664    -1.6598    49.1514    -0.1249
      6    -50.3024    0.5625    47.3505    -0.7965
      7    114.2403    14.6168 -122.6911    2.0553
      8       -0.0502    -0.0004    0.0124    -0.0127
      9    -165.5556    -2.3197 172.9369    1.6872
      10       11.5813    19.8680 -30.2736    0.3919

按照Rotatory Strengths(length)的计算公式

R=|μ||m|cosθ θ为电跃迁偶极和磁跃迁偶极的夹角

可以计算：

> e=c(0.0084,0.0002,0.0132)
> m=c(-0.1599,0.1211,0.0549)
> sum(e*m)/(sqrt(sum(e*e))*sqrt(sum(m*m)))
[1] -0.1826233

说明cosθ=-0.1826233， 可以计算出R= -0.1404 in cgs (10**-40 erg-esu-cm/Gauss)

刚好与打印出的结果相反，我不清楚是否是由于单位转换时存在符号变换还是Gaussian在计算中有无特殊处理，非常希望能得到及时回复！

作者
Author: sobereva 时间: 2022-4-12 01:16
看http://bbs.keinsci.com/thread-28865-1-1.html里我的回复

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