CyRhmU.jpeg
查看: 6043  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

春华秋实2020

金虫 (正式写手)


[交流] gauss 计算电偶极矩

优化了一个团簇,计算了电偶极矩,导师希望得到各个分子的电偶极矩,我就用matlab把各个分子的电荷与坐标相乘去计算各个分子的电偶极矩。现在又这样一个问题,我认为用matlab把所有的电荷与坐标相乘得到的应该是总的电偶极矩,但算完后发现与gauss输出文件中的dipole有误差,已经考虑了guass中dipole的单位是德拜。不知怎么回事。另外,我计算用的电荷是ATP charge。
回复此楼

» 收录本帖的淘帖专辑推荐

激发态计算 QC Learning

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家


★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
gmy1990: 金币+3 2013-09-17 03:02:10
春华秋实2020: 金币+1 2013-09-17 21:38:49
春华秋实2020: 金币+2 2013-09-17 21:39:32
所谓的ATP charge有两种,一个是GAPT,Gaussian在做频率计算的时候会打印此电荷,但是GATP charge的计算方法没有理论依据,结果也不太好。另一个是意大利人提出的“IR charge”,不过这个名字误导人,其实并未用到红外强度信息。IR charge理论上比较严谨,能够重复电偶极矩,但是仅适用于平面分子。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2013-09-16 20:18:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上


★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
春华秋实2020: 金币+1 2013-09-16 19:18:14
gmy1990: 金币+2 2013-09-17 03:02:03
高斯中计算dipole的估计是Mulliken电荷.
另, 为啥不直接用高斯中的偶极矩了?
前途光明,出路难觅!
2楼2013-09-13 15:45:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)


我用mulliken 电荷计算差的更多,gauss计算的结果是整个分子团簇的电偶极矩,我想把每个分子的电偶极矩都算出来
3楼2013-09-16 19:19:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)


★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
gmy1990: 金币+2 2013-09-17 03:02:20
春华秋实2020: 金币+3 2013-09-17 21:39:24
团簇内一个分子的电偶极矩不是一个有明确定义的物理量,除非相互作用很弱,否则本来就是没有确切答案的,其实差不多算一下就可以了。

电荷与坐标相乘是得不到总偶极矩的,输出的偶极矩是用电子密度积分算出来的,把原子视为点电荷来算的话结果就不一样了。
5楼2013-09-17 02:57:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见