| 3 | 1/1 | 返回列表 |
| 查看: 1985 | 回復: 2 | |||
[求助]
求教:關于Partical charge density的做法?做過的前輩請進 已有2人參與
|
|
如題,按照這個帖子里的做法: http://www.gaoyang168.com/bbs/viewthread.php?tid=7355162&fpage=1 非常費解的是,為什么要保留自洽的WAVECAR?自洽不就是弛豫嗎?那么弛豫的WAVECAR有意義嗎? 個人感覺應該用非自洽計算,也就是靜態(tài)計算的WAVECAR才有意義…… 請做過的前輩指教一下…… 先謝過了 ![]() 以下為轉帖: http://www.gaoyang168.com/bbs/viewthread.php?tid=7355162&fpage=1 [交流]分享“Partial charge density計算-VASP” 自己初學,受教很多。故抄過來,望能幫助到更多的新人。 文章出處:http://plianghust.blog.163.com/b ... 090822012892478779/ Partial charge density計算-VASP 2012-09-09 14:47:08| 分類: 量子計算 |舉報|字號 訂閱 Partial charge density計算或稱為Band decomposed charge density計算,即計算特定的某個(或某些)k點和本征值(這些k點和本征值是相互對應的)所對應的本征波函數(shù)的平方(也就是電荷密度)。特別是用在 STM的計算中,以及分析特定能量范圍內(nèi)或能量點的化學鍵特征(或atomic characteristic)。 LPARD:可賦予的值為.TRUE.或.FALSE.,它的默認值是.FALSE.,當為.TRUE.時,表示讀入自洽收斂的CHGCAR和WAVECAR并進行Band decomposed charge density計算。 IBAND:它的賦值就是你想要計算的第幾條能帶或哪幾條能帶(比如要計算第4、5、6條能帶,那么就設置IBAND = 4 5 6)。此時NBMOD的值就是所要計算的能帶的條數(shù)(本例中,只計算一條能帶,那么設置為1)。它和EINT不能一起用。 EINT:另一種指定所要計算能帶的方式,它是指定計算某能量范圍的Band decomposed charge density,一般是設置為兩個實數(shù),比如EINT= 4.00 5.00,此時NBOMD的值設置為-2。如果只設置了一個數(shù),那么表示計算從EINT到費米能級這個范圍內(nèi)的,此時NBOMD的值為-3. NBOMD: 當設置0時,表示計算全部的電荷密度(并包括了未占據(jù)的能帶);如果設置-1,就是通常的電荷密度計算,這也是它的默認值。 KPUSE:指定所要計算的k點(哪個或哪幾個),要計算一個,并相應的在KPOINTS設置所要計算的k點。如果要計算好幾個,那么設置KPUSE的值為這些k點在KPOINTS中的序號。 LSEPB:指定是不是要把計算的partial charge density按每個帶分別寫到各自對應的文件PARCHG.nb.(設置為.TRUE.)中,還是把它們合并寫到一個文件中(相當于把各個帶對應的partial charge density加起來,設置為.FALSE.)。默認值為.FALSE.。 LSEPK:指定是不是把要計算的partial charge density按每個k分別寫到各自對應的文件PARCHG.nk.(設置為.TRUE.)中,還是把它們合并寫到一個文件中(相當于把各個帶對應的partial charge density加起來,設置為.FALSE.)。默認值為.FALSE.。 決定畫哪條能帶,或者那些感興趣的K點之后,有如下幾種方法計算不同的Partial Charge。如果你希望計算價帶頂端的Partial Charge,則需要首先通過能帶結構圖確定價帶的能帶標號。需要注意,進行Partial Charge分析必須要保留有自洽計算的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點所對應的Partial Charge。分析導帶、價帶等的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。這種模式適合于分析某個能量區(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。這種模式最利于分析費米面附近的波函數(shù)的性質(zhì)。 |
vasp |

|
本帖內(nèi)容被屏蔽 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +15 | 太容易1018 2026-04-01 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 362求調(diào)劑 +14 | 西南交材料專碩3 2026-03-31 | 14/700 |
|
|
[考研] 413求調(diào)劑 +3 | 柯某某 2026-03-31 | 3/150 |
|
|
[考研] 求調(diào)劑!生物與醫(yī)藥專碩 +4 | 逆轉陸先生 2026-04-01 | 4/200 |
|
|
[考研] 289求調(diào)劑 +23 | 新時代材料 2026-03-27 | 26/1300 |
|
|
[考研]
|
shzhou_ 2026-04-01 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學 ,080500材料科學與工程學碩,總分289分 +11 | @taotao 2026-03-29 | 11/550 |
|
|
[考研] 291求調(diào)劑 +20 | Y-cap 2026-03-29 | 25/1250 |
|
|
[考研] 求調(diào)劑 +4 | DADA怪 2026-03-31 | 4/200 |
|
|
[考研] 一志愿華南師范361分,化學求調(diào)劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +18 | zhubinhao 2026-03-27 | 18/900 |
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
[考研] 0856 335分 +9 | cccchenso 2026-03-29 | 9/450 |
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 本科新能源科學與工程,一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 330一志愿中國海洋大學 化學工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|