| 查看: 9580 | 回复: 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算出来的电子数，也就是自旋向上的电子数却等于总的电子数？？？ 所以我觉得这种方法可能不对，想请教一下大家，谢谢！ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
行侠仗义
铁虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2590.9
- 散金: 33
- 红花: 2
- 帖子: 752
- 在线: 595.5小时
- 虫号: 1452634
- 注册: 2011-10-20
- 专业: 凝聚态物性 II ：电子结构
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
> 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" 请问这是什么原因？ |
27楼2014-12-07 16:18:11
huangyc
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2242.1
- 散金: 12
- 红花: 2
- 帖子: 411
- 在线: 312.5小时
- 虫号: 624241
- 注册: 2008-10-12
- 专业: 理论和计算化学
穷成这样也不容易啊，金币免了吧2楼2009-10-25 09:44:31
wangyujia
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1704.3
- 散金: 1042
- 红花: 5
- 帖子: 369
- 在线: 79.7小时
- 虫号: 459079
- 注册: 2007-11-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
aylayl08(金币+3,VIP+0):感谢分享你的经验，鼓励新虫 10-25 15:00
aylayl08(金币+0,VIP+0):提示楼主，未求助到书的话可以向所在版的版主反应，可以返还金币，没有打水漂一说：）欢迎讨论 10-25 15:03
aylayl08(金币+3,VIP+0):感谢分享你的经验，鼓励新虫 10-25 15:00
aylayl08(金币+0,VIP+0):提示楼主，未求助到书的话可以向所在版的版主反应，可以返还金币，没有打水漂一说：）欢迎讨论 10-25 15:03
|
谢谢你啊！但是用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
ice_rain
铁杆木虫 (正式写手)
- 应助: 34 (小学生)
- 贵宾: 0.871
- 金币: 5590.5
- 散金: 60
- 红花: 6
- 帖子: 848
- 在线: 963.6小时
- 虫号: 283377
- 注册: 2006-10-08
- 专业: 凝聚态物性 II ：电子结构
4楼2009-11-17 11:55:36
投票:
10