| 查看: 2430 | 回復(fù): 16 | |||
mazuju028木蟲 (正式寫手)
|
[求助]
GW計(jì)算流程問題
|
|
GW計(jì)算在vasp網(wǎng)站上有簡單說明,還是有些不明確的地方想和大家交流一下。 先簡單羅列一下步驟, 1,較少空帶的自洽計(jì)算。 System = Si groundstate occupied orbitals ISMEAR = 0 ; SIGMA = 0.05 ! small sigma is required to avoid partial occupancies EDIFF = 1E-8 ! required tight tolerance for groundstate orbitals 2,較多能帶的自洽計(jì)算。 System = Si unoccupied orbitals ALGO = Exact ! use exact diagonalization of the Hamiltonian NELM = 1 ! since we are already converged stop after one step NBANDS = 96 ISMEAR = 0 ; SIGMA = 0.05 ! small sigma is required to avoid partial occupancies LOPTICS = .TRUE. 3,G0W0計(jì)算。 System = Si NBANDS = 96 ISMEAR = 0 ; SIGMA = 0.05 LOPTICS = .TRUE. ALGO = GW0 ; NOMEGA = 50 ---疑問如下--------------------------------------------------------------------------- 1,這三步都應(yīng)該是在同一目錄下進(jìn)行,所以第二步應(yīng)該是默認(rèn)讀取第一步的波函數(shù),第三步默認(rèn)讀取第二步的波函數(shù)? 2,有關(guān)K點(diǎn)設(shè)置,是否可以從第一步自洽計(jì)算開始就用特殊K點(diǎn)? 3,第二步中NELM=1如果不設(shè)置,對(duì)最后GW計(jì)算的結(jié)果會(huì)有多大影響? |

木蟲 (正式寫手)
|
NELM = 1 ! since we are already converged stop after one step 看注釋,是因?yàn)閟tep1進(jìn)行了自洽運(yùn)算,故這里設(shè)置NELM=1是為了減少計(jì)算量 2,有關(guān)K點(diǎn)設(shè)置,是否可以從第一步自洽計(jì)算開始就用特殊K點(diǎn)? 樓主是為了算能帶嗎? 我曾經(jīng)試過,特殊的k點(diǎn)一開始可以使用的,但是最后GW的計(jì)算仍然得不到能帶圖。 |
榮譽(yù)版主 (著名寫手)
木蟲 (正式寫手)

木蟲 (正式寫手)

木蟲 (正式寫手)
木蟲 (正式寫手)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +83 |
榮譽(yù)版主 (著名寫手)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +83 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|