G03 输出文件中dipole moment 表现形式问题

nuke

考察一下HF的电偶极矩。但是遇到两个值，不知有什么关联 非常感谢。
看了gaussian03用户参考手册和Exploring Chemistry with Electronic Structure Methods 第二版也没得要领。
请大神指点。
hf.gjf:
%chk=hf.chk
%mem=6MW
%nproc=1
# opt=(calcall,tight) freq=raman b3lyp/6-311++g(d,p) geom=connectivity

Title Card Required

0 1
F
H                  1              B1

   B1             0.88000000

1 2 1.0
2

输出文件中
Dipole moment (field-independent basis, Debye):
    X=     0.0000    Y=     0.0000    Z=    -1.9818 Tot=     1.9818
还有另一个带有Dipole的在输出文件末尾
.........Dipole=0.,0.,0.7796949|DipoleDeriv=-0.4474047,0.,0.,0.
,-0.4474047,0.,0.,0.,-0.357435,0.4474047,0.,0.,0.,0.4474047,0.,0.,0.,0
.357435|Polar=2.8324529,0.,2.8324529,0.,0.,5.1108361|........
我的理解是-1.9818的dipole是总的dipole只有z方向分量；
另外一个0.7796949的dipole代表什么 还是用了不同的单位。0.7796949的dipole是氢原子的dipole？

zjxitcc

既然怀疑是不同的单位，为什么不先自己计算转化一下呢····一个是Debye，另一个是原子单位。都是H-F分子的电偶极矩，没有错。

milkxx

这个问题 我也不是很明白， 我一直以为是长度用英寸时的偶极矩，但是 真是不知道为什么下面给出的是正值 我很困惑

nuke

zjxitcc 发表于 2017-5-14 11:01
既然怀疑是不同的单位，为什么不先自己计算转化一下呢····一个是Debye，另一个是原子单位。都是H-F分子 ...

1debye=0.20819434 eÅ;
-1.9818x0.20819434 eÅ=-0.4126eÅ 与0.7796949差一些

zjxitcc

nuke 发表于 2017-5-14 11:53
1debye=0.20819434 eÅ;
-1.9818x0.20819434 eÅ=-0.4126eÅ 与0.7796949差一些

1D≈ 0.393430307 ea0，注意是原子单位，是a0，玻尔半径；
1D≈ 0.20819434 eÅ，这是Angstrom，不是原子单位。

zjxitcc

milkxx 发表于 2017-5-14 11:32
这个问题 我也不是很明白， 我一直以为是长度用英寸时的偶极矩，但是 真是不知道为什么下面给出的是正值 我 ...

为什么会有正负？因为不加nosymm关键词的话，高斯根据对称性自动调整了分子的取向，这恰好是个直线分子，按照高斯的标准，高斯把它反180度了。加nosymm关键词之后就可以看到两处正负是吻合的。
PS1：这也是为什么如果要做多片段偶极矩矢量加和的话，最好加nosymm，不然加起来得到的什么都不是。
PS2：建议阅读sob这篇：谈谈Gaussian中的对称性与nosymm关键词的使用 http://sobereva.com/297

nuke

zjxitcc 发表于 2017-5-14 12:19
1D≈ 0.393430307 ea0，注意是原子单位，是a0，玻尔半径；
1D≈ 0.20819434 eÅ，这是Angstrom，不 ...

最近一直在抠后面的那部分dipole 和polarizability还有他们的derivatives   真是太赞了 谢谢指点

milkxx

zjxitcc 发表于 2017-5-14 12:30
为什么会有正负？因为不加nosymm关键词的话，高斯根据对称性自动调整了分子的取向，这恰好是个直线分子， ...

多谢 多谢