| 5 | 1/1 | 返回列表 |
| 查看: 1395 | 回復: 4 | ||
jk射手座新蟲 (著名寫手)
|
[求助]
求教LCALCPOL計算極化的結果怎么處理的,多出來一行
|
|
我先貼下我的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)格是自動 6 6 4 的 POTCAR用的是PAW的 POSCAR有兩個, 其一: 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 然后計算出的結果我提取的 '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 請問 electrons x Angstroem那行的dipolemoment 和下邊兩行電子與離子的偶極矩是啥關系啊?多設置了一個IDIPOLE=3之后就多了 electrons x Angstroem的那一行; 搜了些相關的說明與討論,如果我理解為: electrons x Angstroem那行dipolemoment是 總P(也就是p[ion]+p[elc])在Z方向上的分量,這樣理解對嗎?我計算了下,分別是3.1027和4.20832,那么其與z軸的夾角為arccos(3.1027/4.20832),大概比45°小一點;我算的是兩個四方相組合的超晶格,總P的方向應該離體對角線不遠,不知這樣理解對不? 請各路大神幫忙指點下。。! (另外還想問下,LBERRY大家都知道是BERRY PHASE方法,那么設置LCALCPOL算極化,這個方法有名字嗎??) [ Last edited by jk射手座 on 2014-7-2 at 10:12 ] |

新蟲 (著名寫手)

新蟲 (著名寫手)

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

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 招08考數(shù)學 +4 | laoshidan 2026-03-20 | 8/400 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
[考研] 308求調劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 291 求調劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 318求調劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 求調劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 工科0856求調劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 材料學學碩080502 337求調劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 346求調劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 301求調劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 296求調劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 321求調劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 考研調劑求學校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調劑) +3 | 123阿標 2026-03-17 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 085601求調劑 +4 | Du.11 2026-03-16 | 4/200 |
|