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

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

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

3、答案很简单，看每一步的能量，然后看变化幅度。这些数据OUTCAR和OSZI觊皱筠桡CAR里面都有，问题是，怎么快速的看到能量收敛呢？难道要每一步能量都复制到origin里么（oh my gold, are you kidding me?）

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、这一节就到这里了，有需要，小伙伴可以留言哦，如果觉得还可以，就鼓励下小编吧（小编是个逗比，哈哈）你们的鼓励是小编最大的动力