24小时热门版块排行榜    

查看: 496  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息,有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 8/400 2026-08-02 15:01 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 14:58 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 18/900 2026-08-02 14:52 by jHcrHMhbKCFV
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +17 低垂的野花 2026-07-31 23/1150 2026-08-02 13:41 by RRRKKK
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
信息提示
请填处理意见