| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2948 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
jk射手座新蟲(chóng) (著名寫(xiě)手)
|
[求助]
VASP下用LCALCPOL算極化強(qiáng)度的POSCAR怎么選取
|
||
|
發(fā)此貼是想認(rèn)認(rèn)真真地合適下我計(jì)算極化強(qiáng)度的方法是否理解的正確。無(wú)奈算了下最最簡(jiǎn)單的PbTiO3的極化,但是結(jié)果出來(lái)出來(lái)一下是0 。。。。我想來(lái)想去都是我POSCAR不太對(duì),想請(qǐng)教前輩們,LCALCPOL它是把LBerry 的工作做了哪些?我想核實(shí)一下我理解的vasp英文手冊(cè)上那一節(jié)是不是對(duì)的。說(shuō)明書(shū)里應(yīng)該介紹的是算完一次極化后手動(dòng)移動(dòng)某個(gè)原子再算極化以此來(lái)算BEC。我想請(qǐng)問(wèn)計(jì)算整個(gè)體系的原子要怎么移動(dòng)坐標(biāo)來(lái)選取參考相?(不好意思可能問(wèn)的有些較真,而且我不是第一次在這個(gè)問(wèn)題上犯迷糊,但是必須要弄得清清楚楚了)。是不是算極化的方法都要有兩個(gè)POSCAR ?我貼一下我算PTO時(shí)候的輸入文件: INCAR: SYSTEM=PTO ENCUT=400 PREC=Accurate ISMEAR=-5 SIGMA=0.05 EDIFF=1E-6 ISTART=0 ICHARG=2 IBRION=2 LCALCPOL=.TRUE. DIPOL=0.3 0.3 0.3 K文件是6 6 6的 POSCAR: data_PTO 1.0 3.89 0.00 0.00 0.00 3.89 0.00 0.00 0.00 3.89 1 1 3 Selective dynamics Direct 0.00 0.00 0.00 F F T 0.50 0.50 0.50 F F T 0.50 0.00 0.50 F F T 0.00 0.50 0.50 F F T 0.50 0.50 0.00 F F T 另一個(gè)POSCAR我也就是所有原子的c坐標(biāo)移動(dòng)了0.01的距離,模仿LCALCPOL-tag里的例子。 輸出的OUTCAR里提取 moment 得到: grep 'moment' OUTCAR2(移動(dòng)原子之后算的結(jié)果) Ionic dipole moment: p[ion]=( -7.00200 -7.00200 -8.32460 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 1.32260 ) electrons Angst grep 'moment' OUTCAR1(移動(dòng)原子之前算的結(jié)果) Ionic dipole moment: p[ion]=( -7.00200 -7.00200 -7.00200 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.00000 ) electrons Angst 請(qǐng)問(wèn)怎么算這個(gè)四方相的極化,這兩組偶極矩要怎么處理?勞煩各位前輩同仁了!。!小弟跪求,跪謝,使勁兒跪。。。! |

|
木蟲(chóng) (著名寫(xiě)手)

新蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|