Gaussian 计算Rotatory Strengths(length)符号问题

ice

师兄您好！

最近我使用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在计算中有无特殊处理，非常希望能得到及时回复！

sobereva

看http://bbs.keinsci.com/thread-28865-1-1.html里我的回复