| 查看: 2060 | 回復(fù): 10 | ||||
jugengfans金蟲 (正式寫手)
|
[交流]
vasp中畫Partial Charge的問(wèn)題!感覺(jué)比較嚴(yán)重 已有6人參與
|
|||
|
partial charge有三種畫法,其中一種是設(shè)置一個(gè)能量范圍,畫出這個(gè)能量范圍內(nèi)的電荷密度。這個(gè)能量是從你設(shè)置的值到費(fèi)米能級(jí)的,如: ISTART = 1 (job : 0-new 1-cont 2-samecut ICHARG = 1 charge: 1-file 2-atom 10-const) LPARD=.TRUE. NBMOD=-3 EINT = -1 LSEPB=.FALSE. LSEPK=.FALSE. 上面的init參數(shù)就表示是畫從費(fèi)米能下-1eV開始,到費(fèi)米能的全部電子態(tài)的電荷密度。 問(wèn)題來(lái)了! 1.因?yàn)槲覀兯鉨and的時(shí)候,K點(diǎn)并沒(méi)有均勻撒在全空間,所以算能帶時(shí)得到的費(fèi)米能級(jí),和自洽計(jì)算的時(shí)候是不同的! 比如費(fèi)米能在自洽計(jì)算的時(shí)候是3eV,在計(jì)算band的時(shí)候就變成了2eV 那么,在這種畫能帶的方法中,程序肯定是按照錯(cuò)誤的費(fèi)米能級(jí)來(lái)給能量范圍的,這樣畫出來(lái)不就沒(méi)意義了嗎。 那這個(gè)方法有什么存在價(jià)值呢。 2.另外,突然覺(jué)得某條能帶所對(duì)應(yīng)的電荷密度也是有問(wèn)題的。因?yàn)樗褪前涯阍贙points里面列的那幾個(gè)點(diǎn)所對(duì)應(yīng)的電荷密度加起來(lái)了! 那如果你的line mode路徑偏重于某一個(gè)方向,重復(fù)走了好幾次,而某些方向走了一次,那這個(gè)時(shí)候總電荷密度就會(huì)偏向體現(xiàn)重復(fù)走的那些路徑! 這樣一來(lái),感覺(jué)只有分析能帶N上的單個(gè)K點(diǎn)所對(duì)應(yīng)的電荷密度才有意義。 也就是說(shuō),LSEPB=.FALSE, LSEPK=.FALSE 讓程序把每個(gè)K點(diǎn)單獨(dú)存放成CHGCAR比較有意義。 大家覺(jué)得呢。 當(dāng)然,如果事先確定一條帶是由某個(gè)原子周圍電子態(tài)產(chǎn)生的,把整個(gè)能帶的合電荷密度畫出來(lái)以證明這點(diǎn),這是可以的。 |
專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
1.?????????band?????K????о??????????????????????????????????????????????????????????????????????????????3eV???????band??????????2eV ?????????????????????У????????????????????????????????Χ??????????????????????????? ?????????????????????? > ??????????????????f?? ???^?@?? NBMOD= -3 ????????????????? band structure ?? partial charge decomposed ?? ?? band structure ?? partial charge ?? NBMOD = -2 ???^??? ???P? NBMOD = -3 ?????????e? ????~ ?????????? STM simulation |
版主 (著名寫手)

金蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
金蟲 (正式寫手)
銅蟲 (初入文壇)
|
那請(qǐng)問(wèn),K點(diǎn)怎樣才能在K空間采樣完全?我最近計(jì)算的占據(jù)數(shù)就出現(xiàn)了誤差,麻煩請(qǐng)樓主給看看,請(qǐng)不吝賜教!http://www.gaoyang168.com/bbs/viewthread.php?tid=10442146&fpage=1 |
木蟲 (小有名氣)
|
算能帶的時(shí)候不是把均勻k點(diǎn)和特殊路徑k點(diǎn)都要加進(jìn)去嗎?而且后者的權(quán)重是0,對(duì)能帶的計(jì)算沒(méi)有影響,只在繪圖的時(shí)候有用 發(fā)自小木蟲Android客戶端 |
金蟲 (正式寫手)
|
是這樣的。是我理解錯(cuò)了。但是由于算能帶的時(shí)候要Ismear=0,對(duì)于一些原本適合Ismear=-5的材料,其展寬會(huì)影響到占據(jù)數(shù),從而影響費(fèi)米面的位置。 發(fā)自小木蟲IOS客戶端 |
金蟲 (正式寫手)
|
謝謝回復(fù)。有人說(shuō)計(jì)算band的時(shí)候即使是line mode,也會(huì)進(jìn)行一個(gè)看不見(jiàn)的均勻撒點(diǎn)的。line mode里面設(shè)的K點(diǎn)不占權(quán)重。所以對(duì)費(fèi)米能的影響僅僅是smearing=0和smearing=-5的區(qū)別 發(fā)自小木蟲IOS客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 314求調(diào)劑 +8 | 無(wú)懈可擊的巨人 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|