Insights on molecular P implantation for scalable spin-qubit arrays

技术摘要：关于用于可扩展自旋量子比特阵列的分子 P 注入研究之见解

问题陈述
利用硅衬底中磷施主（$^{31}$P）自旋的固态自旋量子比特是量子计算的主要候选方案之一。尽管自旋相干性已得到显著改善，但一个关键挑战仍然存在：即如何精确且可扩展地在硅晶格中放置施主原子。对于优化的量子比特操作，施主必须被定位在距离表面 7–20 nm 的深度。目前的离子束注入技术虽然具备可扩展性，但其位置不确定性比扫描隧道显微镜（STM）刻蚀技术高出数个数量级。

为了缓解这一问题，近期的提议建议使用分子离子（具体为 $PF_2$）而非单原子磷离子。其假设是，分子离子可以产生用于检测的高电子信号，同时保持与 9 keV 单原子 P 离子相当的放置不确定性。然而，现有的二元碰撞近似（BCA）方法（该方法假设分子在撞击瞬间立即破碎）无法捕捉多体效应、分子在穿透过程中的完整性、电子信号的具体性质，以及分子注入与单原子注入所导致的特定辐射损伤剖面。

方法论
作者采用分子动力学（MD）模拟，利用大规模原子/分子并行模拟器（LAMPS），研究了单原子与分子（$PF_2$）离子对硅进行辐照的差异。

• 模拟设置： 使用大型模拟单元（$43.5 \times 43.5 \times 65.7$ nm$^3$）及结晶硅（c-Si）结构来追踪穿透动力学，包括沟道效应。一个带有 5 nm 非晶二氧化硅（$a$-$SiO_2$）层的较小单元用于模拟晶圆上存在的天然氧化层。
• 势函数： Si-Si 相互作用采用 Stillinger-Weber (SW) 势（因其计算效率高且准确度与机器学习势相当），ZBL 势用于处理离子-晶格间的相互作用，并使用 Lennard-Jones 势处理分子内相互作用。
• 条件： 模拟在 300 K 下进行。$PF_2$ 分子被分配总能量为 20.041 keV，以确保其组成部分的 P 原子能量为 9 keV。入射角设定为 7 度，以最小化沟道效应概率；同时通过改变 $a$-$SiO_2$ 层是否存在，来评估其对能量损失和分子完整性的影响。
• 分析： 研究追踪了能量损失、电子阻止本领以及缺陷形成（使用 Wigner-Seitz 分析）和缺陷团簇。电子阻止本领被建模以模拟用于检测的电子-空穴对生成。

关键结果

1. 分子完整性与能量损失：
   与 BCA 假设分子在表面立即破碎相反，MD 模拟显示，$PF_2$ 分子在进入衬底后通常保持完整或部分完整（表现为 $PF$ 复合物）数纳米。在纯 c-Si 中，超过 40% 的案例显示在一个氟原子在 6 nm 深度前脱离，而第二个氟原子仍紧随磷原子附近。存在 $a$-$SiO_2$ 层会增加碰撞频率，导致更早的解离。

   • 能量分布： 当穿过 $a$-$SiO_2$ 层时，P 离子的能量损失比 F 离子多约 1 keV。
   • 注入剖面： 与单原子注入相比，分子注入产生的深度剖面更平坦，且 P 原子被推向更深的深度。这归因于分子在解离前作为一个单一的、更重的投射体发挥作用。值得注意的是，分子注入产生了更强的“长尾”沟道离子，这可能是由于 P 原子在穿过氧化层后进入 c-Si 时的有效能量更高。

2. 电子信号与沟道检测：
   研究分析了电子阻止信号（用于检测注入深度）与 P 原子最终位置之间的相关性。

   • 信号相关性： 在电子信号强度与 P 原子单独的穿透深度之间发现弱相关性（皮尔逊系数 $r \approx 0.5$）。然而，该信号与所有三个分子组分（P + 2F）的平均穿透深度强相关。
   • 沟道歧义性： 强电子信号并不一定意味着 P 原子发生了沟道效应，因为信号可以由 $PF_2$ 分子中的任何一个原子产生。相反，作者指出，低信号意味着任何粒子均未发生沟道效应。因此，虽然强信号具有歧义性（可能源于 F 原子发生沟道效应而 P 原子位置正确，反之亦然），但低信号能可靠地表明没有任何粒子发生了沟道效应。
   • 检测极限： 研究发现，使用信号阈值来区分沟道事件会导致显著的不确定性；大约一半的事件无法仅根据信号强度被正确标记为沟道或非沟道事件。然而，具有强信号的事件可以被安全地丢弃以规避沟道风险，即便这可能会丢弃正确的 P 原子注入。

3. 辐射损伤与缺陷演化：

   • 协同损伤： 与三个独立单原子离子的损伤总和相比，分子注入产生了更多的缺陷和更宽的缺陷计数分布。这表明存在一种“协同损伤效应”，即分子中紧密排列的原子产生的重叠碰撞级联创造了更广泛的损伤。
   • 团簇化： 分子注入中的缺陷倾向于形成较大的团簇（部分超过 500 个原子），这些团簇以 P 原子为中心。这些大团簇通常是富含空位的非晶口袋。
   • 氧化层的作用： $a$-$SiO_2$ 层降低了原子的进入能量，从而限制了最大缺陷团簇的尺寸，可能有助于减少晶格恢复所需的退火处理。
   • 退火启示： 大规模、富含空位的团簇以及特定的缺陷对（如 E-center）的存在，对电学激活和量子特性构成了风险。这些缺陷需要高温退火才能恢复，但分子注入的特定损伤剖面表明，标准的单原子注入退火协议可能需要调整，以抑制有害的原子扩散。

意义与主张
本文旨在提供一份基于分子动力学的综合评估，挑战了当前基于 BCA 模型及实验解释中的若干假设。

• 驳斥立即解离论： 研究证明，假设分子在表面立即破碎是不准确的；分子通常在解离前会穿透几纳米，从而导致重叠的级联效应和协同损伤。
• 信号解释： 作者警告不要将高电子信号直接与 P 原子沟道效应联系起来。他们认为，在分子注入中，高信号可能源于任何组成部分，使其成为判断特定 P 施主深度或沟道状态的不可靠指标。相反，低信号是所有粒子均未发生沟道的可靠指标。
• 损伤剖面： 研究强调，分子注入创造了一个独特的损伤景观，其特征是比单原子注入具有更大、更集中的缺陷团簇。这要求重新评估后注入退火策略，以确保施主实现电学激活，同时避免残留应变或缺陷破坏量子比特性能。

总之，虽然分子离子为可扩展量子比特生产提供了路径，但研究结论指出，其使用引入了在放置检测和损伤恢复方面必须仔细管理的复杂性，特别是在解释电子信号和优化热退火工艺方面。