当前位置：首页 > news >正文

第一性原理计算在半导体缺陷研究中的应用：以氢掺杂氧化镓为例

news 2026/5/23 6:24:40

1. 项目概述：从“掺杂”与“缺陷”说起

在半导体材料的研究与开发中，我们常常听到“掺杂”这个词。简单来说，就像在炒菜时撒入不同的调料来改变风味，掺杂就是在纯净的半导体材料（本征材料）中，有目的地引入微量的其他元素（杂质原子），从而精确调控其电学、光学等物理性质。氧化镓（Ga₂O₃）作为一种新兴的超宽禁带半导体材料，因其在高压、高温、高功率电子器件以及深紫外光电器件领域的巨大潜力而备受关注。然而，材料的本征性能往往无法满足所有应用需求，这时，氢（H）掺杂就成了一种关键的技术手段。

你可能会问，为什么要研究氢掺杂？氢是宇宙中最轻、最常见的元素，在材料生长和加工过程中极易以各种形式（如H₂、H₂O、OH⁻）被引入晶体内部。这些氢原子并非总是“捣蛋鬼”，它们有时会扮演“调解员”的角色：钝化有害的缺陷、改变载流子浓度、甚至作为施主或受主直接参与导电。因此，理解氢在Ga₂O₃中的存在形式、占据位置以及对材料性能的影响，对于制备高性能的Ga₂O₃基器件至关重要。

“H掺杂Ga₂O₃的缺陷计算”这个项目，其核心就是运用基于密度泛函理论（DFT）的第一性原理计算，在原子尺度上模拟和探究氢原子进入β-Ga₂O₃晶体后可能形成的各种缺陷结构，并计算它们的形成能、电荷转移能级、以及引起的晶格畸变等关键物理量。这相当于在计算机里搭建一个虚拟的“材料实验室”，让我们能够以极低的成本和极高的分辨率，去观察和理解那些在真实实验中难以直接捕捉的微观过程。本次的“准备计算PREPARE02”阶段，正是这个虚拟实验的基石——构建准确、合理的计算模型，为后续复杂的缺陷性质计算铺平道路。

2. 计算整体设计与思路拆解

进行缺陷计算，尤其是涉及轻元素如氢的掺杂计算，绝非简单地在一个晶胞里替换或插入一个原子然后开始算。它需要一个严谨、系统的设计流程，以确保计算结果的物理合理性和可比性。整个项目可以拆解为几个逻辑严密的阶段，而PREPARE02正处于承上启下的关键位置。

2.1 计算流程全景图

一个完整的缺陷计算研究通常遵循以下路径：

结构优化与验证（PREPARE01）：首先获得完美的、无缺陷的β-Ga₂O₃原胞晶体结构。这需要从可靠的实验数据或高质量的理论计算中获取初始结构，然后进行充分的几何优化，直至晶格常数、原子位置、键长键角等与已知结果吻合。这一步是后续所有计算的参考基准。
超胞构建与参数测试（PREPARE02 - 本次核心）：由于缺陷会破坏晶体的周期性，我们需要构建一个足够大的超胞来容纳缺陷，并确保缺陷之间的相互作用可以忽略不计。同时，需要系统测试并确定一套最优的计算参数（截断能、K点网格、赝势等），在计算精度和耗时之间取得最佳平衡。
缺陷模型构建与初步优化：在优化好的超胞中，根据化学直觉和对称性，构建氢原子可能占据的各种位点模型，例如间隙位、替位位（取代氧或镓？这需要根据电荷平衡判断）、以及与空位复合的缺陷对等。对每个初始模型进行初步的几何优化。
缺陷形成能与电荷态计算：这是缺陷物理的核心。计算不同电荷态下缺陷的形成能，绘制形成能随费米能级变化的曲线。由此可以判断在何种生长条件下（富氧或富镓）哪种缺陷最易形成，以及缺陷的电荷转移能级（即缺陷在禁带中的能级位置），这直接决定了它是深能级陷阱还是浅能级杂质。
电子结构及性质分析：分析缺陷态的电子局域函数、态密度、能带结构等，从电子层面理解氢掺杂如何影响Ga₂O₃的导电类型、光学吸收等性质。

2.2 PREPARE02阶段的核心任务与决策逻辑

PREPARE02阶段聚焦于超胞和计算参数的准备，其每一个决策都直接影响后续计算的可靠性和效率。

1. 超胞尺寸的选择：在精度与代价间权衡

为什么需要超胞？在周期性边界条件下，一个单胞中的缺陷会与其在无限空间中的“镜像”发生相互作用，这种人为的“缺陷-缺陷相互作用”会严重歪曲计算结果，尤其是对于带电缺陷。构建超胞相当于将缺陷“稀释”，增大缺陷间的距离。
尺寸如何定？常见的做法是构建2×2×1， 2×2×2， 3×3×1等不同尺寸的超胞。选择标准是：进行一系列测试计算，观察目标性质（如缺陷形成能、缺陷态的能级位置）随超胞尺寸的收敛情况。通常需要确保超胞的边长大于缺陷波函数衰减的长度（一般需要>10 Å）。对于β-Ga₂O₃，其单胞包含12个Ga原子和18个O原子，一个2×2×2的超胞就含有96个Ga原子和144个O原子，这已经是一个240个原子的体系，对于DFT计算来说负担不小。因此，我们可能从2×2×1开始测试。
决策逻辑：我们会先用一个中等大小的超胞（如2×2×1）完成全套参数测试和部分缺陷模型的初步探索。如果发现某些关键性质（如带电缺陷的形成能）对尺寸敏感，再考虑升级到更大的超胞进行验证计算。这符合“先探索，后深化”的务实策略。

2. 计算参数的系统测试：找到“甜蜜点”DFT计算包含多个可调参数，错误的选择会导致结果谬以千里或计算资源浪费。

平面波截断能（ENCUT）：决定平面波基组的精度。测试方法：计算完美超胞的总能量随ENCUT的变化。当能量变化小于某个阈值（如1 meV/atom）时，认为能量收敛。选择刚好收敛的值为佳，过高会指数级增加计算量。
K点网格密度：用于布里渊区积分，影响电子态的描述。测试方法：计算完美超胞的总能量随K点网格密度的变化（如从Γ点-only，到2×2×2， 3×3×3…）。同样观察能量收敛情况。对于超胞，由于实空间变大，倒空间缩小，K点密度可以比原胞计算时取得更稀疏一些。
赝势的选择：这是基石。对于Ga₂O₃，Ga的3d电子是半芯态，对化学环境敏感，必须使用将Ga 3d电子作为价电子处理的赝势（如PAW_PBE Ga_d）。对于O，标准PAW_PBE即可。对于H，需要注意有些赝势对氢键的描述有特别优化。通常，从一个被广泛验证的、同一套赝势库中的选择是最稳妥的。
电子自洽收敛标准（EDIFF）和离子弛豫收敛标准（EDIFFG）：前者控制电子步迭代，通常设为1e-5到1e-6 eV；后者控制原子受力，当所有原子受力小于某个值（如0.01 eV/Å）时停止优化。设得太松，结构不准确；设得太紧，迭代步数激增。

实操心得：参数测试看似繁琐，但绝不能跳过。我习惯为每个新材料体系建立一个“参数测试日志”，记录下不同ENCUT和K点下的能量、计算时间、以及最终选定的值。这不仅是当前项目的依据，未来处理类似体系时也能快速上手，事半功倍。

3. 核心细节解析与实操要点

在PREPARE02阶段，除了宏观的参数选择，一些微观的细节处理同样决定了计算的成败。这些细节往往隐藏在软件的操作流程和输入文件的设置中。

3.1 氢掺杂初始模型的构建化学

将氢原子放入Ga₂O₃超胞，不是随便找个空位塞进去就行。我们需要基于化学知识和晶体结构，构建物理上合理的初始猜测模型。

间隙氢（H_interstitial）：氢原子最可能占据晶体中较大的间隙空间。需要分析β-Ga₂O₃的晶体结构（通常属于单斜晶系），找出由氧原子八面体或四面体配位形成的空腔。可以使用VESTA等可视化软件，或者通过计算“空隙”来辅助定位。将氢原子放置在这些空腔中心附近。
替位氢（H_substitutional）：氢原子取代晶格原子。这需要考虑电荷补偿。
- H替代O（H_O）：一个O原子是-2价，H通常是+1价。直接替换会导致局部正电荷过剩，因此这个模型通常需要伴随一个额外的电子（即缺陷带负电，H_O^(-1)）或者附近有一个带正电的缺陷（如镓空位V_Ga）来补偿，才能稳定存在。
- H替代Ga（H_Ga）：一个Ga原子是+3价。直接替换会导致局部负电荷过剩，因此这个模型通常带正电（H_Ga^(+1)）或需要其他补偿。从原子尺寸看，H（共价半径~30 pm）远小于Ga（~126 pm）和O（~66 pm），直接替位的可能性较低，但作为理论探索仍需验证。
氢-空位复合体（H-Vacancy Complex）：氢极易与材料中固有的点缺陷（空位）结合。例如：
- H与氧空位（V_O）结合：V_O在Ga₂O₃中通常是一个双电子施主。一个或多个H原子可以填入这个空位，形成(H)_nV_O复合体，这可能会钝化V_O的施主特性。
- H与镓空位（V_Ga）结合：V_Ga是受主。H的引入可能与它形成复合体，改变其受主能级。
- 构建技巧：先创建空位缺陷超胞，然后将H原子放置在空位周围能量有利的位置（如沿着原来断裂的化学键方向），进行多种初始取向的测试。

3.2 自洽计算与结构优化的关键设置

在VASP（一个常用的第一性原理计算软件）中，INCAR文件的设置是核心。

# INCAR 关键参数示例 (用于PREPARE02阶段的测试和优化) SYSTEM = beta-Ga2O3 supercell test # 系统名称 PREC = Accurate # 高精度模式，影响FFT网格等 ENCUT = 520 # 截断能，根据测试结果设定 (例如520 eV) ISMEAR = 0; SIGMA = 0.05 # 对于半导体/绝缘体，使用Gaussian smearing，SIGMA要小 IVDW = 11 # 可选，加入DFT-D3色散修正，对含氢体系尤为重要 ISIF = 3 # 优化晶胞形状和体积（对于完美超胞测试）或=2（固定晶胞优化原子位置，对于缺陷模型） IBRION = 2 # 使用共轭梯度法进行离子弛豫 NSW = 200 # 最大离子弛豫步数 EDIFF = 1E-6 # 电子步收敛标准 EDIFFG = -0.01 # 离子步收敛标准（当所有力 < 0.01 eV/A时停止） LORBIT = 11 # 输出详细的态密度和局域电荷信息，便于后续分析

对于缺陷计算的两个特殊设置：

LCALCPOL = .TRUE.：计算Berry相位电极化。对于非共价键结合的离子晶体中的带电缺陷，这有助于修正由于周期性边界条件引起的虚假静电相互作用，但计算量较大。在初步优化时可以不开启，在最终单点能计算时开启。
LSUBROT = .TRUE.和LCHARG = .TRUE.：在优化带电缺陷时，有时需要固定缺陷周围原子的某些对称性以避免结构弛豫到不合理的局域极小值，并保存电荷密度用于分析。

注意事项：在测试K点网格时，对于超胞，经常使用Gamma-centered的网格（KGAMMA = .TRUE.）。因为超胞的倒格子原胞很小，有时仅用Γ点（1x1x1）采样就是够的。但这必须通过测试总能量是否收敛来验证。盲目使用高密度K点网格只会浪费计算资源。

4. 实操过程与核心环节实现

让我们以构建一个2×2×1的β-Ga₂O₃超胞，并测试计算参数为例，走一遍PREPARE02的核心流程。

4.1 步骤一：获取并优化完美原胞结构

首先，我们需要一个可靠的β-Ga₂O₃原胞结构作为起点。

获取初始结构：从Materials Project、COD晶体学数据库或已发表的高质量理论文章中获取β-Ga₂O₃的晶体结构文件（通常是CIF格式）。确保空间群信息正确（β相通常为C2/m）。
原胞优化：使用VASP对该原胞进行充分的几何优化（ISIF=3，允许晶胞形状和体积变化）。优化后，将计算得到的晶格常数（a, b, c, β角）与实验值或公认的理论值对比。误差应在1-2%以内。这一步的CONTCAR文件就是我们后续所有工作的“标准晶体”。

4.2 步骤二：构建超胞并测试截断能（ENCUT）

构建超胞：使用VASP自带的makeposcar脚本（或其他工具如pymatgen、ASE）将优化好的原胞沿三个晶格矢量方向扩胞。例如，执行makeposcar -s 2 2 1来生成2×2×1超胞的POSCAR文件。
ENCUT收敛性测试：
- 创建一系列计算文件夹，如encut_400,encut_450,encut_500,encut_550,encut_600。
- 在每个文件夹中，放置相同的POSCAR（超胞）、POTCAR（赝势）和KPOINTS文件（先使用一个较密的网格，如3x3x3，以确保K点误差不干扰ENCUT测试）。
- 修改INCAR文件，仅设置SYSTEM,PREC=Accurate,ISMEAR,SIGMA,EDIFF=1E-6,NSW=0（只做单点能计算，不优化），以及不同的ENCUT值。
- 提交这些计算任务。完成后，提取每个任务的OSZICAR或OUTCAR文件中的总能量（E0）。
- 绘制总能量（单位：eV/atom）随ENCUT变化的曲线。通常能量会随着ENCUT增大而降低并逐渐趋于平缓。我们选择能量变化进入平台区（如变化< 1 meV/atom）时的最小ENCUT值。假设在520 eV后能量基本不变，我们就选定ENCUT = 520。

4.3 步骤三：测试K点网格密度

在确定了ENCUT后，接下来测试K点网格。

创建一系列文件夹：kpoint_111,kpoint_222,kpoint_333,kpoint_221等。这里的数字代表在倒易空间各方向上的采样点数。
使用上一步确定的ENCUT（520 eV），NSW=0，进行单点能计算。
绘制总能量随K点密度变化的曲线。对于2×2×1的超胞，其倒易空间原胞比原胞小，因此可能不需要很密的K点。可能发现2x2x2网格的能量与3x3x3网格的能量相差无几（< 1 meV/atom）。那么我们就可以选择2x2x2网格作为后续计算的KPOINTS。

KPOINTS文件示例 (2x2x2 Monkhorst-Pack网格)：

Automatic mesh 0 # 0 -> 自动生成网格 Monkhorst-Pack # 网格类型 2 2 2 # 三个方向上的格点数 0 0 0 # 偏移量

4.4 步骤四：优化完美超胞结构

现在，我们有了优化的ENCUT和K点，可以对我们构建的2×2×1完美超胞进行完整的几何优化了。

创建一个新文件夹perfect_supercell_opt。
拷贝POSCAR（2x2x1超胞）、POTCAR、以及测试好的KPOINTS文件（2x2x2）进去。
编写INCAR文件，使用测试好的参数，并设置ISIF=2（固定超胞大小，只优化原子位置，因为超胞尺寸已经由原胞优化决定）、IBRION=2、NSW=200、EDIFFG=-0.01。
提交优化计算。优化完成后，检查CONTCAR中的原子位置是否合理（没有原子异常靠近），并记录下最终的总能量。这个优化后的超胞结构（CONTCAR）和对应的总能量，将成为后续所有缺陷形成能计算的参考基准（即“完美晶体”的能量E_{perfect}）。

5. 常见问题与排查技巧实录

即使按照上述流程操作，在实际计算中仍会遇到各种问题。以下是一些典型问题及排查思路。

5.1 计算不收敛或收敛极慢

问题现象：电子自洽循环（SCF）迭代次数超过默认限制（通常60步）仍未收敛，或者离子弛豫震荡无法达到受力标准。
排查与解决：
1. 检查初始结构：用可视化软件检查POSCAR，确认原子没有重叠（距离过近），氢原子的初始位置是否合理（例如，不要放在离氧原子0.5 Å的地方）。
2. 调整混合参数：在INCAR中，ALGO = All或ALGO = Normal是默认的阻塞Davidson算法。对于难收敛的体系，可以尝试ALGO = Damped（阻尼动力学）结合较小的时间步长TIME = 0.1，或者使用ALGO = Conjugate（共轭梯度）。对于带隙较大的绝缘体如Ga₂O₃，有时ALGO = Normal配合合理的AMIX和BMIX参数效果更好。
3. 检查赝势：确认所有元素的赝势来自同一库且版本兼容。不匹配的赝势是导致收敛问题的常见元凶。
4. 分步优化：对于复杂的缺陷模型，可以先使用较松的收敛标准（EDIFF=1E-5, EDIFFG=-0.05）和较低的ENCUT进行预优化，得到一个相对合理的中间结构，再用严格的参数进行精细优化。

5.2 缺陷优化后结构“跑飞”

问题现象：优化完成后，查看CONTCAR发现氢原子跑到了晶胞边界，或者整个结构发生了意想不到的大畸变。
排查与解决：
1. 对称性破缺：缺陷本身降低了原晶格的对称性。在优化带电缺陷时，初始对称性设置过高可能导致离子弛豫陷入不稳定的路径。可以尝试在创建缺陷时手动将对称性降低（如使用pymatgen的DefectGenerator并设置symprec=1e-3），或者直接在VASP的POSCAR中使用笛卡尔坐标并关闭对称性检测（ISYM=0）。
2. 初始位置不佳：氢原子的初始位置可能位于势能面的鞍点或高点，导致弛豫时滑向远处。尝试多个不同的初始位置（例如，间隙位的中心、偏向某个氧原子、偏向某个镓原子）分别进行优化，比较最终的能量，取能量最低的稳定结构。
3. 检查电荷态：确认你为缺陷模型设置的净电荷（NELECT或在INCAR中通过NUPDOWN间接设置）是否正确。一个错误的电荷态会导致极强的库仑力，将原子“推开”。

5.3 缺陷形成能异常高或为负值

问题现象：按照形成能公式计算出的结果不符合物理预期（例如，在一切条件下形成能都为负，意味着该缺陷会自发无限生成；或者形成能高达几十eV，几乎不可能形成）。
排查与解决：
1. 参考能量错误：确保用于计算的“完美超胞”总能量E_{perfect}是在完全相同的计算参数（ENCUT, K点，赝势，ISIF等）下得到的。任何微小的参数差异都会直接带入形成能误差。
2. 化学势取值：形成能公式依赖于原子化学势 μ_Ga, μ_O, μ_H。这些化学势的取值范围由Ga₂O₃的生成焓和外界条件（富Ga/富O）限定。必须仔细检查你使用的化学势值是否在热力学允许的区间内。μ_H 的参考态通常是 H₂ 分子，需要计算一个孤立的H₂分子在足够大盒子中的能量作为基准。
3. 静电修正遗漏：对于带电缺陷，在周期性边界条件下计算的总能量包含了与超胞尺寸相关的虚假静电相互作用能。必须对此进行修正（如使用LCALCPOL相关的修正，或更常用的，使用专门的后处理工具如pycdt,DefectFormation或VASPKIT中的功能进行 Makov-Payne 或 Freysoldt 修正）。未修正的带电缺陷能量通常是错误的。
4. 结构未完全优化：比较的缺陷结构和完美结构都必须充分弛豫到基态。如果缺陷结构被困在亚稳态，其能量偏高，会导致形成能计算偏高。务必检查优化是否真正收敛（EDIFFG达标，且最后几步能量/受力变化很小）。

实操心得：我强烈建议建立一个计算日志表格，记录每一个缺陷模型的计算信息。表格列可以包括：缺陷名称、电荷态、初始结构描述、使用的INCAR关键参数、计算文件夹路径、最终总能量、是否收敛、备注（遇到的问题）。这不仅能避免混乱，当需要复查或写论文时，这些记录就是无价的宝藏。对于氢掺杂计算，由于氢原子质量轻、迁移势垒可能较低，在优化时可以考虑使用POTIM = 0.1（较小的离子移动步长）来获得更平滑的收敛路径。

查看全文

http://www.jsqmd.com/news/869163/