ATAT快速生成VASP输入:ezvasp八种示例
从一个文本生成INCAR、POSCAR和KPOINTS,演示磁矩、空位、约束及二维网格语法。
阅读本篇目录
ezvasp把计算参数、结构和k网格控制放在一个输入文本中,再生成VASP需要的文件。下载实例包。本篇演示ATAT3.52原生输入生成;下一篇接实际VASP6.4.2计算。
1 先生成文件,不启动计算#
mkdir ATAT-ezvasp
python3 -m zipfile -e atat-ezvasp-*.zip ATAT-ezvasp
cd ATAT-ezvasp
bash prepare.sh
八套输入写在work,此步骤不运行VASP。脚本加载atat/3.52-8173m;ATAT工具使用GCC13.3环境。为了每个实例自带配置,脚本从已安装的ezvasp生成一个仅改变配置文件读取位置的私有副本,使用当前目录的ezvasp.rc,不改用户主目录配置。
赝势由有访问权限的集群PBE64势库读取;下载包不包含POTCAR。换到自己的机器时,要修改prepare.py里的势库路径并使用自己获授权的赝势。
2 单文件输入的主体#
cat inputs/Al-gamma.in
cat work/Al-gamma/INCAR
cat work/Al-gamma/POSCAR
cat work/Al-gamma/KPOINTS
输入的完整内容如下,其中KPPRA控制按原子数折算的网格目标,最终实际网格还取决于晶胞形状;它不是直接写死每个方向的k点数。
[INCAR]
SYSTEM=ezvasp input and native CPU MPI example
ENCUT=500
PREC=Accurate
EDIFF=1e-6
ISMEAR=1
SIGMA=0.1
ISPIN=1
ISYM=0
NSW=0
IBRION=-1
ALGO=Normal
NELM=100
LREAL=.FALSE.
LASPH=.TRUE.
LWAVE=.FALSE.
LCHARG=.FALSE.
USEPOT=PAWPBE
KPPRA=512
KSCHEME=Gamma
[POSCAR]
Al FCC primitive
1.0
0 2.02 2.02
2.02 0 2.02
2.02 2.02 0
Direct
0 0 0 Al
这是ezvasp扩展语法,结构行尾带元素名;不能直接把整个文件当成VASP的POSCAR。该Al结构的常规晶格参数4.04 Å,原胞体积16.484816 ų,实际生成8×8×8 Gamma网格。
3 八套输入分别演示什么#
| 输入目录 | 重点 | 实际生成网格 |
|---|---|---|
| Al-gamma | Gamma中心网格 | 8×8×8 |
| Al-monkh | Monkhorst-Pack网格语法 | 8×8×8 |
| Al-cell-patch | 原生-c晶胞处理选项 |
8×8×8 |
| Al-explicit-mesh | 显式[KPOINT]网格 |
6×6×6 |
| Al-selective-vac | T T F约束、Vac占位去除 |
7×7×7 |
| Al-alias | SUBATOM将标签映射为元素 |
8×8×8 |
| Fe-magnetic-seed | Fe+2.2生成磁矩初值 |
8×8×8 |
| Al-kz1 | KZ1把第三方向设为1 |
23×23×1 |

4 修改网格和磁矩#
显式网格使用本工具的[KPOINT]节名,生成后文件名仍为VASP标准的KPOINTS:
[KPOINT]
Explicit Gamma grid
0
Gamma
6 6 6
0 0 0
cat inputs/Al-explicit-mesh.in
cat work/Fe-magnetic-seed/INCAR
cat work/Al-selective-vac/POSCAR
Fe+2.2会生成ISPIN=2与MAGMOM=2.2,这是磁矩初猜,不是计算后的磁矩;Vac在生成结构时去掉;T T F保留相应坐标约束。Fe和含真空的Al例子仅演示输入语法,本批没有对它们作DFT材料结论。
生成文件后先核对原子数、元素顺序、约束方向、网格和势库,再按两阶段计算教程提交。自己的体系仍需检查截断能与k网格是否满足精度需要。
需要根据你的体系选择资源?联系讨论计算需求。实际可用软件与权限,以当前账户和运行环境为准。