编辑推荐：
本综述系统探讨了二维硅玻璃（SiO2）与α、β、γ型石墨炔（graphyne）杂化形成的Glaphynes系列材料。研究通过自洽电荷密度泛函紧束缚（SCC-DFTB）方法证实了其结构稳定性，并发现SiO2层可诱导β-glaphyne产生1.16 eV带隙，展现了通过界面工程调控二维材料电子结构（如能带开口、轨道重构）的巨大潜力，为纳米电子器件和半导体应用提供了新平台。
引言
近年来，纳米科学领域呈现出爆发式增长态势，这主要得益于二维纳米结构（如石墨烯）的预测与成功合成。石墨烯因其卓越的机械强度和热导率等特性，激发了人们对新型二维材料的探索热情。在此背景下，石墨炔作为一类由sp和sp2杂化碳原子组成的二维碳同素异形体，通过乙炔链接器插入石墨烯碳原子之间形成，其乙炔链的比例和排列方式形成了α、β和γ-石墨炔等不同类型。随着γ-石墨炔的实验合成，其在超级电容器、电池、光伏器件等领域的应用潜力被广泛研究。与此同时，硅酸盐家族中的二氧化硅纳米结构因其绝缘特性和约6.7 eV的宽禁带特性，在电子器件中扮演重要角色。
异质结构策略已成为增强纳米材料性能的重要手段。例如，将六方氮化硼域嵌入石墨烯片可实现力学性能调控；石墨烯/石墨炔异质结可用于碳基晶体管设计；不同石墨炔之间的异质结可调节热导率。层间堆叠是构建异质结构的另一重要方法，如石墨炔与XSe2堆叠可工程化带隙，石墨炔-氮化石墨烯异质结可用于纳米电容器等。最近，Iyengar等人通过实验制备了名为glaphene的杂化二维材料，该材料将二维硅玻璃单层与石墨烯结合，其层间相互作用超越典型范德华力，导致强层间杂化，使原本零带隙的石墨烯和8.2 eV带隙的硅玻璃组合形成带隙约3.6 eV的半导体材料，证明了通过电子邻近效应可实现能带工程。
本研究受glaphene启发，提出了一类名为glaphynes的异质结构，即通过将SiO2单层堆叠在α、β和γ-石墨炔上形成三种构型。研究旨在评估其能量和结构稳定性，并分析电子性质，探索通过堆叠策略实现带隙调控的潜力。
结果与讨论
通过密度泛函紧束缚方法对孤立单层进行几何优化后，结构参数与文献数据吻合良好：SiO2、α、β和γ-石墨炔的晶格参数偏差分别为2.5%、0.57%、0.74%和0.73%，γ角偏差可忽略不计。优化后的石墨炔乙炔链键长分别为：α-石墨炔（1.411, 1.234, 1.411 ?）、β-石墨炔（1.427, 1.227, 1.427 ?）和γ-石墨炔（1.438, 1.221, 1.438 ?）。SiO2单层的平均Si–O键长为1.63 ?，厚度为4.33 ?，与实验值一致。
glaphynes的构建采用3×3超胞的α和γ-石墨炔与β-石墨炔单胞，分别与4×4和2×2超胞的二维SiO2组合，初始层间距约3.5 ?。优化后所有构型均收敛至P6MM对称性，乙炔链键长仅轻微改变：α-glaphyne（1.421, 1.239, 1.421 ?）、β-glaphyne（1.576, 1.247, 1.576 ?）和γ-glaphyne（1.457, 1.228, 1.457 ?）。SiO2与石墨炔层的垂直间距分别为3.16、3.22和3.12 ?，与glaphene的3.20 ?相近。SiO2层厚度约0.44 nm，与实验值0.43 nm吻合。β-glaphyne通过7×7 SiO2与4×4 β-石墨炔超胞组合可进一步将晶格失配降至2.8%。
Mulliken布居分析显示α、β和γ-glaphyne的重叠布居平均增加约0.87、10.88和0.96 m|e|，主要源于碳2pz轨道分布变化。层分辨电荷分析表明界面存在微小但空间相干的电荷重分布（Δq ～ 10 m|e|），主要涉及硅氧层的氧原子和石墨炔的最近邻碳原子。
能量分析表明，γ-石墨炔的内聚能最低（-8.41 eV/atom），而glaphynes的内聚能略低于孤立石墨炔，提示SiO2覆盖层可增强结构稳定性。热力学模拟显示所有结构在300 K下保持完整；力学稳定性满足Born-Huang条件（C11> |C12|, C66> 0）；β-glaphyne的声子谱无虚频，高频模式（>66 THz）与sp1杂化碳原子一致。
电子结构分析显示，α-石墨炔的狄拉克锥位于K点，带隙为零；α-glaphyne的狄拉克锥移至Γ点，带隙仅10 meV，费米能级附近态主要由碳2p轨道贡献。β-石墨炔沿Γ–M路径存在狄拉克锥，带隙为零；而β-glaphyne在Γ点打开1.16 eV直接带隙，伴有能带平坦化，价带主要由碳2p轨道贡献，导带含碳2p和硅3d轨道贡献。DFT基准计算验证了DFTB结果的可靠性：半局域泛函预测β-glaphyne为金属性，杂化泛函给出0.63–1.14 eV带隙，DFTB结果与B3LYP-D3高度一致。γ-石墨炔带隙1.47 eV（M点），γ-glaphyne带隙增至1.51 eV，能带平坦化，电子态分布与α-glaphyne相似。
前沿轨道分布显示，α和γ-glaphyne的HOCO和LUCO相对离域，而β-glaphyne的轨道高度局域化，可能有利于电子传输。电荷密度分析表明石墨炔层向SiO2层存在轻微电荷转移。
声子态密度分析揭示了950–1150 cm–1区间的Si–O–C键合特征：α-glaphyne在980和1080 cm–1处出现强峰；β-glaphyne在1132 cm–1处有特征峰并伴肩峰，提示其Si–O–C键合最强；γ-glaphyne谱图与β-glaphyne相似。这些振动特征与glaphene中1050 cm–1附近的Si–O–C模式相对应。
结论
通过DFTB方法系统研究了glaphynes系列材料的结构与电子特性。α-glaphyne保持近零带隙，电子态重分布；β-glaphyne实现1.16 eV带隙开口，但能带平坦化可能降低电子迁移率；γ-glaphyne带隙小幅增加至1.51 eV。研究表明SiO2/石墨炔界面可通过电子邻近效应调控能带结构，但这种调控效果因石墨炔类型而异，并非均伴随Si–O–C键形成而必然发生带隙开口。该系列材料在纳米电子器件、半导体应用和表面催化等领域展现出通过界面工程精确调控电子性质的潜力。
模拟方法
研究采用DFTB+软件中的SCC-DFTB方法，使用matsci-0-3参数集，包含Lennard-Jones势和UFF力场描述 dispersion作用。几何优化收敛标准为SCC相互作用10–8a.u.，原子力10–5a.u.，布里渊区采样采用10×10×1 Monkhorst-Pack网格。glaphynes构建时通过超胞匹配最小化晶格失配，结构可视化采用OVITO和VESTA软件实现。
