吸附与催化是当前热门研究领域，常见的计算结果包括催化剂结构，磁性，电子性质等，这些性质对催化剂有着决定性作用。本章华算科技朱老师将正式介绍催化领域的表面与低维结构模型构建，磁性体系处理，关键电子性质计算。
二维异质结构是指通过将两种或多种不同性质的二维材料（如石墨烯、过渡金属硫族化合物TMDCs、六方氮化硼h-BN等）在原子尺度上进行集成而形成的复合结构。
这类结构因其独特的物理、化学和电子性质，在光电子器件、自旋电子学、催化、能源存储等领域展现出巨大的应用潜力。构建高质量的二维异质结构模型，无论是从理论计算角度还是从实验制备角度，都是当前凝聚态物理和材料科学领域的重要研究课题。
根据构建方式和材料维度组合，二维异质结构主要可以分为以下几类：
垂直异质结构是通过范德华力将不同材料沿面外方向（c轴）堆叠在一起形成的层状结构。这是最常见也最容易构建的异质结构类型。其特点是界面结合力较弱，主要是范德华相互作用，不会破坏各层材料的本征电子结构。
横向异质结构是在同一平面内通过共价键或配位键将不同二维材料区域并置形成的结构。这种结构通过原子级精准的界面连接，在界面处会产生独特的电子性质。
根据参与构建的二维材料与零维（量子点）、一维（纳米线）或三维（块体材料）材料的组合方式，还可以形成2D/0D、2D/1D、2D/3D等多种复合异质结构。
从理论计算角度构建二维异质结构模型，通常采用密度泛函理论（DFT）结合材料模拟软件如Materials Studio（MS）、VASP、Quantum ESPRESSO等工具。以下是典型的建模流程：
准备原胞文件
首先，需要分别获取构成异质结构的各层二维材料的原始晶体结构文件（通常为POSCAR格式，用于VASP计算；或CIF格式，用于其他软件）。在准备过程中，需要准确记录每种材料的晶格常数、原子位置、空间群等关键参数。
晶格匹配与超胞构建
由于不同二维材料的晶格常数存在差异，直接堆叠会引入较大的晶格失配，导致界面应变过大而不稳定。因此需要进行晶格匹配计算，确定合适的扩胞组合。
常用的方法是采用m-n匹配法，即对两种材料分别进行扩胞，寻找使得两者晶格常数最接近的扩胞倍数组合。具体操作可使用专门的脚本工具，如m_n_supercell.py等。该脚本输入两种材料的晶格常数，会输出可能的m-n扩胞组合和对应的扭转角（twist angle），为后续建模提供指导。
m-n组合的选择需要考虑几个因素：一是晶格失配度（mismatch），通常要求低于5%以保证结构的稳定性；二是超胞大小，过大超胞会导致计算成本急剧增加；三是可能的扭转角，某些特定的扭转角（如魔角）会赋予异质结构特殊的电子性质。
构建多层结构
如果需要构建多层异质结构，可以使用multiplayer.py等脚本工具。在运行脚本时，需要输入以下参数：
层数（number of layers）
层间距（distance between layers），通常设置为2-3.5埃
原子配位类型（coordination type），建议优先选择D类型以获得更好的VESTA显示效果
旋转角度（rotation angle）
脚本会根据输入参数自动生成多层结构的POSCAR文件。
异质结组装
完成各层材料的独立构建后，使用heterojunction.py等工具进行最终的异质结组装。在此过程中，需要设置层间距离（interlayer spacing），通常在3-4埃之间，并选择适当的堆叠方式（stacking mode），如AA堆叠、AB堆叠、AA’堆叠等。
此外，还需要在堆叠方向上添加足够的真空层（vacuum layer），通常设置为15-20埃，以避免周期性边界条件下相邻镜像层之间的相互作用。
结构优化与验证
异质结构模型构建完成后，需要进行结构优化（relaxation）计算，使得原子位置和晶格参数达到能量最低的稳定状态。在优化过程中，可以对水平方向的晶格矢量（a、b轴）进行选择性优化，而对垂直方向（c轴）和真空层保持不变。
优化完成后，需要计算异质结构的结合能（binding energy）、声子谱（phonon spectrum）等参数，验证结构的力学和动力学稳定性。
在Materials Studio软件中构建异质结构的流程相对更为直观：
1.导入结构：首先分别导入构成异质结的两个或多个单体结构文件。
2.晶面选择：对每个单体结构选择合适的晶面（如(001)面或(111)面）并建立表面模型（slab model）。
3.晶格匹配检查：通过软件的匹配工具检查两个晶面是否匹配。系统会显示晶格失配度（mismatch），一般低于5%即可认为匹配良好。
4.构建异质层：使用Build Layers功能，选择layer1（或Average）作为基面，执行构建操作。当晶面不匹配时，系统会弹出提示，可选择接受失配。
5.添加真空层：通过Build → Crystal → Rebuild Crystal功能，将c矢量的长度增大至15埃以上，形成足够的真空层空间。
6.原子坐标调整：如有必要，可以对原子位置进行微调，以寻找最稳定的堆叠方式。
二维异质结构的主要类型：垂直异质结构、横向异质结构、维度组合异质结构
理论建模方法：基于第一性原理的计算建模、Materials Studio建模流程
下一章将正式引入本次教程的核心—vasp无磁结构计算。我们将从无磁结构计算的基本原理、INCAR关键参数设置、计算流程与步骤方面详细介绍无磁结构计算。
