| 24小時(shí)熱門(mén)版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 5983 | 回復(fù): 3 | ||||
[交流]
分享“Partial charge density計(jì)算-VASP” 已有3人參與
|
|
自己初學(xué),受教很多。故抄過(guò)來(lái),望能幫助到更多的新人。 文章出處: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è)帶分別寫(xiě)到各自對(duì)應(yīng)的文件PARCHG.nb.(設(shè)置為.TRUE.)中,還是把它們合并寫(xiě)到一個(gè)文件中(相當(dāng)于把各個(gè)帶對(duì)應(yīng)的partial charge density加起來(lái),設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。 LSEPK:指定是不是把要計(jì)算的partial charge density按每個(gè)k分別寫(xiě)到各自對(duì)應(yīng)的文件PARCHG.nk.(設(shè)置為.TRUE.)中,還是把它們合并寫(xiě)到一個(gè)文件中(相當(dāng)于把各個(gè)帶對(duì)應(yīng)的partial charge density加起來(lái),設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。 決定畫(huà)哪條能帶,或者那些感興趣的K點(diǎn)之后,有如下幾種方法計(jì)算不同的Partial Charge。如果你希望計(jì)算價(jià)帶頂端的Partial Charge,則需要首先通過(guò)能帶結(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ì)。 |
電荷密度分離 |

送紅花一朵
|
新蟲(chóng) (小有名氣)
送紅花一朵 |
本帖內(nèi)容被屏蔽 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +9 | z1z2z3879 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|