《physica status solidi (b)–– basic solid state physics》：Alternative Structures for the Si(111) 7?×?7 Surface: Evidence for a Novel Class of Distorted Honeycomb Reconstructions
编辑推荐：
利用密度泛函理论（DFT）计算创建并优化了多种Si(111) 7×7表面的新型原子模型，并针对多种实验测量结果进行评估，以理解被广泛接受的二聚体-吸附原子-层错（DAS）重构模型中发现的悖论。这些模型包括配位数偏移的吸附层、不同的终止层、额外原子、缺失二聚体或
利用密度泛函理论（DFT）计算创建并优化了多种Si(111) 7×7表面的新型原子模型，并针对多种实验测量结果进行评估，以理解被广泛接受的二聚体-吸附原子-层错（DAS）重构模型中发现的悖论。这些模型包括配位数偏移的吸附层、不同的终止层、额外原子、缺失二聚体或畸变蜂窝（HC）及硅烯样结构中的哑铃原子。研究人员利用DFT及近似的准粒子DFT-1/2方案，在考虑和不考虑自旋极化的条件下评估了这些模型的几何与电子结构，结果表明若干新型同素异形体在多种实验特征上与实验的吻合度优于DAS模型。计算得到的形成能用于构建以硅化学势为函数的相图，该相图在体相原子储库条件下倾向于DAS模型，而在修正的衬底化学势条件下其他几何构型更有利。研究发现大多数及/或少数非磁性同素异形体的磁性相存在表面带隙。与实验总体符合最佳的是一个磁性相结构，该结构不同于DAS，具有显著横向畸变的HC结构，且在晶胞两侧分别存在高翘曲和低翘曲的吸附原子。文中讨论了该结构的本质及驱动这类新型重构的相互作用机制。
论文解读文章
一、研究背景与问题
硅基电子器件是现代电子工业的基石，而Si(111) 7×7表面重构作为半导体表面科学中最经典、研究最深入的体系，其原子结构模型一直被视为表面物理的典范。自1959年该表面重构被发现以来，经过长达26年的探索，科学界最终形成了以二聚体-吸附原子-层错（dimer-adatom-stacking-fault，DAS）模型为核心的共识结构。该模型通过透射电子衍射（TED）分析、扫描隧道显微镜（STM）成像和低能电子衍射（LEED）测量逐步确立，并由后续的第一性原理DFT计算所证实。
然而，近三十年来积累的大量实验数据与基于DAS模型的理论预测之间出现了系统性矛盾。这些悖论主要包括六个方面：第一，实验在低温下观测到约0.1 eV的能隙，而DAS模型预测该表面为金属性（半满带）；第二，实验测量的吸附原子高度比DFT预测值高0.2–0.4 ?，远超DFT的公认精度范围；第三，实验中存在第四个“被遗忘”的表面态，而DFT对DAS的计算未预测到该态；第四，该第四表面态及其特征在实验上表现出自旋极化特征；第五，与其他表面态不同，该态在表面平面内呈极化特征；第六，低温下观测到侧向位移的吸附原子，而DAS模型不存在该特征。这些悖论促使研究人员质疑DAS模型是否真正描述了Si(111) 7×7表面的真实原子结构，以及当前的计算方法是否足以准确描述这一复杂体系。
二、研究内容与结论
针对上述问题，研究人员系统构建并评估了多种替代的7×7表面重构模型，包括配位数偏移的吸附层模型（MBA、MHA系列）、不同终止层模型（TFA）、缺失二聚体或哑铃原子模型（MAS）、以及一类被称为非对称畸变蜂窝（asymmetric distorted honeycomb，ADH）家族的新型重构模型。研究人员从超过1500次VASP计算、约50种不同模型结构中筛选出最具代表性的若干同素异形体进行深入分析。
研究得出以下核心结论：首先，在体相硅化学势条件下，DAS模型仍是最稳定的重构结构，但这一优势仅在体相原子储库假设下成立；当考虑修正的衬底化学势时，其他几何构型（尤其是ADH系列）可以变得更加有利。其次，通过DFT-1/2准粒子方法计算发现，多个同素异形体在磁性相下呈现绝缘性，而DAS结构始终为金属性，这与实验观测到的低温能隙一致。第三，与实验各项测量指标（包括X射线衍射图样、STM形貌、表面态能量位置等）综合比较后，研究人员发现一个磁性ADH相与实验的定性吻合度最佳，该结构在晶胞两侧具有显著不同的横向畸变蜂窝结构，并分别存在高翘曲和低翘曲的吸附原子。这一发现表明，7×7表面的真实结构可能并非传统认为的DAS模型，而是一类具有强对称性破缺特征的新型重构。
三、关键技术方法
研究人员主要采用了以下关键技术方法：基于投影缀加波（PAW）方法和平面波截断能250 eV的密度泛函理论（DFT）计算，使用PBEsol交换关联泛函进行结构优化和总能计算；采用DFT-1/2近似准粒子方案修正带隙计算；构建双双层和四双层Si(111)衬底模型，底部以氢原子钝化；利用形成能相图分析不同化学势条件下的结构稳定性；计算二维结构因子（SF）和帕特森函数（PM）与实验衍射数据对比验证结构模型。
四、研究结果
4.1 DAS结构与新型同素异形体的构建
研究人员首先回顾了DAS模型的基本结构特征，包括吸附原子（AAs）、剩余原子（RAs）、角洞（CHs）、二聚体（Ds）和哑铃（Db）原子等组成要素。通过移除DAS晶胞中特定位置的Db原子，研究人员发现晶胞两侧表现出截然不同的弛豫行为：一侧形成倾斜翘曲的蜂窝结构，另一侧吸附原子被显著抬升并产生约10°的倾斜，这种强对称性破缺结构成为ADH家族的核心特征。
4.2 各同素异形体的结构特征
MHA系列模型考察了无二聚体但有角洞的层错结构，弛豫后左侧蜂窝发生翻转形成更强共价键，同时重新分布了剩余原子位置。DFA模型是基于π键合的“互锁层错吸附原子”模型，弛豫后蜂窝结构严重畸变并出现sp2杂化特征。TFA模型基于三重悬挂键表面构建，弛豫后形成具有“过渡层”的压缩四聚体结构。ADH系列模型则源自早期“对称化三聚体模型”，弛豫后形成不完全的π键合链结构，其中ADH3模型在F侧添加Db原子后获得了定义良好的7×7结构。
4.3 形成能与相图分析
形成能计算表明，在体相化学势条件下DAS模型最稳定，但MHA3和ADH5等模型在弱键合原子储库条件下可以变得更有利。相图分析显示，对于?5.41 eV至约?4.28 eV的化学势范围，DAS保持最稳定；而对于更强键合的储库条件，ADH、ADH2和ADH4模型趋于有利。TFA?模型在考虑过渡层原子后，其形成能显著降低，可在特定化学势条件下优于DAS。
4.4 衍射图样与结构验证
X射线衍射图样对比显示，仅DAS、ADH3和MAS三种同素异形体能够复现实验衍射图样的主要“雪花状”特征。最新的X射线测量显示晶胞两侧的(0 3/7)反射强度存在16%的差异，而DAS模型仅预测0.5%的不对称性，ADH类结构则能更好地解释这种实验观测到的对称性破缺。
4.5 电子结构与磁性
DFT-1/2准粒子计算表明，多个ADH同素异形体在磁性相下呈现绝缘性，带隙约为0.1 eV量级，与实验低温测量一致。自旋极化主要分布在吸附原子上，且该磁性相的表面态具有面内极化特征，与实验观测相符。
五、讨论与结论
讨论部分指出，传统衍射分析中采用的六重对称平均会掩盖晶胞两侧的结构差异，导致DAS模型被错误地认为是唯一可能的解。投影杆（projection rod）建模方法对强畸变结构的适用性也存在局限。研究人员强调，7×7表面的真实结构可能是一类具有强对称性破缺的新型畸变蜂窝重构，而非传统认为的DAS模型。
研究结论指出，磁性ADH相与实验各项指标的综合吻合度最佳，该结构具有显著横向畸变的蜂窝结构和晶胞两侧不对称的吸附原子构型。这一发现不仅对理解Si(111) 7×7表面的真实原子结构具有重要意义，也为其他复杂表面和二维材料体系的结构确定提供了方法论启示：必须综合考虑多种实验测量手段，而非仅依赖单一衍射技术。未来研究建议包括更高精度的GW计算、低温高分辨测量以及更系统的结构搜索方法。
