CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 8305  |  回复: 30
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wangyujia

木虫 (正式写手)

[交流] 【求助】用Bader计算磁矩的困惑已有13人参与

最近在本论坛上看到了用Bader计算磁矩的方法,我试过之后发现有些东西很奇怪。不知道是我的理解有误不是其它什么原因。
下面这个网页就简单的介绍了怎么用Bader计算磁矩。
http://theory.cm.utexas.edu/forum/viewtopic.php?f=1&t=529
我的理解就是把CHGCAR划分为两个文件,一个包括自旋向上的部分CHGCAR-1,一个包括自旋向下的部分CHGCAR-2。再用bader CHGCAR-2 -ref CHGCAR,就算出了各个原子的磁矩(算出来的值比较像磁矩)。按理来说,自旋向上的电子数和自旋向下的电子数之和等于总电子数。可是我发现用bader CHGCAR-1 -ref CHGCAR算出来的电子数,也就是自旋向上的电子数却等于总的电子数???
所以我觉得这种方法可能不对,想请教一下大家,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuqun

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
16179336楼: Originally posted by wangyujia at 2009-10-25 09:25:37:
最近在本论坛上看到了用Bader计算磁矩的方法,我试过之后发现有些东西很奇怪。不知道是我的理解有误不是其它什么原因。
下面这个网页就简单的介绍了怎么用Bader计算磁矩。
http://theory.cm.utexas.edu/forum/v ...

你好,我先用split程序作用CHGCAR,得到cf1和cf2两个文件,然后bader cf1 -ref CHGCAR,就不行了,一直出现

Image              PC                Routine            Line        Source            
bader              00000000004C004D  Unknown               Unknown  Unknown
bader              00000000004BEB55  Unknown               Unknown  Unknown
bader              0000000000477239  Unknown               Unknown  Unknown
bader              000000000043550D  Unknown               Unknown  Unknown
bader              0000000000434D5A  Unknown               Unknown  Unknown
bader              000000000045585B  Unknown               Unknown  Unknown
bader              0000000000452C07  Unknown               Unknown  Unknown
bader              000000000040BCFB  Unknown               Unknown  Unknown
bader              0000000000410B67  Unknown               Unknown  Unknown
bader              000000000040039C  Unknown               Unknown  Unknown
bader              00000000004002EC  Unknown               Unknown  Unknown
bader              00000000004CB3E0  Unknown               Unknown  Unknown
bader              00000000004001B9  Unknown               Unknown  Unknown
。楼主知道怎么回事吗?一定要回我啊,我自己琢磨了好几天也不知道,刚接触VASP。处于一片迷茫中,楼主一定要救救我啊
24楼2012-04-09 17:12:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 31 个回答

huangyc

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★
qasd(金币+2,VIP+0):谢谢~ 10-25 12:47
wangyujia(金币+5,VIP+0):谢谢你的回帖,激发了我的一个想法。 10-25 14:02
这个-ref我也用不好,哪位解释一下
我觉得你分三次不就好了嘛,bader CHGCAR-1, bader CHGCAR-2, bader CHGCAR.这样前两次之和应该基本等于总电子数

很好奇LZ的金币数目 穷成这样也不容易啊,金币免了吧

[ Last edited by huangyc on 2009-10-25 at 09:46 ]
2楼2009-10-25 09:44:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

★ ★ ★
aylayl08(金币+3,VIP+0):感谢分享你的经验,鼓励新虫 10-25 15:00
aylayl08(金币+0,VIP+0):提示楼主,未求助到书的话可以向所在版的版主反应,可以返还金币,没有打水漂一说:)欢迎讨论 10-25 15:03
引用回帖:
Originally posted by huangyc at 2009-10-25 09:44:
这个-ref我也用不好,哪位解释一下
我觉得你分三次不就好了嘛,bader CHGCAR-1, bader CHGCAR-2, bader CHGCAR.这样前两次之和应该基本等于总电子数

很好奇LZ的金币数目 穷成这样也不容易啊,金币免了吧

[ ...

谢谢你啊!但是用bader CHGCAR-1得到的电子数和bader CHGCAR得到的电子数是一样的。
不过你的回帖给了我一个想法。有可能CHGCAR-1,也就是CHGCAR的上半部分其实是up+down,CHGCAR-2其实是up-down。按照磁性的定义,up-down部分才是系统的磁性。只有这样才能解释我用bader CHGCAR-2 -ref CHGCAR算出来的值比较像磁矩,这样也能解释用bader CHGCAR-1得到的电子数和bader CHGCAR得到的电子数是一样的。
至于-ref,我发现用不用都没有关系。
金币我是从来不在乎的,最近求了一本书,花了100,没求到,打了水漂,呵呵!先给你5个,等高人确认了我的想法我再把剩下的5个给那个人。
3楼2009-10-25 14:02:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice_rain

铁杆木虫 (正式写手)

★ ★ ★ ★ ★
wangyujia(金币+5,VIP+0):谢谢指出错误!! 11-18 08:57
你们算bader电荷么,我算了一个用-ref 和不用时有很大差距的。不知道原因是啥
4楼2009-11-17 11:55:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见