序列淬火用于基于形成能与迁移能预测半导体缺陷浓度：以CdTe:As掺杂为例

《Journal of Materials Chemistry C》：Sequential quenching to predict semiconductor defect concentrations from formation & migration energies: the case of CdTe:As doping

　　

半导体中的缺陷浓度受生长与降温过程中的热历史显著影响，然而多数缺陷计算通常假定两种极限情形之一：要么从高温瞬时淬火，要么在整个过程中始终维持完全平衡——而这两种情形在现实中都很少真正达到。该研究提出序列淬火（sequential quenching, SQ）作为第3类缺陷计算方法，利用缺陷形成能与迁移能，对以有限冷却速率降温样品中受扩散限制动力学（diffusion-limited kinetics）支配的缺陷浓度进行建模。在SQ中，每一种缺陷的浓度会在由其扩散速率、至源/汇的距离以及冷却速率共同决定的特征温度处被冻结。由于不同电荷态通过电荷中性条件相互耦合，但又在不同温度下冻结，因此冻结事件的先后序列不满足交换性。关键的是，并非所有室温下的SQ解都能够由完全平衡（equilibrium, EQ）或完全淬火（full quenching, FQ）计算推导得到——若不采用SQ，极有可能得到错误预测。研究人员以As掺杂CdTe为例说明SQ方法；实验数据表明，该体系中的掺杂行为会随冷却速率变化，并且光伏用多晶薄膜与块体晶体之间存在差异。SQ计算揭示，诸如Cd间隙原子（Cd-interstitials, Cd_i）等快速扩散缺陷可在更低温度下保持迁移能力，并在更大的特征距离上发生冻结，从而在快速冷却或块体样品中导致强补偿并呈现n型行为。较慢冷却以及更小特征距离会抑制施主冻结并增强p型激活。这些结果确立了SQ作为一种物理图像清晰且计算高效的分析框架，可用于将冷却条件、样品几何形貌和缺陷动力学与CdTe及相关材料中的掺杂激活联系起来。

该论文发表于《Journal of Materials Chemistry C》，围绕半导体中缺陷浓度对热历史高度敏感这一核心问题展开。点缺陷及缺陷复合体决定了半导体的载流子浓度、补偿效应与输运行为，而现实材料在生长、退火和冷却过程中既受热力学约束，也受扩散动力学、样品尺寸以及位错、晶界、表面和界面等源/汇结构影响。传统缺陷分析多采用两类极限假设：一是完全平衡（equilibrium, EQ），即冷却无限缓慢、各类缺陷始终维持平衡；二是完全淬火（full quenching, FQ），即高温缺陷态在瞬时降温后整体冻结。这两类方法虽便于计算，却难以描述真实加工条件下有限冷却速率、有限扩散长度和多种缺陷依次冻结的非平衡过程。尤其对于CdTe这类技术上重要的半导体，V族掺杂受本征缺陷和复合缺陷强烈补偿，在块体晶体与多晶薄膜中的激活行为差异明显，仅依靠传统EQ或FQ框架难以给出一致解释，因此有必要建立一种兼顾缺陷热力学与扩散动力学的新方法。

研究人员据此提出序列淬火（sequential quenching, SQ）框架，并将其集成到KROGER缺陷浓度计算体系中。该方法的核心思想是：在有限冷却速率γ下，不同缺陷或不同电荷态因迁移势垒不同，其能够继续与表面、晶界、位错等源/汇交换的温度范围并不相同；当其在后续降温过程中的均方根扩散长度不足以跨越特征距离L时，该缺陷总浓度即在对应冻结温度T_freeze处锁定。与此同时，体系中的带电缺陷仍通过全局电荷中性条件彼此耦合，因此某一种缺陷冻结后，将反过来约束费米能级与其他缺陷浓度的后续演化。这意味着不同缺陷的冻结顺序具有非交换性，室温缺陷组合并不一定包含于EQ或FQ两类极限解之中。研究以As掺杂CdTe为案例，重点考察Cd富与Te富边界条件、固定As总浓度与第二相限制μ_As两类掺杂约束，以及冷却速率和特征距离对受主激活与施主补偿的联合作用。

研究人员使用的主要技术方法包括：基于KROGER平台进行序列淬火（SQ）、完全平衡（EQ）和完全淬火（FQ）缺陷浓度计算；采用全局电荷中性与元素守恒约束，通过数值求解确定费米能级与各缺陷浓度；利用Arrhenius扩散模型和连续冷却积分扩散长度判据定义T_freeze；采用第一性原理密度泛函理论（density functional theory, DFT）计算CdTe中缺陷形成能、跃迁能级和迁移势垒，其中使用投影缀加波（projector augmented-wave, PAW）方法、HSE杂化泛函并包含自旋–轨道耦合（spin–orbit coupling, SOC）；迁移势垒通过攀爬图像弹性带（climbing nudged elastic band, cNEB）方法获得。文中未涉及生物学样本队列。

在研究结果部分，作者首先比较了SQ与EQ在Cd平衡条件下的温度依赖缺陷演化。结果表明，在EQ计算中，缺陷浓度可随温度持续重排，甚至在低温区域仍表现出非单调变化，例如As_Te-Cd_i复合缺陷在一定温区内可上升，这本质上源于平衡模型默认长程扩散始终足以维持重构。相比之下，在SQ框架下，高迁移势垒缺陷会在较高温度率先冻结，而更高迁移率物种如Cd_i则可继续演化至更低温度。由此，低温下某些缺陷浓度的变化并非意味着所有缺陷仍在平衡，而是移动性较高的缺陷通过改变电荷平衡和化学势，间接驱动开放缺陷种类的重新分配。这一冻结层级结构正是SQ区别于EQ和FQ的核心。

随后，作者分析了固定[As]与第二相Cd₃As₂限制μ_As两种掺杂情形下的差异。在固定[As]情形中，SQ在室温下预测的空穴浓度远高于EQ结果，说明顺序冻结能够显著改变最终电学活化水平。该差异主要源于Cd_i施主浓度对冻结过程极其敏感，而As_Te-Cd_i复合体在较高温度冻结后，会影响后续费米能级及补偿行为。在第二相限制条件下，SQ虽然同样会使As相关缺陷在较高温度被冻结，但由于热力学上As固溶度受限，仅靠动力学冻结并不足以在所有边界条件下获得p型行为。这表明，动力学效应虽关键，但并不能脱离热力学边界独立决定掺杂结果。

在径向分布与特征距离分析中，作者给出了Cd富和Te富条件下、不同冷却速率对应的SQ预测。结果显示，更快冷却会使缺陷在更小特征距离上冻结；而更大的特征距离则要求更长扩散路径，从而促进更早冻结。在Cd富条件下，Cd_i是主导补偿的关键施主。快速冷却且特征距离较大时，高浓度Cd_i在冷却过程中被保留下来，导致强施主补偿乃至n型行为；较慢冷却可使Cd_i保持迁移到更低温，从而减少冻结保留量并减轻补偿。与之相对，在Te富条件下，补偿主要由As_Cd反位缺陷主导；当特征距离较大时，其浓度较低，补偿较弱，而在较小特征距离时As_Cd浓度明显增加，可与As_Te发生强烈相互补偿并降低空穴浓度。该结果说明，主导补偿缺陷的种类会随热化学边界条件而发生转换。

作者进一步构建了冷却后去除可动Cd_i的理想化情形，用以估计掺杂激活的上限。在Cd富条件下，若低温阶段移动Cd_i能够扩散至表面或晶界并被去除，则补偿显著降低，空穴浓度可接近As_Te浓度，激活率接近约98%；但仍有冻结在体内的As_Te-Cd_i复合体保留一定补偿作用。在Te富条件下，去除Cd_i前后结果几乎不变，因为此时主导补偿并非Cd_i而是As_Cd。这部分结果强调了移动间隙施主在后处理激活中的关键作用，同时指出残余复合缺陷可能构成激活上限。

在As激活率的系统比较中，作者总结了不同热力学和动力学场景下的趋势。对于固定[As]情形，Cd富条件下SQ给出的最大激活率约为30%，且激活率随特征距离和冷却速率系统变化；这一趋势与Nagaoka等人的块体实验结果在定性上相符。Te富条件则整体表现出更高激活率，在较大特征距离和中等掺杂浓度下可接近完全激活，但在高浓度掺杂或小特征距离下仍会下降，表明补偿缺陷在高掺杂条件下重新稳定。对于第二相限制μ_As的情形，Te富条件下SQ可预测接近完全激活，而Cd富条件下激活仍非常低，说明第二相约束带来的热力学极限可压制单纯的动力学获益。总体而言，Cd_i在Cd富条件下是控制补偿和激活的中心因素。

在讨论部分，作者进一步指出，CdTe多晶薄膜中的晶界和接触耗尽区可能形成内建电场，影响带正电的Cd_i²⁺长程迁移，使其难以简单地向晶界和界面完全逸出；光照诱导的能带弯曲减弱可能允许一部分自由Cd_i重新分布，但As_Te-Cd_i复合体具有更高结合能，单靠常见静电势垒不足以使其解离。因此，生长后残余的移动Cd_i可能为后续亚稳态演化提供路径，并与实验中观察到的“唤醒”现象相关。作者同时指出，文中未纳入显式复合反应动力学、团簇沉淀以及晶界专属扩散势垒，这些过程仍有待未来工作补充，但这并不影响SQ在连接冷却条件、样品几何和缺陷动力学方面的框架价值。

论文结论部分可译为：本研究建立了用于描述有限冷却速率下缺陷演化的序列淬火（SQ）框架，并将其应用于As掺杂CdTe。通过引入缺陷特异的迁移势垒、冷却速率和至有效源/汇的特征距离，SQ提供了一个具有物理依据的非平衡缺陷布居描述，超越了完全平衡与完全淬火模型的极限假设。该方法能够在实验相关的加工条件下预测空间变化的缺陷与载流子浓度。SQ对CdTe的应用表明，动力学限制在决定掺杂激活方面起核心作用。快速扩散的施主缺陷，尤其是Cd_i，可在低温下仍保持迁移能力，并在较大特征距离下冻结，从而在Cd富平衡条件下导致强补偿和n型行为；相反，较慢冷却与较小特征距离会限制施主保留并促进更高的p型激活。这些结果说明，CdTe中的掺杂激活不仅受缺陷形成能控制，也取决于冷却的时间–温度轨迹以及冷却过程中移动缺陷可扩散的长度尺度。研究还表明，生长后移动Cd_i的再分布或去除决定了可实现激活的上限，而残余缺陷复合体则构成基本补偿极限。在这一背景下，生长后存在的移动间隙施主为亚稳行为提供了自然通道，在光照等引发静电条件变化时，Cd_i可能发生再分布而无需改变缺陷形成过程，这类过程可能有助于解释CdTe中实验观察到的“唤醒”现象。总体而言，这些结果确立了序列淬火作为一种实用且具有预测能力的框架，可用于解释非平衡缺陷布居及其在变化热环境与静电环境下的演化，并为CdTe及相关半导体材料中冷却条件、样品几何、缺陷动力学与掺杂激活、亚稳性之间的联系提供一般性方法。