ATAT 专题

ATAT快速生成VASP输入:ezvasp八种示例

从一个文本生成INCAR、POSCAR和KPOINTS,演示磁矩、空位、约束及二维网格语法。

更新于 2026-11-11约 3 分钟阅读按任务与环境选择适用配置
阅读本篇目录

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网格是否满足精度需要。

需要根据你的体系选择资源?联系讨论计算需求。实际可用软件与权限,以当前账户和运行环境为准。

查看配图

放大后可滚动查看细节,按 Esc 关闭。

100%查看原图 ↗