吸附与催化是当前热门研究领域，常见的计算结果包括催化剂结构，磁性，电子性质等，这些性质对催化剂有着决定性作用。本章将正式介绍催化领域的表面与低维结构模型构建，磁性体系处理，关键电子性质计算。
掺杂是指在纯净的晶体材料中有目的地引入外来原子（杂质原子）的过程。根据杂质原子在晶格中的位置，掺杂主要分为以下两类：
替位掺杂：杂质原子占据晶格中原有原子的位置。例如在硅晶体中掺入磷原子，磷原子占据硅原子的位置，形成n型半导体。这种掺杂方式会显著改变材料的电子结构，是最常见的掺杂类型。
间隙掺杂：杂质原子位于晶格空隙位置，不取代原有原子。例如在钛晶格中掺入碳原子，碳原子位于钛原子之间的间隙中。间隙掺杂通常会导致晶格畸变，影响材料的力学性能。
根据电荷性质，掺杂还可分为n型掺杂（提供电子）和p型掺杂（提供空穴），在半导体器件设计中具有重要应用价值。
在理论计算中，构建掺杂模型主要采用以下方法：
超胞法（Supercell Method）
超胞法是研究掺杂体系最常用的方法，其基本原理是：将原始晶胞按照特定倍数进行周期性扩展，形成超胞，然后在超胞中替换特定位置的原子为杂质原子，从而模拟掺杂效应。
例如，构建(2×2×2)的超胞后替换其中一个原子，相当于掺杂浓度为1/16=6.25%。超胞尺寸越大，掺杂浓度越低，计算精度越高，但计算成本也相应增加。
虚晶近似法（VCA）
虚晶近似法将掺杂原子视为平均化后的等效原子，每个格点位置上的原子是母体原子和杂质原子的加权平均。这种方法计算效率高，但无法捕捉杂质原子的局域效应。
格林函数法
基于密度泛函理论结合格林函数方法，适用于处理点缺陷和杂质问题，能够精确计算缺陷形成能和能级位置。
Materials Studio（MS）
Materials Studio是业界广泛使用的材料建模与模拟软件，在掺杂结构构建方面具有强大的功能：
核心功能模块：
Build菜单：提供”Build > Supercell”功能，可构建不同尺寸的超胞
Build菜单：提供”Build > Modify > Replace Atoms”功能，可进行原子替换和掺杂建模
CASTEP模块：内置DFT计算功能，可直接进行掺杂体系的几何优化和性质计算
具体操作流程：
1.导入原始晶体的CIF文件或POSCAR文件
2.选择”Build > Supercell”，设置扩胞矩阵（如2×2×2）
3.进入”Edit > Unit Cell Transformation”，去除对称性限制
4.使用”Build > Modify > Replace Atoms”功能，选择待替换原子并指定目标掺杂元素
5.进行几何结构优化后导出计算文件
VESTA
VESTA是一款免费开源的晶体结构可视化软件，适合进行结构编辑和掺杂建模：
主要功能特点：
可导入多种格式的结构文件（CIF、POSCAR、XYZ等）
支持从在线数据库（Materials Project、ICSD）直接下载结构数据
可通过”Edit > Edit Data”手动替换原子位置
支持原子坐标的直接编辑和原子删除操作
具体操作步骤：
1.从晶体结构数据库下载目标材料的CIF文件
2.导入到VESTA中，通过”Edit > Edit Data > Unit Cell”取消对称性
3.根据掺杂比例选择目标原子进行替换
4.调整晶胞参数后导出为VASP等计算软件兼容的POSCAR格式
获取初始晶体结构
构建掺杂模型的第一步是获取原始晶体结构。常用的结构数据库包括：
1. Materials Project
网址：https://materialsproject.org
包含超过14万种已知化合物结构
可直接下载CIF格式文件，支持在线搜索特定材料
2. ICSD（无机晶体结构数据库）
包含数十万种无机化合物结构
需通过机构订阅访问，是权威的结构数据来源
超胞尺寸选择与计算
步骤一：确定掺杂浓度
掺杂浓度与超胞尺寸的关系为：
掺杂浓度 = 1 / (超胞总原子数)
常见的掺杂浓度选择：
(2×2×2)超胞：6.25%掺杂浓度（适用于大多数研究）
(3×3×2)超胞：约3.47%掺杂浓度（适用于精确计算）
(4×4×4)超胞：约1.56%掺杂浓度（适用于高浓度模拟）
步骤二：构建超胞
在建模软件中选择超胞构建功能，设置扩胞矩阵。扩胞矩阵通常为正整数，如(2,2,2)表示在三个方向各扩展两倍。
关键注意事项：
扩胞后需检查原子数是否正确
注意超胞的对称性变化，必要时需取消对称性限制
超胞尺寸越大，计算时间越长，需平衡计算成本与精度
步骤三：原子替换操作
替换目标原子为掺杂元素时，需要考虑以下因素：
1.替位位置选择：选择合适的原子位置进行替换，通常选择等价位置
2.局域环境保持：尽量保持替代原子的配位环境不变
3.多重替位构型：对于较高掺杂浓度，可能存在多种替位构型，需考虑能量最低的构型
几何结构优化
构建完成的掺杂模型需要进行几何结构优化，优化参数包括：
1. 原子弛豫策略：
全原子弛豫：所有原子均可自由移动，计算精度高但计算量大
部分弛豫：仅弛豫掺杂区域附近的原子，远离区域原子固定
固定晶胞：保持晶胞参数不变，仅优化原子位置
2. 晶胞参数设置：
对于低掺杂浓度体系，通常固定晶胞参数为原始晶体结构
对于高掺杂浓度体系，允许晶胞参数弛豫以反映晶格畸变
3. 计算参数：
泛函选择：PBE-GGA是最常用的选择
平面波截断能：一般设置为400-600 eV
k点密度：随超胞尺寸增大而降低
收敛标准：能量收敛精度10⁻⁴-10⁻⁵ eV，力收敛精度0.02 eV/Å
超胞尺寸的选择直接影响计算结果的准确性。需要进行超胞尺寸收敛性测试：
测试方法：
1.保持计算参数一致，依次构建不同尺寸的超胞（如2×2×2、3×3×3、4×4×4）
2.计算形成能、能带结构和态密度等性质
3.当结果随超胞尺寸变化小于阈值时，认为已收敛
对于电子结构研究，一般建议至少使用(2×2×2)超胞，对于精确计算建议使用更大尺寸的超胞。
k点密度的选择需与超胞尺寸相匹配：
(2×2×2)超胞：建议使用4×4×4或6×6×6的k点网格
(3×3×2)超胞：建议使用3×3×3或4×4×2的k点网格
(4×4×4)超胞：建议使用2×2×2的k点网格
对于大超胞，可以使用Γ点近似（Gamma-only）来降低计算成本。
对于二维材料或表面掺杂模型，需要设置足够的真空层以避免上下表面相互作用。真空层厚度一般建议为15-25 Å，具体需通过测试确定。
掺杂结构模型的基本概念：掺杂的定义与分类、掺杂模型的计算方法
掺杂结构模型构建的具体步骤：获取初始晶体结构、超胞尺寸选择与计算、几何结构优化
关键参数设置与收敛性测试：超胞尺寸收敛性测试、k点采样设置、真空层设置
下一章将正式引入本次教程的核心—VASP一维纳米结构建模。我们将从一维纳米结构的基本概念与分类、一维纳米结构模型构建的具体步骤、关键参数设置与收敛性测试方面详细介绍一维纳米结构建模。
