光电热力机械耦合机制在光电材料中的揭示：从非线性静力学到动力学

《Materials Today Physics》：Uncovering Opto-Electro-Thermo-Mechanical Couplings for Optoelectronic Materials: From Nonlinear Statics to Dynamics

　　

金属卤化物钙钛矿是光伏与光电子学中极具前景的候选材料，然而其结构行为，尤其是在光学、电学、热学与力学耦合场作用下的变形与稳定性，仍未明晰。因此，捕捉这种场—结构反馈对于预测非线性响应并评估真实服役条件下的机械可靠性至关重要。本文针对 MAPbI₃ 型铅卤钙钛矿，构建了一个全耦合光-电-热-力本构框架，通过统一多场诱导应力、电位移、光压敏感性以及热熵密度来实现。基于该耦合模型，研究人员采用矩形钙钛矿板作为代表性平台，结合 Mindlin-Reissner 运动学与 von Kármán 非线性进行分析。耦合系统总能量被纳入 Hamilton 原理，以推导控制方程；随后对方程进行无量纲化，并采用微分求积法（DQ，differential quadrature）与迭代特征值格式求解。研究人员对 MAPbI₃ 的电学与光学性质进行了对比验证，并进一步对多物理场作用下的力学响应进行基准测试，从而确保了耦合模型与数值方法的准确性。结果揭示出两条相互竞争的光学作用路径。其一是由电位移主导的、受光照驱动的光电耦合路径，它控制高场区间，并通过放大挠度、降低屈曲承载能力以及诱发动力学软化来起作用。其二是与光压敏感性相关的光机械路径，其主要贡献是一个内部光学反馈项，可提高基线切线刚度。因此，忽略该路径将会高估场致失稳与机械退化程度。

该文发表于《Materials Today Physics》，聚焦于金属卤化物钙钛矿，尤其是甲基铵铅碘（MAPbI₃）在真实服役环境中的多场耦合力学行为。此类材料因软晶格、离子—电子混合输运、强极化性以及优异的吸光与光电转换能力，在光伏和光电器件中展现出巨大潜力。随着器件效率不断逼近应用目标，研究重点已从单纯追求转换效率转向器件运行稳定性与机械可靠性。在实际工作条件下，光照、电偏压与温度波动往往同时存在，它们不仅是外部载荷，还会通过离子迁移、极化重构、光致介电增强、热致熵变等机制改变材料内部场分布，并反向影响变形、刚度、屈曲与振动特征。现有研究虽已涉及钙钛矿板和膜结构的非线性弯曲、振动、屈曲等问题，也开始将光、电、热因素纳入模型，但多数工作仍局限于单场分析，或仅停留于多场载荷下的应力—应变关系，尚未系统刻画光、电、热、力之间相互嵌套的反馈机制。特别是光致介电响应、光伏极化、热释电（pyroelectricity，温度变化诱发极化）、应变调制光压敏感性以及熵相关热耦合等关键机制常被忽略，导致对结构软化、失稳边界和动态响应的认识不完整。因此，开展该研究的根本目的在于建立一个统一、热力学自洽的多物理场耦合理论，用于更准确地预测钙钛矿光电材料在复杂工况下的非线性静力与动力学行为，并为器件机械稳健性设计提供理论基础。

在方法上，研究人员首先建立了适用于 MAPbI₃ 的全耦合光电热力机械本构模型，将压电（piezoelectricity，电场与应变耦合）、光致伸缩（photostriction，光诱导应变）、光致介电增强、光伏极化、热释电、光压敏感性及热熵密度统一纳入广义焓密度框架。随后，以矩形钙钛矿板为代表结构，采用 Mindlin-Reissner 厚板理论与 von Kármán 几何非线性建立控制方程，并将总能量代入 Hamilton 原理完成推导。数值上，对无量纲方程使用微分求积法（DQ）离散，并借助直接迭代方案求解非线性弯曲、临界屈曲、后屈曲和平板自由振动问题；同时通过文献中的压电、光伏、光致电学、光致应变与屈曲/振动结果进行交叉验证。本文未涉及生物样本或队列来源。

论文主体结果可概括如下。

1. Model validation
研究人员首先对所建模型的物理合理性和数值精度进行了系统验证。在压电响应方面，模型得到的纵向压电应力常数 e₃₃ 与既有理论和实验数据高度一致，说明电—机械耦合关系能够正确刻画 MAPbI₃ 晶体在应变作用下的极化变化。在光伏响应方面，由开路电压 V_oc 与光照强度之间的理论关系得到的结果与已有实验吻合，验证了光生载流子生成、分离与内建电场贡献的合理性。在光致电响应方面，模型预测的电容随光照增强的变化规律与参考实验高度一致，说明光致介电与光极化项的引入是有效的。在光机械耦合方面，光致应变—光强关系以及由光致伸缩系数换算得到的光机械应力尺度均与文献报道一致。这些验证共同表明，所提出的耦合本构模型不仅在单项物性上可靠，而且能够一致地整合多类场致响应。

2. Convergence studies and numerical methods validation
为保证计算结果的稳定性，研究人员对微分求积离散的网格密度进行了收敛分析。结果表明，当两个方向均取 17 个网格点时，非线性挠度、单轴和双轴临界屈曲载荷、非线性频率及频率比等关键量均已达到网格无关精度。在强耦合、高光强和高电场条件下，非线性弯曲、屈曲、后屈曲与自由振动的迭代解均可在较小误差容限内收敛，说明该数值策略对于求解强多场耦合、几何非线性问题是稳定且高效的。进一步与已有文献中关于挠度、屈曲边界、线性频率和非线性频率的结果比较，也均表现出良好一致性，从而支撑了本文理论与算法的可信度。

3. Nonlinear bending and stability
在非线性弯曲、屈曲与后屈曲分析中，研究人员比较了三类情形：忽略多物理耦合、只考虑电机械耦合、以及引入完整光电热力机械耦合。结果显示，若忽略耦合，结构挠度主要受光致伸缩驱动，随光强近似线性增长，电场影响很弱；一旦引入电机械耦合，光照诱导的介电增强会放大极化与静电能，使电场对刚度调制作用明显增强；而在完整耦合框架下，挠度对外场的依赖呈现更显著非线性，高场下出现明显软化。研究特别指出，完整耦合模型在零外场基线附近反而可能表现出更高的等效切线刚度，这并非材料本征“变硬”，而是内部电学和光学广义变量在边界约束下被机械变形激活后形成的反馈恢复作用。随着外加光场、电场和温升增强，软化效应逐渐占主导，导致挠度增大、屈曲能力下降。单轴和双轴屈曲载荷均随光强增大而单调降低，且在高电场下退化更为明显，双轴压缩情形最不稳定。后屈曲路径分析进一步表明，多场耦合显著降低后屈曲承载能力。

4. Effects of mechanically related coupling terms
通过参数消除分析，研究人员区分了压电项、热机械项和光机械项对弯曲与稳定性的贡献。结果表明，忽略压电效应只引起较小变化，说明在本文考察的载荷水平下，压电系数对宏观挠度和屈曲的影响相对次要。相比之下，忽略热致应力会低估非线性挠度并高估临界屈曲载荷，而且这种误差随光强增加而增大，说明光热转换及热膨胀相关项在强照明条件下更加重要。最显著的是光机械耦合：若去除由内部光学变量介导的双向反馈，基线挠度明显上升，临界屈曲载荷和后屈曲路径整体下移，说明该光学反馈通道为结构提供了额外恢复作用，是控制基线刚度与稳定性的关键机制。

5. Effects of optoelectrical coupling
关于光电耦合，研究人员采用“单因素关闭”方式考察多场诱导电位移和光压敏感性的作用。结果显示，若抑制电位移关系中的耦合项，电场对挠度和屈曲的调制显著减弱，说明电位移通道是主要的场致软化来源。若仅忽略光压敏感性，结构仍对光电场高度敏感，但会表现出更大的挠度和更低的屈曲载荷，说明缺失内部光反馈将高估失稳风险。进一步分析表明，照明通过光致介电增强与光极化显著强化电场相关内能与 Maxwell 应力效应，使高场区域的挠曲增长更快、屈曲与后屈曲能力衰减更明显。该结果直接支持摘要中提出的第一条竞争路径：由电位移主导的光电耦合控制高场区的变形放大与动力学软化。

6. Effects of thermo-electrical and thermo-optical couplings
在热相关耦合方面，忽略热电耦合和热光耦合会低估弯曲挠度并高估结构稳定性，但二者影响程度并不相同。热电耦合由热释电系数主导，因其数值较小，对刚度调制作用有限；而热光耦合则同时受到光热升温与光压热敏性影响，对弯曲、屈曲和后屈曲路径均有更显著影响，尤其在强照明下表现突出。这说明在多场服役条件下，温度不仅通过热膨胀进入力学问题，更通过熵相关项和光学载荷热敏性改变能量分配与稳定边界。

7. Field distributions and mode characteristics in bending/buckling
研究人员还分析了弯曲与屈曲状态下的位移、转角、电势、光势、中面应力、电位移 D_z 以及光压敏感性 Φ_z 的空间分布。对于简支边界（SSSS），板的柔顺性较高，挠度及内部电光变量幅值通常大于四边固支（CCCC）情形；固支虽抑制位移和转角，但会增强曲率、应变梯度及边界层效应，从而改变电势与光势局域化形态。D_z 与 Φ_z 在静态问题中表现出近似相同的中心主导拓扑，说明二者均主要受共同的变形模态控制，但幅值响应机制不同：D_z 更受介电与光极化项驱动而表现出整体抬升，Φ_z 则更表现为局部峰值锐化。该结果支持摘要中第二条竞争路径：光机械路径主要通过内部光反馈项提高基线切线刚度，而非像自由空间辐射压那样直接施加强大表面载荷。

8. Nonlinear vibration response
在动力学方面，若忽略耦合，非线性频率随光强增加而下降，频率比上升，整体近似线性；引入电机械耦合后，高光强下电场对频率的调制显著增强；在完整耦合模型中，频率和频率比对多场变化呈高度非线性响应。研究人员指出，完整耦合会提高零外场参考态下的基线频率，这是因为振动中的曲率与转角会激活内部电学和光学变量，储存额外场能并形成恢复反馈，从而抬升等效刚度。该机制在固支边界下更明显。随着光场和电场增大，光电耦合导致频率更快下降、频率比更快上升，体现出明显动力学软化。相比之下，去除光机械反馈会降低非线性频率、增大频率比并缩短振动周期，说明光学内部反馈对动力学稳定性同样至关重要。热相关耦合中，热光效应仍显著强于热电效应，对频率降低和频率比升高的贡献更大。

9. Vibration-mode, stress, and field distributions
在自由振动模态与场分布上，SSSS 板表现出更平滑、更广泛的位移和转角分布，W、电势 Ψ_E 与光势 Ψ_I 呈中心单峰、无节点形态；CCCC 板虽然整体位移较小，但由于零转角边界导致应变梯度局域化，其 Ψ_E 与 Ψ_I 峰值可能反而更高。振动状态下的中面应力分布较静弯曲更复杂，法向应力呈压缩主导并伴随局部集中，剪应力 τ_XY 呈多胞状交替分布。D_z 与 Φ_z 在振动工况下总体接近均匀，仅有弱空间调制，且两种边界条件的整体水平接近，说明在动力学场景中，这些输出量更多受模态控制，而外加光电输入主要通过改变平均偏置放大 D_z，对增量型 Φ_z 的影响相对有限。

总体而言，本文最核心的结论是：MAPbI₃ 类钙钛矿的结构响应不能被简化为单纯的外场加载问题，而应被理解为一个由内部诱导电场、光场与热场变量共同参与的能量反馈系统。弱耦合时，材料表现为以光照主导、近似线性的力学和振动响应；全耦合时，系统在高场区进入显著非线性软化状态。压电与热释电效应影响相对有限，而热机械、热光及尤其是光电、光机械通道对弯曲、屈曲、后屈曲及振动都有实质性控制作用。其中，电位移主导的光电路径会在高场下放大挠度、降低屈曲与后屈曲承载力，并加剧动态软化；而与光压敏感性相关的光机械路径则提供内部光学反馈刚化，忽略该路径会系统性高估场致失稳和机械退化。论文因此为理解钙钛矿器件在复杂服役条件下的多场结构行为提供了统一理论框架，也为高可靠性光电器件的结构设计与稳定性评估提供了重要理论依据。

研究结论部分可译为：本文提出了一种适用于铅卤钙钛矿的全耦合光电热力机械本构模型，并据此研究了多物理场条件下钙钛矿板的机械变形、电学与光学响应、稳定性及振动特征。主要结论包括：当耦合效应较弱或被忽略时，钙钛矿板表现出以光场主导、近似线性的力学与振动响应；而在引入完整耦合后，控制机制转变为内部诱导场变量共同参与的能量反馈过程。忽略压电和热释电效应仅带来轻微变化；相反，忽略热机械和热光效应会低估非线性挠度与频率比，并高估屈曲稳定性与非线性频率。光电效应主要由电位移通道控制，且在光照下受到光致介电增强与光极化强化，使响应在高场区高度敏感；该耦合会导致更大的挠度、更低的屈曲与后屈曲承载力、更快的频率下降以及更陡的频率比上升。光机械耦合及光压敏感性通道提供了耦合诱导的基线光反馈刚化，若忽略这一反馈，则会放大基线挠度、削弱屈曲/后屈曲能力并软化动力学响应。电势和光势总体上跟随挠度变化，在 SSSS 条件下表现为无节点中心峰；CCCC 则抑制位移与转角，但会增强曲率/应变梯度与边界层，并可能在自由振动中提高 Ψ_E 与 Ψ_I 的峰值。电位移与光压敏感性在空间上呈现近乎相同且受边界条件影响的拓扑分布；静态下，增大场强会显著抬升 D_z 并扩大 SSSS 与 CCCC 的差异，而 Φ_z 的增强主要表现为局部峰值锐化；在振动状态下，二者则近乎均匀、对边界不敏感，并主要受模态控制。