计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 15509 | 评论数: 21 | 收藏 Add to favorites 8
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2019-8-18 14:17

正文摘要:

本帖最后由 linqiaosong 于 2021-11-14 23:43 编辑 Fermi-Softness-for-VASP:VASP计算费米软度的脚本 最近更新了一个处理VASP结果得到总费米软度,简缩费米软度和局域费米软度的脚本,基本上可以实现费米软度的 ...

回复 Reply

abite 发表于 Post on 2021-11-15 10:07:30
linqiaosong 发表于 2021-11-14 23:25
我想应该是python3版本不同导致pathlib库出现了点问题,我用的python 3.8.5
如果正常生成了LFS.cube文件 ...

好的,谢谢
linqiaosong 发表于 Post on 2021-11-14 23:25:29
本帖最后由 linqiaosong 于 2021-11-14 23:27 编辑
abite 发表于 2021-11-14 12:13
谢谢你的回复,intermediate_file_options改了以后可以运行了,最后还有一个错误,我不知道这样[/backcol ...

我想应该是python3版本不同导致pathlib库出现了点问题,我用的python 3.8.5
如果正常生成了LFS.cube文件,那就不用管它了
abite 发表于 Post on 2021-11-14 12:13:20
linqiaosong 发表于 2021-11-12 10:23
我这里测试Pt3Y的例子是正常的,试一下把脚本更新到最新版本,如果不行试试把intermediate_file_options ...

谢谢你的回复,intermediate_file_options改了以后可以运行了,最后还有一个错误,我不知道这样有没有跑完计算分析过程。
linqiaosong 发表于 Post on 2021-11-12 10:23:57
abite 发表于 2021-11-10 17:41
试了一下PtY的例子,结果报错如下是什么原因呢?vaspkit、bader和ase都装了。
Traceback (most recent c ...

我这里测试Pt3Y的例子是正常的,试一下把脚本更新到最新版本,如果不行试试把intermediate_file_options改为True
abite 发表于 Post on 2021-11-10 17:41:30

试了一下PtY的例子,结果报错如下是什么原因呢?vaspkit、bader和ase都装了。
Traceback (most recent call last):
  File "runfs.py", line 454, in <module>
    run_fs(kbT,dfdd_threshold,band_gap,intermediate_file_options,bader_dir,vaspkit_dir)
  File "runfs.py", line 395, in run_fs
    fs,atoms=calc_lfs(para,kbT,dfdd_threshold,intermediate_file_options,vaspkit_dir)
  File "runfs.py", line 265, in calc_lfs
    q.unlink(True)
TypeError: unlink() takes 1 positional argument but 2 were given
linqiaosong 发表于 Post on 2021-11-9 20:59:19
abite 发表于 2021-11-9 17:03
你好,想问一下非自洽计算的时候就是读取CHGCAR然后写WAVECAR就可以吗,需要计算其它性质吗?

是的
abite 发表于 Post on 2021-11-9 17:03:05
你好,想问一下非自洽计算的时候就是读取CHGCAR然后写WAVECAR就可以吗,需要计算其它性质吗?
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.  
但是用通配符*的话 就报错了。
Para1lel 发表于 Post on 2021-10-13 15:17:53
linqiaosong 发表于 2021-10-12 10:51
get_paraments函数里面
# get ispin下面的
spin=...

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

QQ截图20211013151315.png (4.12 KB, 下载次数 Times of downloads: 61)

QQ截图20211013151315.png

QQ截图20211013151508.png (13.5 KB, 下载次数 Times of downloads: 23)

QQ截图20211013151508.png
Para1lel 发表于 Post on 2021-10-12 09:08:48

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

QQ截图20211012090714.png (24.33 KB, 下载次数 Times of downloads: 35)

QQ截图20211012090714.png
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-9-16 21:04:40
我试了一下,设置ISPIN=2的话,脚本显示的是ISPIN=80, 设置ISPIN=1的话,脚本显示的是11.

202109162103138981..png (12.79 KB, 下载次数 Times of downloads: 31)

202109162103138981..png

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 15:45 , Processed in 0.176260 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list