第一性原理研究Ag(111)上硼烯的结构稳定性与STM成像

摘要

近期实验证明分子束外延法可在Ag(111)上生长原子级薄硼片——硼烯，观测到两种不同表面相，但其精确原子排列和热力学稳定性仍有争议。本研究采用密度泛函理论考察Ag(111)上最可能的硼烯多晶型——屈曲三角、β12和χ3。计算显示三种独立单层均热力学不稳定且呈金属性，而Ag(111)衬底可将其稳定。模拟扫描隧道显微镜（STM）图像复现了实验观测的条纹相和均匀相，明确识别了底层晶格结构。

背景

自石墨烯发现以来，二维（2D）材料因其独特的电子、力学和催化性能备受关注[1-7]。过去十年，多项理论预测了各种2D硼片（统称硼烯），呈现多样的六方空位排列和屈曲模式[8-21]。大量第一性原理研究表明，带周期性六方空位的平面三角硼晶格在能量上最受青睐，但2D形式仍高于块体硼能量，意味着存在向3D结晶的热力学驱动力[22-24]。因此，需要足够“粘性”的衬底来抑制3D成核，促进原子层形成[25-27]。

近期实验突破在金属表面实现了硼烯。Mannix等在Ag(111)沉积硼，利用高分辨STM观测到条纹相和均匀相两种表面相[28]。Feng等报道类似相，将均匀结构指认为χ3晶格，条纹相解释为屈曲三角或矩形β12晶格[29-31]。这些分歧凸显了系统理论研究Ag(111)上可能原子构型及其稳定性的必要性。

除结构鉴定外，硼烯多晶型展现卓越性能。屈曲三角片预测为高各向异性金属，臂椅方向杨氏模量超越石墨烯，晶格热导率强各向异性[32,33]。β12和χ3硼烯预测超导性，潜力可比肩MgB₂薄膜[34-36]。但β12和χ3热力学稳定性存争议，有研究报告Γ点附近存在虚频[37,38]。

计算方法

所有计算采用VASP包进行密度泛函理论（DFT）计算[38]。投影增强波法描述电子-离子相互作用[39]，采用PBE泛函的广义梯度近似（GGA）处理交换相关[40]。平面波截断能500 eV，布里渊区采用25×15×1、15×9×1、11×11×1 k点网格分别对应屈曲三角、β12、χ3相。z向至少20 Å真空层消除周期像相互作用。几何优化收敛至力<0.02 eV Å⁻¹，底层两层Ag固定。PHONOPY有限位移法计算声子色散[41]。

STM图像在Tersoff-Hamann框架下模拟，隧道电流与样品局域态密度（LDOS）相关[42]。采用恒流模式在不同偏压下生成图像，直接对比实验STM数据。

结果与讨论

图1展示三种硼烯多晶型优化几何：屈曲三角(a)、β12(b)、χ3(c)。屈曲三角沿b轴单屈曲，β12和χ3为平面但六方空位排列不同。晶格常数——屈曲三角1.613 Å×2.866 Å，β12 2.916 Å×5.075 Å，χ3 4.448 Å——与前人理论和实验吻合[23,25,27]。

空位浓度η：β12为1/6，χ3为1/5，表明后者缺陷有序度更高。通过平均每原子结合能E_FB=E_boron/n评估相对稳定性（表2），独立β12片最稳定，χ3最不稳定，差异0.08 eV/原子。但声子色散计算（图2）显示三者Γ点附近均有虚频，确认独立状态下热力学不稳定。

电子能带结构（图3）确认三者均为金属。屈曲三角表现强各向异性，S-Y和G-X方向带隙分别为9.63 eV和4.32 eV，暗示导电性局限于未屈曲的a轴。

置于Ag(111)上，β12晶格失配极小（~1%），屈曲三角中等（~3%）。χ3片在衬底上采取两种不同排列（χ3和χ3'），晶格常数分别为8.67 Å和2.89 Å×25.02 Å。硼烯层与Ag表面垂直距离2.4-3.3 Å，表明弱物理吸附，保持β12、χ3、χ3'平面特征，与实验报告的Ag(111)上硼烯厚度一致[23,28]。

模拟STM图像（图5）复现实验观测的条纹相（Mannix等）和均匀相（Feng等）。屈曲三角条纹呈纺锤状亮脊，β12产生椭圆状凸起被亮角包围，χ3显示菱形哑铃图案匹配Feng等报道的S2相。偏压依赖模拟（图6）进一步区分两种条纹构型：屈曲三角条纹在正负偏压间反差反转，β12条纹椭圆形状不随偏压变化，支持将Mannix条纹归属为屈曲三角，Feng条纹归属为β12。

总之，第一性原理计算表明自由站立屈曲三角、β12、χ3硼烯热力学不稳定且呈金属性。Ag(111)衬底稳定这些单层，每硼原子降低能量0.1-0.2 eV。模拟STM图像在观测表面相与底层晶格结构间建立了确定性映射，解决了硼烯在银上原子构型的长期争论。这些见解为硼烯基器件合理设计和指导金属衬底上实验合成奠定基础。

缩写

2D

二维

3D

三维

STM

扫描隧道显微镜