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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于原子轨道电子占据band的问题
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 2432  |  回复: 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的回帖

jpchou

专家顾问 (著名写手)
引用回帖:
5楼: Originally posted by mika at 2013-11-19 17:14:22
还想请教一下，在vasp里面有没有可能指定电子在具体轨道上的占据情况呢？谢谢...

LORBIT = 11
產生 PROCAR
這個檔就是你要的
赞一下 回复此楼
6楼2013-11-19 18:40:58
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂：085600材料与化工，考材科基，总分319 +13 678lucky 2026-03-31 14/700 2026-03-31 11:43 by JourneyLucky
[考研] 311求调剂一志愿合肥工业大学 +7 秋二十二 2026-03-30 7/350 2026-03-31 10:51 by Jaylen.
[考研] 274求调剂 +6 xiao爱同学 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 083000环境科学与工程调剂，总分281 +4 橙子（胜意） 2026-03-30 4/200 2026-03-31 00:44 by Linzejun
[考研] 277跪求调剂 +8 1915668 2026-03-27 12/600 2026-03-30 21:01 by dophin1985
[考研] 085601材料工程找调剂 +17 oatmealR 2026-03-29 18/900 2026-03-30 19:21 by Wang200018
[考研] 26考研-291分-厦门大学（085601）-柔性电子学院材料工程专业求调剂 +5 min3 2026-03-24 6/300 2026-03-30 18:42 by 544594351
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 318求调剂 +7 陈晨79 2026-03-30 7/350 2026-03-30 10:49 by 探123
[考研] 311求调剂 +6 冬十三 2026-03-24 6/300 2026-03-29 20:45 by 无际的草原
[考研] 299求调剂 +10 15188958825 2026-03-25 10/500 2026-03-29 17:51 by 王亮_大连医科大
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 数一英一271专硕（085401）求调剂，可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 求调剂推荐 材料 304 +15 荷包蛋hyj 2026-03-26 15/750 2026-03-28 04:13 by fmesaito
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培 +5 YangTyu1 2026-03-26 5/250 2026-03-26 18:27 by 猫咪猫咪呀
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
信息提示
关闭
请填处理意见
关闭