将VASP中的结构优化文件中的INCAR修改成SCF的样子(Python脚本实现)
代码描述
已经有了SCF需要的5个文件,但是INCAR跟POSCAR中的INCAR文件一样(其实就是上个脚本整理好的文件)。此代码可以自动修改INCAR文件,具体修改内容如下:
- 将IBRION设置成-1
- 在文件的第15行,加上ICORELEVEL=1
- 如果LORBIT=11被注释掉,将其打开
- 如果LVHAR=.TRUE.被注释掉,将其打开
- 设置NSW=0
对文件的要求
实现过程
代码有3个函数
函数INCAR(address):是将此地址下的INCAR文件按照要求修改,生成INCAR_new文件
函数action1(address):遍历当前地址下的所有INCAR文件,每找到一个执行一次INCAR(addressd)函数
函数action2(address):遍历当前地址下的所有INCAR_new文件,每找到一个,将其名字修改为INCAR
使用方法
- 将第二行的address修改成想要操作的地址
- 将最后一行修改成action1(address)
- 运行,看输出是否正确
- 将最后一行修改成action2(address)
- 运行即可
import os
address="/home-gg/users/nsgg276_SJ/makaichuang/metal-2DGaO-1/simlulation/SCF"
def INCAR(address):
f=open(address+"/INCAR",encoding="UTF-8")
ff=open(address+"/INCAR_new","w",encoding="UTF-8"