吸附与催化是当前热门研究领域，常见的计算结果包括催化剂结构，磁性，电子性质等，这些性质对催化剂有着决定性作用。一个可靠的计算模型是获得可信结果的前提。本章聚焦模型构建与参数测试，这正是初学者最容易忽视、却又最关键的环节。电子性质分析（态密度、差分电荷、d带中心等）将在第二章详细讲解。
晶体表面模型是用于模拟和分析晶体表面物理化学性质的计算模型。由于真实晶体具有三维周期性结构，而表面则是三维空间中的二维界面，因此构建表面模型时需要将三维周期性转换为二维周期性结构。
在计算方法上，最常用的模型是slab模型（薄片模型），即通过将三维晶胞沿特定方向切割，在垂直于表面方向上添加真空层，从而形成二维周期性的表面结构。这种模型广泛应用于密度泛函理论（DFT）计算中，能够有效模拟无限大表面的物理和化学行为，同时通过真空层隔离上下表面之间的相互作用。
晶体表面模型的构建涉及多种专业软件，目前常用的工具包括：
Materials Studio（MS）
Materials Studio是业界广泛使用的晶体建模软件，其内置的”cleave surfaces”功能可以轻松构建任意晶面的表面模型。具体操作流程为：导入晶体结构文件，选择”Build > Surfaces > Cleave Surfaces”，输入目标晶面指数（如110、100、111等），设置层数参数，然后添加真空层，最后导出为计算软件所需的格式。
VESTA
VESTA是一款优秀的晶体结构可视化软件，可用于查看、编辑和导出晶体结构。通过VESTA可以进行表面切割操作，并将结果导出为CIF、POSCAR等格式，便于后续VASP等计算软件使用。
Avogadro
Avogadro是一款开源的分子建模软件，内置了丰富的晶体结构库。通过”Crystallography > Build > Cleave”功能，可以生成指定Miller平面的表面模型，支持2D和3D超胞的构建。
ASE（Atomic Simulation Environment）
ASE是Python编写的一套原子模拟环境工具包，提供丰富的表面构建函数。例如，ase.build.fcc111可以快速构建面心立方晶体的(111)表面，参数包括尺寸大小、真空层厚度等。
Packmol
对于复杂体系或特殊结构，Packmol可以解决约90%的复杂建模问题，关键在于合理编写输入脚本。
步骤一：获取晶体结构
从晶体结构数据库获取初始晶胞数据，常用的数据库包括：
ICSD（无机晶体结构数据库）
Materials Project
American Mineralogist Crystal Structure Database
需要获取的信息包括空间群、晶格参数、原子坐标等 。
步骤二：表面切割操作
在建模软件中选择“cleave surface”或”切片”功能，设置以下关键参数：
晶面指数（hkl）：根据研究目标选择合适的晶面，如(100)、(110)、(111)等。不同晶面的表面原子排列和性质存在显著差异，直接影响计算结果 1 1 。
层数参数（depth/Thickness）：决定模型中包含的原子层数。一般建议设置3-7层，具体取决于材料体系和计算精度要求。较厚的模型可以更准确地模拟体相性质，但会增加计算成本 1 1 。
位置参数（position）：控制切割的起始位置，可以选择不同的切割面以得到所需的表面结构 1 1 。
步骤三：添加真空层
真空层的作用是将三维周期性结构在垂直于表面的方向上断开，形成二维周期性结构。真空层厚度的典型值为10-20 Å，具体需根据体系进行收敛性测试 1 1 1 。
步骤四：构建超胞
对于需要更大表面面积或特定对称性的体系，可通过“supercell”功能构建超晶胞。超胞尺寸（如2×2、3×3）的选择影响k点采样和计算精度 1 2 。
步骤五：结构优化与导出
构建完成的模型需要进行结构优化，优化参数包括：
是否固定底层原子（固定底层可以加快计算速度）
是否全原子弛豫
晶胞形状和体积是否允许变化
优化完成后，将模型导出为计算软件所需的格式，如VASP的POSCAR文件、Quantum ESPRESSO的输入文件等。
真空层厚度
真空层厚度的选择直接影响计算结果的准确性。通常需要进行真空层收敛性测试：保持基底层数不变，依次增加真空层厚度（如8 Å、10 Å、15 Å、20 Å），计算体系总能，当总能随真空层厚度变化小于某一阈值（如0.01 eV）时，认为已收敛。
对于吸附体系，真空层厚度还需考虑吸附物的高度，一般建议设置为吸附物最高点到表面距离的2-3倍。
基底厚度
基底厚度（即slab层数）的确定需要进行收敛性测试。常见的测试方法是：
1.固定真空层厚度
2.逐步增加基底层数（如3层、5层、7层、9层）
3.计算表面能或吸附能
4.当结果随层数变化小于设定阈值时，认为已收敛
对于金属表面，一般建议至少7-9层原子；对于氧化物表面，可能需要更多层以保证体相性质的准确模拟。
原子弛豫设置
表面原子的弛豫对计算结果影响很大，常见的弛豫策略包括：
1.全原子弛豫：所有原子层均可自由移动，计算精度高但计算量大
2.部分弛豫：仅弛豫表面2-3层原子，底层原子固定
3.单面弛豫：仅弛豫上层表面，底层完全固定
部分弛豫策略在保证计算精度的同时显著提高效率，是目前常用的做法。
k点采样
表面模型是二维周期性结构，k点采样需要在平行于表面的方向充分，垂直于表面的方向只需Γ点。k点密度的选择需要进行收敛性测试，一般从较密的k网格开始，逐步降低密度直到结果变化小于阈值。
偶极子修正
当表面模型上下不对称时（如仅单层表面吸附），体系会产生偶极矩，可能导致计算误差。此时需要启用偶极子修正（dipole correction），在VASP中通过设置LDIPOL = .TRUE.和IDIPOL = 3。
构建表面模型时需注意：
切割后可能失去原晶体的全部对称性，需重新确定空间群
某些晶面的切割可能导致表面重构，需要考虑这种可能性
超胞构建时需保持必要的对称性以简化计算
对于化合物表面，构建的slab模型应尽量保持与体相相同的化学计量比。例如TiO₂(110)表面，模型中Ti与O的比例应为1:2，这样电子结构更稳定，表面能计算也更简便。
进行收敛性测试时，需保持计算级别一致（相同的泛函、截断能、收敛标准等），否则测试失去意义。特别是不同层数测试时，需保证真空层方向的晶胞长度一致3。
表面能的计算是验证模型合理性的重要手段，计算公式为：
γ = (Eslab – n*Ebulk) / (2*A)
其中Eslab是表面模型能量，Ebulk是单位晶胞能量，n是原子数比例，A是表面积。
晶体表面模型的构建方法：介绍了软件工具选择、具体操作步骤
关键参数设置与收敛性测试：真空层厚度、基底厚度、原子弛豫设置
常见问题与注意事项：模型对称性、化学计量比、计算参数一致性、表面能计算
下一章将正式引入本次教程的核心—VASP二维材料建模。我们将从二维材料模型的构建方法、二维材料模型构建的具体步骤、异质结与堆垛模型构建方面详细介绍二维材料建模。
