应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】用Bader计算磁矩的困惑
131/2返回列表12下一页
查看: 8350  |  回复: 30
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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算出来的电子数，也就是自旋向上的电子数却等于总的电子数？？？
所以我觉得这种方法可能不对，想请教一下大家，谢谢！
回复此楼

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

VASP 第一性原理相关文档 表界面计算 VASP
ms and vasp学习 第一性 凝聚态物理 siesta & vasp

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2009-10-25 09:25:37
已阅   回复此楼   关注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个给那个人。
赞一下(12人) 回复此楼
3楼2009-10-25 14:02:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
引用回帖:
Originally posted by ice_rain at 2009-11-17 11:55:
你们算bader电荷么，我算了一个用-ref 和不用时有很大差距的。不知道原因是啥

呵呵！你说的对。我以前测试-ref的时候就算了bader CHGCAR-1 -ref CHGCAR，没有算bader CHGCAR-2 -ref CHGCAR。导致了我得出错误的结论。谢谢你给我指出！剩下5个金币都给你了。
赞一下 回复此楼
5楼2009-11-18 08:57:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
★ ★
yjmaxpayne(金币+2,VIP+0):谢谢 11-18 14:56
引用回帖:
Originally posted by ice_rain at 2009-11-17 11:55:
你们算bader电荷么，我算了一个用-ref 和不用时有很大差距的。不知道原因是啥

我现在觉得-ref的作用是先按照总的电子密度对每个原子划分区域，再对磁矩部分的电子密度按照这种划分区域来求每个原子的磁矩。如果不用-ref，直接按照磁矩部分的电子密度来对原子进行划分的话，则会产生误差。原因可能是bader这个划分原子区域的方法对总的电子密度还比较准确，但是原子的磁矩电子密度的分布就不那么规则，不是大体以球状分布在原子核的周围，bader在这种情况下就处理不好了。
个人意见，欢迎大家批评。
赞一下(6人) 回复此楼
6楼2009-11-18 10:51:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

总结

★ ★
wuchenwf(金币+2,VIP+0):谢谢回帖交流 12-3 21:44
总结：
用vtstscripts中的split_chg.sh脚本（split_chg.sh CHGCAR）把CHGCAR分成Total spin部分（file cf1）和Magnetization部分（file cf2），然后用bader cf2 -ref cf1（或CHGCAR），得到的ACF.dat中的charge栏就是每个原子的磁矩。



[ Last edited by wangyujia on 2011-3-17 at 15:59 ]
赞一下(13人) 回复此楼
7楼2009-12-03 15:17:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
引用回帖:
Originally posted by cenwanglai at 2011-03-15 12:27:24:
bader的用处还是不少呢，看来需要做一下总结。谁来呢？

我就是有一段时间对磁矩比较感兴趣，研究了一下bader，目前知道的仍然是这么多，呵呵！
赞一下 回复此楼
10楼2011-03-17 14:09:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
★ ★
zzy870720z(金币+2): 谢谢提供 2011-03-17 20:29:18
以前写总结时没有把脚本传上去，现在把脚本传上去。

用vtstscripts中的split_chg.sh脚本（split_chg.sh CHGCAR）把CHGCAR分成Total spin部分（file cf1）和Magnetization部分（file cf2），然后用bader cf2 -ref cf1（或bader cf2 -ref CHGCAR），得到的ACF.dat文件中的charge一栏就是每个原子的磁矩。
赞一下(7人) 回复此楼
11楼2011-03-17 16:20:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
引用回帖:
Originally posted by ddang100 at 2011-03-17 22:23:21:
为什么不用OUTCAR里面得到的磁矩呢？？用bader得到的磁矩更准确吗？？

OUTCAR里面有吗？我怎么没看到啊！能否告知在哪里找？谢谢！
赞一下 回复此楼
13楼2011-03-21 10:17:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
引用回帖:
Originally posted by uuv2010 at 2011-03-21 13:56:33:
设置LORBIT=11，在OUTCAR的最后面有磁性数据，不过每个原子的磁矩相加与总磁矩往往对应不起来。

谢谢！
回复此楼
15楼2011-03-23 10:13:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
引用回帖:
Originally posted by maliangcai at 2011-05-20 11:42:40:
file cf1给出了原子总的电子分布，我想问一下，怎么才能算出来每个原子分别的up 和 down的电子数。

我想bader cf1可以得到每个原子的Nup+Ndown，bader cf2 -ref cf1可以得到每个原子的Nup-Ndown，从这两个量就可以得到你要的数据。
赞一下(2人) 回复此楼
17楼2011-05-21 10:08:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangyujia 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭