24小时热门版块排行榜    

查看: 6397  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

nkleof

银虫 (正式写手)

[交流] 【求助】表面静电势怎么计算?

看到有的文献中做出了从表面深入到体材料的各层的静电势,如下图所示:

请问这个静电势数据从vasp计算的结果中如何提取?
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理和电化学 第一性原理 不错的帖子

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
nkleof(金币+1):谢谢参与
cenwanglai(金币+2):谢谢回复! 2010-09-30 09:18:36
nkleof(金币+7):非常感谢! 2010-09-30 10:24:19
网上可以找到相关的工具,vasp组弄的workshop的材料里面好像就有类似的工具,台湾学术圈办的workshop的材料里面也有处理LOCPOT的工具。自己写一个也很简单。

读入LOCPOT,然后写一个求和的子过程就可以得到了。

pot (i,jk) 这是local potential的数组,它的维数在LOPOT文件里面有给出,也一般是NGX, NGY, NGZ。假设表面或界面的法向方向是在z方向上,那么
do k =1, NGZ
    summ = 0.0
    do j= 1, NGY
       do i = 1, NGX
           summ = pot(i,j,k) +summ
        end do
    end do
     summ =  表面或界面的面积 * summ/ (NGX * NGY)
     write(*,*)   k*原胞c轴的长度,   summ
  end do
要注意的就是各个量的单位。
4楼2010-09-30 08:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

liangab_234620

木虫 (小有名气)

★ ★
nkleof(金币+1):谢谢参与
cenwanglai(金币+1):你是说那个计算功函数的离子? 2010-09-29 11:36:03
nkleof(金币+4):谢谢。 2010-09-29 17:07:04
VASP老家的例子文件里面有
2楼2010-09-29 10:40:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)

引用回帖:
Originally posted by liangab_234620 at 2010-09-29 10:40:53:
VASP老家的例子文件里面有

谢谢回复,看起来好像是在INCAR里面设置LVTOT=TRUE,然后在LOCPOT文件里面能够得到the planar-averaged potential 。
不过还是不太明白如何处理,能不能提供详细一点的计算并处理相关数据的资料?
3楼2010-09-29 17:06:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangpengru

铜虫 (小有名气)

★ ★
nkleof(金币+1):谢谢参与
cenwanglai(金币+1):谢谢回复! 2010-10-10 09:03:21
看帖http://muchong.com/bbs/viewthread.php?tid=2357737&fpage=1
复制过来的时候记住:
WRITE(*,*) 'Error opening or reading file LOCPOT.'
之前少了个空格
5楼2010-10-09 09:31:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
信息提示
请填处理意见