| 查看: 6930 | 回复: 15 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
wbt11129木虫 (正式写手)
|
[交流]
【求助】VASP计算出的电极化强度数据，怎么做后处理
|
|||||
|
用VASP中的Berry相方法计算出极化数据以后，不会做后处理， 比如我计算所得 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。 那么，我所求体系的电极化强度是多少？沿哪个方向？离子贡献多少？电子贡献又为多少？ 本人实在不会后面的处理，请侯博和各位大虾指教。 注意，我现在计算的材料的prototype structure根本就不知道。 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp计算资料 | 第一性原理计算经验 | 第一性原理 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有85人回复
-
感觉是下周放榜了
已经有9人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算PDOS，怎么分辨分波态属于哪个具体轨道
已经有7人回复
- 用vasp脚本算极化强度，脚本出错，急死了
已经有6人回复
- 在方法中直接使用G3关键词算出的结果与用G3组合方法计算的结果有什么区别
已经有8人回复
- 关于G09计算的极化率的比较
已经有3人回复
- vasp自旋极化计算时，初始设置自旋能否不被优化
已经有6人回复
- 利用VASP算压电极化率
已经有10人回复
- 求问，VASP计算born有效电荷后估算极化强度的方法
已经有19人回复
- 压电陶瓷极化前后，剩余极化强度为什么差别这么大？
已经有7人回复
- VASP Berry phase计算极化强度
已经有18人回复
- vasp计算中的dE是那两个数值之差
已经有5人回复
- vasp计算振动强度
已经有15人回复
- 计算出极化强度矢量后，怎么判断它是沿哪个晶向？（如[111]，[001]）
已经有5人回复
- 【求助】vasp Berry phase计算自发极化
已经有4人回复
- 【求助】VASP怎么考虑自旋极化？
已经有14人回复
- 【求助】凝胶强度计算数据分析问题
已经有15人回复
wbt11129
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2378.5
- 散金: 1444
- 帖子: 349
- 在线: 122.9小时
- 虫号: 834541
- 注册: 2009-08-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2010-09-30 10:19:14
wbt11129
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2378.5
- 散金: 1444
- 帖子: 349
- 在线: 122.9小时
- 虫号: 834541
- 注册: 2009-08-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2010-09-28 16:32:26
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★
wuli8(金币+1):谢谢 2010-10-01 22:57:08
wbt11129(金币+3):谢谢侯博了，不过还是不会我这种体系的电极化强度计算，因为我这种体系没有原型结构。 2010-10-20 09:14:04
wuli8(金币+1):谢谢 2010-10-01 22:57:08
wbt11129(金币+3):谢谢侯博了，不过还是不会我这种体系的电极化强度计算，因为我这种体系没有原型结构。 2010-10-20 09:14:04
|
只有这样的数据是得到不到的。 可以看看我blog上的一个例子： http://valenhou.blog.edu.cn/2008/133226.html |
4楼2010-09-30 12:09:30
bnliu
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1032.6
- 散金: 21
- 帖子: 149
- 在线: 30.5小时
- 虫号: 536536
- 注册: 2008-03-31
- 专业: 凝聚态物性I:结构、力学和
结果处理
★ ★ ★
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
|
本人之前计算了钛酸钡的自发极化，最后结果处理一下就可以得到自发极化。最好用脚本程序，把最后要处理的相关数据放到一个文件里面，再用FORTRAN对文件数据进行处理就可以了。 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
