| 5 | 1/1 | 返回列表 |
| 查看: 1396 | 回復(fù): 4 | ||
jk射手座新蟲 (著名寫手)
|
[求助]
求教LCALCPOL計(jì)算極化的結(jié)果怎么處理的,多出來一行
|
|
我先貼下我的INCAR: SYSTEM=PTOLAO ENCUT=450 ISMEAR=-5 SIGMA=0.02 EDIFF=1E-3 ISTART=0 ICHARG=2 IBRION=2 PREC=Accurate LCALCPOL=.TRUE. DIPOL=0.3 0.3 0.3 IDIPOL=3 LDAU=.TRUE. LDAUTYPE=2 LDAUL=-1 -1 3 -1 -1 LDAUU=0 0 6 0 0 LDAUJ=0 0 0 0 0 LDAUPRINT=2 K網(wǎng)格是自動(dòng) 6 6 4 的 POTCAR用的是PAW的 POSCAR有兩個(gè), 其一: data_PTOLAO 1.00000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 8.1400000000000000 Pb Ti La Al O 1 1 1 1 6 Selective dynamics Direct 0.00 0.00 0.00 F F T 0.50 0.50 0.25 F F T 0.00 0.00 0.50 F F T 0.50 0.50 0.75 F F T 0.50 0.00 0.25 F F T 0.00 0.50 0.25 F F T 0.50 0.50 0.00 F F T 0.50 0.00 0.75 F F T 0.00 0.50 0.75 F F T 0.50 0.50 0.50 F F T 其二: data_PTOLAO 1.00000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 7.9772000000000000 Pb Ti La Al O 1 1 1 1 6 Selective dynamics Direct 0.0000000000000000 0.0000000000000000 -0.0444396041984977 F F T 0.5000000000000000 0.5000000000000000 0.2327315837951796 F F T 0.0000000000000000 0.0000000000000000 0.4971976926301170 F F T 0.5000000000000000 0.5000000000000000 0.7423875204210301 F F T 0.5000000000000000 0.0000000000000000 0.2778793310366013 F F T 0.0000000000000000 0.5000000000000000 0.2778793310366013 F F T 0.5000000000000000 0.5000000000000000 0.0092994307667890 F F T 0.5000000000000000 0.0000000000000000 0.7482258333968913 F F T 0.0000000000000000 0.5000000000000000 0.7482258333968913 F F T 0.5000000000000000 0.5000000000000000 0.5106130477183970 F F T 然后計(jì)算出的結(jié)果我提取的 'moment'.提取出來是醬嬸兒的: angular momentum for each species LDAUL = -1 -1 3 -1 -1 dipolmoment 0.000000 0.000000 -21.359368 electrons x Angstroem Ionic dipole moment: p[ion]=( 12.57300 12.57300 -20.04233 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.28231 ) electrons Angst [root@node16 -1.5calc]# grep 'moment' OUTCAR1 angular momentum for each species LDAUL = -1 -1 3 -1 -1 dipolmoment 0.000000 0.000000 -24.451270 electrons x Angstroem Ionic dipole moment: p[ion]=( 12.57300 12.57300 -24.01300 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.15750 ) electrons Angst [root@node16 -1.5calc]# grep 'volume' OUTCAR2 volume/ion in A,a.u. = 11.64 78.54 volume of cell : 116.39 volume of cell : 116.39 請(qǐng)問 electrons x Angstroem那行的dipolemoment 和下邊兩行電子與離子的偶極矩是啥關(guān)系?多設(shè)置了一個(gè)IDIPOLE=3之后就多了 electrons x Angstroem的那一行; 搜了些相關(guān)的說明與討論,如果我理解為: electrons x Angstroem那行dipolemoment是 總P(也就是p[ion]+p[elc])在Z方向上的分量,這樣理解對(duì)嗎?我計(jì)算了下,分別是3.1027和4.20832,那么其與z軸的夾角為arccos(3.1027/4.20832),大概比45°小一點(diǎn);我算的是兩個(gè)四方相組合的超晶格,總P的方向應(yīng)該離體對(duì)角線不遠(yuǎn),不知這樣理解對(duì)不? 請(qǐng)各路大神幫忙指點(diǎn)下。。! (另外還想問下,LBERRY大家都知道是BERRY PHASE方法,那么設(shè)置LCALCPOL算極化,這個(gè)方法有名字嗎??) [ Last edited by jk射手座 on 2014-7-2 at 10:12 ] |

新蟲 (著名寫手)

新蟲 (著名寫手)

金蟲 (正式寫手)
新蟲 (著名寫手)

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +5 | 苑豆豆 2026-03-20 | 5/250 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無言 2026-03-21 | 6/300 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|