本教程演示如何使用 VASPKIT 的 803 功能模块，基于体相结构快速生成特定 Miller 指数表面的 slab 模型。以 (101) 面为例，展示完整的交互式操作流程。
一、前置准备
确保已正确安装 VASPKIT（本教程基于 Standard Edition 1.5.1）。
在工作目录下准备好体相结构的 POSCAR 文件（或 CONTCAR）。
确认已加载 VASP 环境变量，命令行可正常调用 vaspkit。
注意： 输入文件必须为标准的 VASP POSCAR 格式，且晶胞坐标建议采用笛卡尔坐标（Cartesian）或分数坐标（Direct）均可，但需保证结构合理、无重叠原子。
二、启动 VASPKIT
在包含 POSCAR 的目录下打开终端，输入以下命令启动：
$ vaspkit
启动后将显示 VASPKIT 主菜单，包含结构工具、电子工具及其他辅助功能：
\/// / _ _ \ Hey, you must know what you are doing. (| (o)(o) |) Otherwise you might get wrong results. o-----.OOOo--()--oOOO.------------------------------------------o | VASPKIT Standard Edition 1.5.1 (27 Jan. 2024) | | Lead Developer: Vei WANG (wangvei@icloud.com) | | Main Contributors: Gang TANG, Nan XU & Jin-Cheng LIU | | Online Tutorials Available on Website: https://vaspkit.com | o-----.oooO-----------------------------------------------------o ( ) Oooo. VASPKIT Made Simple \ ( ( ) _) ) / (_/ ===================== Structural Utilities ====================== 01) VASP Input-Files Generator 02) Mechanical Properties 03) K-Path for Band-Structure 04) Structure Editor 05) Catalysis-ElectroChem Kit 06) Symmetry Analysis 07) Materials Databases 08) Advanced Structure Models ===================== Electronic Utilities ====================== 11) Density-of-States 21) Band-Structure 23) 3D Band-Structure 25) Hybrid-DFT Band-Structure 26) Fermi-Surface 28) Band-Structure Unfolding 31) Charge-Density Analysis 42) Potential Analysis 44) Piezoelectric Properties 51) Wave-Function Analysis 62) Magnetic Analysis 65) Spin-Texture 68) Transport Properties ======================== Misc Utilities ========================= 71) Optical Properties 72) Molecular-Dynamics Kit 74) User Interface 78) VASP2other Interface 84) ABACUS Interface 91) Semiconductor Kit 92) 2D-Material Kit 95) Phonon Analysis 0) Quit
三、进入 Advanced Structure Models 模块
切面功能位于 08) Advanced Structure Models 子菜单中。在主菜单提示符后输入 803 直接调用切面生成工具：
------------>>
803
程序将提示输入 Miller 指数：
Input the Miller Indices (hkl, e.g. 1 1 0):
------------>>
四、输入 Miller 指数
以构建 (101) 表面为例，输入三个整数并用空格分隔：
------------>>
1 0 1
提示： Miller 指数 h k l 对应切面的法向方向。常见低指数面如 1 0 0、1 1 0、1 1 1 等，根据研究体系选择。
五、设置 Slab 重复层数
程序继续提示输入沿表面法向重复的晶胞层数（厚度控制）：
Input the Number of Repeated Unit Slab Layer (e.g., 1, 2, ...):
------------>>
1
此处输入 1 表示仅重复 1 层体相晶胞。若需更厚的 slab 以消除周期性镜像相互作用，可输入更大的整数（如 3 或 4）。
六、设置真空层厚度
接下来输入真空层厚度（单位：Å），用于隔离上下表面的周期性相互作用：
Input the Vacuuum Thickness (e.g. 10):
------------>>
15
注意： 真空层厚度通常建议不小于 15 Å，以确保表面偶极矩及静电势在真空区域充分衰减。若体系存在较强表面极化，可适当增加至 20–25 Å。
七、输出结果
参数输入完成后，VASPKIT 将自动枚举所有不等价的表面终止方式（surface terminations），并依次写入独立的 SLAB 文件：
-->> (01) written SLAB_101_0001.vasp file!
-->> (02) written SLAB_101_0002.vasp file!
-->> (03) written SLAB_101_0003.vasp file!
-->> (04) written SLAB_101_0004.vasp file!
-->> (05) written SLAB_101_0005.vasp file!
-->> (06) written SLAB_101_0006.vasp file!
-->> (07) written SLAB_101_0007.vasp file!
-->> (08) written SLAB_101_0008.vasp file!
最终程序输出结构摘要信息：
+-------------------------- Summary ----------------------------+ Possible Configurations: 8
Planar Lattice Angle: 73.132
Lattice Constants: 7.247 4.206
Slab Thickness: 5.046
Total Atoms: 8
+---------------------------------------------------------------+
八、结果解读
参数项
数值
含义
Possible Configurations
8
该 Miller 指数下共有 8 种不等价的表面终止构型
Planar Lattice Angle
73.132°
表面二维晶格的夹角
Lattice Constants
7.247, 4.206 Å
表面单胞的二维晶格常数 a 和 b
Slab Thickness
5.046 Å
原子层的实际厚度（不含真空）
Total Atoms
8
每个 slab 模型中的总原子数
九、后续处理建议
构型筛选：
根据表面化学环境（如暴露的原子种类、配位数、极性等）从 8 个文件中选择最合理的终止面。
结构优化：
将选定的 SLAB_101_000X.vasp 重命名为 POSCAR，并设置合理的 INCAR、KPOINTS、POTCAR 进行 VASP 结构优化。
固定底层原子：
若仅需研究表面性质，建议在 POSCAR 中固定底部半数原子（添加 Selective dynamics 标签），以模拟体相支撑效应并减少计算量。
偶极校正：
对于非对称 slab（如两种表面终止不同），建议在 INCAR 中设置 LDIPOL = .TRUE. 与 IDIPOL = 3 进行偶极矩校正。
