 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6382 | 回復(fù): 15 | ||||||
wbt11129木蟲 (正式寫手)
|
[交流]
【求助】VASP計算出的電極化強(qiáng)度數(shù)據(jù)，怎么做后處理
|
|||||
|
用VASP中的Berry相方法計算出極化數(shù)據(jù)以后，不會做后處理， 比如我計算所得 G1方向P[elc]=(-0.00001, -0.00006, 0.01615), G2方向P[elc]=(0.00000, 0.00000, 0.01615)， G3方向P[elc]=(0.00000, 0.00000, -9.59885), 三個方向的P[ion]=(-61.97390,133.49391,-194.56332)是一致的，單位都是e*Angst。 那么，我所求體系的電極化強(qiáng)度是多少？沿哪個方向？離子貢獻(xiàn)多少？電子貢獻(xiàn)又為多少？ 本人實(shí)在不會后面的處理，請侯博和各位大蝦指教。 注意，我現(xiàn)在計算的材料的prototype structure根本就不知道。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 vasp計算資料 | 第一性原理計算經(jīng)驗 | 第一性原理 |
» 猜你喜歡
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有5人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研求調(diào)劑
已經(jīng)有4人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有4人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- vasp計算PDOS，怎么分辨分波態(tài)屬于哪個具體軌道
已經(jīng)有7人回復(fù)
- 用vasp腳本算極化強(qiáng)度，腳本出錯，急死了
已經(jīng)有6人回復(fù)
- 在方法中直接使用G3關(guān)鍵詞算出的結(jié)果與用G3組合方法計算的結(jié)果有什么區(qū)別
已經(jīng)有8人回復(fù)
- 關(guān)于G09計算的極化率的比較
已經(jīng)有3人回復(fù)
- vasp自旋極化計算時，初始設(shè)置自旋能否不被優(yōu)化
已經(jīng)有6人回復(fù)
- 利用VASP算壓電極化率
已經(jīng)有10人回復(fù)
- 求問，VASP計算born有效電荷后估算極化強(qiáng)度的方法
已經(jīng)有19人回復(fù)
- 壓電陶瓷極化前后，剩余極化強(qiáng)度為什么差別這么大？
已經(jīng)有7人回復(fù)
- VASP Berry phase計算極化強(qiáng)度
已經(jīng)有18人回復(fù)
- vasp計算中的dE是那兩個數(shù)值之差
已經(jīng)有5人回復(fù)
- vasp計算振動強(qiáng)度
已經(jīng)有15人回復(fù)
- 計算出極化強(qiáng)度矢量后，怎么判斷它是沿哪個晶向？（如[111]，[001]）
已經(jīng)有5人回復(fù)
- 【求助】vasp Berry phase計算自發(fā)極化
已經(jīng)有4人回復(fù)
- 【求助】VASP怎么考慮自旋極化？
已經(jīng)有14人回復(fù)
- 【求助】凝膠強(qiáng)度計算數(shù)據(jù)分析問題
已經(jīng)有15人回復(fù)
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★
wuli8(金幣+1):謝謝 2010-10-01 22:57:08
wbt11129(金幣+3):謝謝侯博了，不過還是不會我這種體系的電極化強(qiáng)度計算，因為我這種體系沒有原型結(jié)構(gòu)。 2010-10-20 09:14:04
wuli8(金幣+1):謝謝 2010-10-01 22:57:08
wbt11129(金幣+3):謝謝侯博了，不過還是不會我這種體系的電極化強(qiáng)度計算，因為我這種體系沒有原型結(jié)構(gòu)。 2010-10-20 09:14:04
|
只有這樣的數(shù)據(jù)是得到不到的。 可以看看我blog上的一個例子： http://valenhou.blog.edu.cn/2008/133226.html |
4樓2010-09-30 12:09:30
8樓2013-01-05 10:18:36
wbt11129
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 2378.5
- 散金: 1444
- 帖子: 349
- 在線: 122.9小時
- 蟲號: 834541
- 注冊: 2009-08-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2010-09-28 16:32:26
wbt11129
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 2378.5
- 散金: 1444
- 帖子: 349
- 在線: 122.9小時
- 蟲號: 834541
- 注冊: 2009-08-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
3樓2010-09-30 10:19:14
bnliu
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1032.6
- 散金: 21
- 帖子: 149
- 在線: 30.5小時
- 蟲號: 536536
- 注冊: 2008-03-31
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
結(jié)果處理
★ ★ ★
zzy870720z(金幣+3):謝謝分享 2010-10-14 19:40:43
wbt11129(金幣+3):謝謝分享。你能否把你這個程序使用的輸入文件給我一份，我好參考。我QQ：158604782 2010-10-19 10:35:16
zzy870720z(金幣+3):謝謝分享 2010-10-14 19:40:43
wbt11129(金幣+3):謝謝分享。你能否把你這個程序使用的輸入文件給我一份，我好參考。我QQ：158604782 2010-10-19 10:35:16
|
本人之前計算了鈦酸鋇的自發(fā)極化，最后結(jié)果處理一下就可以得到自發(fā)極化。最好用腳本程序，把最后要處理的相關(guān)數(shù)據(jù)放到一個文件里面，再用FORTRAN對文件數(shù)據(jù)進(jìn)行處理就可以了。 program polar !this program is write for calculated the s_p of the solid after the !scrip "berry" has been done and the file data developed by the scrip !"berry" must be needed implicit real*8 (a-h,o-z) dimension gpar(9,3),pelec(2,3),berry(2,3),pionic(2,3) dimension p(3) character*86 nosen,nosen1 character aa open (9,file='data') open (10,file= 'result') read(9,'(a86)')nosen write(10,'(a86)')nosen do i=1,9 read(9,'(a7,f12.5,a,f12.5,a,f12.5)')aa,gpar(i,1),aa,gpar(i,2),aa,gpar(i,3) ! write(10,110)gpar(i,1),gpar(i,2),gpar(i,3) enddo do j=1,3 pelec(1,j)=(gpar(1,j)+gpar(4,j)+gpar(7,j))/3 berry(1,j)=gpar(2,j)+gpar(5,j)+gpar(8,j) pionic(1,j)=gpar(3,j) enddo ! write(10,*)'the elec, Berry, ionic is' write(10,110)pelec(1,1),pelec(1,2),pelec(1,3) write(10,110)berry(1,1),berry(1,2),berry(1,3) write(10,110)pionic(1,1),pionic(1,2),pionic(1,3) 110 format(f12.5,f12.5,f12.5) read(9,'(a86)')nosen1 write(10,'(a86)')nosen1 do i=1,9 read(9,'(a7,f12.5,a,f12.5,a,f12.5)')aa,gpar(i,1),aa,gpar(i,2),aa,gpar(i,3) ! write(10,110)gpar(i,1),gpar(i,2),gpar(i,3) enddo do j=1,3 pelec(2,j)=(gpar(1,j)+gpar(4,j)+gpar(7,j))/3 berry(2,j)=gpar(2,j)+gpar(5,j)+gpar(8,j) pionic(2,j)=gpar(3,j) enddo write(10,110)pelec(2,1),pelec(2,2),pelec(2,3) write(10,110)berry(2,1),berry(2,2),berry(2,3) write(10,110)pionic(2,1),pionic(2,2),pionic(2,3) read (9,*)nosen read (9,*)nosen read (9,'(a20,f10.2)')nosen,volum read (9,'(a13,f18.10)')nosen,alat read (9,'(a13,f18.10)')nosen,ratio correct= ((1.000/alat)/alat)*1600 !to calculate the polaraztion in x, y, z direction! do k=1,3 p(k)=(pelec(2,k)+berry(2,k)+pionic(2,k)-pelec(1,k)-berry(1,k)-pionic(1,k)) !*1600/volum write(10,*)'the diff is',p(k),'the polar is',(p(k)*(-1600))/volum enddo write(10,*)'The free correction is ', correct end 這是程序，可以看一下。 |
5樓2010-10-14 16:30:55
byin
銅蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 245
- 散金: 332
- 紅花: 3
- 帖子: 395
- 在線: 174.4小時
- 蟲號: 948127
- 注冊: 2010-01-25
- 性別: GG
- 專業(yè): 固體力學(xué)
6樓2011-03-23 08:39:35
future_wl
木蟲 (著名寫手)
- 應(yīng)助: 87 (初中生)
- 金幣: 3295.4
- 散金: 1881
- 紅花: 5
- 帖子: 1404
- 在線: 676.6小時
- 蟲號: 950891
- 注冊: 2010-01-30
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
7樓2013-01-04 23:49:58
9樓2013-01-05 10:19:25
10樓2014-05-28 22:32:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wbt11129 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
|
[考研] 工科，求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|