小木虫

当前位置：首页 > 第一原理 >【求助】表面静电势怎么计算？

【求助】表面静电势怎么计算？

作者 nkleof: 来源: 小木虫 400 8 举报帖子

+关注

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

请问这个静电势数据从vasp计算的结果中如何提取？返回小木虫查看更多

今日热帖

精华评论

liangab_234620

VASP老家的例子文件里面有
nkleof

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

谢谢回复，看起来好像是在INCAR里面设置LVTOT=TRUE，然后在LOCPOT文件里面能够得到the planar-averaged potential 。
不过还是不太明白如何处理，能不能提供详细一点的计算并处理相关数据的资料？
valenhou001

网上可以找到相关的工具，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
要注意的就是各个量的单位，
huangpengru

看帖http://muchong.com/bbs/viewthread.php?tid=2357737&fpage=1
复制过来的时候记住：
WRITE(*,*) 'Error opening or reading file LOCPOT.'
之前少了个空格
huangpengru

lz能否告诉我你截图的那篇文章的题目？多些！！！！
nkleof

引用回帖:
Originally posted by huangpengru at 2010-10-09 10:05:52:
lz能否告诉我你截图的那篇文章的题目？多些！！！！

Reuter, K. and M. Scheffler (2001). "Composition, structure, and stability of RuO2(110) as a function of oxygen pressure." Physical Review B 65(3): 035406.
这篇好像是做DFT表面稳定性计算比较早而且成熟的一个工作，后来有不少文献都引用它。
huangpengru

多谢lz！！！