| 5 | 1/1 | 返回列表 |
| 查看: 2953 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
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ě)手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +7 | 阿斯蒂芬2004 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +6 | Li李魚(yú) 2026-03-26 | 6/300 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書(shū) 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|