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

24小时热门版块排行榜    

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

金虫 (正式写手)

纳米宝贝
引用回帖:
9楼: Originally posted by jpchou at 2013-11-20 15:07:38
指定一个电子占据指定的轨道...
你是想做 excited state calculation ?...

恩，是啊，不知道能不能指导一下~
赞一下 回复此楼
10楼2013-11-23 15:33:06
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭