| 10 | 1/1 | 返回列表 |
| 查看: 2941 | 回復(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)該介紹的是算完一次極化后手動移動某個原子再算極化以此來算BEC。我想請問計(jì)算整個體系的原子要怎么移動坐標(biāo)來選取參考相?(不好意思可能問的有些較真,而且我不是第一次在這個問題上犯迷糊,但是必須要弄得清清楚楚了)。是不是算極化的方法都要有兩個POSCAR ?我貼一下我算PTO時候的輸入文件: 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 另一個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 請問怎么算這個四方相的極化,這兩組偶極矩要怎么處理?勞煩各位前輩同仁了!。!小弟跪求,跪謝,使勁兒跪啊。。! |

木蟲 (著名寫手)

新蟲 (著名寫手)

木蟲 (著名寫手)

新蟲 (著名寫手)

木蟲 (著名寫手)

新蟲 (著名寫手)

版主 (著名寫手)

|
新蟲 (著名寫手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +4 | adaie 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-19 | 8/400 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +6 | Doleres 2026-03-19 | 6/300 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么! 2026-03-15 | 19/950 |
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|