| 查看: 2927 | 回复: 6 | |||
[求助]
求教各位前辈，如何利用LOCPOT做求出真空能级？
|
|
如题，LOCPOT导入origin的数据如图中所示，应该以哪一列作为x，哪一列作为y呢？这三列各列的含义是什么？origin数据中的那三个图，放大后就是123与234两幅图，感觉怪怪的。 请各位大侠不吝赐教啊   LOCPOT.jpg  123.jpg  234.jpg |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有12人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位有处理LOCPOT 文件的脚本呢？想计算真空能级和功函
已经有39人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+3, EDIFFG 2013-10-27 15:06:24
奔跑的爷们: 金币+3, ★★★★★最佳答案 2013-11-02 12:26:15
感谢参与，应助指数 +1
liliangfang: 金币+3, EDIFFG 2013-10-27 15:06:24
奔跑的爷们: 金币+3, ★★★★★最佳答案 2013-11-02 12:26:15
|
LOCPOT只有三列数吗？不应该是n1xn2xn3吗 (n1,n2,n3是格点数) 文件前几行有这个信息 这个网格在你的超晶胞分成n1xn2xn3个点 每个点上有一个charge density的值（就是chargcar的数据） 分网格点的时候沿着a方向变化最快 然后是b 再是c 比如你的超晶胞是正交的 10埃x9埃x8埃 原点在（0，0，0） 格点是10x9x8 那么你一共有720个数据在chargcar里 前10个数据对应超晶胞中的格点位置是（0，0，0），（1， 0，0），（2，0，0）。。。 后十个数据对应格点位置是（0，1，0），（1，1，0），（2，1，0）。。。。 一般沿着z方向是真空吧 你就要每90个数据做一个平均（x,y平面上）这样一共得到8个值 就行了 不知道说明白没有 自己再看看vasp说明书吧 |
2楼2013-10-26 21:28:32
3楼2013-11-01 20:35:01
【答案】应助回帖
★ ★ ★
sunyang1988: 金币+1, 谢谢交流 2013-11-01 22:48:30
奔跑的爷们: 金币+2, ★★★很有帮助 2013-11-02 12:26:23
sunyang1988: 金币+1, 谢谢交流 2013-11-01 22:48:30
奔跑的爷们: 金币+2, ★★★很有帮助 2013-11-02 12:26:23
|
肯定要写程序了 我用网上的一个程序 https://github.com/compphys/ase_tools asetools文件夹里有个文件vtotave.py |
4楼2013-11-01 22:23:09
5楼2013-11-02 12:25:55
|
这个是python script 不用编译 有一些python的package要安装 最主要的是ase https://wiki.fysik.dtu.dk/ase/ 如果没接触过python 可能要花点时间了解 使用很简单 直接python vtotave.py LOCPOT z |
6楼2013-11-02 12:59:50
7楼2013-11-02 13:05:09
50