| 4 | 1/1 | 返回列表 |
| 查看: 5980 | 回復(fù): 3 | ||||
[交流]
分享“Partial charge density計(jì)算-VASP” 已有3人參與
|
|
自己初學(xué),受教很多。故抄過來,望能幫助到更多的新人。 文章出處:http://plianghust.blog.163.com/b ... 090822012892478779/ Partial charge density計(jì)算-VASP 2012-09-09 14:47:08| 分類: 量子計(jì)算 |舉報(bào)|字號(hào) 訂閱 Partial charge density計(jì)算或稱為Band decomposed charge density計(jì)算,即計(jì)算特定的某個(gè)(或某些)k點(diǎn)和本征值(這些k點(diǎn)和本征值是相互對(duì)應(yīng)的)所對(duì)應(yīng)的本征波函數(shù)的平方(也就是電荷密度)。特別是用在 STM的計(jì)算中,以及分析特定能量范圍內(nèi)或能量點(diǎn)的化學(xué)鍵特征(或atomic characteristic)。 LPARD:可賦予的值為.TRUE.或.FALSE.,它的默認(rèn)值是.FALSE.,當(dāng)為.TRUE.時(shí),表示讀入自洽收斂的CHGCAR和WAVECAR并進(jìn)行Band decomposed charge density計(jì)算。 IBAND:它的賦值就是你想要計(jì)算的第幾條能帶或哪幾條能帶(比如要計(jì)算第4、5、6條能帶,那么就設(shè)置IBAND = 4 5 6)。此時(shí)NBMOD的值就是所要計(jì)算的能帶的條數(shù)(本例中,只計(jì)算一條能帶,那么設(shè)置為1)。它和EINT不能一起用。 EINT:另一種指定所要計(jì)算能帶的方式,它是指定計(jì)算某能量范圍的Band decomposed charge density,一般是設(shè)置為兩個(gè)實(shí)數(shù),比如EINT= 4.00 5.00,此時(shí)NBOMD的值設(shè)置為-2。如果只設(shè)置了一個(gè)數(shù),那么表示計(jì)算從EINT到費(fèi)米能級(jí)這個(gè)范圍內(nèi)的,此時(shí)NBOMD的值為-3. NBOMD: 當(dāng)設(shè)置0時(shí),表示計(jì)算全部的電荷密度(并包括了未占據(jù)的能帶);如果設(shè)置-1,就是通常的電荷密度計(jì)算,這也是它的默認(rèn)值。 KPUSE:指定所要計(jì)算的k點(diǎn)(哪個(gè)或哪幾個(gè)),要計(jì)算一個(gè),并相應(yīng)的在KPOINTS設(shè)置所要計(jì)算的k點(diǎn)。如果要計(jì)算好幾個(gè),那么設(shè)置KPUSE的值為這些k點(diǎn)在KPOINTS中的序號(hào)。 LSEPB:指定是不是要把計(jì)算的partial charge density按每個(gè)帶分別寫到各自對(duì)應(yīng)的文件PARCHG.nb.(設(shè)置為.TRUE.)中,還是把它們合并寫到一個(gè)文件中(相當(dāng)于把各個(gè)帶對(duì)應(yīng)的partial charge density加起來,設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。 LSEPK:指定是不是把要計(jì)算的partial charge density按每個(gè)k分別寫到各自對(duì)應(yīng)的文件PARCHG.nk.(設(shè)置為.TRUE.)中,還是把它們合并寫到一個(gè)文件中(相當(dāng)于把各個(gè)帶對(duì)應(yīng)的partial charge density加起來,設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。 決定畫哪條能帶,或者那些感興趣的K點(diǎn)之后,有如下幾種方法計(jì)算不同的Partial Charge。如果你希望計(jì)算價(jià)帶頂端的Partial Charge,則需要首先通過能帶結(jié)構(gòu)圖確定價(jià)帶的能帶標(biāo)號(hào)。需要注意,進(jìn)行Partial Charge分析必須要保留有自洽計(jì)算的WAVECAR才可以。 第一種Partial Charge分析的INCAR ISTART = 1 job : 0-new 1-cont 2-samecut ICHARG = 1 charge: 1-file 2-atom 10-const LPARD=.TRUE. IBAND= 20 21 22 23 KPUSE= 1 2 3 4 LSEPB=.TRUE. LSEPK=.TRUE. 這樣的INCAR給出的是指定能帶,指定K點(diǎn)所對(duì)應(yīng)的Partial Charge。分析導(dǎo)帶、價(jià)帶等的Partial Charge特性,通常采用的都是這種模式。 第二種Partial Charge分析的INCAR ISTART = 1 job : 0-new 1-cont 2-samecut ICHARG = 1 charge: 1-file 2-atom 10-const LPARD=.TRUE. EINT = -10.3 -5.1 LSEPB=.FALSE. LSEPK=.FALSE. 這樣的INCAR給出的是在[-10.3 -5.1]能量之間的Partial Charge。這種模式適合于分析某個(gè)能量區(qū)間內(nèi)的波函數(shù)的性質(zhì)。 第三種Partial Charge分析的INCAR 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. 這樣的INCAR給出的是從[Ef-1.0 Ef]能量之間的Partial Charge。這種模式最利于分析費(fèi)米面附近的波函數(shù)的性質(zhì)。 |
電荷密度分離 |

送紅花一朵
|
新蟲 (小有名氣)
送紅花一朵 |
本帖內(nèi)容被屏蔽 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|