多语言展示
当前在线:1657今日阅读:176今日分享:34

vasp入门到精通[35]查看收敛情况(脚本)

这一节小编教大家写一个脚本来查看收敛情况
工具/原料

vasp,linux shell

方法/步骤
1

有时候我们想看一个比较大的体系的收敛情况,往往要跑很多离子步,你怎么才能看到能量的收敛情况呢?

2

比如小编的例子,有69个原子。而且小编的原子位置猜测的可能不合理,这样的话vasp会跑很多步,但是我怎么才能知道收敛的情况呢?

4

嘻嘻,小编强行配图了,这里小编给大家写一个命令来快速的实现查看能量收敛其实很简单,懂shell的不要喷小编,小编shell也很强的,如果小伙伴留言,小编会给大家写一些高大上的脚本或者程序(废话有点多)

5

命令如下:grep F= OSZICAR |awk '{print $1,$5}'(可以把命令复制到你的terminal窗口运行)

6

举个例子,如上图,小编算完了,看了下OSZICAR尾部,发现有可能没收敛,因为小编设置的EDIFF=100,跑了100步才停下来,一般就是没有收敛。小编想把上图红色部分提取出来,就是步数和能量。

7

运行命令:grep F= OSZICAR |awk '{print $1,$5}'之后的结果如下图

8

这就把数据提取出来了。我们需要把他保存成dat文件那就运行命令:grep F= OSZICAR |awk '{print $1,$5}' >conv.dat

9

会产生一个数据文件conv.dat,把它复制出来,直接拉到origin里就可以直观地看能量收敛了

10

这一节就到这里了,有需要,小伙伴可以留言哦,如果觉得还可以,就鼓励下小编吧(小编是个逗比,哈哈)你们的鼓励是小编最大的动力

推荐信息