| 24小時(shí)熱門(mén)版塊排行榜 |
| 3 | 1/1 | 返回列表 |
| 查看: 1709 | 回復(fù): 2 | ||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)2次,作者maomao7910增加金幣 1.6 個(gè) | ||
| 本帖產(chǎn)生 1 個(gè) 1ST強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | ||
maomao7910金蟲(chóng) (小有名氣)
|
[資源]
自發(fā)極化設(shè)置LCALCPOL時(shí)存在的問(wèn)題
|
|
|
如題,用LBERRY = .TRUE. 的方法候博已經(jīng)給出詳細(xì)例子, http://valenhou.blog.edu.cn/2008/133226.html ,我不再贅述。vasp5.2后加入了LCALCPOL=.TRUE.開(kāi)關(guān),可以一次性算出電子的極化貢獻(xiàn)。不過(guò)目前這樣設(shè)置后在OUTCAR中并沒(méi)有出來(lái)如手冊(cè)上所講的 Ionic dipole moment: p[ion]=( 2.25510 2.25510 2.25510 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.00000 ) electrons Angst 兩項(xiàng),而只是給出 CALCP: cpu time 334.26: real time 341.92 Total macroscopic dipole moment: p[mac]=( 0.00000 0.00000 0.19288 ) electrons Angst 這一項(xiàng),是電子極化的貢獻(xiàn),那么離子極化的貢獻(xiàn)如何解決?給兩個(gè)辦法,如上面候博講的求原子移動(dòng)和未移動(dòng)時(shí)的berry phase,但這里只求一次,即 IGPAR = 1 or2 or3就可以了。另一種方法。如果你看關(guān)于berry phase 的手冊(cè) http://cms.mpi.univie.ac.at/vasp/Berry_phase/node4.html 你會(huì)發(fā)現(xiàn)離子極化項(xiàng)其實(shí)就是 (價(jià)電荷*原子位移)所以,自己手動(dòng)一算就出來(lái)了。先在勢(shì)函數(shù)里面找到你要移動(dòng)的原子的價(jià)電子數(shù) ZVAL*(原子移動(dòng)前的坐標(biāo)-原子移動(dòng)后的坐標(biāo))就等于離子極化的貢獻(xiàn)。 歡迎指正! [ Last edited by youzhizhe on 2011-7-6 at 12:59 ] |
銀蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫(xiě)手)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|