| 5 | 1/1 | 返回列表 |
| 查看: 2063 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
jugengfans金蟲 (正式寫手)
|
[交流]
vasp中畫Partial Charge的問題!感覺比較嚴(yán)重 已有6人參與
|
||
|
partial charge有三種畫法,其中一種是設(shè)置一個能量范圍,畫出這個能量范圍內(nèi)的電荷密度。這個能量是從你設(shè)置的值到費米能級的,如: 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ù)就表示是畫從費米能下-1eV開始,到費米能的全部電子態(tài)的電荷密度。 問題來了! 1.因為我們算band的時候,K點并沒有均勻撒在全空間,所以算能帶時得到的費米能級,和自洽計算的時候是不同的! 比如費米能在自洽計算的時候是3eV,在計算band的時候就變成了2eV 那么,在這種畫能帶的方法中,程序肯定是按照錯誤的費米能級來給能量范圍的,這樣畫出來不就沒意義了嗎。 那這個方法有什么存在價值呢。 2.另外,突然覺得某條能帶所對應(yīng)的電荷密度也是有問題的。因為它就是把你在Kpoints里面列的那幾個點所對應(yīng)的電荷密度加起來了! 那如果你的line mode路徑偏重于某一個方向,重復(fù)走了好幾次,而某些方向走了一次,那這個時候總電荷密度就會偏向體現(xiàn)重復(fù)走的那些路徑! 這樣一來,感覺只有分析能帶N上的單個K點所對應(yīng)的電荷密度才有意義。 也就是說,LSEPB=.FALSE, LSEPK=.FALSE 讓程序把每個K點單獨存放成CHGCAR比較有意義。 大家覺得呢。 當(dāng)然,如果事先確定一條帶是由某個原子周圍電子態(tài)產(chǎn)生的,把整個能帶的合電荷密度畫出來以證明這點,這是可以的。 |
金蟲 (正式寫手)
|
謝謝回復(fù)。有人說計算band的時候即使是line mode,也會進(jìn)行一個看不見的均勻撒點的。line mode里面設(shè)的K點不占權(quán)重。所以對費米能的影響僅僅是smearing=0和smearing=-5的區(qū)別 發(fā)自小木蟲IOS客戶端 |
版主 (著名寫手)

金蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 314求調(diào)劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|