 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3643 | 回復(fù): 10 | ||||||
utimasa新蟲 (初入文壇)
|
[求助]
求助VASP 處理帶電缺陷的問題！能量變化太大
|
|||||
|
本人正在算一個半導(dǎo)體氧化物中引入N雜質(zhì)，N替換O，計算其形成能。 首先計算了中性缺陷，體系的總能： ETOT= -648.481449 eV 然后計算N（-1）,體系的總能，將NELECT設(shè)為NELECT(原來的中性體系）+1 得到ETOT=-643.658924 eV 竟然相差4.82eV。 也就是說如果N得到一個電子，體系的能量竟然增加了4.82eV. 其他參數(shù)都是一樣，沒有變化。 根據(jù)實(shí)驗(yàn)，N是很好的受主，應(yīng)該是容易得到電子。 百思不得其解。 望有高手指點(diǎn)一二！ 不勝感激！ 以下是輸入文件INCAR： NWRITE = 2 # verbosity write-tag(how much is written) PREC = Normal ISTART = 0 # startjob 0=new 1=cont 2=samecut ISPIN = 2 # spin polarized calculation (2=yes,1=no) INIWAV = 1 # initial electr wf 0=lowe 1=rand ENCUT = 450.0 eV EDIFF = 1.0E-06 NSW = 0 # number of steps for ionic upd IBRION = -1 # ionic relaxation 0=MD 1=quasi-New 2=CG ISIF = 2 ISMEAR = 0 # partial occupancies -5=Blochl(energy) -4=tet -1=fermi 0=gaussian >0 MP(force) SIGMA = 0.05 LORBIT = 11 # 0 simple, 1 ext, 2 COOP (PROOUT) # NELECT = 702.0000 #plus U LDAU=.TRUE. LDAUTYPE=2 LDAUL= -1 -1 -1 2 LDAUU=0.00 0.00 0.00 5.2 LDAUPRINT=2 LMAXMIX=4 [ Last edited by utimasa on 2012-2-12 at 21:31 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一原理資料匯編 | 帶電摻雜體系形成能 | first principle | 雜七雜八 |
» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 表面帶電缺陷體系VASP計算參數(shù)設(shè)置和模型搭建
已經(jīng)有8人回復(fù)
- 【求助】VASP收斂問題
已經(jīng)有15人回復(fù)
- 【求助】vasp 算出的0k下的能量如何擴(kuò)展到其他溫度壓力
已經(jīng)有12人回復(fù)
- 【求助】vasp靜態(tài)計算能量值錯誤
已經(jīng)有21人回復(fù)
- 【求助】新手求助vasp計算氧分子結(jié)合能的問題
已經(jīng)有12人回復(fù)
- 【求助】vasp中關(guān)于電場單位的問題
已經(jīng)有12人回復(fù)
acridine
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 11
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 4901.6
- 紅花: 13
- 帖子: 706
- 在線: 389.5小時
- 蟲號: 305478
- 注冊: 2006-12-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
賀儀(金幣+2, 專家考核): 多謝指教 2012-02-15 22:22:17
utimasa(金幣+8): ★★★★★最佳答案 2012-02-16 20:03:46
賀儀(金幣+2, 專家考核): 多謝指教 2012-02-15 22:22:17
utimasa(金幣+8): ★★★★★最佳答案 2012-02-16 20:03:46
|
兩個體系電子數(shù)目不等，直接比較總能量是沒有意義的。 如果真的要比帶電前后的能量差別，那也應(yīng)該是比較：（"中性體系+一個額外電子"的總能量)與（"帶一個電子體系"的總能量) 后一項(xiàng)自然是你計算得到的Etot 前一項(xiàng)還需要加上一個電子的能量，但這個能量的定義方式需要結(jié)合具體情況。比如你這里，如果認(rèn)為這個電子是從半導(dǎo)體氧化物中轉(zhuǎn)移來的，那么這個電子的能量就是半導(dǎo)體氧化物的費(fèi)米能級的能量。不過實(shí)際計算的時候不能直接把bulk計算得到的Ef直接拿來用，因?yàn)檫�要考慮不同計算時候的band align。 關(guān)于缺陷形成能的計算和討論可以看看這篇綜述，里面介紹的比較詳細(xì) J. Appl. Phys. 95, 3851 (2004) |
6樓2012-02-15 17:46:46
uuv2010
榮譽(yù)版主 (職業(yè)作家)
- 1ST強(qiáng)帖: 6
- 應(yīng)助: 160 (高中生)
- 貴賓: 3.14
- 金幣: 3681.7
- 散金: 5802
- 紅花: 79
- 沙發(fā): 35
- 帖子: 4560
- 在線: 916小時
- 蟲號: 966216
- 注冊: 2010-03-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
2樓2012-02-13 15:01:39
utimasa
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3
- 帖子: 10
- 在線: 16.4小時
- 蟲號: 1243869
- 注冊: 2011-03-24
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2012-02-14 09:10:03
buct2010
金蟲 (正式寫手)
- 1ST強(qiáng)帖: 10
- 應(yīng)助: 27 (小學(xué)生)
- 金幣: 2328.6
- 散金: 306
- 紅花: 11
- 帖子: 351
- 在線: 231.2小時
- 蟲號: 1129027
- 注冊: 2010-10-22
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
4樓2012-02-14 14:36:15
WDD880227
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 25 (小學(xué)生)
- 貴賓: 0.098
- 金幣: 1348.5
- 散金: 194
- 紅花: 16
- 帖子: 761
- 在線: 371.6小時
- 蟲號: 1129580
- 注冊: 2010-10-22
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2012-02-14 19:34:19
WDD880227
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 25 (小學(xué)生)
- 貴賓: 0.098
- 金幣: 1348.5
- 散金: 194
- 紅花: 16
- 帖子: 761
- 在線: 371.6小時
- 蟲號: 1129580
- 注冊: 2010-10-22
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
7樓2012-02-16 08:46:07
utimasa
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3
- 帖子: 10
- 在線: 16.4小時
- 蟲號: 1243869
- 注冊: 2011-03-24
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
8樓2012-02-16 19:47:03
9樓2012-02-19 19:24:28
linghu.jj
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 298
- 帖子: 6
- 在線: 4.3小時
- 蟲號: 1346336
- 注冊: 2011-07-15
- 專業(yè): 材料物理化學(xué)
10樓2012-06-26 11:43:57
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 utimasa 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿中國石油大學(xué)（華東） 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|