本文档是一篇题为“Wide-Band Gap Binary Semiconductor P3N5 with Highly Anisotropic Optical Linearity and Nonlinearity”的原创性研究论文，发表于Inorganic Chemistry期刊，于2024年3月8日在线发布。通讯作者为中国科学院理化技术研究所功能晶体与激光技术重点实验室的康磊研究员和复旦大学微电子学院的严晓兰研究员。第一作者为中国科学院理化技术研究所的李世航。其他合作作者包括理化技术研究所的林哲帅研究员。该研究获得国家自然科学基金资助。
这项研究属于无机非金属材料、计算材料科学与非线性光学交叉领域。宽禁带二元半导体（如SiC、GaN、AlN、Ga2O3、Si3N4）因其优异的电子、光学和缺陷特性，在高压器件、日盲光电探测器和波导器件中具有重要应用。它们在100-400 nm紫外光谱范围内的透明性和可调谐特性，使其成为紫外非线性光学频率转换（如倍频）的潜在候选者。然而，实现高效紫外非线性光学频率转换不仅需要宽禁带和大的倍频效应，还需要足够大的双折射来实现相位匹配。例如，AlN虽然具有宽禁带和显著的倍频效应，但其双折射不足，限制了其实现高效相位匹配倍频输出的能力。这种不足源于其四面体[AlN4]结构单元的结构各向异性不足。因此，寻找兼具宽禁带、强倍频效应和大双折射的新型宽禁带二元半导体材料是该领域的一个重要目标。磷-氮键具有更多样的杂化形式，可以形成不同配位数的[Pnx]结构基元，这为设计具有更高结构各向异性的新型材料提供了可能。基于此背景，本研究旨在通过第一性原理计算，系统研究二元化合物P3N5两种结构（α相和γ相）的线性、非线性光学及缺陷性质，评估其对外部压力调控的响应，并探索其作为新型紫外非线性光学材料的潜力。研究的核心目标是阐明P3N5中结构、压力、光电性质与缺陷行为之间的内在关联。
详细研究流程：
研究对象确立与结构分析：研究以二元化合物P3N5的两种晶体结构为研究对象，分别是常压稳定相α-P3N5和高压相γ-P3N5。首先，研究团队从无机晶体数据库中获取了这两种结构的晶体学信息。α-P3N5属于单斜晶系，空间群为C2/c，是中心对称结构，其基本结构单元是四配位的[PN4]四面体，这些四面体通过共顶点或共边连接形成致密的三维网络。γ-P3N5是高压相，属于正交晶系，空间群为Imm2，是非中心对称结构，其结构单元包含四配位的[PN4]四面体和五配位的[PN5]金字塔，且两者呈平行排列。这一结构上的差异是导致其光学性质显著不同的根本原因。研究没有涉及实验样品的合成与制备，所有分析均基于已知的晶体结构数据。
计算方法与参数设置：本研究全程采用第一性原理计算方法。这是一种基于量子力学密度泛函理论的先进计算模拟手段，能够高精度预测材料的电子结构、光学性质及缺陷形成能等。研究特别指出，该方法在预测氮化物系统（如GaN和Si3N4）的光学性质方面已被证实具有高准确性，从而为本研究的可靠性提供了依据。具体计算细节见支撑信息，主要包括：使用平面波赝势方法进行结构优化和电子结构计算；采用杂化泛函或GW近似等方法以获得更精确的禁带宽度；利用密度泛函微扰理论计算线性介电函数和二阶非线性光学系数（即倍频系数）；通过分析电子能带结构和态密度，确定材料的带隙及价带顶、导带底的原子轨道贡献；最后，使用专门的点缺陷模拟软件包DASP计算了两种P3N5相中各种本征缺陷（包括空位Vn、Vp，反位缺陷Pn、Np，以及间隙原子Ni、Pi）在不同化学势条件下的形成能，并分析了缺陷能级在带隙中的位置。整个工作流程严格遵循计算材料学的标准范式。
高压相变与性质演化模拟：为了探究压力对材料性质的影响，研究系统计算了在0至11 GPa压力范围内，α-P3N5和γ-P3N5的相对焓变（δH）、带隙（Eg）、双折射（δn）和倍频系数（dij）随压力的变化关系。这是一个关键步骤，用于确定相变压力点以及压力对光电性能的调制作用。计算过程中，通过对施加不同压力的晶体结构进行弛豫，然后计算其能量和光学性质，从而得到一系列压力依赖性的数据点。
对比分析与性能评估：为了评估P3N5的性能水平，研究将其关键参数（带隙Eg、最大倍频系数dij、双折射δn）与目前已知的主流宽禁带二元半导体材料，包括4H-SiC、GaN、AlN、α-Si3N4和β-Ga2O3，进行了横向对比。这种对比旨在突出γ-P3N5在综合性能上的优势或特点。
相位匹配能力与紫外吸收边分析：对于非线性光学材料，实现相位匹配至关重要。研究计算了γ-P3N5的折射率色散曲线，并根据相位匹配条件，确定了其可实现相位匹配倍频输出的最短波长（λpm）。此外，还计算了其光学吸收谱，确定了材料的紫外吸收边（λuv），以评估其潜在的工作光谱范围。
缺陷性质系统研究：采用第一性原理计算软件包DASP，系统地研究了α-P3N5和γ-P3N5中所有可能本征缺陷的形成能。计算考虑了不同电荷态和不同化学环境（富N条件和贫N条件）。通过比较不同缺陷的形成能高低，可以预测在实验生长条件下哪种缺陷最容易形成、浓度最高。同时，计算了缺陷能级在带隙中的位置，并评估了最容易形成的缺陷（如Vn+1）对材料光学吸收谱的影响，从而判断缺陷是否会引入有害的光吸收，进而影响其在紫外波段的倍频输出效率。
主要研究结果：
结构、电子性质与高压相变：计算确认α-P3N5具有约5.92 eV的超宽禁带，与实验值（~5.87 eV）吻合良好；γ-P3N5的禁带宽度为~4.92 eV，属于日盲紫外波段。两者的导带底主要由P原子贡献，价带顶主要由N原子贡献。γ-P3N5带隙较小的原因是[PN5]中P-N键长比[PN4]中的更长。通过相对焓变分析，确定了α相到γ相的相变压力约为7.7 GPa，与先前计算结果一致。在压力作用下，两种相的带隙均因P-N键长缩短而增大，但γ相的带隙始终小于α相。
光学性质的各向异性与压力响应：γ-P3N5展现出高度各向异性的光学性质。其双折射（δn）在常压下高达~0.12，远大于α-P3N5的0.03，也优于对比材料AlN（~0.02）、Si3N4等。这归因于[PN5]金字塔基元比[PN4]四面体具有更高的结构各向异性（极化各向异性因子分别为32.84和6.71），且两种基元在结构中平行排列。γ-P3N5还表现出显著的倍频效应，其d24系数约为1.8 pm/V，是经典紫外非线性光学材料KDP（d36=0.39 pm/V）的约4.6倍。由于结构的中心对称性，α-P3N5没有倍频响应。在0-10 GPa压力下，γ-P3N5的Eg和dij变化幅度分别小于2%和10%，表现出优异的稳定性。这归因于其独特的晶格不可压缩性：计算发现其沿b轴的弹性常数远大于a轴和c轴，导致a-c面内更易压缩，而b方向几乎不可压缩，这种各向异性压缩反而增强了其结构各向异性和倍频效应。
综合性能对比优势：横向对比显示，γ-P3N5在关键性能指标上具有独特优势：其带隙（~4.92 eV）与β-Ga2O3相当；倍频系数（~1.85 pm/V）强于AlN和Si3N4（但弱于SiC和GaN）；而其双折射（~0.12）是所有对比材料中最大的。大双折射是实现相位匹配的关键。计算表明，γ-P3N5可以实现相位匹配倍频输出的最短波长（λpm）为286 nm，远低于其他已知的宽禁带二元半导体材料（如AlN、GaN等通常只能实现更长波长的输出）。这意味着γ-P3N5能将激光转换到更深的紫外波段。
缺陷性质及其对光学性能的影响：缺陷计算结果表明，在两种P3N5相中，氮空位（Vn）在大多数条件下形成能最低，是占主导地位的本征缺陷。在贫N条件下，Vn和反位缺陷Pn的形成能较低；在富N条件下，间隙N原子（Ni）的形成能显著降低。总体而言，γ-P3N5中各种缺陷的形成能略高于α-P3N5，意味着在相同生长条件下，γ-P3N5的缺陷密度可能更低。最重要的发现是，计算了最容易形成的缺陷（Vn+1）对γ-P3N5光学吸收的影响，发现在~300 nm附近引入了一个小的吸收峰，这可能对280-320 nm（UVB）波段的激光倍频输出产生轻微负面影响，但对320-400 nm（UVA）波段的影响基本可以忽略。研究强调，低浓度的缺陷基本不影响材料的双折射和非线性光学系数本身，但可能通过引入寄生吸收影响实际器件的转换效率。
结论： 本研究通过系统的第一性原理计算，首次深入揭示了高压相γ-P3N5作为一种新型宽禁带二元半导体在紫外非线性光学应用方面的巨大潜力。研究表明，从α-P3N5到γ-P3N5的结构相变（~7.7 GPa）伴随着光学性质的显著跃变：带隙从~6.03 eV降至~5.03 eV，双折射从0.029增至0.099，并从无倍频效应转变为具有显著倍频效应（~5.2倍KDP）。γ-P3N5集宽紫外日盲禁带（~4.9 eV）、强倍频效应（d24 ~1.8 pm/V）和大双折射（δn ~0.12）于一身，并能实现短至286 nm的相位匹配相干光输出，其综合性能超越了目前已知的SiC、GaN、AlN、Ga2O3、Si3N4等传统宽禁带二元半导体。此外，其关键光学性质在高压下变化微小，展现出良好的稳定性。对缺陷性质的深入分析表明，其本征缺陷对紫外光转换的影响有限，且γ相相对α相可能具有更低的缺陷密度。这些结果表明，高质量的γ-P3N5单晶有望成为应用于深紫外波段的高性能非线性光学材料。
研究亮点： 1. 重要发现：首次系统预测并论证了二元化合物P3N5的高压相（γ-P3N5）是一种极具潜力的紫外非线性光学材料，其最大的亮点在于实现了大双折射与大倍频系数在紫外材料中的结合，并因此能将相位匹配输出波长推进到286 nm的深紫外区域。 2. 研究方法的创新性与系统性：研究采用了从晶体结构分析到电子结构、线性/非线性光学性质、高压响应再到缺陷物理的完整第一性原理计算链条，对材料进行了多角度、多层次的全方位评估，逻辑严密，论证充分。特别地将高压调控、性能对比和缺陷分析有机结合，为材料的设计与性能评估提供了范本。 3. 研究对象的特殊性：聚焦于二元氮化物P3N5这一相对新颖的体系，并特别关注其高压相。通过揭示[PN5]非四面体结构基元在增强结构各向异性和光学非线性方面相对于传统[PN4]、[AlN4]等四面体基元的优势，为未来设计新型非线性光学材料提供了新的思路，即通过增加配位数和结构基元多样性来打破四面体结构的对称性限制。
其他有价值的内容： 研究中对“晶格不可压缩性”与“性质稳定性”关联的阐述颇具价值。通过计算弹性常数，发现γ-P3N5在b方向（[PN5]通过共边连接的方向）极其难以压缩，而在a-c面内相对容易压缩。这种独特的各向异性压缩行为不仅解释了其光学性质在高压下变化较小的原因，也暗示了该材料在极端条件下可能具有稳健的性能，拓展了其在高压光学器件方面的应用想象空间。同时，研究中对缺陷形成能随化学势（富N/贫N条件）变化的详细分析，为实验上通过控制生长条件来调控缺陷类型和浓度提供了重要的理论指导。