| 查看: 6314 | 回复: 8 | ||||||
nkleof银虫 (正式写手)
|
[交流]
【求助】表面静电势怎么计算？
|
看到有的文献中做出了从表面深入到体材料的各层的静电势，如下图所示： 请问这个静电势数据从vasp计算的结果中如何提取？ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理和电化学 | 第一性原理 | 不错的帖子 |
» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- Multiwfn在AIM分析中的准确性如何？
已经有7人回复
- 利用Multiwfn计算Gaussian输出的激发态之间的跃迁偶极矩
已经有5人回复
- Multiwfn2013年暑期培训班Sobereva总结的相互作用研究文章的写作方法的ppt
已经有5人回复
- multiwfn绘制的ESP等值面图里静电势的默认单位是什么？
已经有8人回复
- 使用Multiwfn结合VMD分析和绘制分子表面静电势分布
已经有14人回复
- 将Zeta 电势换算成纳米颗粒自身表面存在的净电荷
已经有6人回复
- 如何计算静电势ESP以及作图
已经有12人回复
- 求助Cl-， Br-和I-（离子）的表面静电势或是电子密度
已经有2人回复
- 电化学计算中有关电极电势的问题
已经有7人回复
- 表面能计算公式中的Ebulk怎么计算？
已经有4人回复
- Gaussian计算和其他方法得到的热力学数据哪个更准确呢
已经有11人回复
- 【求助】计算表面的电势分布或者电荷分布
已经有18人回复
- 【求助】求一个基本物理知识！ 导体接电源，电源提供电势已知。导体表面电势怎么分布�
已经有10人回复
- 【讨论】关于静电势图的描述
已经有11人回复
- 【求助】关于晶体能级与水分子氧化还原电极电势的相对位置
已经有20人回复
liangab_234620
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1341.2
- 红花: 2
- 帖子: 131
- 在线: 145.5小时
- 虫号: 461188
- 注册: 2007-11-18
- 性别: GG
- 专业: 金属材料的力学行为
2楼2010-09-29 10:40:53
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
3楼2010-09-29 17:06:37
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
nkleof(金币+1):谢谢参与
cenwanglai(金币+2):谢谢回复！ 2010-09-30 09:18:36
nkleof(金币+7):非常感谢！ 2010-09-30 10:24:19
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
huangpengru
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 67.8
- 散金: 10
- 红花: 1
- 帖子: 67
- 在线: 77.3小时
- 虫号: 989003
- 注册: 2010-04-04
- 专业: 凝聚态物性 II ：电子结构
★ ★
nkleof(金币+1):谢谢参与
cenwanglai(金币+1):谢谢回复！ 2010-10-10 09:03:21
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
huangpengru
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 67.8
- 散金: 10
- 红花: 1
- 帖子: 67
- 在线: 77.3小时
- 虫号: 989003
- 注册: 2010-04-04
- 专业: 凝聚态物性 II ：电子结构
6楼2010-10-09 10:05:52
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
7楼2010-10-10 09:46:37
huangpengru
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 67.8
- 散金: 10
- 红花: 1
- 帖子: 67
- 在线: 77.3小时
- 虫号: 989003
- 注册: 2010-04-04
- 专业: 凝聚态物性 II ：电子结构
8楼2010-10-10 18:55:07
huangpengru
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 67.8
- 散金: 10
- 红花: 1
- 帖子: 67
- 在线: 77.3小时
- 虫号: 989003
- 注册: 2010-04-04
- 专业: 凝聚态物性 II ：电子结构
9楼2010-10-10 18:59:08
40