计算化学公社

标题: G03 输出文件中dipole moment 表现形式问题 [打印本页]

作者
Author: nuke 时间: 2017-5-14 10:21
标题: G03 输出文件中dipole moment 表现形式问题
考察一下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？

作者
Author: zjxitcc 时间: 2017-5-14 11:01
既然怀疑是不同的单位，为什么不先自己计算转化一下呢····一个是Debye，另一个是原子单位。都是H-F分子的电偶极矩，没有错。

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

作者
Author: nuke 时间: 2017-5-14 11:53

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

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

作者
Author: zjxitcc 时间: 2017-5-14 12:19

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，不是原子单位。

作者
Author: zjxitcc 时间: 2017-5-14 12:30

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

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

作者
Author: nuke 时间: 2017-5-14 13:20

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

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

作者
Author: milkxx 时间: 2017-5-14 13:51

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

多谢多谢

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