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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】用Bader计算磁矩的困惑
北京石油化工学院2026年研究生招生接收调剂公告
311/4返回列表1234下一页
查看: 8680  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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

木虫 (正式写手)

总结

★ ★
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 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的回帖

行侠仗义

铁虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
25楼: Originally posted by wangyujia at 2012-04-11 08:46:18
bader是需要编译的。另外应该使用bader cf2 -ref CHGCAR...

> split_chg.sh CHGCAR
提示：

Input file :: "CHGCAR"
(standard_in) 1: illegal character: N
(standard_in) 1: syntax error
(standard_in) 1: syntax error
(standard_in) 1: illegal character: M
Total spin (file "cf1"
sed: -e expression #1, char 4: unexpected `,'
Magneticzatition (file "cf2"


请问这是什么原因？
赞一下(2人) 回复此楼
27楼2014-12-07 16:18:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 谢谢指教 2011-03-21 16:37:42
引用回帖:
Originally posted by wangyujia at 2011-03-21 10:17:14:
OUTCAR里面有吗？我怎么没看到啊！能否告知在哪里找？谢谢！

设置LORBIT=11，在OUTCAR的最后面有磁性数据，不过每个原子的磁矩相加与总磁矩往往对应不起来。
赞一下(3人) 回复此楼

» 本帖已获得的红花（最新10朵）

啊洪江
14楼2011-03-21 13:56:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)
引用回帖:
17613010楼: Originally posted by liuqun at 2012-04-09 17:12:12:
你好，我先用split程序作用CHGCAR，得到cf1和cf2两个文件，然后bader cf1 -ref CHGCAR，就不行了，一直出现

Image              PC                Routine            Line        Source            
bad ...

bader是需要编译的。另外应该使用bader cf2 -ref CHGCAR
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

liuqun
25楼2012-04-11 08:46:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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 ]
赞一下(12人) 回复此楼
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个给那个人。
赞一下(12人) 回复此楼
3楼2009-10-25 14:02:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice_rain

铁杆木虫 (正式写手)
★ ★ ★ ★ ★
wangyujia(金币+5,VIP+0):谢谢指出错误！！ 11-18 08:57
你们算bader电荷么，我算了一个用-ref 和不用时有很大差距的。不知道原因是啥
赞一下(1人) 回复此楼
4楼2009-11-17 11:55:36
已阅   回复此楼   关注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的回帖

cenwanglai

荣誉版主 (知名作家)
8楼2011-03-15 11:35:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bader的用处还是不少呢，看来需要做一下总结。谁来呢？
赞一下(1人) 回复此楼
9楼2011-03-15 12:27:24
已阅   回复此楼   关注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 的主题更新
311/4返回列表1234下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 315分 085602 求调剂 +8 26考研上岸版26 2026-04-02 8/400 2026-04-02 20:00 by chran16
[考研] 材料调剂 +9 一样YWY 2026-04-02 9/450 2026-04-02 18:19 by sprinining
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 311求调剂 +9 勇敢的小吴 2026-04-02 9/450 2026-04-02 11:37 by Sammy2
[考研] 273求调剂 +19 李芷新1 2026-03-31 19/950 2026-04-01 21:49 by chyhaha
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 材料专硕调剂 +17 椰椰。 2026-03-29 17/850 2026-04-01 20:44 by cq2548
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 086000生物与医药298调剂求助 +4 元元青青 2026-03-31 6/300 2026-04-01 11:13 by syh9288
[考研] 318求调剂 +8 七忆77 2026-04-01 8/400 2026-04-01 10:37 by Jaylen.
[考研] 0710生物学求调剂！ +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 085600，专业课化工原理，320分求调剂 +6 大馋小子 2026-03-29 6/300 2026-03-31 10:03 by 氯化亚硝酰
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by 保护地球你我做�
信息提示
关闭
请填处理意见
关闭