武汉理工大学赵焱课题组开发脚本 MS建模一键获取VASP输入文件POSCAR


本文给大家推介并分享一个通过Materials Studio一键获取VASP输入文件POSCAR的脚本,xsd2pos.pl。此脚本是武汉理工大学首席教授赵焱在帮助学生解决VASP计算问题时编写,解决了其他间接非手动获取POSCAR存在的问题,同时也相较于其他方法更为便捷,相信会对从事VASP计算的同行们带来便利。赵焱老师怀着哺育天下学子的情怀,希望给大家分享一下。

下文介绍一下它的使用方法以及它所解决的问题,在此就不给大家介绍脚本的内容了,会编程的朋友们可自行学习。

使用方法

此脚本适用于用Materials Studio 建模获取POSCAR,操作非常简单。

(1)把xsd2pos.pl脚本直接拖至需要转化成POSCAR的xsd格式文件所在目录。示例中,我们将TiO2_single.xsd转化成POSCAR。

(2)将xsd2pos.pl打开,将第12行引号中的内容改写成需转化的xsd文件的名字

(3)点击上图紫红色框中的三角形按钮运行脚本便生成了POSCAR文件,如下图。

解决的问题

此脚本的诞生源于赵焱老师的学生运用cel2pos.pl 脚本和Vesta生成POSCAR时不能便利地选择性固定原子。赵焱老师为了解决这个问题,花了半天时间为学生们创造了该脚本。不得不说,尽管赵焱老师作为计算化学专家,编个脚本不费吹灰之力,但在繁重的教学任务、科研任务之下,半天的时间于他而言实为珍贵。

一般大家构造表面模型时,需要选择性固定几层原子以模拟体相结构。

如下图是一个TiO2(101)表面结构模型,通过Materials Studio我们可以很轻松地按需要固定原子。该示例中选择性固定了图中标黄的原子。

运用先前大家比较熟悉的cel2pos.pl将cell文件转化为POSCAR得到的原子结构信息如图,所有的原子都不是固定的,而要手动固定将会有很大的工作量。

另外,还有一个常用的方法就是用Vesta生成POSCAR,这种方法得到的POSCAR如下图所示,这种就没有原子固定与否的信息,等同于cel2pos.pl得到的POSCAR。

而我们运用xsd2pos.pl脚本能够直接将原子固定信息写进POSCAR。如图所示,在这里需要注意一下,选择性固定原子做动力学的时候“Selective Dynamics”必须有。

现在大家知道这个脚本的厉害了吧。

该脚本诞生仅有几天时间,应广大VASP计算同行的需求,将它分享出来。可能还有错误之处没有来得及发现和纠正,如果大家在之后的应用中发现错误,请反馈给我们(邮箱:yan60@hotmail.com),不胜感激。如陈述内容有错误之处,也请大家指正,我们会虚心接受并改正。

请点击本链接下载脚本

本文系赵老师的学生何同学供稿。

欢迎大家到材料人宣传科技成果并对文献进行深入解读,投稿邮箱tougao@cailiaoren.com。

内容合作及投稿请加微信cailiaorenvip.

分享到