版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

mika

金蟲 (正式寫手)

納米寶貝

應(yīng)助: 6 (幼兒園)
金幣: 992.4
散金: 76
帖子: 400
在線: 120.9小時
蟲號: 836617
注冊: 2009-08-29
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 關(guān)于原子軌道電子占據(jù)band的問題

對于電子占據(jù)哪個能帶的問題很迷惑，特發(fā)帖來請教一下，希望大家多多指教。
我是將一個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軌道只有一個單電子占據(jù)，那mag應(yīng)該約等于1，可這里是確實約等于2。請問是怎么回事呢？
而在outcar結(jié)尾部分有這么個輸出：

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不一樣，請問到底應(yīng)該用哪個值呢？
接下來我查看了一下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軌道，但是從上面的結(jié)果看來在f的spin up約有1.7個電子，而spin down有1個，這不是跟f軌道只有一個電子的事實相違背嗎？請大家多多指教，謝謝~

回復(fù)此樓

1樓 2013-11-11 13:29:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

jpchou

專家顧問 (著名寫手)

專家經(jīng)驗: +20
1ST強帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時
蟲號: 585915
注冊: 2008-08-08
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+3, 謝謝指導(dǎo) 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軌道只有一個單電子占據(jù)，那mag應(yīng)該約等于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不一樣，請問到底應(yīng)該用哪個值呢？"
用 OSZICAR 的
OUTCAR 內(nèi)的是透過 wigner radius
畫一個圓(球?)
然後去求圓內(nèi)的總 moment
所以分?巖歡ú粫竭_(dá) 100%
因此 OUTCAR 內(nèi)的會小於 OSZICAR

贊一下(1人)

回復(fù)此樓

3樓2013-11-12 21:41:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

mika

金蟲 (正式寫手)

納米寶貝

應(yīng)助: 6 (幼兒園)
金幣: 992.4
散金: 76
帖子: 400
在線: 120.9小時
蟲號: 836617
注冊: 2009-08-29
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

自己頂一下，求幫助啊

回復(fù)此樓

2樓2013-11-12 09:42:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mika

金蟲 (正式寫手)

納米寶貝

應(yīng)助: 6 (幼兒園)
金幣: 992.4
散金: 76
帖子: 400
在線: 120.9小時
蟲號: 836617
注冊: 2009-08-29
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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軌道只有一個單電子占據(jù)，那mag應(yīng)該約等于1，可這里是確實約等于2 ...

謝謝你，我具體查看了一下那個磁性來源，發(fā)現(xiàn)有一個uB是來自于Ce的d軌道。所以有2個。
現(xiàn)在知道為什么OUTCAR里面的值偏小了，謝謝~

贊一下

回復(fù)此樓

4樓2013-11-19 17:13:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mika

金蟲 (正式寫手)

納米寶貝

應(yīng)助: 6 (幼兒園)
金幣: 992.4
散金: 76
帖子: 400
在線: 120.9小時
蟲號: 836617
注冊: 2009-08-29
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

還想請教一下，在vasp里面有沒有可能指定電子在具體軌道上的占據(jù)情況呢？謝謝

贊一下

回復(fù)此樓

5樓2013-11-19 17:14:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專家顧問 (著名寫手)

專家經(jīng)驗: +20
1ST強帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時
蟲號: 585915
注冊: 2008-08-08
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

引用回帖:

5樓: Originally posted by mika at 2013-11-19 17:14:22
還想請教一下，在vasp里面有沒有可能指定電子在具體軌道上的占據(jù)情況呢？謝謝...

LORBIT = 11
產(chǎn)生 PROCAR
這個檔就是你要的

贊一下

回復(fù)此樓

6樓2013-11-19 18:40:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專家顧問 (著名寫手)

專家經(jīng)驗: +20
1ST強帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時
蟲號: 585915
注冊: 2008-08-08
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

引用回帖:

5樓: Originally posted by mika at 2013-11-19 17:14:22
還想請教一下，在vasp里面有沒有可能指定電子在具體軌道上的占據(jù)情況呢？謝謝...

不對~等等
PROCAR 是 "顯示電子在軌道上的占據(jù)情形"
你是要 "顯示" 還是要 "指定" ?

贊一下

回復(fù)此樓

7樓2013-11-19 18:43:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mika

金蟲 (正式寫手)

納米寶貝

應(yīng)助: 6 (幼兒園)
金幣: 992.4
散金: 76
帖子: 400
在線: 120.9小時
蟲號: 836617
注冊: 2009-08-29
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

7樓: Originally posted by jpchou at 2013-11-19 18:43:37
不對~等等
PROCAR 是 "顯示電子在軌道上的占據(jù)情形"
你是要 "顯示" 還是要 "指定" ?...

要指定，因為從PROCAR里面可以分析出具體哪個原子的哪條軌道占據(jù)了哪條band，現(xiàn)在我想指定一個電子占據(jù)指定的軌道

贊一下

回復(fù)此樓

8樓2013-11-20 13:00:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專家顧問 (著名寫手)

專家經(jīng)驗: +20
1ST強帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時
蟲號: 585915
注冊: 2008-08-08
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

引用回帖:

8樓: Originally posted by mika at 2013-11-20 13:00:02
要指定，因為從PROCAR里面可以分析出具體哪個原子的哪條軌道占據(jù)了哪條band，現(xiàn)在我想指定一個電子占據(jù)指定的軌道...

指定一個電子占據(jù)指定的軌道...
你是想做 excited state calculation ?

贊一下

回復(fù)此樓

9樓2013-11-20 15:07:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mika

金蟲 (正式寫手)

納米寶貝

應(yīng)助: 6 (幼兒園)
金幣: 992.4
散金: 76
帖子: 400
在線: 120.9小時
蟲號: 836617
注冊: 2009-08-29
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

9樓: Originally posted by jpchou at 2013-11-20 15:07:38
指定一個電子占據(jù)指定的軌道...
你是想做 excited state calculation ?...

恩，是啊，不知道能不能指導(dǎo)一下~

贊一下

回復(fù)此樓

10樓2013-11-23 15:33:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 mika 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3	脫穎而出 2026-03-16	6/300	2026-03-17 21:19 by z1z2z3879
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 279求調(diào)劑 +3	Dizzy123@ 2026-03-10	3/150	2026-03-13 23:02 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6	步川酷紫123 2026-03-13	6/300	2026-03-13 21:59 by 星空星月
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 求調(diào)劑資源與環(huán)境 285 +3	未名考生 2026-03-10	3/150	2026-03-13 10:31 by houyaoxu

24小時熱門版塊排行榜

[求助] 關(guān)于原子軌道電子占據(jù)band的問題

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: