ATAT 多元簇展开与统计

ATAT 专题导读 · 查看全部文章

三元 mmaps 和 memc2#

三元例采用固定 FCC 晶格常数 3.9 Å，物种为 Cu、Au、Al，使用 ASE EMT 势计算小晶胞结构的静态能量。它用于演示多元输入和数据组织，不作为该三元系的定量相图模型。

首先把 lat.in 的物种行写成 Cu,Au,Al，然后枚举和计算结构能量。拟合使用 mmaps，不能用仅支持二元的 maps 替代。多元基函数、组分坐标、化学势维数都会随物种数量变化。

本例 control.in 为：

300 0 0 0
1500 0 0 0 7

第一行给初始温度和各物种化学势；后续行规定扫描方向及采样数量。本例使用 -il 包含终点，得到 300、500、…、1500 K 共七点。

memc2 -gs=-1 -er=6 -eq=100 -n=500 -keV \
      -il -sd=2410 -tstat=0 -o=mc.out

物种顺序以 atoms.out 和程序输出为准，不凭 lat.in 中的书写顺序猜测。读取结果先打开 mcheader.out，它逐列标明本次输出的含义；不要套用二元 emc2 的固定列号。-gs=-1 是随机初态，不表示整个过程中组分固定。三个化学势都为零也不意味着三种元素各占 1/3；不同物种的相对能量仍会影响平衡组分。

需要守恒组分约束时，研究 conccons.in 和对应的初始构型设置；任意小数拼出的不相容约束可能使交换操作无法完成。

数据和配套算例#

下载完整实算资料包，保留 examples/common，进入所需算例目录执行 sbatch submit.sh。输入、参考结果、科学设置和重跑方式均在各目录的 README 中。