| 10 | 1/1 | 返回列表 |
| 查看: 2945 | 回復(fù): 9 | ||
jk射手座新蟲 (著名寫手)
|
[求助]
VASP下用LCALCPOL算極化強(qiáng)度的POSCAR怎么選取
|
|
發(fā)此貼是想認(rèn)認(rèn)真真地合適下我計(jì)算極化強(qiáng)度的方法是否理解的正確。無奈算了下最最簡單的PbTiO3的極化,但是結(jié)果出來出來一下是0 。。。。我想來想去都是我POSCAR不太對,想請教前輩們,LCALCPOL它是把LBerry 的工作做了哪些?我想核實(shí)一下我理解的vasp英文手冊上那一節(jié)是不是對的。說明書里應(yīng)該介紹的是算完一次極化后手動移動某個(gè)原子再算極化以此來算BEC。我想請問計(jì)算整個(gè)體系的原子要怎么移動坐標(biāo)來選取參考相?(不好意思可能問的有些較真,而且我不是第一次在這個(gè)問題上犯迷糊,但是必須要弄得清清楚楚了)。是不是算極化的方法都要有兩個(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)移動了0.01的距離,模仿LCALCPOL-tag里的例子。 輸出的OUTCAR里提取 moment 得到: grep 'moment' OUTCAR2(移動原子之后算的結(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(移動原子之前算的結(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 請問怎么算這個(gè)四方相的極化,這兩組偶極矩要怎么處理?勞煩各位前輩同仁了。。!小弟跪求,跪謝,使勁兒跪啊。。! |

木蟲 (著名寫手)

新蟲 (著名寫手)

木蟲 (著名寫手)

新蟲 (著名寫手)

木蟲 (著名寫手)

新蟲 (著名寫手)

版主 (著名寫手)

|
新蟲 (著名寫手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|