Fermi-Softness-for-VASP：VASP计算费米软度的脚本

linqiaosong · 发表于 Post on 2019-8-18 14:17:09

本帖最后由 linqiaosong 于 2021-11-14 23:43 编辑

Fermi-Softness-for-VASP：VASP计算费米软度的脚本

最近更新了一个处理VASP结果得到总费米软度，简缩费米软度和局域费米软度的脚本，基本上可以实现费米软度的计算了。
这个脚本和使用教程可以从以下网站获得：
Fermi-Softness-for-VASP: A script dealing with VASP result for calculating fermi-softness (github.com)

关于什么是费米软度，实际上是用于衡量表面活性的一个反应性指标，具体的内容可以参看下面的参考文献：
B. Huang, L. Xiao, J. Lu, L. Zhuang, Angew. Chem. Int. Ed. 2016, 55, 6239 –6243 出版商处原文

使用这个脚本需要预先安装:

Anaconda3 (Python >= 3.8.5)：下载地址
ASE：安装方法
bader：官方主页
vaspkit 1.3：安装方法

暂时只测试了VASP 5.4.4版本和vaspkit 1.3版本可以兼容，其它版本是否兼容尚不明确。

附件是参考教程，仅供参考。

QE和CP2K的费米软度脚本最近也会推出。

sobereva · 发表于 Post on 2019-8-18 14:42:23

可以参考参考http://www.andrew.cmu.edu/user/feenstra/wavetrans/

linqiaosong · 发表于 Post on 2019-8-18 14:58:30

sobereva 发表于 2019-8-18 14:42
可以参考参考http://www.andrew.cmu.edu/user/feenstra/wavetrans/

感谢sob老师

ggdh · 发表于 Post on 2019-8-19 12:47:01

你快写一个傻瓜都能用的自动分析费米软度的程序，加一个猴子都能学会的教程+手册，把庄林牌方法推广出去

卡开发发 · 发表于 Post on 2019-8-20 15:54:28

本帖最后由卡开发发于 2019-8-20 15:56 编辑

ggdh 发表于 2019-8-19 12:47
你快写一个傻瓜都能用的自动分析费米软度的程序，加一个猴子都能学会的教程+手册，把庄林牌方法推广出去

建议把那个弟弟易吸和屁嗯必哦的鬼东西集成到一块去2333333，搞个好用点的.

ggdh · 发表于 Post on 2019-8-20 19:08:20

卡开发发发表于 2019-8-20 15:54
建议把那个弟弟易吸和屁嗯必哦的鬼东西集成到一块去2333333，搞个好用点的.

这就是周期性领域的Multiwfn的节奏了。要么再出一个sob这样的妖孽，要么就等sob有空了。。。
可是这两个的概率同样低。。。

linqiaosong · 发表于 Post on 2021-8-29 18:53:25

计算费米软度之前，VASP需要进行Slab模型的弛豫计算、自洽计算，然后读取自洽计算的电荷密度进行大k点的非自洽计算。

Para1lel · 发表于 Post on 2021-9-16 17:24:25

为什么我使用这个脚本总是提示ISPIN error。用的是VASP 5.4.4和vaspkit1.3，ASE也安装了。 INCAR里面设置的ISPIN=2。

Para1lel · 发表于 Post on 2021-9-16 21:04:40

我试了一下，设置ISPIN=2的话，脚本显示的是ISPIN=80，设置ISPIN=1的话，脚本显示的是11.

linqiaosong · 发表于 Post on 2021-10-6 23:09:46

本帖最后由 linqiaosong 于 2021-10-6 23:26 编辑

Para1lel 发表于 2021-9-16 21:04
我试了一下，设置ISPIN=2的话，脚本显示的是ISPIN=80，设置ISPIN=1的话，脚本显示的是11.

新脚本已更新

Para1lel · 发表于 Post on 2021-10-12 09:08:48

linqiaosong 发表于 2021-10-6 23:09
新脚本已更新

新脚本好像还是有问题，重复您的例子出现NameError: name 'ispin' is not defined

linqiaosong · 发表于 Post on 2021-10-12 10:51:49

Para1lel 发表于 2021-10-12 09:08
新脚本好像还是有问题，重复您的例子出现NameError: name 'ispin' is not defined

get_paraments函数里面
# get ispin下面的
spin=...
应该是ispin=...
应该是我不小心把变量名前面的i给搞掉了，改过来就行了

Para1lel · 发表于 Post on 2021-10-13 15:17:53

linqiaosong 发表于 2021-10-12 10:51
get_paraments函数里面
# get ispin下面的
spin=...

261行这个地方的*通配符好像读取不了，python总是把它当具体字符了。

Para1lel · 发表于 Post on 2021-10-13 15:20:20

linqiaosong 发表于 2021-10-12 10:51
get_paraments函数里面
# get ispin下面的
spin=...

如果把这个地方改成具体文件名，脚本可以正常运行了。
比如ISPIN=1的时候把*去掉。
ISPIN=2的时候，改成_UP, _DW.
但是用通配符*的话就报错了。

linqiaosong · 发表于 Post on 2021-10-13 23:09:03

Para1lel 发表于 2021-10-13 15:20
如果把这个地方改成具体文件名，脚本可以正常运行了。
比如ISPIN=1的时候把*去掉。
ISPIN=2的时候， ...

做了如下更新，最近网不好，GitHub同步不上去，暂时还没有同步代码

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] Fermi-Softness-for-VASP：VASP计算费米软度的脚本

评分 Rate