24小时热门版块排行榜    

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

nkleof

银虫 (正式写手)

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

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

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

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

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

» 猜你喜欢

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

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

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

huangpengru

铜虫 (小有名气)

我现在也在做类似的,可以借鉴借鉴!!
9楼2010-10-10 18:59:08
已阅   回复此楼   关注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的回帖

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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
请填处理意见