24小时热门版块排行榜    

查看: 416  |  回复: 1

yongleli

木虫 (正式写手)


[交流] 【求助】Gaussian中的坐标系统

最近做了一个dipole的计算,

因为计算的是CN-这个带电体系,所以需要坐标转换。

Nosymm优化,坐标,Mulliken charge:
    0.000000    0.000000    0.944109   -0.669314
    0.000000    0.000000    2.135038   -0.330686

根据公式

Rc=\sum{Ri*qi}/\sum{qi}       (1)

先求出其电荷中心,然后用

D=\sum{(Ri-Rc)*qi}

计算dipole,

跟在此基础上再启用Symmetry计算一个单点(电荷一样)
坐标,Mulliken电荷:
    0.000000    0.000000   -0.641269   -0.669314
    0.000000    0.000000    0.549659   -0.330686

然后直接
D=\sum{Ri*qi}                       (2)

计算得到的dipole不一样。怎么回事呢?
而且,用这个坐标,计算得到的Rc不是坐标原点!
按理说这套坐标应该是转换之后的,我再计算Rc应该是原点才对呀?

P. S. 使用了启用Symmetry的坐标系,直接用(2)计算得到的dipole
跟从头算得到的dipole是一致的。而用NoSymm的结果结合(1)式,
计算得到的dipole就差别很大。
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

ykwang

金虫 (正式写手)


yongleli(金币+1): 2010-12-20 00:39:34
yongleli(金币+1): 好久没有来了,这个贴子一直没有更多的回帖,这个金币给你了! 2011-04-14 18:54:00
对于带电荷的体系。电偶极矩是原点依赖性的,与是否使用Symmetry没有关系。
2楼2010-12-19 22:11:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yongleli 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见