| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 5982 | 回復(fù): 3 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||||
[交流]
分享“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| 分類(lèi): 量子計(jì)算 |舉報(bào)|字號(hào) 訂閱 Partial charge density計(jì)算或稱(chēng)為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ì)。 |
電荷密度分離 |

送紅花一朵 |
本帖內(nèi)容被屏蔽 |
新蟲(chóng) (小有名氣)
|
樓主求助,怎么賦予的電荷,譬如硫酸根離子,那個(gè)電荷怎么賦予的呢 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +3 | ioodiiij 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學(xué)。調(diào)劑到任何專(zhuān)業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 求調(diào)劑(材料與化工327) +4 | 愛(ài)吃香菜啦 2026-03-11 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專(zhuān)碩 +3 | YXCT 2026-03-11 | 3/150 |
|