前面文章:建无序格点参杂&合金模形工具包ATAT的安装与其各模块的详细使用例子示范了ATAT-mcsqs的安装与使用。本章笔者用LAMMPS做高熵合金模拟的过程中发现:对于大体系如果用mcsqs建模,其in文件中每行坐标结尾都需要手动添加元素成分及其含量,对于上千上万的原子数体系工作量特别大。因此本文通过shell操作展示其快速写in文件的建模方法。
第一步:用atomsk建立一个10*10*10的bcc_W晶体盒子,并转换成VASP的POSCAR。
atomsk --create bcc 3.165 W -duplicate 10 10 10 -frac vasp
第二部将:POSCAR写成mcsqs的in文件形式。
先对POSCAR中原子坐标前几行修改一下
修改前:
修改后:
需要在每行原子坐标后面添加元素成分及其含量,但两千行手动添加特别麻烦,所以可以用以下命令添加。
sed -i "s/$/\t&Nb=0.2,Mo=0.2,V=0.2,Ta=0.2,W=0.2/g" POSCAR
处理后POSCAR如下:
只需要删除前4行多加的字符就可以了
按照建无序格点参杂&合金模形工具包ATAT的安装与其各模块的详细使用例子中的建模步骤可以得出以下高熵合晶模型:
PS:由于10*10*10的胞虚拟机运行较慢,笔者使用了5*5*5的胞
祝大家科研愉快!!!