24小时热门版块排行榜    

查看: 561  |  回复: 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 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6PrO1XyZ5URF 2026-11-13 3/150 2026-11-15 16:23 by 0aDbZAghdVq9
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 TZqyNTBa07NP 2026-11-13 4/200 2026-11-15 10:08 by xSrQNSSjcgoX
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:22 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:05 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:20 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 9/450 2026-11-12 11:22 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:27 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 07:09 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
信息提示
请填处理意见