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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于原子轨道电子占据band的问题
51/1返回列表
查看: 2409  |  回复: 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

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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的回帖
查看全部 11 个回答

mika

金虫 (正式写手)

纳米宝贝
自己顶一下，求帮助啊
回复此楼
2楼2013-11-12 09:42:29
已阅   回复此楼   关注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的回帖

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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +4 墨墨漠 2026-03-18 4/200 2026-03-19 21:46 by wang_dand
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 0703化学调剂 ，六级已过，有科研经历 +12 曦熙兮 2026-03-15 12/600 2026-03-19 19:42 by maocaozhuxi
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 12/600 2026-03-19 16:17 by 脱颖而出
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
关闭
请填处理意见
关闭