| 查看: 2437 | 回复: 10 | |||
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轨道只有一个电子的事实相违背吗？请大家多多指教，谢谢~ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性 |
» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- d轨道对应的数值
已经有4人回复
- d band center 计算时d电子数如何确定
已经有15人回复
- 两个原子的局域态密度在费米处重叠
已经有4人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】能带结构图中，最高占据能带与最低未占据能带的专业名称是什么呢？
已经有7人回复
- 【讨论】原子/分子在金属表面的吸附过程中外层轨道初始宽化过程
已经有3人回复
- 电子结构和态密度以及能带之间有什么关系?
已经有6人回复
- 【求助】怎么看HOMO LUMO
已经有11人回复
- 【求助】关于HOMO和LUMO的问题
已经有5人回复
- 【重点讨论】Castep中的Empty Band的用处和设置？（参与讨论就有机会赢取大礼包！）
已经有54人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 谢谢指导 2013-11-13 09:21:01
mika: 金币+10, ★★★很有帮助 2013-11-19 17:12:14
感谢参与，应助指数 +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 |
3楼2013-11-12 21:41:59
mika
金虫 (正式写手)
纳米宝贝
- 应助: 6 (幼儿园)
- 金币: 992.4
- 散金: 76
- 帖子: 401
- 在线: 121小时
- 虫号: 836617
- 注册: 2009-08-29
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
2楼2013-11-12 09:42:29
mika
金虫 (正式写手)
纳米宝贝
- 应助: 6 (幼儿园)
- 金币: 992.4
- 散金: 76
- 帖子: 401
- 在线: 121小时
- 虫号: 836617
- 注册: 2009-08-29
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
4楼2013-11-19 17:13:35
mika
金虫 (正式写手)
纳米宝贝
- 应助: 6 (幼儿园)
- 金币: 992.4
- 散金: 76
- 帖子: 401
- 在线: 121小时
- 虫号: 836617
- 注册: 2009-08-29
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
5楼2013-11-19 17:14:22
6楼2013-11-19 18:40:58
7楼2013-11-19 18:43:37
mika
金虫 (正式写手)
纳米宝贝
- 应助: 6 (幼儿园)
- 金币: 992.4
- 散金: 76
- 帖子: 401
- 在线: 121小时
- 虫号: 836617
- 注册: 2009-08-29
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
8楼2013-11-20 13:00:02
9楼2013-11-20 15:07:38
mika
金虫 (正式写手)
纳米宝贝
- 应助: 6 (幼儿园)
- 金币: 992.4
- 散金: 76
- 帖子: 401
- 在线: 121小时
- 虫号: 836617
- 注册: 2009-08-29
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
10楼2013-11-23 15:33:06
