晶体无机材料中的力致发光：局部无序与弹性畸变假说

《Optical Materials: X》：Mechanoluminescence in crystalline inorganic materials: local disorder and the elastic distortion hypothesis

　　

在这项探索性工作中，研究人员旨在理解各种无机化合物在弹性加载范围（elastico-mechanoluminescence, EML）内表现出的力致发光（mechanoluminescence）现象，重点关注机械加载对活性位点局部畸变和对称性丧失的影响。为此，研究人员分析了若干相关晶相的弹性变形，并通过基于弹性常数的运动学分析表明，加载诱导的畸变在幅度上小于但接近于由Baur描述符（Baur descriptor）量化的内在结构畸变。尽管结构畸变先前已被提出作为化合物EML潜力的结构指纹，但其本质上是静态参数，必须通过机械加载下产生的动态畸变贡献来补充。通过将晶体变形近似为简单多边形图案的变形，研究人员能够为其他令人困惑的观察结果提供合理解释：特别是对静水压力和剪切敏感性的显著差异，在某些化合物中卸载时出现强烈的EML峰，以及初步紫外线照射是在加载下还是静止状态下进行的行为对比。

力致发光（mechanoluminescence, ML）是材料在机械加载下发光的能力，其中弹性力致发光（elastico-mechanoluminescence, EML）发生在弹性加载范围内，通常需紫外光（UV）预照射。尽管已有研究聚焦于电子跃迁、点缺陷陷阱等机制，但关于宏观应力场如何引起晶体结构畸变和局部无序变化，进而影响EML性能的报道很少。特别是，不同晶体结构即使化学成分相同，EML行为差异显著，纯剪切或纯静水压力下EML可能消失，且加载和卸载时均可能出现发光峰，这些现象缺乏统一解释。因此，研究人员旨在通过分析一系列EML活性化合物，结合弹性常数运动学分析，探讨局部畸变和对称性丧失在EML中的作用，为理解其机制提供新视角。该论文发表在《Optical Materials: X》。

研究人员主要采用了以下关键技术方法：1）使用PBEsol泛函的密度泛函理论（DFT）计算多种EML活性化合物（包括SrAl₂O₄、BaAl₂O₄、Ca₂MgSi₂O₇、ZrO₂、ZnS、BaZnOS、CaZnOS、Ba₄Si₆O₁₆和Ba₅Si₈O₂₁等）的弹性常数，获得多晶平均剪切模量（G）和体积模量（B）；2）利用Baur结构畸变指数（Ds）基于CIF文件通过Python代码量化配位多面体的畸变程度；3）通过运动学分析比较剪切和静水压力对弹性储能密度的贡献，引入剪切-体积能量比（R_S/V=3B/G）评估材料对剪切或压力的敏感倾向；4）采用简单多边形模型（三角形、正方形、六边形、八边形）模拟剪切诱导畸变与Ds的关系。

**Ten key observations**：通过文献综述总结出十个关键观察：UV预照射必要性（KO1）、EML信号弱于持久发光（PersL）（KO2）、共掺杂增强EML（KO3）、卸载时出现EML峰（KO4）、强度与机械功率成正比（KO5）、对静水压力和剪切的不同敏感性（KO6）、加载下UV照射影响EML响应（KO7）、卸载中间阶段出现强度极小值（KO8）、存在阈值应力（KO9）、活性位点局域对称性破缺有利于EML（KO10）。这些观察表明EML机制超出经典陷阱-去陷阱模型，需考虑晶体畸变和应力场传递。

**Tendencies through different chemical systems**：比较多个化学体系（SrAl₂O₄、BaAl₂O₄、Ca₂MgSi₂O₇/Ba₂MgSi₂O₇、ZrO₂、ZnS、BaZnOS/CaZnOS、Ba₄Si₆O₁₆/Ba₅Si₈O₂₁）后，发现较大结构畸变（高Ds）有利于EML，但Ds是静态参数，需补充动态畸变。例如，单斜SrAl₂O₄（Ds较大）相比六方相表现出EML，而BaAl₂O₄的低温柔性相（P6₃）因非中心对称和较大畸变而具有EML。在Ca₂MgSi₂O₇中，畸变较大（Ds=0.058）且EML活性强，而Ba₂MgSi₂O₇畸变较小（Ds=0.018）且EML弱。

**The driving force: Stress field**：宏观尺度上，通过弹性储能密度分析，剪切和静水压力对畸变和体积变化贡献不同。R_S/V大于7表示剪切敏感，小于6表示体积敏感，计算显示大多数EML化合物的R_S/V在5-10之间。结构尺度上，弹性各向异性导致畸变优先沿特定晶面发生。在SrAl₂O₄中，剪切沿（bc）面产生畸变，与内在畸变相当；在层状结构如Ca₂MgSi₂O₇中，剪切集中在易滑移层，对活性位点影响小。原子尺度上，Ds虽有用，但需考虑动态加载贡献，且Ds应加权并归一化至体积。

**Consequences: Kinematics of the crystal deformation**：通过简单多边形模型证明Ds与畸变角近似线性关系。在典型EML测试应力（几十MPa）下，动态畸变与内在畸变大小相当。这解释了为何剪切敏感化合物（如SrAl₂O₄）在加载和卸载时均有EML峰，而压力敏感化合物（如Ba₄Si₆O₁₆）卸载时无峰。还解释了加载下UV照射导致不同EML响应（KO7）的原因：陷阱填充依赖于加载状态，且卸载中间出现强度极小值（KO8）与畸变引起的陷阱深度变化相关。

讨论部分总结：EML行为难以仅由R_S/V解释，需考虑宏观畸变向活性位点的传递方式。结构畸变指数（Ds）的静态性质不能完全反映实际EML启动情况，应补充加载诱导的动态畸变估计。未来研究应系统比较剪切和静水压力行为，确定弹性各向异性张量，并关注应力场特征、局部畸变与晶体场强度之间的关联。研究结论部分翻译如下：科学进步源于提出问题。本文旨在鼓励进一步聚焦于EML基本理解的研究，最终形成基于应力场如何影响活性位点局部晶体场强度的清晰描述的物理图像和定量模型。未来研究应特别关注（i）应力场特征（方向、静水压力与剪切贡献）、（ii）由此产生的局部畸变以及（iii）局部晶体场强度之间的相关性。更好地理解这些相关性需通过系统比较剪切（或剪切与压力组合，如单轴测试）与纯静水压力加载下的行为，并结合原位结构监测，同时确定各向异性弹性张量并分析外部应力场如何在原子尺度传递至活性位点。结构畸变指数应理想地由加载诱导的动态畸变估计来补充，该动态畸变在许多情况下与静止时的内在结构畸变大小相当。但需注意，EML主要受陷阱群体特性、载流子俘获/释放过程、掺杂剂价态、缺陷复合物、压电场、位错和界面效应控制，畸变有助于陷阱达到适当深度并控制释放动力学，同时可能促进压电场进一步推动去陷阱。