应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】表面静电势怎么计算？
51/1返回列表
查看: 5974  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

nkleof

银虫 (正式写手)

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

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

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-09-29 09:17: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.'
之前少了个空格
赞一下(2人) 回复此楼
5楼2010-10-09 09:31:17
已阅   回复此楼   关注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人) 回复此楼
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
要注意的就是各个量的单位。
赞一下(2人) 回复此楼
4楼2010-09-30 08:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭