1、利用atat软件包的mcsqs产生随机点阵;并说明使用的过程和注意事项。
2、具体步骤
2.1、输入文件的格式;名字叫做rndstr.in 格式如下
1 1 1 90 90 90 #晶胞,也可以说是坐标系
0 0.5 0.5 #原胞
0.5 0 0.5 #原胞
0.5 0.5 0 #原胞
0 0 0 Fe=0.5 Cr=0.25 Al=0.25 #等效点位,以及比例,原胞中只有一个所以只有一个点位。
2.2 输入 corrdump -l=rndstr.in -ro -noe -nop -clus -2=1.1 对于hcp -2=1.2是可以的,且一个原胞包含俩个原子。所以是用54比较合理
-l 指定输入文件
-ro 读取占用比例
-noe 不要空的群
-nop 不要点群
-clus 仅查找群
-2或-3 指定第二近邻和第三近邻的距离,可以和自己建立的晶胞进行比较。
2.3 完成计算后,使用mcsqs -n xx 其中xx 是指定超胞中的原子的个数。注意这个个数要能被rndstr.in中的比例所整除,不然会报错。
2.4 但是有时后我们会要求晶胞是立方的而不是其他的形状。所以先暂时中止计算,然后运行sqscell.py脚本,查找sqscell.out中的等边的超胞结构。然后打开sqscell.out找到立方超胞,删除多余的结构。如果出现索引报错时可能是sqscel.out还没输出完,多等一会就行
2.5 输入mcsqs -rc 重新计算。
2.6 完成后可以用str2cif <bestsqs.out> 输出文件名.cif 转换成cif格式。
3、参考文献如下 其他教程
3.2 mcsqs使用教程