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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于原子轨道电子占据band的问题
51/1返回列表
查看: 2407  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

mika

金虫 (正式写手)

纳米宝贝

[求助] 关于原子轨道电子占据band的问题

对于电子占据哪个能带的问题很迷惑，特发帖来请教一下，希望大家多多指教。
我是将一个Ce原子放在一个30*31*32的box里面，K点只用了gamma点，用的potcar是把它的4f电子作为价电子来处理的，所以从potcar里面读出来有个12个价电子，如： POMASS =  140.115; ZVAL   =   12.000    mass and valenz
我的INCAR是这样的：
SYSTEM = ce
ENCUT = 500
ISPIN = 2
MAGMOM = 1
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.05
NELM = 200
NSW = 0
IBRION = -1
ISIF = 2
EDIFF = 1.0E-4
EDIFFG = -0.02
LMAXMIX = 6
AMIX = 0.2
BMIX = 0.0001
LORBIT = 11
===================================
计算出来的OSZICAR最后显示：   1 F= -.11425886E+01 E0= -.10967022E+01  d E =-.917729E-01  mag=     1.9995
这里的mag=1.9995，而Ce的4f轨道只有一个单电子占据，那mag应该约等于1，可这里是确实约等于2。请问是怎么回事呢？
而在outcar结尾部分有这么个输出：

magnetization (x)

# of ion     s       p       d       f       tot
------------------------------------------------
  1        0.005   0.002   0.197   1.455   1.659
而这里的tot magnetization是1.659，跟OSZICAR里面的mag值1.9995不一样，请问到底应该用哪个值呢？
接下来我查看了一下OUTCAR里面的band分布：
spin component 1

k-point   1 :       0.0000    0.0000    0.0000
  band No.  band energies     occupation
      1     -39.1057      1.00000
      2     -22.0951      1.00000
      3     -22.0886      1.00000
      4     -22.0206      1.00000
      5      -3.6462      1.00000
      6      -2.8049      0.36649
      7      -2.7957      0.27353
      8      -2.7938      0.25579
      9      -2.7921      0.24132
     10      -2.7910      0.23128
     11      -2.7908      0.22965
     12      -2.7800      0.14824
     13      -2.7800      0.14786
     14      -2.7728      0.10583
     15      -2.4980      0.00000
     16      -2.4456      0.00000

spin component 2

k-point   1 :       0.0000    0.0000    0.0000
  band No.  band energies     occupation
      1     -38.2853      1.00000
      2     -21.4831      1.00000
      3     -21.2979      1.00000
      4     -21.2696      1.00000
      5      -3.4209      1.00000
      6      -2.3084      0.00000
      7      -2.2845      0.00000
      8      -2.2805      0.00000
      9      -1.9670      0.00000
     10      -1.9379      0.00000
     11      -1.8569      0.00000
     12      -1.8343      0.00000
     13      -1.8235      0.00000
     14      -1.8167      0.00000
     15      -1.7639      0.00000
     16      -1.7638      0.00000

从PROCAR里面分析得band5-11是Ce原子的f轨道，但是从上面的结果看来在f的spin up约有1.7个电子，而spin down有1个，这不是跟f轨道只有一个电子的事实相违背吗？请大家多多指教，谢谢~
回复此楼

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

第一性

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-11 13:29:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mika

金虫 (正式写手)

纳米宝贝
引用回帖:
3楼: Originally posted by jpchou at 2013-11-12 21:41:59
"计算出来的OSZICAR最后显示：   1 F= -.11425886E+01 E0= -.10967022E+01  d E =-.917729E-01  mag=     1.9995
这里的mag=1.9995，而Ce的4f轨道只有一个单电子占据，那mag应该约等于1，可这里是确实约等于2 ...

还想请教一下，在vasp里面有没有可能指定电子在具体轨道上的占据情况呢？谢谢
赞一下 回复此楼
5楼2013-11-19 17:14:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

mika

金虫 (正式写手)

纳米宝贝
自己顶一下，求帮助啊
回复此楼
2楼2013-11-12 09:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 谢谢指导 2013-11-13 09:21:01
mika: 金币+10, ★★★很有帮助 2013-11-19 17:12:14
"计算出来的OSZICAR最后显示：   1 F= -.11425886E+01 E0= -.10967022E+01  d E =-.917729E-01  mag=     1.9995
这里的mag=1.9995，而Ce的4f轨道只有一个单电子占据，那mag应该约等于1，可这里是确实约等于2。"

可能是 trap 在 moment=2 的 local minimum state
可手動給與 mag moment (MAGMOM)
或是鎖定 mag moment (NUPDOWN)
?泶_定 moment=2 是 global minimum or local minimum


"而这里的tot magnetization是1.659，跟OSZICAR里面的mag值1.9995不一样，请问到底应该用哪个值呢？"
用 OSZICAR 的
OUTCAR 內的是透過 wigner radius
畫一個圓(球?)
然後去求圓內的總 moment
所以分?岩欢ú粫竭_ 100%
因此 OUTCAR 內的會小於 OSZICAR
赞一下(1人) 回复此楼
3楼2013-11-12 21:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mika

金虫 (正式写手)

纳米宝贝
引用回帖:
3楼: Originally posted by jpchou at 2013-11-12 21:41:59
"计算出来的OSZICAR最后显示：   1 F= -.11425886E+01 E0= -.10967022E+01  d E =-.917729E-01  mag=     1.9995
这里的mag=1.9995，而Ce的4f轨道只有一个单电子占据，那mag应该约等于1，可这里是确实约等于2 ...

谢谢你，我具体查看了一下那个磁性来源，发现有一个uB是来自于Ce的d轨道。所以有2个。
现在知道为什么OUTCAR里面的值偏小了，谢谢~
赞一下 回复此楼
4楼2013-11-19 17:13:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +7 Liwangman 2026-03-15 7/350 2026-03-18 20:08 by walc
[考研] 0703化学调剂，求各位老师收留 +9 秋有木北 2026-03-14 9/450 2026-03-18 19:49 by macy2011
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 285化工学硕求调剂（081700） +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
关闭
请填处理意见
关闭