版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

wbt11129

木蟲 (正式寫手)

應助: 3 (幼兒園)
金幣: 2378.5
散金: 1444
帖子: 349
在線: 122.9小時
蟲號: 834541
注冊: 2009-08-26
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

[交流] 【求助】VASP計算出的電極化強度數(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。

那么，我所求體系的電極化強度是多少？沿哪個方向？離子貢獻多少？電子貢獻又為多少？

本人實在不會后面的處理，請侯博和各位大蝦指教。

注意，我現(xiàn)在計算的材料的prototype structure根本就不知道。

回復此樓

» 收錄本帖的淘帖專輯推薦

vasp計算資料

第一性原理計算經(jīng)驗

第一性原理

» 猜你喜歡

26考研求調(diào)劑已經(jīng)有6人回復
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以已經(jīng)有5人回復
085600材料與化工求調(diào)劑已經(jīng)有3人回復
311求調(diào)劑已經(jīng)有7人回復
303求調(diào)劑已經(jīng)有3人回復
有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有7人回復
梁成偉老師課題組歡迎你的加入已經(jīng)有10人回復
一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項已經(jīng)有8人回復
070300化學319求調(diào)劑已經(jīng)有3人回復
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有4人回復

» 本主題相關商家推薦: (我也要在這里推廣)

» 本主題相關價值貼推薦，對您同樣有幫助:

vasp計算PDOS，怎么分辨分波態(tài)屬于哪個具體軌道已經(jīng)有7人回復
用vasp腳本算極化強度，腳本出錯，急死了已經(jīng)有6人回復
在方法中直接使用G3關鍵詞算出的結果與用G3組合方法計算的結果有什么區(qū)別已經(jīng)有8人回復
關于G09計算的極化率的比較已經(jīng)有3人回復
vasp自旋極化計算時，初始設置自旋能否不被優(yōu)化已經(jīng)有6人回復
利用VASP算壓電極化率已經(jīng)有10人回復
求問，VASP計算born有效電荷后估算極化強度的方法已經(jīng)有19人回復
壓電陶瓷極化前后，剩余極化強度為什么差別這么大？已經(jīng)有7人回復
VASP Berry phase計算極化強度已經(jīng)有18人回復
vasp計算中的dE是那兩個數(shù)值之差已經(jīng)有5人回復
vasp計算振動強度已經(jīng)有15人回復
計算出極化強度矢量后，怎么判斷它是沿哪個晶向？（如[111]，[001]）已經(jīng)有5人回復
【求助】vasp Berry phase計算自發(fā)極化已經(jīng)有4人回復
【求助】VASP怎么考慮自旋極化？已經(jīng)有14人回復
【求助】凝膠強度計算數(shù)據(jù)分析問題已經(jīng)有15人回復

BaotianWang

1樓 2010-09-26 09:28:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

valenhou001

至尊木蟲 (職業(yè)作家)

1ST強帖: 13
應助: 241 (大學生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時
蟲號: 1007127
注冊: 2010-04-27
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
wuli8(金幣+1):謝謝 2010-10-01 22:57:08
wbt11129(金幣+3):謝謝侯博了，不過還是不會我這種體系的電極化強度計算，因為我這種體系沒有原型結構。 2010-10-20 09:14:04

只有這樣的數(shù)據(jù)是得到不到的。
可以看看我blog上的一個例子：
http://valenhou.blog.edu.cn/2008/133226.html

贊一下(6人)

回復此樓

4樓2010-09-30 12:09:30

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

闞敏

新蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 178.3
散金: 8
紅花: 1
帖子: 122
在線: 69小時
蟲號: 942507
注冊: 2010-01-14
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

731528樓: Originally posted by valenhou001 at 2010-09-30 12:09:30
只有這樣的數(shù)據(jù)是得到不到的。
可以看看我blog上的一個例子：
http://valenhou.blog.edu.cn/2008/133226.html

侯老師，我就是用你博客里面的方法計算的，但是我發(fā)現(xiàn)只要DIPOL那一項稍微變一下，算出來的極化就差別特別大

贊一下(1人)

回復此樓

8樓2013-01-05 10:18:36

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

wbt11129

木蟲 (正式寫手)

應助: 3 (幼兒園)
金幣: 2378.5
散金: 1444
帖子: 349
在線: 122.9小時
蟲號: 834541
注冊: 2009-08-26
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

有誰算過多鐵材料的，請過來看看

贊一下

回復此樓

BaotianWang

2樓2010-09-28 16:32:26

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wbt11129

木蟲 (正式寫手)

應助: 3 (幼兒園)
金幣: 2378.5
散金: 1444
帖子: 349
在線: 122.9小時
蟲號: 834541
注冊: 2009-08-26
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

頂起來，希望有人關注。

贊一下

回復此樓

BaotianWang

3樓2010-09-30 10:19:14

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

bnliu

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1032.6
散金: 21
帖子: 149
在線: 30.5小時
蟲號: 536536
注冊: 2008-03-31
專業(yè): 凝聚態(tài)物性I:結構、力學和

結果處理

★ ★ ★
zzy870720z(金幣+3):謝謝分享 2010-10-14 19:40:43
wbt11129(金幣+3):謝謝分享。你能否把你這個程序使用的輸入文件給我一份，我好參考。我QQ：158604782 2010-10-19 10:35:16

本人之前計算了鈦酸鋇的自發(fā)極化，最后結果處理一下就可以得到自發(fā)極化。最好用腳本程序，把最后要處理的相關數(shù)據(jù)放到一個文件里面，再用FORTRAN對文件數(shù)據(jù)進行處理就可以了。       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
這是程序，可以看一下。

贊一下(1人)

回復此樓

人生只有一次

5樓2010-10-14 16:30:55

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

byin

銅蟲 (正式寫手)

1ST強帖: 1
應助: 0 (幼兒園)
金幣: 245
散金: 332
紅花: 3
帖子: 395
在線: 174.4小時
蟲號: 948127
注冊: 2010-01-25
性別: GG
專業(yè): 固體力學

引用回帖:

Originally posted by bnliu at 2010-10-14 16:30:55:
本人之前計算了鈦酸鋇的自發(fā)極化，最后結果處理一下就可以得到自發(fā)極化。最好用腳本程序，把最后要處理的相關數(shù)據(jù)放到一個文件里面，再用FORTRAN對文件數(shù)據(jù)進行處理就可以了。 program polar
!this progr ...

求輸入文件格式。。

回復此樓

6樓2011-03-23 08:39:35

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

future_wl

木蟲 (著名寫手)

應助: 87 (初中生)
金幣: 3295.4
散金: 1881
紅花: 5
帖子: 1404
在線: 676.6小時
蟲號: 950891
注冊: 2010-01-30
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

請問侯老師，用Berry phase算出的BEC同 DFPT算出的BEC有何區(qū)別？如果算極化強度的話，兩個值是否都可以用？
謝謝!

贊一下

回復此樓

未來就是現(xiàn)在

7樓2013-01-04 23:49:58

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

闞敏

新蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 178.3
散金: 8
紅花: 1
帖子: 122
在線: 69小時
蟲號: 942507
注冊: 2010-01-14
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

2416727樓: Originally posted by future_wl at 2013-01-04 23:49:58
請問侯老師，用Berry phase算出的BEC同 DFPT算出的BEC有何區(qū)別？如果算極化強度的話，兩個值是否都可以用？
謝謝!...

我是算的壓電張量，搭車同問

贊一下

回復此樓

9樓2013-01-05 10:19:25

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyqh

鐵桿木蟲 (著名寫手)

應助: 33 (小學生)
金幣: 6698
散金: 1454
紅花: 22
帖子: 2375
在線: 593.4小時
蟲號: 2762045
注冊: 2013-10-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

5樓: Originally posted by bnliu at 2010-10-14 16:30:55
本人之前計算了鈦酸鋇的自發(fā)極化，最后結果處理一下就可以得到自發(fā)極化。最好用腳本程序，把最后要處理的相關數(shù)據(jù)放到一個文件里面，再用FORTRAN對文件數(shù)據(jù)進行處理就可以了。 program polar
!this program ...

大蝦，你好，請問你可以把你的處理極化的腳本傳一份給我嗎？謝謝了

贊一下

回復此樓

雨后方能見彩虹

10樓2014-05-28 22:32:13

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 wbt11129 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 化學工程321分求調(diào)劑 +10	大米飯！ 2026-03-15	13/650	2026-03-17 06:57 by laoshidan
[文學芳草園] 伙伴們，祝我生日快樂吧 +17	myrtle 2026-03-10	26/1300	2026-03-16 18:32 by 青橙Ln
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 0703化學調(diào)劑，六級已過，有科研經(jīng)歷 +7	曦熙兮 2026-03-15	7/350	2026-03-16 16:34 by houyaoxu
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 復試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 四川大學085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[考研] 0703化學求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao
[基金申請] 提交后的基金本子，已讓學校撤回了，可否換口子提交 +3	dut_pfx 2026-03-10	3/150	2026-03-11 08:38 by kudofaye
[考研] 279求調(diào)劑 +3	莫xiao 2026-03-10	4/200	2026-03-11 08:06 by 斬魂滴兔子！

24小時熱門版塊排行榜

[交流] 【求助】VASP計算出的電極化強度數(shù)據(jù)，怎么做后處理

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關商家推薦: (我也要在這里推廣)

» 本主題相關價值貼推薦，對您同樣有幫助:

結果處理

[交流] 【求助】VASP計算出的電極化強度數(shù)據(jù)，怎么做后處理